SEARCH · Search NASA

Results for “Shape 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1996

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗

NASA NTRS2011

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗

DOE OSTI2019

Shape memory alloy engine for high efficiency low-temperature gradient thermal to electrical conversion

More than half of the energy generated worldwide is lost as unused thermal energy because of the lack of efficient methodology for harnessing the low-grade heat. In this paper, we demonstrate that shape-memory alloy can be an effective mechanism for recovering low-grade heat. Shape memory alloys exhibit thermally induced martensite to austenite phase transformation and super-elasticity (stress-induced martensitic transformation). Employing these two characteristics, we demonstrate a thermal engine for harnessing waste energy through all modes of heat transfer: convection, conduction, and radiation. In this work, we performed material and heat transfer analysis for achieving high frequency, sustainable and efficient operation of our engine. An optimized shape memory alloy engine generated 36 W per kilogram or 234 kW of electricity per cubic meter of active material. A continuous three-day operation of several SMA engines could generate 7.2 kWh of electricity when installed on a 500 m long hot pipe network. This generated power can potentially reduce the carbon footprint by 5.1 kg of CO 2 illustrating the promise of this technology for addressing climate change.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS2013

Geometric Effects on the Amplification of First Mode Instability Waves

The effects of geometric changes on the amplification of first mode instability waves in an external supersonic boundary layer were investigated using numerical techniques. Boundary layer stability was analyzed at Mach 6 conditions similar to freestream conditions obtained in quiet ground test facilities so that results obtained in this study may be applied to future test article design to measure first mode instability waves. The DAKOTA optimization software package was used to optimize an axisymmetric geometry to maximize the amplification of the waves at first mode frequencies as computed by the 2D STABL hypersonic boundary layer stability analysis tool. First, geometric parameters such as nose radius, cone half angle, vehicle length, and surface curvature were examined separately to determine the individual effects on the first mode amplification. Finally, all geometric parameters were allowed to vary to produce a shape optimized to maximize the amplification of first mode instability waves while minimizing the amplification of second mode instability waves. Since first mode waves are known to be most unstable in the form of oblique wave, the geometries were optimized using a broad range of wave frequencies as well as a wide range of oblique wave angles to determine the geometry that most amplifies the first mode waves. Since first mode waves are seen most often in flows with low Mach numbers at the edge of the boundary layer, the edge Mach number for each geometry was recorded to determine any relationship between edge Mach number and the stability of first mode waves. Results indicate that an axisymmetric cone with a sharp nose and a slight flare at the aft end under the Mach 6 freestream conditions used here will lower the Mach number at the edge of the boundary layer to less than 4, and the corresponding stability analysis showed maximum first mode N factors of 3.

Kirk, Lindsay C.↗

NASA NTRS1997

Variational Methods in Design Optimization and Sensitivity Analysis for Two-Dimensional Euler Equations

Variational methods (VM) sensitivity analysis employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗

NASA NTRS2017

Mitigation of Engine Inlet Distortion Through Adjoint-Based Design

The adjoint-based design capability in FUN3D is extended to allow efficient gradient- based optimization and design of concepts with highly integrated aero-propulsive systems. A circumferential distortion calculation, along with the derivatives needed to perform adjoint-based design, have been implemented in FUN3D. This newly implemented distortion calculation can be used not only for design but also to drive the existing mesh adaptation process and reduce the error associated with the fan distortion calculation. The design capability is demonstrated by the shape optimization of an in-house aircraft concept equipped with an aft fuselage propulsor. The optimization objective is the minimization of flow distortion at the aerodynamic interface plane of this aft fuselage propulsor.

Ordaz, Irian↗

DOE OSTI2023

STRESS AND FATIGUE ANALYSIS OF AIR-TO-REFRIGERANT HEAT EXCHANGERS WITH NON-ROUND TUBE SHAPES

Air-to-refrigerant heat exchangers are fundamental components in Heating, Ventilation, Air Conditioning, and Refrigeration (HVAC&R) systems. Recent research has shown that heat exchangers utilizing small hydraulic diameter, non-round, shape-optimized tubes to be promising next design update. Such tubes can improve performance, reduce material and lower refrigerant charge. However, such tube designs may be more susceptible to fatigue issues. Therefore, stress and fatigue analyses of such heat exchangers are critical to ensure structural integrity, reliability, and manufacturability. In the present work, a mechanical model is developed to obtain the stress distributions in heat exchangers with round, ellipse, and airfoil-shaped tubes. The results indicate that under 3.45 MPa internal pressure, the maximum stresses don't surpass 190 MPa.

Zhang, Mingkan↗

DOE OSTI2019

CFD-Guided Combustion System Optimization of a Gasoline Range Fuel in a Heavy-Duty Compression Ignition Engine Using Automatic Piston Geometry Generation and a Supercomputer

A computational fluid dynamics (CFD) guided combustion system optimization was conducted for a heavy-duty diesel engine running with a gasoline fuel that has a research octane number (RON) of 80. The goal was to optimize the gasoline compression ignition (GCI) combustion recipe (piston bowl geometry, injector spray pattern, in-cylinder swirl motion, and thermal boundary conditions) for improved fuel efficiency while maintaining engine-out NO x within a 1-1.5 g/kW-hr window. The numerical model was developed using the multi-dimensional CFD software CONVERGE. A two-stage design of experiments (DoE) approach was employed with the first stage focusing on the piston bowl shape optimization and the second addressing refinement of the combustion recipe. For optimizing the piston bowl geometry, a software tool, CAESES, was utilized to automatically perturb key bowl design parameters. This led to the generation of 256 combustion chamber designs evaluated at several engine operating conditions. The second DoE campaign was conducted to optimize injector spray patterns, fuel injection strategies and in-cylinder swirl motion for the best performing piston bowl designs from the first DoE campaign. This comprehensive optimization study was performed on a supercomputer, Mira, to accelerate the development of an optimized fuel-efficiency focused design. Compared to the production combustion system in the baseline engine, the new combustion recipe from this study showed significantly improved closed-cycle fuel efficiency across key engine operating points while meeting the engine-out NO x targets. Optimized piston bowl designs and injector spray patterns were predicted to provide enhanced in-cylinder air utilization and more rapid mixing-controlled combustion, thereby leading to a fuel efficiency improvement. In addition, shifting the engine thermal boundary conditions toward leaner operation was also key to the improved fuel efficiency.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1999

Viscous Design of TCA Configuration

The goal in this effort is to redesign the baseline TCA configuration for improved performance at both supersonic and transonic cruise. Viscous analyses are conducted with OVERFLOW, a Navier-Stokes code for overset grids, using PEGSUS to compute the interpolations between overset grids. Viscous designs are conducted with OVERDISC, a script which couples OVERFLOW with the Constrained Direct Iterative Surface Curvature (CDISC) inverse design method. The successful execution of any computational fluid dynamics (CFD) based aerodynamic design method for complex configurations requires an efficient method for regenerating the computational grids to account for modifications to the configuration shape. The first section of this presentation deals with the automated regridding procedure used to generate overset grids for the fuselage/wing/diverter/nacelle configurations analysed in this effort. The second section outlines the procedures utilized to conduct OVERDISC inverse designs. The third section briefly covers the work conducted by Dick Campbell, in which a dual-point design at Mach 2.4 and 0.9 was attempted using OVERDISC; the initial configuration from which this design effort was started is an early version of the optimized shape for the TCA configuration developed by the Boeing Commercial Airplane Group (BCAG), which eventually evolved into the NCV design. The final section presents results from application of the Natural Flow Wing design philosophy to the TCA configuration.

Krist, Steven E.↗

NASA NTRSDate not supplied

High-Fidelity Aeropropulsive Optimization of a Mail-Slot Distributed Electric Propulsion System for the SUSAN Electrofan

Hybrid- and all-electric aircraft concepts use electric motors for power rather than a conventional jet engine. Electric propulsors open the door to new ways to synergistically integrate the propulsion system with the airframe. For example, many small electric propulsors can be distributed along the wing to increase the effective bypass ratio for better overall efficiency. Furthermore, these propulsors can be attached to the wing surface for boundary layer ingestion(BLI) to further the efficiency gains. However, these novel methods of aeropropulsive integration also create challenges such as nonuniform inflow and complex nacelle geometries. Here we use gradient-based aerodynamic shape optimization to address the design challenges of the wing-mounted distributed electric propulsion system of the Subsonic Single Aft Engine (SUSAN)concept. In doing so, we aim to more accurately benchmark the flow power of SUSAN’s mail slot propulsors relative to a conventional propulsion system in both an isolated and BLI configuration. Our preliminary results found relative to an optimized podded propulsor the optimized mailslot and BLI mailslot design required 8% and 17% more flow power respectively.The methods and key design insights also apply to other aircraft concepts that utilize distributed electric propulsion and boundary layer ingestion.

CAS↗

NASA NTRS2012

Parametric Deformation of Discrete Geometry for Aerodynamic Shape Design

We present a versatile discrete geometry manipulation platform for aerospace vehicle shape optimization. The platform is based on the geometry kernel of an open-source modeling tool called Blender and offers access to four parametric deformation techniques: lattice, cage-based, skeletal, and direct manipulation. Custom deformation methods are implemented as plugins, and the kernel is controlled through a scripting interface. Surface sensitivities are provided to support gradient-based optimization. The platform architecture allows the use of geometry pipelines, where multiple modelers are used in sequence, enabling manipulation difficult or impossible to achieve with a constructive modeler or deformer alone. We implement an intuitive custom deformation method in which a set of surface points serve as the design variables and user-specified constraints are intrinsically satisfied. We test our geometry platform on several design examples using an aerodynamic design framework based on Cartesian grids. We examine inverse airfoil design and shape matching and perform lift-constrained drag minimization on an airfoil with thickness constraints. A transport wing-fuselage integration problem demonstrates the approach in 3D. In a final example, our platform is pipelined with a constructive modeler to parabolically sweep a wingtip while applying a 1-G loading deformation across the wingspan. This work is an important first step towards the larger goal of leveraging the investment of the graphics industry to improve the state-of-the-art in aerospace geometry tools.

Anderson, George R.↗

NASA NTRS1999

Approximation of the Newton Step by a Defect Correction Process

In this paper, an optimal control problem governed by a partial differential equation is considered. The Newton step for this system can be computed by solving a coupled system of equations. To do this efficiently with an iterative defect correction process, a modifying operator is introduced into the system. This operator is motivated by local mode analysis. The operator can be used also for preconditioning in Generalized Minimum Residual (GMRES). We give a detailed convergence analysis for the defect correction process and show the derivation of the modifying operator. Numerical tests are done on the small disturbance shape optimization problem in two dimensions for the defect correction process and for GMRES.

Arian, E.↗

NASA NTRS2005

Design and Optimization of the SPOT Primary Mirror Segment

The 3m Spherical Primary Optical Telescope (SPOT) will utilize a single ring of 0.86111 point-to-point hexagonal mirror segments. The f2.85 spherical mirror blanks will be fabricated by the same replication process used for mass-produced commercial telescope mirrors. Diffraction-limited phasing will require segment-to-segment radius of curvature (ROC) variation of approx.1 micron. Low-cost, replicated segment ROC variations are estimated to be almost 1 mm, necessitating a method for segment ROC adjustment & matching. A mechanical architecture has been designed that allows segment ROC to be adjusted up to 400 microns while introducing a minimum figure error, allowing segment-to-segment ROC matching. A key feature of the architecture is the unique back profile of the mirror segments. The back profile of the mirror was developed with shape optimization in MSC.Nastran(TradeMark) using optical performance response equations written with SigFit. A candidate back profile was generated which minimized ROC-adjustment-induced surface error while meeting the constraints imposed by the fabrication method. Keywords: optimization, radius of curvature, Pyrex spherical mirror, Sigfit

Budinoff, Jason G.↗

DOE OSTI2020

Gradient-based constrained optimization using a database of linear reduced-order models

A methodology grounded in model reduction is presented for accelerating the gradient-based solution of a family of linear or nonlinear constrained optimization problems where the constraints include at least one linear Partial Differential Equation (PDE). A key component of this methodology is the construction, during an offline phase, of a database of pointwise, linear, Projection-based Reduced-Order Models (PROM)s associated with a design parameter space and the linear PDE(s). A parameter sampling procedure based on an appropriate saturation assumption is proposed to maximize the efficiency of such a database of PROMs. A real-time method is also presented for interpolating at any queried but unsampled parameter vector in the design parameter space the relevant sensitivities of a PROM. The practical feasibility, computational advantages, and performance of the proposed methodology are demonstrated for several realistic, nonlinear, aerodynamic shape optimization problems governed by linear aeroelastic constraints.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1997

Sensitivity Analysis and Optimization of Aerodynamic Configurations with Blend Surfaces

A novel (geometrical) parametrization procedure using solutions to a suitably chosen fourth order partial differential equation is used to define a class of airplane configurations. Inclusive in this definition are surface grids, volume grids, and grid sensitivity. The general airplane configuration has wing, fuselage, vertical tail and horizontal tail. The design variables are incorporated into the boundary conditions, and the solution is expressed as a Fourier series. The fuselage has circular cross section, and the radius is an algebraic function of four design parameters and an independent computational variable. Volume grids are obtained through an application of the Control Point Form method. A graphic interface software is developed which dynamically changes the surface of the airplane configuration with the change in input design variable. The software is made user friendly and is targeted towards the initial conceptual development of any aerodynamic configurations. Grid sensitivity with respect to surface design parameters and aerodynamic sensitivity coefficients based on potential flow is obtained using an Automatic Differentiation precompiler software tool ADIFOR. Aerodynamic shape optimization of the complete aircraft with twenty four design variables is performed. Unstructured and structured volume grids and Euler solutions are obtained with standard software to demonstrate the feasibility of the new surface definition.

Thomas, A. M.↗

DOE OSTI2021

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗

NASA NTRSDate not supplied

High-Fidelity Aerodynamic Analysis and Optimization of the SUSAN Electrofan Concept

The Subsonic Single Aft Engine (SUSAN) Electrofan is a new single-aisle regional jet class transport aircraft being studied by NASA, which benefits from reduced fuel consumption through a reduction in the number of hydrocarbon fuel-burning engines from two to one. This is achieved through a hybrid-electric system architecture, where a series of generators are used to extract power from a fuel-consuming tail cone thruster to drive electric propulsors mounted on the wings. To support the development of these propulsion systems, which involve high levels of propulsion-airframe integration, this paper presents high-fidelity analyses through computational fluid dynamics (CFD) simulations provided by the Launch Ascent and Vehicle Aerodynamics (LAVA) framework, which also includes capabilities for modeling the effect of active propulsors. A trade space exploration is performed to investigate the advantages of several wing mounted distributed electric propulsion system configurations, including over-wing, under-wing, and trailing-edge arrangements. For the aft fuselage propulsor, focus is on the sizing of the inlet, core, and bypass ducts. High-fidelity aerodynamic shape optimization is also used to help minimize the inlet distortion intensity experienced by the aft propulsor fan through a reshaping of the fuselage.

CAS↗

DOE OSTI2024

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗