Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field 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,549 records · Page 86

Lax-Oleinik-Type Formulas and Efficient Algorithms for Certain High-Dimensional Optimal Control Problems

Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimension. In this paper, we provide analytical solutions to certain optimal control problems whose running cost depends on the state variable and with constraints on the control. We also provide Lax-Oleinik-type representation formulas for the corresponding Hamilton-Jacobi partial differential equations with state-dependent Hamiltonians. Additionally, we present an efficient, grid-free numerical solver based on our representation formulas, which is shown to scale linearly with the state dimension, and thus, to overcome the curse of dimensionality. Using existing optimization methods and the min-plus technique, we extend our numerical solvers to address more general classes of convex and nonconvex initial costs. We demonstrate the capabilities of our numerical solvers using implementations on a central processing unit (CPU) and a field-programmable gate array (FPGA). In several cases, our FPGA implementation obtains over a 10 times speedup compared to the CPU, which demonstrates the promising performance boosts FPGAs can achieve. Furthermore, our numerical results show that our solvers have the potential to serve as a building block for solving broader classes of high-dimensional optimal control problems in real-time.

97 MATHEMATICS AND COMPUTING↗

Systematic Operating Temperature Differences Between Al-BSF, PERC, and PERT-With-Optimized-Rear-Reflector Solar Mini-Modules Due to Rear Reflectance

Reflecting sub-bandgap light from photovoltaic modules has the potential to improve lifetime energy generation of fielded systems by reducing operating temperature. In this article, the temperature of fielded aluminum back-surface field (Al-BSF) and passivated emitter and rear contact (PERC) mini-modules was monitored every 5 minutes for 75 days along with corresponding meteorological data. Additionally, passivated emitter rear totally diffused (PERT) mini-modules with high-performance sub-bandgap rear reflectors were tested and compared to the state-of-the-art industrial modules. These reflectors consisted of a >300-nm-thick silicon dioxide nanoparticle film with a low refractive index. The impact of reflectance on measured operating temperature was isolated with a previously developed thermal model and quantified as the reflectance-induced median temperature difference between each tested module at representative outdoor conditions (1000 W m -2 , 25 degrees C ambient temperature, and 1.43 m s -1 wind speed). We found that, because of their reflectance differences, PERC modules ran systematically cooler than Al-BSF modules by 1.0 degrees C, whereas the PERT-with-optimized-rear-reflector systematically operated 1.4 degrees C cooler than the Al-BSF module and 0.4 degrees C than the PERC module. We also found that the rear reflector provided the greatest temperature benefit during periods of highest irradiance.

14 SOLAR ENERGY↗

Dephasing of ion beams as magnetic vortex acceleration regime transitions into a bubble-like field structure

The interaction of an ultra-intense laser pulse with a near critical density target results in the formation of a plasma channel, a strong azimuthal magnetic field and moving vortices. An application of this is the generation of energetic and collimated ion beams via magnetic vortex acceleration. The optimized regime of magnetic vortex acceleration is becoming experimentally accessible with new high intensity laser beamlines coming online and advances made in near critical density target fabrication. The robustness of the acceleration mechanism with realistic experimental conditions is examined with three-dimensional simulations. Of particular interest is the acceleration performance with different laser temporal contrast conditions, in some cases leading to pre-expanded target profiles prior to the arrival of the main pulse. Preplasma effects on the structure of the accelerating fields are explored, including a detailed analysis of the ion beam properties and the efficiency of the process. Improved scaling laws for the magnetic vortex acceleration mechanism, including the laser focal spot size effects, are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tunable Superconducting Gravity Gradiometer for Mars Climate, Atmosphere, and Gravity Field Investigation

We are developing a compact tensor superconducting gravity gradiometer (SGG) for obtaining gravimetric measurements from planetary orbits. A new and innovative design gives a potential sensitivity of approximately 10(sup -4) E Hz(sup - 1/2)( 1 E = 10(sup -9 S(sup -2) in the measurement band up to 0.1 Hz (suitale for short wavelength static gravity) and of approximately 10(sup -4) E Hz(sup - 1/2) in the frequency band less than 1 mHz (for long wavelength time-variable gravity) from the same device with a baseline just over 10 cm. The measurement band and sensitiy can be optimally tuned in-flight during the mission by changing resonance frequencies, which allows meaurements of both static and time-variable gravity fields from the same mission. Significant advances in the technologies needed for space-based cryogenic instruments have been made in the last decade. In particular, the use of cryocoolers will alleviate the previously severe constraint on mission lifetime imposed by the use of liquid helium, enabling mission durations in the 5 - 10 year range.

Mars↗

Fabrication of a Hybrid Transition Edge Sensor Array for Lynx

Lynx is a proposed NASA X-Ray telescope flight mission aimed at achieving state-of-the-art angular and energy resolution with a 100 kilopixel array to probe the hot energetic young universe in unprecedented detail. To achieve these goals, our team plans on leveraging our current work in development of the focal plane for the Athena X-Ray Integral Field Unit (X-IFU) while advancing the state-of-the-art in transition edge sensor (TES) X-ray detector technology. The TES is an optimal technology for achieving both high energy and fine angular resolution at the same time because pixel features can be made extremely small and the absorber which dominates the heat capacity can be tuned to meet resolution requirements. Specifically, the proposed mission concept calls for a hybrid detector of three different arrays fabricated in the same planar process in one focal plane and optimized for different science goals. The main arrays consist of 5x5 hydras, 25 pixels of 4 micron thick Au absorbers each with a different thermal link to one common TES. The outer array has absorbers on a 50-micron pitch for most of the 5 arc-minute field-of-view, and the inner array has 25-micron absorbers for the central 1 arc-minute region. A high resolution array consisting of single pixel 1 micron thick Au absorbers on 50-micron pitch will lie off to the side. Reading out an array of this magnitude will likely require improvements in indium bump bonding to superconducting flexible wiring. Fabrication of absorbers of two different sizes requires electroplating through a photoresist mold by careful tuning of the current density to achieve uniform flat absorbers on a fine pitch scale, followed by ion milling to yield narrow streets separating the pixels while preserving high quantum efficiency. We report on progress made at fabricating the hybrid array with different absorber sizes and thicknesses. Further, we also report on ongoing work to adequately heat sink the pixels with backside wire bonding and copper coating. We also report on work to improve detector pixel yield and top side indium bump bonding to flexible wiring.

Focal plane array↗

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗

EUV debris mitigation using magnetic nulls

Next generation EUV sources for photolithography use light produced by laser-produced plasmas (LPP) from ablated tin droplets. A major challenge for extending the lifetime of these devices is mitigating damage caused by deposition of tin debris on the sensitive collection mirror. Especially difficult to stop are high energy (up to 10 keV) highly charged tin ions created in the plasma. Existing solutions include the use of stopping gas, electric fields, and magnetic fields. One common configuration consists of a magnetic field perpendicular to the EUV emission direction, but such a system can result in ion populations that are trapped rather than removed. Here we investigate a previously unconsidered mitigation geometry consisting of a magnetic null by performing full-orbit integration of the ion trajectories in an EUV system with realistic dimensions and optimize the coil locations for the null configuration. The magnetic null prevents a fraction of ions from hitting the mirror comparable to that of the perpendicular field, but does not trap any ions due to the chaotic nature of ion trajectories that pass close to the null. This technology can potentially improve LPP-based EUV photolithography system efficiency and lifetime and may allow for a different, more efficient formulation of buffer gas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Efficient Annual-Performance Model of a Geothermal Network for Improved System Design, Operation, and Control: Preprint

Geothermal district energy systems (DES), and specifically geothermal networks, provide a viable solution for decarbonizing residential and commercial heating and cooling loads. District energy systems of all kinds enable a thermal resource with a relatively high capital cost (such as a geothermal borehole field) to be shared among a large number of users. While district heating and cooling has been studied for many years, geothermal networks, fifth generation DES that utilize water-source heat pumps and an ambient temperature loop to meet heating and cooling loads, have not been implemented extensively and thus require additional technical and economic optimization to obtain maximum benefits. This paper presents a newly developed reduced-order model that captures the flow of energy within the network, including the commercial and residential users' electrical usage, at an hourly rate over a year. The model includes building loads, heat pumps, borehole fields, and auxiliary heat/cool input, all connected with an ambient-temperature thermal loop model. In the model, operational control is possible for the borehole fields, circulation pump, and auxiliary system. For a given system, the model can output the complete state parameters for each component, the thermal loop, and the collective system, such as flow rate over time, average thermal loop temperature over time, and total electricity usage. The model can also be used to optimize the system control for maximizing system efficiency or minimizing system operational cost. For example, one initial assessment of the borehole controller for an example system showed that a controller with on/off operation of the borehole field reduces annual electrical usage by 33%, compared with continuous operation mode. Hence, the model can assist in optimizing a given system's operation to get the most value out of a geothermal network installation. Future work will consider the model's application to a demonstration project, including the model validation against operational data and system operation optimization.

ambient-temperature loop↗

3D In-situ Stop Action Study of Recrystallization in Additively Manufactured 316L Stainless Steel: Reconstruction Optimization and Observations

A volume of an additively manufactured 316L stainless steel sample has been tracked during its recrystallization using near-and far-field High Energy Diffraction Microscopy (HEDM) and absorption tomography at Advanced Photon Source beamline 1-ID. A nearfield compatible in situ furnace allows monitoring of Bragg diffraction signals as they evolve out of a weak and diffuse background while the sample temperature is approximate to 1250 degrees C. The sample is rapidly cooled to room temperature after observation of significant signal evolution and ti similar to 0.035 mm(3) is mapped by the near-field method. Four cycles of heat treatment follow the structure from a state of small, isolated grains through impingement of domains to near completion of recyrstallization. Here, the experiment and reconstructions are described, and recrystallized fractions, twin domains, and distributions of grain boundary types are discussed.

Liu, He↗

Development of Analysis Methods that Integrate Numeric and Textual Equipment Reliability Data

Within the Light Water Reactor Sustainability (LWRS) program, the Risk-Informed Systems Analysis (RISA) Pathway is performing collaborative research on the development and deployment of technologies designed to assist operating nuclear power plants (NPPs) to reduce operating costs improve plant reliability and availability. One of the RISA research areas is focusing on the development of methods and tools designed to optimize plant operations (e.g., maintenance/replacement schedules, optimal maintenance postures for plant structures, systems, and components [SSCs]) in a manner that is more cost effective than current approaches and makes better use of available SSC health data. The Risk-Informed Asset Management (RIAM) project targets this research area by creating a direct bridge between component equipment reliability (ER) data and system engineer decision making regarding maintenance activity scheduling and component aging management. In this respect, one challenge that NPP system engineers are facing is that the amount of ER data being continuously generated is not only extremely large in size, but it comes in different forms: textual (e.g., condition or maintenance reports) and numeric (e.g., generated by monitoring systems). All these data elements provide them with valuable insights and information regarding: 1) the discovery of anomalous behaviors or degradation trends, 2) the identification of the possible causes behind such behaviors/trends, and 3) the prediction of their direct consequences. However, several challenges have proved to be roadblocks to this process. While some of these challenges are technical in nature (i.e., data are often distributed over several physical servers/databases), others are conceptual in nature: data elements come in different formats (e.g., numeric or textual), and measured values have different scales (e.g., vibration spectra and oil temperature). The activities performed by the RIAM project during FY23 directly tackles the need to simultaneously integrate the analysis of ER data in all its forms, numeric and textual. Note that such task has never been performed before due to the complexity of the systems under consideration but, most importantly, because of the technical challenges behind the harmonization of ER data formats and the lack of adequate computational methods to analyze them. Our approach borrows ideas and concepts from the medical field where integration of several data sources is vital to assist medical practitioners to perform correct diagnosis and indicate optimal treatments. In our view a NPP asset is equivalent to a patient in a medical context. The main difference is the complexity of a human body is a magnitude more complex when compared to typical assets commonly present in NPPs (e.g., centrifugal pumps, or motor operated valves). This simplifies our first requirement when analyzing heterogenous ER data formats: to put data into “context”. Context is here intended as the additional piece of information that is needed by ER data analysis tools to understand what these data elements are referring to, i.e., which king of knowledge they are generating. In our context, this knowledge can be translated into models that capture the form and functional architecture of assets/systems, their dependencies, and how they interact. These models actually emulate the knowledge that that NPP system engineers possess about assets and systems; this is their key of success when analyzing ER data, their challenge is ability to handle large amount of data. Here, we employ model-based system engineering (MBSE) models of systems and assets to represent and capture their architecture and functional, i.e. cause-effect, relations. Then, ER data elements are processed by identifying first of all which elements of the developed MBSE elements they are referring to. For numeric ER data this task is fairly easy since it is possible to precisely pinpoint what MBSE elements the corresponding sensor are observing (e.g., bearing temperature of a centrifugal pump). Task is much harder for textual data since the information contained in issue or maintenance reports needs to “be understood” by a computational tool. Here we called this process as “knowledge extraction”. Once again, we borrow the experience in the medical field where methods to extract knowledge from textual data have been developed in the past decade. The missing element for us is the availability of a complete dictionary of NPP related concepts (in addition to the MBSE models presented earlier) that can put “text into context”. In FY23, such dictionary has been developed along with all the computational elements required for knowledge extraction. Lastly, once numeric and textual ER data elements have been processed and “understood”, then the last step is the discovery of possible cause-effect relations among them. This is performed by observing if a logical connection through the MBSE models exists, and if the

97 MATHEMATICS AND COMPUTING↗

Error Estimation of An Ensemble Statistical Seasonal Precipitation Prediction Model

This NASA Technical Memorandum describes an optimal ensemble canonical correlation forecasting model for seasonal precipitation. Each individual forecast is based on the canonical correlation analysis (CCA) in the spectral spaces whose bases are empirical orthogonal functions (EOF). The optimal weights in the ensemble forecasting crucially depend on the mean square error of each individual forecast. An estimate of the mean square error of a CCA prediction is made also using the spectral method. The error is decomposed onto EOFs of the predictand and decreases linearly according to the correlation between the predictor and predictand. Since new CCA scheme is derived for continuous fields of predictor and predictand, an area-factor is automatically included. Thus our model is an improvement of the spectral CCA scheme of Barnett and Preisendorfer. The improvements include (1) the use of area-factor, (2) the estimation of prediction error, and (3) the optimal ensemble of multiple forecasts. The new CCA model is applied to the seasonal forecasting of the United States (US) precipitation field. The predictor is the sea surface temperature (SST). The US Climate Prediction Center's reconstructed SST is used as the predictor's historical data. The US National Center for Environmental Prediction's optimally interpolated precipitation (1951-2000) is used as the predictand's historical data. Our forecast experiments show that the new ensemble canonical correlation scheme renders a reasonable forecasting skill. For example, when using September-October-November SST to predict the next season December-January-February precipitation, the spatial pattern correlation between the observed and predicted are positive in 46 years among the 50 years of experiments. The positive correlations are close to or greater than 0.4 in 29 years, which indicates excellent performance of the forecasting model. The forecasting skill can be further enhanced when several predictors are used.

Shen, Samuel S. P.↗

Study of Cathode Gas Diffusion Architecture for Improved Oxygen Transport in Hydroxide Exchange Membrane Fuel Cells

The high pH environment in hydroxide exchange membrane fuel cells (HEMFCs) has the potential to reach lower costs than the current proton exchange membrane fuel cells (PEMFCs), the incumbent technology. A significant difference between HEMFCs and PEMFCs is the location of water production within the cell. In PEMFCs, the water is produced on the cathode, limiting oxygen transport. In HEMFCs, the water is produced on the anode where the fuel is pure hydrogen. This allows the cathode to be optimized for oxygen transport without the presence of excess liquid water. Limiting current analysis, a technique previously used in PEMFCs, is adopted in HEMFCs to evaluate the oxygen mass transport resistances for different sections of the cathode. Through elimination of the microporous layer (MPL), gas diffusion layer (GDL), and traditional flow field and using porous nickel foam for gas distribution, the transport resistance at an operating condition of 150 kPa(g) and with the cell temperature at 80 °C was decreased from 112 s m −1 to 48 s m −1 , effectively halved. The optimal configuration for performance was found with Ni foam and a GDL, eliminating the MPL and traditional flow field, which vastly improved oxygen transport while maintaining adequate electrical contact with the cathode catalyst layer.

25 ENERGY STORAGE↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

HPFTP flow diverter analysis

A computational fluid flow analysis on the flow diverter system under consideration for the Space Shuttle main engine high pressure fuel turbopump (SSME HPFTP) is proposed. A three dimensional viscous flow environment is computed to optimize the geometric configuration and location of the flow diverter system. The analysis consists of a fully turbulent cold flow calculation by Navier-Stokes equations and a Baldwin-Lomax turbulence model. The equations are numerically by a finite difference/element procedure. The results will provide the steady and unsteady pressure field and thermal environment required to assess the usefulness of the flow diverter system in deflecting the cold flow away from the hot turbine components. A geometry optimization study determines the best diverter shape and location to avoid larger thermal gradients on the rotor/stator components.

Spadley, L. W.↗

An investigation of design optimization using a 2-D viscous flow code with multigrid

Computational fluid dynamics (CFD) codes have advanced to the point where they are effective analytical tools for solving flow fields around complex geometries. There is also a need for their use as a design tool to find optimum aerodynamic shapes. In the area of design, however, a difficulty arises due to the large amount of computer resources required by these codes. It is desired to streamline the design process so that a large number of design options and constraints can be investigated without overloading the system. There are several techniques which have been proposed to help streamline the design process. The feasibility of one of these techniques is investigated. The technique under consideration is the interaction of the geometry change with the flow calculation. The problem of finding the value of camber which maximizes the ratio of lift over drag for a particular airfoil is considered. In order to test out this technique, a particular optimization problem was tried. A NACA 0012 airfoil was considered at free stream Mach number of 0.5 with a zero angle of attack. Camber was added to the mean line of the airfoil. The goal was to find the value of camber for which the ratio of lift over drag is a maximum. The flow code used was FLOMGE which is a two dimensional viscous flow solver which uses multigrid to speed up convergence. A hyperbolic grid generation program was used to construct the grid for each value of camber.

Doria, Michael L.↗

Mass drivers. 3: Engineering

This paper is the last of a series of three by the Mass-Driver Group of the 1977 Ames Summer Study. It develops the engineering principles required to implement the basic mass-driver design defined in paper I and as presented previously in the literature. Optimum component mass trade-offs are derived from a set of four input parameters, and the program used to design a lunar launcher. The mass optimization procedure is then incorporated into a more comprehensive mission optimization program called OPT-4, which evaluates an optimized mass-driver reaction engine and its performance in a range of specified missions. Finally, this paper discusses, to the extent that time permitted, certain peripheral problems: heating effects in buckets due to magnetic field ripple; an approximate derivation of guide force profiles; the mechanics of inserting and releasing payloads; the reaction mass orbits; and a proposed research and development plan for implementing mass drivers. This paper also considers further optimization that depends on advances in component technology beyond the present state. It is concluded that significant advantage can be gained in the case of small reaction engines with high emission velocity by developing lightweight packaging technology for silicon-controlled rectifiers (SCR's).

William H. Arnold↗