Results for “shape optimization” — Page 68

NASA NTRS2002

Three-Dimensional High-Order Spectral Finite Volume Method for Unstructured Grids

Many areas require a very high-order accurate numerical solution of conservation laws for complex shapes. This paper deals with the extension to three dimensions of the Spectral Finite Volume (SV) method for unstructured grids, which was developed to solve such problems. We first summarize the limitations of traditional methods such as finite-difference, and finite-volume for both structured and unstructured grids. We then describe the basic formulation of the spectral finite volume method. What distinguishes the SV method from conventional high-order finite-volume methods for unstructured triangular or tetrahedral grids is the data reconstruction. Instead of using a large stencil of neighboring cells to perform a high-order reconstruction, the stencil is constructed by partitioning each grid cell, called a spectral volume (SV), into 'structured' sub-cells, called control volumes (CVs). One can show that if all the SV cells are partitioned into polygonal or polyhedral CV sub-cells in a geometrically similar manner, the reconstructions for all the SVs become universal, irrespective of their shapes, sizes, orientations, or locations. It follows that the reconstruction is reduced to a weighted sum of unknowns involving just a few simple adds and multiplies, and those weights are universal and can be pre-determined once for all. The method is thus very efficient, accurate, and yet geometrically flexible. The most critical part of the SV method is the partitioning of the SV into CVs. In this paper we present the partitioning of a tetrahedral SV into polyhedral CVs with one free parameter for polynomial reconstructions up to degree of precision five. (Note that the order of accuracy of the method is one order higher than the reconstruction degree of precision.) The free parameter will be determined by minimizing the Lebesgue constant of the reconstruction matrix or similar criteria to obtain optimized partitions. The details of an efficient, parallelizable code to solve three-dimensional problems for any order of accuracy are then presented. Important aspects of the data structure are discussed. Comparisons with the Discontinuous Galerkin (DG) method are made. Numerical examples for wave propagation problems are presented.

Liu, Yen↗

NASA NTRS1990

Optical neural network system for pose determination of spinning satellites

An optical neural network architecture and algorithm based on a Hopfield optimization network are presented for multitarget tracking. This tracker utilizes a neuron for every possible target track, and a quadratic energy function of neural activities which is minimized using gradient descent neural evolution. The neural net tracker is demonstrated as part of a system for determining position and orientation (pose) of spinning satellites with respect to a robotic spacecraft. The input to the system is time sequence video from a single camera. Novelty detection and filtering are utilized to locate and segment novel regions from the input images. The neural net multitarget tracker determines the correspondences (or tracks) of the novel regions as a function of time, and hence the paths of object (satellite) parts. The path traced out by a given part or region is approximately elliptical in image space, and the position, shape and orientation of the ellipse are functions of the satellite geometry and its pose. Having a geometric model of the satellite, and the elliptical path of a part in image space, the three-dimensional pose of the satellite is determined. Digital simulation results using this algorithm are presented for various satellite poses and lighting conditions.

Lee, Andrew↗

NASA NTRS1998

Acoustics and Aeroperformance of Nozzles With Screwdriver Shaped and Axisymmetric Plugs

The recent experimental and numerical tests of corrugated nozzles have shown some acoustic and thrust benefits relative to traditional round nozzles. For example, a Bluebell nozzle which was obtained by 3D nozzle design incorporating a corrugated cross section nozzle shape with a sinusoidal lip line nozzle edge, can provide an acoustic benefit up to 4dB with about a 1% thrust augmentation. In references, this effect was explained as being the result of the corrugated design producing more efficient mixing of the exhausted jet with ambient air. Based on this argument, the authors have proposed the application of this concept for a centerbody (plug) which can form several vortices downstream from the centerbody. Several different corrugated designs are proposed and described in detail in this paper. The main design is a Screwdriver shaped centerbody or plug (SCR) which was tested experimentally and numerically. The acoustic tests were conducted in the anechoic chamber of the Central AeroHydrodynamics Institute (TsAGI, Moscow) under Civilian Research and Development Foundation (CRDF) grant. These experiments have shown an essential acoustic benefit of about 10-13% with the application of the co-annular nozzles by comparison with the reference round nozzle with the same mass flow rate. However, the expected acoustic benefits with the application of the 4-petal Screwdriver shaped centerbody were not obtained by comparison with the reference axisymmetric centerbody (CON) having the same length and the same cross section areas at the same distance from the nozzle throat. Moreover, for some angles (Theta = 60 deg and 90 deg) noise increase was observed (about 1-3%). These tests will be continued with the goal of obtaining better acoustic results. In particular, acoustic characteristics are hoped to be improved by moving t lie centerbody into the nozzle and using penetrable walls for the SCR and/or for the main nozzle. Preliminary results for such approach are very promising. Aeroperformance effects were analyzed numerically. The nozzle thrust, calculations were based on a full Navier--Stokes equations solver (NSE), and both full and marching Euler codes: CFL3D, CRAFT, and Krayko-Godunov . Grid preparation and its optimization were conducted using GRIDGEN and our own codes. The general conclusion of this numerical analysis is some thrust loss with the application of SCR design (about 1 - 1.5%). But, again, some constructive features of SCR design give some promising perspectives for its application in aviation and domestic industries.

Gilinsky, M.↗

NASA NTRS2024

Mineral Dust Optical Properties for Remote Sensing and Global Modeling: A Review

Dust plays a key role in many Earth system processes and is ubiquitous in the Martian atmosphere. Various intensive field campaigns, laboratory analyses, space-based remote sensing missions, and global modeling efforts aim to characterize dust optical properties. This is a bountiful time for dust scientists, and yet the interpretation of retrievals and comparison to models remains complicated by various conflicting assumptions that are part of each algorithm. For example, the conversion of satellite radiance measurements into products like aerosol optical depth for model evaluation depends upon aerosol properties like particle size and shape that are often prescribed and not part of the retrieval. Conversely, the model calculation of aerosol optical depth often uses different assumptions. The goal of this review is to first document algorithmic assumptions by various satellite retrieval products and models, and identify where there is consistency and where there are differences. In general, the differences documented in this paper reflect uncertainties resulting from incomplete observational characterization of dust aerosols and limitations in our understanding. Second, we note what observations might reduce uncertainties in our knowledge and bring greater consistency to retrievals and models, allowing for a more rigorous and harmonious comparison. The lack of comprehensive and realistic shape models for dust is an outstanding issue, such that closure between forward modeling from particle refractive index, shape, and size and observed optical properties cannot be achieved. Limitations in the computational methods that must be applied to model scattering from complex shapes also makes accurate optical modeling for dust challenging. Field observations indicate the persistence of coarse and giant dust particles at higher altitudes and farther downwind from their source than previously expected. Remote sensing retrieval algorithms based on observations at visible wavelengths have limited sensitivity to these particles and generally do not consider them, although a recent product based on longwave radiances is encouraging. Current measurements of the refractive index of bulk dust and fundamental dust minerology components such as hematite vary widely, inhibiting attempts to represent the variability in dust optical properties and forcing, as expected from different major dust source regions on Earth that have varying mineralogical composition. Some remote sensing retrieval algorithms allow for limited refractive index variability in their inversion solutions through mixing with other fine mode aerosol models, or optimizing the single scattering albedo, but Earth system models surveyed for this paper assume a globally uniform, size-invariant refractive index. Although no Martian dust samples have yet been returned to Earth, remote sensing observations indicate that Martian dust is globally homogenous in composition, and a single spectral refractive index assumption has been widely adopted to represent Martian dust. The lack of comprehensive, statistically representative measurements of dust particle microphysical properties (size distribution, morphology, complex index of refraction spectra, internal structure heterogeneity), and the resulting optical properties, limits our ability to verify the fidelity of these assumptions. A chain of measurements is needed, ranging from characterizing individual dust mineralogy components (e.g., pure hematite and goethite) to in situ sampling of complex atmospheric aerosol mixtures. Such results could be applied to both remote sensing retrievals that characterize the optical properties of the total aerosol burden in the atmosphere from total radiance measurements, and to global models that represent the total aerosol burden in the atmosphere by building it up from the balance of individual aerosol sources and sinks.

Aeolian dust↗

NASA NTRS2000

Radar Observations of Asteroid 10115 (1992 SK)

The approach of 1992 SK to within 0.056 AU of Earth on 1999 March 26 provided an excellent opportunity for radar observations. We observed it daily between March 22-26 at Goldstone at a transmitter frequency of 8560 MHz (3.50 cm). Weighted, optimally smoothed sums of cw (Doppler-only) echoes achieve a signal-to-noise ratio of about 300. Our highest resolution cw setup places up to approx. 30 0.5-Hz Doppler cells on the target at echo powers exceeding two standard deviations of the noise. Our highest resolution imaging setup places about 20 0.25-us (37.5 m) x 2.0 Hz pixels on 1992 SK at echo powers above the same threshold. We completed nearly 300 runs with each setup and obtained thorough rotational phase coverage (period = 7.320 h, P. Pravec, pers. comm.). The dispersion of the echoes in time delay indicates a lower bound on 1992 SK's maximum pole-on breadth of approx. 0.4 km that is consistent with the echo bandwidth and rotation period. Variations in the cw spectral bandwidths and shapes and in the delay-Doppler images as a function of rotation phase are evident on each day, are consistent in bandwidth with each other and with the rotation period, and indicate that the asteroid is asymmetric and modesty elongated.

Benner, L. A. M.↗

NASA NTRSDate not supplied

Progress of 2.05 uM Fiber Laser Development for A Martian Co2 Dial

We have proposed a new concept for differential absorption lidar (DIAL) operating in the 2050 nm CO 2 absorption band for atmospheric CO 2 and pressure observations on Mars. This concept has earned us funding from NASA's PICASSO Program to advance 2050 nm fiber laser technology for future space applications. The laser design is an all-fiber master oscillator and power amplifier (MOPA) system. The master oscillator with pulse shaping and several stages of fiber preamplifiers has been built. Preliminary tests show the output pulse energy can reach 1 mJ, meeting our baseline goal. However, the current power amplifier design has an issue with undesired parasitic lasing at a different frequency, which exhausts pump energy and limits the laser energy. We have modified the power amplifier design to suppress parasitic lasing, and higher laser power energy is expected. To achieve high measurement accuracy, the laser frequency must be stabilized. A fraction of the master oscillator is split to lock the laser frequency to the center of the selected absorption line. The laser can then be shifted to a frequency far from the absorption line center, where extinction due to CO 2 and other trace gases is minimal, serving as the offline reference. The online wavelength is optimized at 2050.44156 nm, ensuring a CO 2 absorption optical depth (AOD) of approximately 1.1 at 3 km, which maximizes the signal-to-noise ratio (SNR) for measurements in the lower Martian atmosphere. At the conference, we will provide more details about this project and report on the progress of instrument development.

DIAL↗

NASA NTRS1991

On the theory of the synthesis of offset dual-shaped reflectors - Case examples

In an earlier paper by Galindo-Israel et al. (1987), the geometrical optics (GO) principles, constraints, and requirements of the dual- and single-offset-shaped reflector synthesis problem were collected and developed into a set of nonlinear first-order PDEs. An extension of the methods by which solutions to the PDEs can be obtained is presented, together with several case examples. These examples are independently analyzed by GO and physical optics diffraction methods. The starting point for the integration over each reflector can be taken on the outer rim, at the center, or at an intermediate point-the intermediate starting point being the more general case. The utility of the speed of this synthesis method is demonstrated. For example, the program utilized completes a synthesis in less than 45 s on a 386/20 personal computer and within a fraction of a second on a mainframe Cray. This makes practical the incorporation of the synthesis into a search algorithm that can optimize one or more parameters of the reflector system. As an example, the optimization of the mapping equations for low cross polarization is discussed.

Galindo-Israel, Victor↗

NASA NTRS1997

Fiber Contraction Approaches for Improving CMC Proportional Limit

The fact that the service life of ceramic matrix composites (CMC) decreases dramatically for stresses above the CMC proportional limit has triggered a variety of research activities to develop microstructural approaches that can significantly improve this limit. As discussed in a previous report, both local and global approaches exist for hindering the propagation of cracks through the CMC matrix, the physical source for the proportional limit. Local approaches include: (1) minimizing fiber diameter and matrix modulus; (2) maximizing fiber volume fraction, fiber modulus, and matrix toughness; and (3) optimizing fiber-matrix interfacial shear strength; all of which should reduce the stress concentration at the tip of cracks pre existing or created in the matrix during CMC service. Global approaches, as with pre-stressed concrete, center on seeking mechanisms for utilizing the reinforcing fiber to subject the matrix to in-situ compressive stresses which will remain stable during CMC service. Demonstrated CMC examples for the viability of this residual stress approach are based on strain mismatches between the fiber and matrix in their free states, such as, thermal expansion mismatch and creep mismatch. However, these particular mismatch approaches are application limited in that the residual stresses from expansion mismatch are optimum only at low CMC service temperatures and the residual stresses from creep mismatch are typically unidirectional and difficult to implement in complex-shaped CMC.

DiCarlo, James A.↗

NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the Overflow Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent System Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and unique grid systems optimized for any given testing condition. Solution independence from grid resolution was determined by successively increasing grid refinement until asymptotic convergence of the mean aerodynamic loads was observed. Overflow Reynolds-averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinusoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle. Initial investigations into the effect of turbulence modeling fidelity demonstrated significant differences between the aerodynamic loads simulated with RANS and Detached Eddy Simulation (DES). Multiple flow mechanisms were identified as root causes of these differences, including bow shock shape augmentation and reduced entrainment due to reduced turbulence in the DES simulations. Initial comparisons of pitching moment control authority and aerodynamic drag performance between the tested nozzle configurations demonstrated that the 1E configuration may be more advantageous for the Mars Entry, Descent, and Landing (EDL) task.

Supersonic Retropropulsion↗

NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the Overflow Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent System Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and unique grid systems optimized for any given testing condition. Solution independence from grid resolution was determined by successively increasing grid refinement until asymptotic convergence of the mean aerodynamic loads was observed. Overflow Reynolds-averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinusoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle. Initial investigations into the effect of turbulence modeling fidelity demonstrated significant differences between the aerodynamic loads simulated with RANS and Detached Eddy Simulation (DES). Multiple flow mechanisms were identified as root causes of these differences, including bow shock shape augmentation and reduced entrainment due to reduced turbulence in the DES simulations. Initial comparisons of pitching moment control authority and aerodynamic drag performance between the tested nozzle configurations demonstrated that the 1E configuration may be more advantageous for the Mars Entry, Descent, and Landing (EDL) task.

Supersonic Retropropulsion↗

NASA NTRSDate not supplied

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

NASA NTRS2021

A Summary of Normal Incidence Tube Tests of NASA'S Bio-inspired Broadband Acoustic Absorber

Motivated by the need to reduce noise pollution from aircraft engines, NASA has continued to design, manufacture and test recently patented structures that mimic the geometry and the broadband and low frequency acoustic absorption of assemblies of natural reeds. This report documents a screening test, which is a precursor to more controlled modelling experiments. The purpose of this experiment was to identify variables and prototype design features that might play a dominant role in the acoustic performance of these structures. To date, twenty-four prototypes, each with an overall shape of a cube (nominally 50 mm x 50 mm x 50 mm), were designed and additively manufactured from thermoplastics using a Fused Filament Fabrication technique. The prototypes were tested in the NASA Glenn and Langley Normal Incidence Tubes to experimentally determine the acoustic absorption as a function of frequency from 400 Hz to 3000 Hz. Results indicate that a variety of structures exhibit substantial acoustic absorption in that frequency range, with an average absorption coefficient greater than 0.6, especially in the frequency range of 400 Hz to 1000 Hz. Six basic prototype designs were chosen to inform future modeling research. Results of these experiments can be used to evaluate existing physics-based models of the interaction of sound waves with these structures or to develop new models. Validated physics-based models of these structures can help engineers optimize the designs for different commercial noise control applications.

Acoustics↗

NASA NTRS2011

Impact of Measurement System Characteristics on Advanced Sounder Information Content

Advanced satellite sensors are tasked with improving global observations of the Earth's atmosphere, clouds, and surface to enable enhancements in weather prediction, climate monitoring capability, and environmental change detection. Achieving such an improvement in geophysical information inferred from these observations requires optimal usage of data from current systems as well as instrument system enhancements for future sensors. This presentation addresses results of tradeoff studies evaluating the impact of spectral resolution, spectral coverage, instrument noise, and a priori knowledge on remote sensing system information content, with a specific emphasis on thermodynamic state and trace species information obtainable from advanced atmospheric sounders. Particular attention will be devoted toward information achievable from the Atmospheric InfraRed Sounder (AIRS) on the NASA EOS Aqua satellite in orbit since 2002, the Infrared Atmospheric Sounding Interferometer (IASI) aboard MetOp-A since 2006, and the Cross-track Infrared Sounder (CrIS) instrument to fly aboard the NPP and JPSS series of satellites expected to begin in late 2011. While all of these systems cover nearly the same infrared spectral extent, they have very different number of channels, instrument line shapes, coverage continuity, and instrument noise. AIRS is a grating spectrometer having 2378 discrete spectral channels ranging from about 0.4 to 2.2/cm resolution; IASI is a Michelson interferometer with 8461 uniformly-spaced spectral channels of 0.5/cm (apodized) resolution; and CrIS is a Michelson interferometer having 1305 spectral channels of 0.625, 1.250, and 2.50/cm (unapodized) spectral resolution, respectively, over its three continuous but non-overlapping bands. Results of tradeoff studies showing information content sensitivity to assumed measurement system characteristics will be presented.

Larar, Allen M.↗

NASA NTRS2014

Development and Characterization of Embedded Sensory Particles Using Multi-Scale 3D Digital Image Correlation

A method for detecting fatigue cracks has been explored at NASA Langley Research Center. Microscopic NiTi shape memory alloy (sensory) particles were embedded in a 7050 aluminum alloy matrix to detect the presence of fatigue cracks. Cracks exhibit an elevated stress field near their tip inducing a martensitic phase transformation in nearby sensory particles. Detectable levels of acoustic energy are emitted upon particle phase transformation such that the existence and location of fatigue cracks can be detected. To test this concept, a fatigue crack was grown in a mode-I single-edge notch fatigue crack growth specimen containing sensory particles. As the crack approached the sensory particles, measurements of particle strain, matrix-particle debonding, and phase transformation behavior of the sensory particles were performed. Full-field deformation measurements were performed using a novel multi-scale optical 3D digital image correlation (DIC) system. This information will be used in a finite element-based study to determine optimal sensory material behavior and density.

Cornell, Stephen R.↗

NASA NTRS2005

The Influence of Local Geometric Effects on Mars Polar Processes

Using simple, qualitative heat balance models, this paper addresses textures and structures that will result from the evolution of volatile layers by accretion and by ablation. Such phenomena may have global implications that are not apparent when only flat or sloped surfaces are modeled. In general, structures such as mounds or depressions formed out of volatile materials will evolve in shape such that the growth or retreat of any particular surface will be maximized. It can be shown that the local radius of curvature is proportional to the growth or retreat rate. For example, icy surfaces will tend to form facets that face the dominant sun direction. Two such cases are evaluated: a) Features associated with condensation of volatiles, include cold-trapping and redistribution, such as the concentration of frost around the Viking 2 lander [1]. Here I will focus on textures that likely result from the formation of seasonal CO2 deposits. b) Features associated with sublimation of volatiles, such as those described by Ingersoll et. al. [2] result in textured surfaces that affect both the apparent emissivity and albedo. Similar calculations have been performed with respect to the "Swiss cheese" features on the South Polar Cap [3]. Here, I evaluate the likely sublimation rates from optimal ice scarp structures and their implications for the long-term evolution of the polar caps and formation of layered terrain.

Hecht, M. H.↗

NASA NTRS1994

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained (i.e., all complicated source code for the derivative calculations is constructed quickly with accuracy). The basic equations for second-order sensitivity derivatives are presented; four methods are compared. Each scheme requires that large systems are solved first for the first-order derivatives and, in all but one method, for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent flow example; both geometric-shape and flow-condition design variables are considered. Several methods are tested; results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method; for six independent variables, it is at least two to four times faster than central finite differences and requires only 60 percent more memory than the original code; the performance is expected to improve further in the future.

Taylor, Arthur C., III↗

NASA NTRS2005

Fractal Reference Signals in Pulse-Width Modulation

A report proposes the use of waveforms having fractal shapes reminiscent of sawteeth (in contradistinction to conventional regular sawtooth waveforms) as reference signals for pulse-width modulation in control systems for thrusters of spacecraft flying in formation. Fractal reference signals may also be attractive in some terrestrial control systems - especially those in which pulse-width modulation is used for precise control of electric motors. The report asserts that the use of fractal reference signals would enable the synchronous control of several variables of a spacecraft formation, such that consumption of propellant would be minimized, intervals between thruster firings would be long (as preferred for performing scientific observations), and delays in controlling large-thrust maneuvers for retargeting would be minimized. The report further asserts that whereas different controllers would be needed for different modes of operation if conventional pulsewidth modulation were used, the use of fractal reference signals would enable the same controller to function nearly optimally in all regimes of operation, so that only this one controller would be needed.

Lurie, Boris↗