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,063 records · Page 59

Engineering Overview of a Multidisciplinary HSCT Design Framework Using Medium-Fidelity Analysis Codes

An objective of the HPCC Program at NASA Langley has been to promote the use of advanced computing techniques to more rapidly solve the problem of multidisciplinary optimization of a supersonic transport configuration. As a result, a software system has been designed and is being implemented to integrate a set of existing discipline analysis codes, some of them CPU-intensive, into a distributed computational framework for the design of a High Speed Civil Transport (HSCT) configuration. The proposed paper will describe the engineering aspects of integrating these analysis codes and additional interface codes into an automated design system. The objective of the design problem is to optimize the aircraft weight for given mission conditions, range, and payload requirements, subject to aerodynamic, structural, and performance constraints. The design variables include both thicknesses of structural elements and geometric parameters that define the external aircraft shape. An optimization model has been adopted that uses the multidisciplinary analysis results and the derivatives of the solution with respect to the design variables to formulate a linearized model that provides input to the CONMIN optimization code, which outputs new values for the design variables. The analysis process begins by deriving the updated geometries and grids from the baseline geometries and grids using the new values for the design variables. This free-form deformation approach provides internal FEM (finite element method) grids that are consistent with aerodynamic surface grids. The next step involves using the derived FEM and section properties in a weights process to calculate detailed weights and the center of gravity location for specified flight conditions. The weights process computes the as-built weight, weight distribution, and weight sensitivities for given aircraft configurations at various mass cases. Currently, two mass cases are considered: cruise and gross take-off weight (GTOW). Weights information is obtained from correlations of data from three sources: 1) as-built initial structural and non-structural weights from an existing database, 2) theoretical FEM structural weights and sensitivities from Genesis, and 3) empirical as-built weight increments, non-structural weights, and weight sensitivities from FLOPS. For the aeroelastic analysis, a variable-fidelity aerodynamic analysis has been adopted. This approach uses infrequent CPU-intensive non-linear CFD to calculate a non-linear correction relative to a linear aero calculation for the same aerodynamic surface at an angle of attack that results in the same configuration lift. For efficiency, this nonlinear correction is applied after each subsequent linear aero solution during the iterations between the aerodynamic and structural analyses. Convergence is achieved when the vehicle shape being used for the aerodynamic calculations is consistent with the structural deformations caused by the aerodynamic loads. To make the structural analyses more efficient, a linearized structural deformation model has been adopted, in which a single stiffness matrix can be used to solve for the deformations under all the load conditions. Using the converged aerodynamic loads, a final set of structural analyses are performed to determine the stress distributions and the buckling conditions for constraint calculation. Performance constraints are obtained by running FLOPS using drag polars that are computed using results from non-linear corrections to the linear aero code plus several codes to provide drag increments due to skin friction, wave drag, and other miscellaneous drag contributions. The status of the integration effort will be presented in the proposed paper, and results will be provided that illustrate the degree of accuracy in the linearizations that have been employed.

Weston, R. P.↗

National Combustion Code: A Multidisciplinary Combustor Design System

The Internal Fluid Mechanics Division conducts both basic research and technology, and system technology research for aerospace propulsion systems components. The research within the division, which is both computational and experimental, is aimed at improving fundamental understanding of flow physics in inlets, ducts, nozzles, turbomachinery, and combustors. This article and the following three articles highlight some of the work accomplished in 1996. A multidisciplinary combustor design system is critical for optimizing the combustor design process. Such a system should include sophisticated computer-aided design (CAD) tools for geometry creation, advanced mesh generators for creating solid model representations, a common framework for fluid flow and structural analyses, modern postprocessing tools, and parallel processing. The goal of the present effort is to develop some of the enabling technologies and to demonstrate their overall performance in an integrated system called the National Combustion Code.

Stubbs, Robert M.↗

Parameter Trade Studies For Coherent Lidar Wind Measurements of Wind from Space

The design of an orbiting wind profiling lidar requires selection of dozens of lidar, measurement scenario, and mission geometry parameters; in addition to prediction of atmospheric parameters. Typical mission designs do not include a thorough trade optimization of all of these parameters. We report here the integration of a recently published parameterization of coherent lidar wind velocity measurement performance with an orbiting coherent wind lidar computer simulation; and the use of these combined tools to perform some preliminary parameter trades. We use the 2006 NASA Global Wind Observing Sounder mission design as the starting point for the trades.

Kavaya, Michael J.↗

Computations of Aerodynamic Performance Databases Using Output-Based Refinement

Objectives: Handle complex geometry problems; Control discretization errors via solution-adaptive mesh refinement; Focus on aerodynamic databases of parametric and optimization studies: 1. Accuracy: satisfy prescribed error bounds 2. Robustness and speed: may require over 105 mesh generations 3. Automation: avoid user supervision Obtain "expert meshes" independent of user skill; and Run every case adaptively in production settings.

Nemec, Marian↗

In-situ Mechanical Neutron Diffraction Loading Characteristics of GRCop-84 Fabricated by SLM

GRCop-84 is a precipitation strengthened alloy composed of Cu-8Cr-4Nb at % with Cr2Nb precipitates that provide dispersion and precipitation strengthening characteristics and limited solubility in the Cu matrix. The particle role of Cr2Nb(C15Laves) is unusual only contributing 1/3 of strengthening at high temperatures while the matrix provides the remainder. The particles mechanically and thermally stabilize the matrix retaining purity and preventing coarsening and loss of strength. At high temperatures (50-85%TmCu), GRCop-84 provides the best thermal and mechanical properties of available alloys. GRCop-84 is currently in development for reusable launch vehicles including the Space Launch System (SLS) with a focus on fabrication with additive manufacturing (AM) techniques. GRCop-84 is an optimal material for consolidating with AM. The base material is costly, the production times are long, and more geometry control can considerably improve cooling efficiency. Development of AMGR Cop-84 with selective laser melting (SLM) has rapidly progressed due to ease of printing and limited operator adjustment between builds, but the necessary knowledge-base of thermal history and stress state during builds is still under development.

Minneci, R.↗

Extreme Environment Hot Fire Durability of Post-Processed Additively Manufactured GRCop-Alloy Combustion Chambers in LOX/Hydrogen and LOX/Methane

Extreme environment survivability of metal additive manufactured (AM) GRCop-alloy thrust chambers has been demonstrated in different bi-propellants at near stoichiometric and even oxygen rich combustion. GRCop-alloy chambers tested at NASA Marshall Space Flight Center (MSFC) have accumulated over 26,000 seconds of hot fire duration and over 500 starts. These chambers are produced using an AM process called laser powder bed fusion (L-PBF). A major feature of this process is high wall surface roughness which can be customized in post processing to leverage various performance advantages. Post-processing can include heat treatment, final machining, polishing, and welding and is key to hardware survivability in extreme environments. Surface finish enhancement techniques were applied to the hot wall and coolant channels to reduce the overall total heat load to the chamber walls. Performance optimization of various thrust class TCA’s is a strategic technology goal of NASA MSFC. Three different chamber geometries using cryogenic methane and de-ionized water as coolants were hot fire tested to obtain their life cycle, pressure drop, and heat load performances. Several 1.2K lbf LOX/H2 chambers, 1K lbf LOX/CH4 chambers, and 7K lbf LOX/CH4 chambers were tested. All post-processed configurations performed remarkably well when subjected to extreme hot fire test conditions. Streaking and blanching due to localized oxygen rich conditions were observed on some test articles in their as-built surface finish state. However, this is actually a function of the injector mixing and only serves to further establish the durability of L-PBF GRCop-alloy thrust chambers. Several different polishing techniques were applied to the hot wall and integrated coolant channels prior to hot fire testing and their performances assessed. Overall, AM produced GRCop chambers are extremely reliable, durable, and customizable to the desired performance metrics.

grcop↗

Computer program optimizes design of nuclear radiation shields

Computer program, OPEX 2, determines minimum weight, volume, or cost for shields. Program incorporates improved coding, simplified data input, spherical geometry, and an expanded output. Method is capable of altering dose-thickness relationship when a shield layer has been removed.

Lahti, G. P.↗

Optimization of circular orifice jets mixing into a heated crossflow in a cylindrical duct

To examine the mixing characteristics of circular jets in an axisymmetric can geometry, temperature measurements were obtained downstream of a row of cold jet injected into a heated cross stream. The objective was to obtain uniform mixing within one duct radius downstream of the leading edge of the jet orifices. An area weighted standard deviation of the mixture fraction was used to help quantify the degree of mixedness at a given plane. Non-reacting experiments were conducted to determine the influence of the number of jets on the mixedness in a cylindrical configuration. Results show that the number of orifices significantly impacts the mixing characteristics of jets injected from round hole orifices in a can geometry. Optimum mixing occurs when the mean jet trajectory aligns with the radius which divides the cross sectional area of the can into two equal parts at one mixer radius downstream of the leading edge of the orifice. The optimum number of holes at momentum-flux ratios of 25 and 52 is 10 and 15 respectively.

Kroll, J. T.↗

Efficient Implementation of an Optimal Interpolator for Large Spatial Data Sets

Interpolating scattered data points is a problem of wide ranging interest. A number of approaches for interpolation have been proposed both from theoretical domains such as computational geometry and in applications' fields such as geostatistics. Our motivation arises from geological and mining applications. In many instances data can be costly to compute and are available only at nonuniformly scattered positions. Because of the high cost of collecting measurements, high accuracy is required in the interpolants. One of the most popular interpolation methods in this field is called ordinary kriging. It is popular because it is a best linear unbiased estimator. The price for its statistical optimality is that the estimator is computationally very expensive. This is because the value of each interpolant is given by the solution of a large dense linear system. In practice, kriging problems have been solved approximately by restricting the domain to a small local neighborhood of points that lie near the query point. Determining the proper size for this neighborhood is a solved by ad hoc methods, and it has been shown that this approach leads to undesirable discontinuities in the interpolant. Recently a more principled approach to approximating kriging has been proposed based on a technique called covariance tapering. This process achieves its efficiency by replacing the large dense kriging system with a much sparser linear system. This technique has been applied to a restriction of our problem, called simple kriging, which is not unbiased for general data sets. In this paper we generalize these results by showing how to apply covariance tapering to the more general problem of ordinary kriging. Through experimentation we demonstrate the space and time efficiency and accuracy of approximating ordinary kriging through the use of covariance tapering combined with iterative methods for solving large sparse systems. We demonstrate our approach on large data sizes arising both from synthetic sources and from real applications.

Memarsadeghi, Nargess↗

Hot-Fire Testing Results of Additively Manufactured NASA HR-1 towards Large Scale Implementation of Channel Wall Nozzles

Additive manufacturing (AM) has provided new design and manufacturing opportunities to reduce cost and schedules, consolidate parts, and optimize performance. One technique being evaluated is Laser Powder Directed Energy Deposition (LP-DED), which provides a significant increase in scale compared to Laser Powder Bed Fusion (L-PBF). NASA along with industry partners have been developing the LP-DED process to demonstrate internal channel geometry and development components for use in liquid rocket engine channel cooled nozzles. One material that has been the focus of this development is NASA HR-1, an Iron-Nickel based superalloy optimized for use in high pressure hydrogen environments. Material properties and process characterization has been completed on the high density thin-wall material in addition to advancements of the supply chain. NASA has also completed fabrication of several subscale and full scale channel wall nozzles in LP-DED NASA HR-1 and completed hot-fire testing. This includes refinement of the process to produce thin-walls and various channel geometries to meet the requirements for channel wall nozzle applications. This paper will provide an overview of the LP-DED process development, material characterization and properties, component manufacturing, and hot-fire testing. Hot-fire testing was completed for a lander-class 7K-lbf thrust chamber using Liquid Oxygen (LOX)/Methane. The design overview and results from hot-fire testing will be presented in addition to hardware development for future testing on 2K-lbf and 35k-lbf thrust chambers and large-scale manufacturing technology demonstrators.

Paul Gradl↗

Large Array Channel Capacity in the Presence of Interference

We develop a model for a large array ground receiver system for use in deep-space communications, and analyze the resulting array channel capacity. The model includes effects of array geometry, time-dependent spacecraft orbital trajectory, point and extended interference sources, and elevation-dependent noise and tropospheric channel variations. Channel capacity is expressed as the ratio of determinants of covariance matrices characterizing source, interference, and additive noise, and then reduced to a simpler quadratic form more amenable to analysis and numerical computation. This formulation facilitates inclusion of array and channel characteristics into the model, as well as comparison of optimal, suboptimal, and equivalent single antenna configurations on achievable throughput. Realistic examples of ground array channel capacity calculations are presented, demonstrating the impact of array geometry, planetary interference sources, and array combining algorithm design upon the achievable data throughput.

Vilnrotter, V.↗

Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗

Static and Dynamic Structural Response of an Aircraft Wing with Damage Using Equivalent Plate Analysis

A process to generate an equivalent plate based on an optimization approach to predict the static and dynamic response of flight vehicle wing structures is proposed. Geometric-scale and frequency-scale factors are defined to construct an equivalent plate with any desired scale to use in simulation and wind tunnel experiments. It is shown that the stiffness and the displacements are scaled linearly with the geometric-scale factor, whereas the load is scaled as the square of the geometric-scale factor. The scaled stiffness of the reference flight vehicle is matched first to construct the equivalent plate. Then the frequency-scale factor is defined to scale the flight vehicle frequencies. The scaled flight vehicle frequencies are matched by placing arbitrary point masses along the equivalent plate geometry. Two simple stiffened-plate examples, one with damage and another without damage, were used to demonstrate the accuracy of the optimization procedure proposed. Geometric-scale factors ranging from 0.2 to 1.0 were used in the analyses. In both examples, the static and dynamic response of the reference stiffened-panel solution is matched accurately. The scaled equivalent plate predicted the first five frequencies of the stiffened panel very accurately. Finally, the proposed equivalent plate procedure was demonstrated in a more realistic typical aircraft wing structure. Two scale equivalent plate models were generated using the geometric-scale factors 1.0 and 0.2. Both equivalent plate models predicted the static response of the wing structure accurately. The equivalent plate models reproduced the first five frequencies of the wing structure accurately.

Krishnamurthy, T.↗

Calculation of optimum airfoils using direct solutions of the Navier-Stokes equations

A method has been developed which calculates the optimum airfoil shape at a specified operating condition. The Reynolds-averaged Navier-Stokes equations with a turbulence model are used as the governing equations for the fluid flow, and are explicitly introduced into the optimization analysis as constraint equations. The locations of all grid points which define the surface of the airfoil are design variables, eliminating the need for shape functions to describe changes in the airfoil geometry. A generalized reduced-gradient method is used to find the optimum airfoil shape. The use of a direct solution technique (Newton's method) to solve the Navier-Stokes equations allows for the efficient solution of the problem. This paper describes the optimization method and presents sample results.

Felker, Fort F.↗

Optimization of a cyclone using MFIX and Nodeworks

Video depicting the optimization process of a cyclone on NETL's chemical looping reactor (CLR) using MFIX and Nodeworks. MFIX is used to model the cyclone using PIC. Nodeworks is then used to generate proposed geometry changes using a Latin hypercube. Each design is simulated, with an objective value being computed based on the cyclone efficiency and pressure drop. A Gaussian Process surrogate model is then constructed from the objective values. This surrogate model is then used by a differential evolution optimization algorithm to identify the optimal cyclone design. Details published here: Weber, J., Fullmer, W., Gel, A., and Musser, J. (February 4, 2020). "Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics." ASME. J. Fluids Eng. March 2020; 142(3): 031111. https://doi.org/10.1115/1.4045952 OSTI: https://www.osti.gov/pages/servlets/purl/1763893

cyclone↗

Composite panel development at JPL

Parametric computer studies can be use in a cost effective manner to determine optimized composite mirror panel designs. An InterDisciplinary computer Model (IDM) was created to aid in the development of high precision reflector panels for LDR. The materials properties, thermal responses, structural geometries, and radio/optical precision are synergistically analyzed for specific panel designs. Promising panels designs are fabricated and tested so that comparison with panel test results can be used to verify performance prediction models and accommodate design refinement. The iterative approach of computer design and model refinement with performance testing and materials optimization has shown good results for LDR panels.

Mcelroy, Paul↗

Optimization of circular orifice jets mixing into a heated cross flow in a cylindrical duct

To examine the mixing characteristics of circular jets in an axisymmetric can geometry, temperature measurements were obtained downstream of a row of cold jet injected into a heated cross stream. The objective was to obtain uniform mixing within one duct radius downstream of the leading edge of the jet orifices. An area weighted standard deviation of the mixture fraction was used to help quantify the degree of mixedness at a given plane. Non-reacting experiments were conducted to determine the influence of the number of jets on the mixedness in a cylindrical configuration. Results show that the number of orifices significantly impacts the mixing characteristics of jets injected from round hole orifices in a can geometry. Optimum mixing occurs when the mean jet trajectory aligns with the radius which divides the cross sectional area of the can into two equal parts at one mixer radius downstream of the leading edge of the orifice. The optimum number of holes at momentum-flux ratios of 25 and 52 is 10 and 15 respectively.

Kroll, J. T.↗

Efficient Implementation of an Optimal Interpolator for Large Spatial Data Sets

Scattered data interpolation is a problem of interest in numerous areas such as electronic imaging, smooth surface modeling, and computational geometry. Our motivation arises from applications in geology and mining, which often involve large scattered data sets and a demand for high accuracy. The method of choice is ordinary kriging. This is because it is a best unbiased estimator. Unfortunately, this interpolant is computationally very expensive to compute exactly. For n scattered data points, computing the value of a single interpolant involves solving a dense linear system of size roughly n x n. This is infeasible for large n. In practice, kriging is solved approximately by local approaches that are based on considering only a relatively small'number of points that lie close to the query point. There are many problems with this local approach, however. The first is that determining the proper neighborhood size is tricky, and is usually solved by ad hoc methods such as selecting a fixed number of nearest neighbors or all the points lying within a fixed radius. Such fixed neighborhood sizes may not work well for all query points, depending on local density of the point distribution. Local methods also suffer from the problem that the resulting interpolant is not continuous. Meyer showed that while kriging produces smooth continues surfaces, it has zero order continuity along its borders. Thus, at interface boundaries where the neighborhood changes, the interpolant behaves discontinuously. Therefore, it is important to consider and solve the global system for each interpolant. However, solving such large dense systems for each query point is impractical. Recently a more principled approach to approximating kriging has been proposed based on a technique called covariance tapering. The problems arise from the fact that the covariance functions that are used in kriging have global support. Our implementations combine, utilize, and enhance a number of different approaches that have been introduced in literature for solving large linear systems for interpolation of scattered data points. For very large systems, exact methods such as Gaussian elimination are impractical since they require 0(n(exp 3)) time and 0(n(exp 2)) storage. As Billings et al. suggested, we use an iterative approach. In particular, we use the SYMMLQ method, for solving the large but sparse ordinary kriging systems that result from tapering. The main technical issue that need to be overcome in our algorithmic solution is that the points' covariance matrix for kriging should be symmetric positive definite. The goal of tapering is to obtain a sparse approximate representation of the covariance matrix while maintaining its positive definiteness. Furrer et al. used tapering to obtain a sparse linear system of the form Ax = b, where A is the tapered symmetric positive definite covariance matrix. Thus, Cholesky factorization could be used to solve their linear systems. They implemented an efficient sparse Cholesky decomposition method. They also showed if these tapers are used for a limited class of covariance models, the solution of the system converges to the solution of the original system. Matrix A in the ordinary kriging system, while symmetric, is not positive definite. Thus, their approach is not applicable to the ordinary kriging system. Therefore, we use tapering only to obtain a sparse linear system. Then, we use SYMMLQ to solve the ordinary kriging system. We show that solving large kriging systems becomes practical via tapering and iterative methods, and results in lower estimation errors compared to traditional local approaches, and significant memory savings compared to the original global system. We also developed a more efficient variant of the sparse SYMMLQ method for large ordinary kriging systems. This approach adaptively finds the correct local neighborhood for each query point in the interpolation process.

Memarsadeghi, Nargess↗