Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape approximation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Implementation of Combined Feather and Surface-Normal Ice Growth Models in LEWICE/X

Experimental observations have shown that discrete rime ice growths called feathers, which grow in approximately the direction of water droplet impingement, play an important role in the growth of ice on accreting surfaces for some thermodynamic conditions. An improved physical model of ice accretion has been implemented in the LEWICE 2D panel-based ice accretion code maintained by the NASA Lewis Research Center. The LEWICE/X model of ice accretion explicitly simulates regions of feather growth within the framework of the LEWICE model. Water droplets impinging on an accreting surface are withheld from the normal LEWICE mass/energy balance and handled in a separate routine; ice growth resulting from these droplets is performed with enhanced convective heat transfer approximately along droplet impingement directions. An independent underlying ice shape is grown along surface normals using the unmodified LEWICE method. The resulting dual-surface ice shape models roughness-induced feather growth observed in icing wind tunnel tests. Experiments indicate that the exact direction of feather growth is dependent on external conditions. Data is presented to support a linear variation of growth direction with temperature and cloud water content. Test runs of LEWICE/X indicate that the sizes of surface regions containing feathers are influenced by initial roughness element height. This suggests that a previous argument that feather region size is determined by boundary layer transition may be incorrect. Simulation results for two typical test cases give improved shape agreement over unmodified LEWICE.

Matthew T Velazquez↗
NASA NTRS1999

Aeroelastic Analysis of a Trimmed Generic Hypersonic Vehicle

The aeroelastic equations of motion governing a hypersonic vehicle in free flight are derived. The equations of motion for a translating and rotating flexible body using Lagrange's equations in terms of quasi-coordinates are presented. These equations are simplified for the case of a vehicle with pitch and plunge rigid body degrees of freedom and small elastic displacements. The displacements are approximated by a truncated series of the unrestrained mode shapes, which are obtained using equivalent plate theory. Subsequently, the nonlinear equations of motion are linearized about the trim state, which is obtained using a rigid body trim model and steady hypersonic aerodynamics. The appropriate flutter derivatives are calculated from piston theory. Results describing mode shapes, trim behavior, and aeroelastic stability of a generic hypersonic vehicle are presented.

Nydick, I.↗
NASA NTRS2000

Innovative Research Program: Supershields for Gamma-Ray Astronomy

The supershield project evaluated the importance of novel shield configurations for suppressing neutron induced background in new classes of gamma-ray detectors such as CZT. The basic concept was to use a two-part shield. The outer shield material heavily moderates the incoming neutron spectrum. This moderated neutron beam is then more easily absorbed by the inner material, which is an efficient neutron absorber. This approach is, in principle, more efficient than that in previous attempts to make neutron shields. These previous attempts involved biatomic, monlithic shields (eg. LiH) in which the shield consisted of a single material but with two types of atoms - one for moderating and one for absorbing. The problem with this type of monolithic shield is that moderating neutrons, without the efficient absorption of them, leads to the leakage into the detector of neutrons with a low energy component (approx. 10-100 KeV). These energy neutrons are particularly problematic for many types of detectors. The project was roughly divided into phases. In the first phase we attempted to carefully define the neutron source function incident on any space instrument. This is essential since the design of any shield depends on the shape of the incident neutron spectrum. We found that approximations commonly used in gamma-ray astronomy for photon background is inadequate. In addition, we found that secondary neutrons produced in any passive shield, and dominated by inelastic neutron scattering, are far more important than background due to neutron activation. The second phase of our work involved design of supershield geometries (one and three dimensional) in order to compare different shield configurations and materials for their effectiveness as neutron shields. Moreover we wanted to compare these supershields with previous neutron shields to confirm the performance differences between the supershield (two material) and monolithic (one material) designs and to understand the physics origins of these differences more clearly. The third phase of the supershield program involved the benchmarking of the supershield designs through direct experimental verification. This required fabricating various supershields and exposing them to beams of neutrons to directly characterize their performance. With explicit verification that our modeling procedures can be used with confidence, we are now in a position to design shields for realistic space geometries. Using the supershield modeling capacity developed as part of this program we are attempting to evaluate their utility for a specific proposed mission--the Energetic X-ray Imaging Survey Telescope (EXIST). It is anticipated that this experiment, which is limited by internal background at high energies, might benefit from a neutron shield.

Hailey, Charles J.↗
NASA NTRS2002

Scale Closure in Upper Ocean Optical Properties: From Single Particles to Ocean Color

Predictions of chlorophyll concentration from satellite ocean color are an indicator of primary productivity, with implications for foodwebs, fisheries, and the global carbon cycle. Models describing the relationship between optical properties and chlorophyll do not account for much of the optical variability observed in natural waters, because of the presence of seawater constituents that do not covary with phytoplankton pigments. in order to understand variability in these models, the optical contributions of seawater constituents were investigated. A combination of Mie theory and flow cytometry was used to determine the diameter, complex refractive index, and optical cross-sections of individual particles. In New England continental shelf waters, eukaryotic phytoplankton were the main particle contributors to absorption and scaftering. Minerals were the main contributor to backscattering (bb) in the spring, whereas in the summer both minerals and detritus contributed to bb. Synechococcus and heterotrophic bacteria were relatively unimportant optically. Seasonal differences in the spectral shape of remote sensing reflectance, Rrs, were contributed to approximately equally by eukaryotic phytoplankton absorption, dissolved absorption, and non-phytoplankton bb. Differences between measurements of bb and Prs and modeled values based on chlorophyll concentration were caused by higher dissolved absorption and non-phytoplankton bb than were assumed by the model.

Green, Rebecca E.↗
NASA NTRS1991

Lake Hoare, Antarctica: sedimentation through a thick perennial ice cover

Lake Hoare in the Dry Valleys of Antarctica is covered with a perennial ice cover more than 3 m thick, yet there is a complex record of sedimentation and of growth of microbial mats on the lake bottom. Rough topography on the ice covering the lake surface traps sand that is transported by the wind. In late summer, vertical conduits form by melting and fracturing, making the ice permeable to both liquid water and gases. Cross-sections of the ice cover show that sand is able to penetrate into and apparently through it by descending through these conduits. This is the primary sedimentation mechanism in the lake. Sediment traps retrieved from the lake bottom indicate that rates of deposition can vary by large amounts over lateral scales as small as 1 m. This conclusion is supported by cores taken in a 3 x 3 grid with a spacing of 1.5 m. Despite the close spacing of the cores, the poor stratigraphic correlation that is observed indicates substantial lateral variability in sedimentation rate. Apparently, sand descends into the lake from discrete, highly localized sources in the ice that may in some cases deposit a large amount of sand into the lake in a very short time. In some locations on the lake bottom, distinctive sand mounds have been formed by this process. They are primary sedimentary structures and appear unique to the perennially ice-covered lacustrine environment. In some locations they are tens of centimetres high and gently rounded with stable slopes; in others they reach approximately 1 m in height and have a conical shape with slopes at angle of repose. A simple formation model suggests that these differences can be explained by local variations in water depth and sedimentation rate. Rapid colonization of fresh sand surfaces by microbial mats composed of cyanobacteria, eukaryotic algae, and heterotrophic bacteria produces a complex intercalation of organic and sandy layers that are a distinctive form of modern stromatolites.

NASA Discipline Exobiology↗
NASA NTRS2003

Probe Without Moving Parts Measures Flow Angle

The measurement of local flow angle is critical in many fluid-dynamic applications, including the aerodynamic flight testing of new aircraft and flight systems. Flight researchers at NASA Dryden Flight Research Center have recently developed, flight-tested, and patented the force-based flow-angle probe (FLAP), a novel, force-based instrument for the measurement of local flow direction. Containing no moving parts, the FLAP may provide greater simplicity, improved accuracy, and increased measurement access, relative to conventional moving vane-type flow-angle probes. Forces in the FLAP can be measured by various techniques, including those that involve conventional strain gauges (based on electrical resistance) and those that involve more advanced strain gauges (based on optical fibers). A correlation is used to convert force-measurement data to the local flow angle. The use of fiber optics will enable the construction of a miniature FLAP, leading to the possibility of flow measurement in very small or confined regions. This may also enable the tufting of a surface with miniature FLAPs, capable of quantitative flow-angle measurements, similar to attaching yarn tufts for qualitative measurements. The prototype FLAP was a small, aerodynamically shaped, low-aspect-ratio fin about 2 in. (approximately equal to 5 cm) long, 1 in. (approximately equal to 2.5 cm) wide, and 0.125 in. (approximately equal to 0.3 cm) thick (see Figure 1). The prototype FLAP included simple electrical-resistance strain gauges for measuring forces. Four strain gauges were mounted on the FLAP; two on the upper surface and two on the lower surface. The gauges were connected to form a full Wheatstone bridge, configured as a bending bridge. In preparation for a flight test, the prototype FLAP was mounted on the airdata boom of a flight-test fixture (FTF) on the NASA Dryden F-15B flight research airplane.

Corda, Stephen↗
NASA NTRS2017

Design and Fabrication of a Composite Morphing Radiator Panel Using High Conductivity Fibers

Upcoming crewed space missions will involve large internal and external heat loads and require advanced thermal control systems to maintain a desired internal environment temperature. Radiators with at least 12:1 turndown ratios (the ratio between the maximum and minimum heat rejection rates) will be needed. However, current technologies are only able to achieve turndown ratios of approximately 3:1. A morphing radiator capable of altering shape could significantly increase turndown capabilities. Shape memory alloys offer qualities that may be well suited for this endeavor; their temperature-dependent phase changes could offer radiators the ability to passively control heat rejection. In 2015, a morphing radiator prototype was constructed and tested in a thermal vacuum environment, where it successfully demonstrated the morphing behavior and variable heat rejection. Newer composite prototypes have since been designed and manufactured using two distinct types of SMA materials. These models underwent temperature cycling tests in a thermal vacuum chamber and a series of fatigue tests to characterize the lifespan of these designs. The focus of this paper is to present the design approach and testing of the morphing composite facesheet. The discussion includes: an overall description of the project background, definition of performance requirements, composite materials selection, use of analytic and numerical design tools, facesheet fabrication, and finally fatigue testing with accompanying results.

Wescott, Matthew T.↗
NASA NTRSDate not supplied

Fixed- and Free-Mode Flutter Derivatives during Aeroelastic Optimization

Aeroelastic flutter design derivatives may be approximated by ignoring the sensitivity of the structural mode shapes. This “fixed-mode” derivative is less-expensive to compute than the exact “free-mode” derivative (which accommodates the mode shape dependency), but also may provide inaccurate sensitivities for optimization. This work formulates both fixed- and free-mode flutter derivatives, and demonstrates the conditions under which they are equivalent. Two different types of fixed-mode derivatives are also developed in this work. Using a cantilevered-plate demonstration problem with both shape and sizing design parameters, the accuracy of the free-mode derivative and the two fixed-mode derivatives, are all demonstrated via gradient-based optimization.

Bret K. Stanford↗
NASA NTRS2016

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗
DOE OSTI2022

Fracture-based shape optimization built upon the topological derivative

In Silva et al. (2011) and Alidoost et al. (2020), the authors developed an approximation of the energy release rate field associated with a small edge or surface crack at any boundary location and with any orientation using the topological derivative. The approximation is computationally attractive because it requires only a single analysis on the non-cracked domain in contrast with conventional boundary-element and finite-element-based methods, which require a separate and costlier analysis for each crack length-location-orientation combination. Here, a shape optimization scheme for fracture-resistant structures is developed using the energy release rate approximation. In the gradient-based optimization scheme, the domain and its boundary are defined implicitly using level-set functions. The level-set functions of arbitrary geometries are constructed using Boolean operations from the level-set functions of simple primitives. This geometrical representation has the dual advantage of (i) allowing shapes to intersect and/or separate during the optimization and (ii) simplifying the computation of the shape sensitivities.

42 ENGINEERING↗
NASA NTRS1954

Summary of Free-Flight Zero-Lift Drag Results from Tests of 1/5-Scale Models of the Convair YF-102 and F-102A Airplanes and Several Related Small Equivalent Bodies at Mach Numbers from 0.70 to 1.46

One-fifth-scale rocket-propelled models of the Convair YF-102 and F-102A airplanes were tested to determine free-flight zero-lift drag coefficients through the transonic speed range at Reynolds numbers near those to be encountered by the full-scale airplane. Trim and duct characteristics were obtained along with measurements of total-, internal-, and base-drag coefficients. Additional zero-lift drag tests involved a series of small equivalent-body-of-revolution models which were launched to low supersonic speeds by means of a helium gun. The several small models tested corresponded to the following full-scale airplanes: basic, YF-102, 2-foot (full-scale) fuselage extension, F-102A, F-102A (relocated inlets), F-102A (faired nose), and F-102A (parabolic nose) . Equivalent-body models corresponding to the normal area distribution (derived for Mach number 1.0) of each of these airplane shapes were flown and, in addition, equivalent-body models designed to represent the YF-102 and F-102A airplanes at Mach number 1.2 were tested. External-drag coefficients obtained from the 115-scale tests ranged from 0.0094 to 0.0273 for the YF-102 model and from 0.0100 to 0.0255 for the F-102A model. Forebody external-pressure-drag coefficients (drag rise) at Mach number 1.05 of 0.0183 and 0.0134 were obtained from the 115-scale models of the YF-102 and F-102A, respectively, a 16-percent reduction for the F-102A model. Values of drag rise at Mach number 1.05 from the small equivalent-body tests were nearly the same for the basic, YF-102, and 2-foot-fuselage-extension airplane shapes. Equivalent-body tests of the YF-102 and F-102A shapes showed the latter to have about 25 percent less drag rise as compared with a 16-percent reduction illustrated by the 1/5-scale tests. Additional equivalent-body tests illustrating effects of modifications to the F-102A airplane shape shared that relocating the inlets on the fuselage or altering the nose shape to provide a smoother cross-sectional area progression reduced the drag rise by approximately 16 percent. Replacing a major portion of the nose of the F-102A equivalent-body model with one of parabolic shape resulted in about a 21-percent reduction in drag rise. The drag-rise data from the equivalent-body tests include base drag.

Wallskog, Harvey A.↗
NASA NTRS1982

The radar cross section of dielectric disks

A solution is presented for the backscatter (nonstatic) radar cross section of dielectric disks of arbitrary shape, thickness and dielectric constant. The result is obtained by employing a Kirchhoff type approximation to obtain the fields inside the disk. The internal fields induce polarization and conduction currents from which the scattered fields and the radar cross section can be computed. The solution for the radar cross section obtained in this manner is shown to agree with known results in the special cases of normal incidence, thin disks and perfect conductivity. The solution can also be written as a product of the reflection coefficient of an identically oriented slab times the physical optics solution for the backscatter cross section of a perfectly conducting disk of the same shape. This result follows directly from the Kirchhoff type approximation without additional assumptions.

D M Le Vine↗
NASA NTRS1994

The Local Supercluster as a test of cosmological models

The Local Supercluster kinematic properties (the Local Group infall toward the Virgo Cluster and galaxy density distribution about the Virgo Cluster) in various cosmological models are examined utilizing large-scale N-body (PM) simulations 500(exp 3) cells, 250(exp 3) particles, and box size of 400 h(exp -1) Mpc) and are compared to observations. Five models are investigated: (1) the standard, Cosmic Background Explorer Satellite (COBE)-normalized cold dark matter (CDM) model with omega = 1, h = 0.5, and sigma(sub 8) = 1.05; (2) the standard Hot Dark Matter (HDM) model with omega = 1, h = 0.75, and sigma(sub 8) = 1; (3) the tilted CDM model with omega = 1, h = 0.5, n = 0.7, and sigma(sub 8) = 0.5; (4) a CDM + lambda model with omega = 0.3, lambda = 0.7, h = 2/3, and sigma(sub 8) = 2/3; (5) the PBI model with omega = 0.2, h = 0.8, x = 0.1, m = -0.5, and sigma(sub 8) = 0.9. Comparison of the five models with the presently available observational measurements v(sub LG) = 85 - 305 km/s (with mean of 250 km/s), delta(n(sub g))/(n(sub g)-bar) = 1.40 + or - 0.35) suggests that an open universe with omega approximately 0.5 (with or without lambda) and sigma(sub 8) approximately 0.8 is preferred, with omega = 0.3-1.0 (with or without lambda) and sigma(sub 8) = 0.7-1.0 being the acceptable range. At variance with some previous claims based on either direct N-body or spherical nonlinear approaches, we find that a flat model with sigma(sub 8) approximately 0.7-1.0 seems to be reasonably consistent with observations. However, if one favors the low limit of v(sub LG) = 85 km/s, then an omega approximately 0.2-0.3 universe seems to provide a better fit, and flat (omega = 1) models are ruled out at approximately 95% confidence level. On the other hand, if the high limit of v(sub LG) = 350 km/s is closer to the truth, then it appears that omega approximately 0.7-0.8 is more consistent. This test is insensitive to the shape of the power spectrum, but rather sensitive to the normalization of the perturbation amplitude on the relevant scale (e.g., sigma(sub 8)) and omega. We find that neither linear nor nonlinear relations (with spherical symmetry) are good approximations for the relation between radial peculiar velocity and density perturbation, i.e., nonspherical effects and gravitational tidal field are important. The derived omega using either of the two relations is underestimated. In some cases, this error is as large as a factor of 2-4.

Cen, Renyue↗
DOE OSTI2020

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗
NASA NTRS2022

Controls on Greenland Moulin Geometry and Evolution From the Moulin Shape Model

Nearly all meltwater from glaciers and ice sheets is routed englacially through moulins, which collectively comprise approximately 10–14% of the efficient englacial–subglacial hydrologic system. Therefore, the geometry and evolution of moulins has the potential to influence subglacial water pressure variations, ice motion, and the runoff hydrograph delivered to the ocean. We develop the Moulin Shape (MouSh) model, a time-evolving model of moulin geometry. MouSh models ice deformation around a moulin using both viscous and elastic rheologies and models melting within the moulin through heat dissipation from turbulent water flow, both above and below the water line. We force MouSh with idealized and realistic surface melt inputs. Our results show that variations in surface melt change the geometry of a moulin by approximately 30% daily and by over 100% seasonally. These size variations cause observable differences in moulin water storage capacity, moulin water levels, and subglacial channel size compared to a static, cylindrical moulin. Our results suggest that moulins are significant storage reservoirs for meltwater, with storage capacity and water levels varying over multiple timescales. Representing moulin geometry within subglacial hydrologic models would therefore improve their accuracy, especially over seasonal periods or in regions where overburden pressures are high.

Moulin↗
DOE OSTI2022

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗
NASA NTRS1984

Measurement of local connective heat transfer coefficients of four ice accretion shapes

In the analytical study of ice accretions that form on aerodynamic surfaces (airfoils, engine inlets, etc.) it is often necessary to be able to calculate convective heat transfer rates. In order to do this, local convective heat transfer coefficients for the ice accretion shapes must be known. In the past, coefficients obtained for circular cylinders were used as an approximation to the actual coefficients since no better information existed. The purpose of this experimental study was to provide local convective heat transfer coefficients for four shapes that represent ice accretions. The shapes were tested with smooth and rough surfaces. The experimental method chosen was the thin-skin heat rate technique. Using this method local Nusselt numbers were determined for the ice shapes. In general it was found that the convective heat transfer was higher in regions where the model's surfaces were convex and lower in regions where the model's surfaces were concave. The effect of roughness was to increase the heat transfer in the high heat transfer regions by approximately 100% while little change was apparent in the low heat transfer regions.

Smith, M. E.↗