Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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 685 records · Page 38

Prominence Mass Supply and the Cavity

A prevalent but untested paradigm is often used to describe the prominence-cavity system; the cavity is under-dense because it it evacuated by supplying mass to the condensed prominence. The thermal non-equilibrium (TNE) model of prominence formation offers a theoretical framework to predict the thermodynamic evolutin of the prominence and the surrounding corona. We examine the evidence for a prominence-cavity connection by comparing the TNE model and diagnostics of dynamic extreme ultraviolet (EUV) emission surrounding the prominence, specifically prominence horns. Horns are correlated extensions of prminence plasma and coronal plasma which appear to connect the prominence and cavity. The TNE model predicts that large-scale brightenings will occur in the Solar Dynamics Observatory Atmospheric Imaging Assembly 171 A badpass near he prominence that are associated with the cooling phase of condensation formation. In our simulations, variations in the magnitude of footpoint heating lead to variations in the duration, spatial scale, and temporal offset between emission enhancements in the other EUV bandpasses. While these predictions match well a subset of the horn observations, the range of variations in the observed structures is not captured by the model. We discuss the implications of one-dimensional loop simulations for the three-dimensional time-averaged equilibrium in the prominence and the cavity. Evidence suggests that horns are likely caused by condensing prominence plasma, but the larger question of whether this process produces a density-depleted cavity requires a more tightly constrained model of heating and better knowledge of the associated magnetic structure.

PROMINENCES↗

Control Design for an Advanced Geared Turbofan Engine

This paper describes the design process for the control system of an advanced geared turbofan engine. This process is applied to a simulation that is representative of a 30,000 lbf thrust class concept engine with two main spools, ultra-high bypass ratio, and a variable area fan nozzle. Control system requirements constrain the non-linear engine model as it operates throughout its flight envelope of sea level to 40,000 ft and from 0 to 0.8 Mach. The control architecture selected for this project was developed from literature and reflects a configuration that utilizes a proportional integral controller integrated with sets of limiters that enable the engine to operate safely throughout its flight envelope. Simulation results show the overall system meets performance requirements without exceeding system operational limits.

aerothermodynamics↗

Control Design for an Advanced Geared Turbofan Engine

This paper describes the design process for the control system of an advanced geared turbofan engine. This process is applied to a simulation that is representative of a 30,000 lbf thrust class concept engine with two main spools, ultra-high bypass ratio, and a variable area fan nozzle. Control system requirements constrain the non-linear engine model as it operates throughout its flight envelope of sea level to 40,000 ft and from 0 to 0.8 Mach. The control architecture selected for this project was developed from literature and reflects a configuration that utilizes a proportional integral controller integrated with sets of limiters that enable the engine to operate safely throughout its flight envelope. Simulation results show the overall system meets performance requirements without exceeding system operational limits.

gas path dynamics↗

Understanding the Puzzling X-Ray Spectrum of the SO Galaxy NGC 4382; NGC 43819: Spectral Analysis of the Prototypical Early Merger and ASCA Observations of a Dynamically Young Elliptical: NGC 4125

We have analyzed the ASCA observations of NGC 4382, NGC 4038/9, NGC 4125 and produced papers for publication. NGC 4382 is one of the E and SO galaxies detected with the lowest X-ray to optical luminosity ratio. These galaxies have a peculiar X-ray (0.1-3 keV) spectrum, with a significant excess of counts in the lowest spectral channels (less than 1 keV) relative to the spectral count distributions of X-ray brihter E and SO galaxies. Analyzing the ROSAT PSPC observation of NGC 4382 it was unclear whether this soft excess was due to a real very soft component in a multi-component spectrum, or reflected an extremely low metal abundance in an isothermal hot gas. Our ASCA observations show that the low-abundance single-temperature model does not fit well to the X-ray spectrum, in agreement with our previous suggestions. A better explanation is a composite spectrum with a very soft component (0.3 keV) in addition to a harder (5 keV) component from X-ray binaries. In this model, the abundance cannot be constrained. More complex spectral models are also possible. The ASCA observations of The Antennae - (NGC4038/9) show that at least two spectral components are required to describe the emission-thermal emission from a plasma at 0.8 keV, and a component at higher energies. The hot gas contributes about half of the flux in the 0.5 to 6 keV band. If the column density to the higher energy component is greater than 2 x 10 (exp 21) per square centimeter, then the fitted abundance in the hot gas component is less than 0.2 solar. This low abundance is not expected for the hot interstellar medium in NGC4038/9 in which supernovae and star formation (expected to enrich and heat the gas) are ongoing. We do not detect any spatial variations in the spectrum. We relate these findings to data obtained by other satellites (Einstein, ROSAT) for this interacting galaxy pair. NGC4125's ASCA data was analyzed jointly with its Beppo-SAX observation. A hard component (kT /sim 4-10 keV) is observed together with a thermal component with line emission.

Fabbiano, Giuseppina↗

Residue detection for real-time removal of paint from metallic surfaces

Paint stripping from large steel ships and other metallic surfaces is a major issue in the maintenance and refurbishing of structures, and environmental concerns are greatly limiting the possible options. As a result, waterjet with water recycling has become the leading form of paint stripping and robotic manipulators with scanning bridges were constructed by various manufacturers to address this need. The application of such scanning bridges is slow and their access is constrained by the complex shape of the ship hull and various features on the surface. To overcome these limitations, a robotic system that is called Ultrastrip (UltraStrip Systems, Inc., Stuart, FL) is developed. This system uses magnetic wheels to attach the stripper to the structure and travel on it while performing paint stripping. To assure efficient paint stripping feedback data is required to control the travel speed by monitoring the paint thickness before and during the stripping process. Efforts at JPL are currently underway to develop the required feedback capability to assure effective paint stripping. Various possible sensors were considered and issues that can affect the sensitivity, reliability and applicability of the sensors are being investigated with emphasis on measuring the initial conditions of the paint. Issues that affect the sensory data in dynamic conditions are addressed while providing real-time real feedback for the control of the paint stripper speed of travel.

NDT↗

Modeling the Spin Motor Current of the International Space Station's Control Moment Gyroscopes

The International Space Station (ISS) attitude control is provided by two means: The Russian Segment uses thrusters and the U.S. Segment uses double-gimbaled control moment gyroscopes (CMG). CMGs are used as momentum exchange devices, providing non propulsive attitude control for the vehicle. The CMGs are very important for the ISS program because, first, they save propellant - which needs to be transferred to the Station in special cargo vehicles - and, second, they provide the microgravity environment on the Station - which is necessary for scientific experiments planned for the ISS mission. Since 2002, when one of the CMG on the ISS failed, all CMGs are closely monitored. High gimbal rates, vibration spikes, unusual variations of spin motor current and bearing temperatures are of great concern, since these parameters are the CMG health indicators. The telemetry analysis of these and some other CMG parameters is used to determine constrains and make changes to the CMGs operation on board. These CMG limitations, in turn, may limit the ISS attitude control capabilities and may be critical to ISS operation. Therefore, it is important to know whether the CMG parameter is nominal or out of family, and why. The goal of this project is to analyze an important CMG parameter - spin motor current. Some operational decisions are made now based on the spin motor current signatures. The spin motor current depends on gimbal rates, ISS rates, and spin bearing friction. The spin bearing friction in turn depends on the bearing temperatures, wheel rates, normal load - which is a function of gimbal and wheel rates - lubrication, etc. The first task of this project is to create a spin motor current mathematical model based on CMG dynamics model and the current knowledge on bearing friction in microgravity.

Pereira, Miguel A.↗

A Ring Accelerator? Unusual Jet Dynamics in the IceCube Candidate PKS 1502+106

On 2019/07/30.86853 UT, IceCube detected a high-energy astrophysical neutrino can-didate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 percent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Ar-ray (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after theIceCube event, reveals evidence for a radio ring structure which develops with time.Several arc-structures evolve perpendicular to the jet ridge line. We find evidence forprecession of a curved jet based on kinematic modelling and a periodicity analysis.An outflowing broad line region (BLR) based on the C IV line emission (Sloan Dig-ital Sky Survey, SDSS) is found. We attribute the atypical ring to an interaction ofthe precessing jet with the outflowing material.We discuss our findings in thecontext of a spine-sheath scenario where the ring reveals the sheath andits interaction with the surroundings (NLR clouds).We find that the radioemission is correlated with theγ-ray emission, with radio lagging theγ-rays. Basedon theγ-ray variability timescale, we constrain theγ-ray emission zone to the BLR(30-200rg) and within the jet launching region. We discuss that the outflowing BLRprovides the external radiation field forγ-ray production via external Compton scat-tering.The neutrino is most likely produced by proton-proton interactionin the blazar zone (beyond the BLR), enabled by episodic encounters ofthe jet with dense clouds, i.e. some molecular cloud in the NLR.

S Britzen↗

Modal Testing of Seven Shuttle Cargo Elements for Space Station

From December 1996 to May 2001, the Modal and Control Dynamics Team at NASA's Marshall Space Flight Center (MSFC) conducted modal tests on seven large elements of the International Space Station. Each of these elements has been or will be launched as a Space Shuttle payload for transport to the International Space Station (ISS). Like other Shuttle payloads, modal testing of these elements was required for verification of the finite element models used in coupled loads analyses for launch and landing. The seven modal tests included three modules - Node, Laboratory, and Airlock, and four truss segments - P6, P3/P4, S1/P1, and P5. Each element was installed and tested in the Shuttle Payload Modal Test Bed at MSFC. This unique facility can accommodate any Shuttle cargo element for modal test qualification. Flexure assemblies were utilized at each Shuttle-to-payload interface to simulate a constrained boundary in the load carrying degrees of freedom. For each element, multiple-input, multiple-output burst random modal testing was the primary approach with controlled input sine sweeps for linearity assessments. The accelerometer channel counts ranged from 252 channels to 1251 channels. An overview of these tests, as well as some lessons learned, will be provided in this paper.

Kappus, Kathy O.↗

Secure Mass Measurements From Transit Timing: 10 Kepler Exoplanets Between 3 and 8M⊕ with Diverse Densities and Incident Fluxes

We infer dynamical masses in eight multiplanet systems using transit times measured from Keplerʼs complete dataset, including short-cadence data where available. Of the 18 dynamical masses that we infer, 10 pass multiple tests for robustness. These are in systems Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler -105 c has a radius of 1.3R⊕ and a density consistent with an Earth-like composition. Strong transit timing variation (TTV) signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently measured transit times, including planets orbiting Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115), and Kepler-177 (KOI-523). Nonetheless, strong upper limits on the mass of Kepler-177 cimply an extremely low density of 0.1 g cm−3. In most cases, individual orbital eccentricities were poorly constrained owing to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate to high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, nonlibrating solutions also match the transit timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have secure mass measurements. Placing these exoplanets on the mass–radius diagram, we find that a wide range of densities is observed among sub-Neptune-mass planets and that the range in observed densities is anticorrelated with incident flux.

Stars individual↗

Light Microscopy Module: An On-Orbit Microscope Planned for the Fluids and Combustion Facility on the International Space Station

The Light Microscopy Module (LMM) is planned as a fully remotely controllable on-orbit microscope subrack facility, allowing flexible scheduling and control of fluids and biology experiments within NASA Glenn Research Center's Fluids and Combustion Facility on the International Space Station. Within the Fluids and Combustion Facility, four fluids physics experiments will utilize an instrument built around a light microscope. These experiments are the Constrained Vapor Bubble experiment (Peter C. Wayner of Rensselaer Polytechnic Institute), the Physics of Hard Spheres Experiment-2 (Paul M. Chaikin of Princeton University), the Physics of Colloids in Space-2 experiment (David A. Weitz of Harvard University), and the Low Volume Fraction Colloidal Assembly experiment (Arjun G. Yodh of the University of Pennsylvania). The first experiment investigates heat conductance in microgravity as a function of liquid volume and heat flow rate to determine, in detail, the transport process characteristics in a curved liquid film. The other three experiments investigate various complementary aspects of the nucleation, growth, structure, and properties of colloidal crystals in microgravity and the effects of micromanipulation upon their properties. Key diagnostic capabilities for meeting the science requirements of the four experiments include video microscopy to observe sample features including basic structures and dynamics, interferometry to measure vapor bubble thin film thickness, laser tweezers for colloidal particle manipulation and patterning, confocal microscopy to provide enhanced three-dimensional visualization of colloidal structures, and spectrophotometry to measure colloidal crystal photonic properties.

Doherty, Michael P.↗

Debris Disk Structure and Morphology as Revealed by Aggressive STIS Multi-Roll Coronagraphy: A New Look at Some Old Friends

We present new imaging results from a well-selected sample of II circumstellar debris disks, all with HST pedigree, using STIS visible-light 6-roll PSF-template subtracted coronagraphy (PSFTSC). These new observations, pushing HST to its highest levels of coronagraphic performance, simultaneously probe both the interior regions of these debris systems, with inner working distances < app 8 AU for half the stars in this sample (corresponding to the giant planet and Kuiper belt regions within our own solar system), and the exterior regions far beyond. These new images enable direct inter-comparison of the architectures of these exoplanetary debris systems in the context of our own Solar System: These observations also permit us, for the first time, to characterize material in these regions at high spatial resolution and identify disk sub-structures that are signposts of planet formation and evolution; in particular, asymmetries and non-uniform debris structures that signal the presence of co-orbiting perturbing planets, and dynamical interactions (e.g., resulting in posited small grain stripping and disk "pollution") with the ISM. We focus here on recently acquired and reduced images of he circumstellar debris systems about: AU Mic (edge-on, and @ 10 pc the closest star in our sample), HD 61005, HD 32297 and HD 15115 (all with morphologies strongly suggestive of ISM wind interactions), HD 181327 & HDI07146 (close to face-on with respectively narrow and broad debris rings), and MP Mus (a "mature" proto-planetary disk hosted by a cTTS). All of our objects were previously observed in the near-IR with inferior spatial resolution and imaging efficacy, but with NICMOS r = 0.3" inner working angle (IWA) comparable to STIS multi-roll coronagraphy. The combination of new optical and existing near-IR imaging can strongly constrain the dust properties, thus enabling an assessment of grain processing and planetesimal populations. These results will directly inform upon the posited planet formation mechanisms that occur after the approximately 10 My epoch of gas depletion, a time in our solar system when giant planets were migrating and terrestrial planets were forming, and directly test theoretical models of these processes. These observations lmiquely probe both into the interior regions of these systems and are sensitive to and spatially resolve low surface-brightness (SB) material at large stellocentric distances with spatial resolution comparable to ACS and with augmenting NICMOS near-IR disk photometry in hand.

Grady, Carol A↗

Picture this SELFI: Submillimeter Enceladus Life Fundamentals Instrument

The Submillimeter Enceladus Life Fundamentals Instrument (SELFI) is a passive remote sensing submillimeter heterodyne spectrometer being developed at NASA GSFC under NASA's Maturation of Instruments for Solar System Exploration (MatISSE) program. SELFI will advance submillimeter receiver technology by 1) investigating the chemical and isotopic compositions and corresponding densities of Enceladus' plume material, their vertical thermal structures, and the transport mechanisms within the plumes, and 2) characterizing the source regions from which the plumes emerge. The Enceladus plumes are important in the context of life and habitability of its subsurface ocean environment. SELFI remote sensing measurements will 1) measure the spatial and temporal variabilities in the plume chemical compositions, 2) provide insight in to Enceladus' subsurface ocean environment by monitoring H2O, HDO, d18O, and d17O, 3) constrain the oxidation state of the subsurface ocean using H2O2 and O3, and 4) utilize the SO2 and H2S spectral signatures to constrain the impact arising from both the sea-floor volcanoes and pre-biotic molecules. Moreover, the detection of the remaining molecular species (14 in total) is vital to improving the current state of knowledge of Enceladus' subsurface ocean habitability this also permits us to explore the chemical alteration processes arising from primordial volatiles that have been observed in comets. Lastly, SELFI's continuum observations enable the correlation between observed variations in plume activity with surface temperatures.SELFI is currently being developed under a technology maturation program that will advance the RF-to-digital electronics of a submillimeter-wave heterodyne spectrometer to simultaneously observe fourteen molecular species with resonances between 530 GHz and 600 GHz. SELFI will have fine radiometric resolution, high spectral resolution (resolving power R > 106), multiple continuum channels and a high dynamical range, necessary to map Enceladus across its 30 K to 250 K temperature range. Under the MatISSE program, SELFI will advance from TRL 4 to 6 four key technologies of the RF-to-digital subsystem, which are: 1) the RF low noise amplifier design; 2) the single-sideband (SSB) mixer and local oscillator; 3) the IF assembly down-converter that maps the fourteen species to 2 x 500 MHz bandwidth; and 4) the digital spectrometer electronics.

Racette, Paul↗

Secure Mass Measurements From Transit Timing: 10 Kepler Exoplanets Between 3 and 8 M⊕ With Diverse Densities and Incident Fluxes

We infer dynamical masses in eight multiplanet systems using transit times measured from Keplerʼs complete data set, including short-cadence data where available. Of the 18 dynamical masses that we infer, 10 pass multiple tests for robustness. These are in systems Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler-105 c has a radius of 1.3 R(sub ⨁) and a density consistent with an Earth-like composition. Strong transit timing variation (TTV) signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently measured transit times, including planets orbiting Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115), and Kepler-177 (KOI-523). Nonetheless, strong upper limits on the mass of Kepler-177 c imply an extremely low density of ∽0.1 g cm(exp −3). In most cases, individual orbital eccentricities were poorly constrained owing to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate to high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, nonlibrating solutions also match the transit timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have secure mass measurements. Placing these exoplanets on the mass–radius diagram, we find that a wide range of densities is observed among sub-Neptune-mass planets and that the range in observed densities is anticorrelated with incident flux.

Individual stars↗

In situ cosmogenic H-3, C-14, and Be-10 for determining the net accumulation and ablation rates of ice sheets

The usefulness of the in situ cosmogenic H-3, C-14, and Be-10 produced by spallation of oxygen nuclei in ice, as tracers to determine net accumulation/ablation rates of ice sheets is explored. The application of the in situ H-3 and Be-10 is severely constrained because at deposition, ice contains appreciable amounts of these isotopes from the atmosphere. The case is much more favorable for C-14, which is not carried with wet precipitations; atmospheric C-14 gets mechanically trapped in the ice during deposition. It is pointed out that cosmogenic C-14 would probably exist as (C-14)O in ice. This seems to be supported by the published results of Fireman and Norris (1981). Considering their inherent amounts in the ice and the expected in situ production rates, conditions under which these isotopes can be used to study net accumulation and ablation rates are discussed. Available data on C-14 and Be-10 on polar ice from accumulation and ablation zones is also discussed. It is concluded that H-3 and C-14 should find wide applications in studying ice dynamics and Be-10 in very special circumstances.

Lal, D.↗

Kalman Filter Constraint Tuning for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints are often neglected because they do not fit easily into the structure of the Kalman filter. Recently published work has shown a new method for incorporating state variable inequality constraints in the Kalman filter, which has been shown to generally improve the filter s estimation accuracy. However, the incorporation of inequality constraints poses some risk to the estimation accuracy as the Kalman filter is theoretically optimal. This paper proposes a way to tune the filter constraints so that the state estimates follow the unconstrained (theoretically optimal) filter when the confidence in the unconstrained filter is high. When confidence in the unconstrained filter is not so high, then we use our heuristic knowledge to constrain the state estimates. The confidence measure is based on the agreement of measurement residuals with their theoretical values. The algorithm is demonstrated on a linearized simulation of a turbofan engine to estimate engine health.

Simon, Dan↗

NASA Langley's Formal Methods Research in Support of the Next Generation Air Transportation System

This talk will provide a brief introduction to the formal methods developed at NASA Langley and the National Institute for Aerospace (NIA) for air traffic management applications. NASA Langley's formal methods research supports the Interagency Joint Planning and Development Office (JPDO) effort to define and develop the 2025 Next Generation Air Transportation System (NGATS). The JPDO was created by the passage of the Vision 100 Century of Aviation Reauthorization Act in Dec 2003. The NGATS vision calls for a major transformation of the nation s air transportation system that will enable growth to 3 times the traffic of the current system. The transformation will require an unprecedented level of safety-critical automation used in complex procedural operations based on 4-dimensional (4D) trajectories that enable dynamic reconfiguration of airspace scalable to geographic and temporal demand. The goal of our formal methods research is to provide verification methods that can be used to insure the safety of the NGATS system. Our work has focused on the safety assessment of concepts of operation and fundamental algorithms for conflict detection and resolution (CD&R) and self- spacing in the terminal area. Formal analysis of a concept of operations is a novel area of application of formal methods. Here one must establish that a system concept involving aircraft, pilots, and ground resources is safe. The formal analysis of algorithms is a more traditional endeavor. However, the formal analysis of ATM algorithms involves reasoning about the interaction of algorithmic logic and aircraft trajectories defined over an airspace. These trajectories are described using 2D and 3D vectors and are often constrained by trigonometric relations. Thus, in many cases it has been necessary to unload the full power of an advanced theorem prover. The verification challenge is to establish that the safety-critical algorithms produce valid solutions that are guaranteed to maintain separation under all possible scenarios. Current research has assumed perfect knowledge of the location of other aircraft in the vicinity so absolute guarantees are possible, but increasingly we are relaxing the assumptions to allow incomplete, inaccurate, and/or faulty information from communication sources.

Butler, Ricky W.↗

Attitude Control Approach for Solar Cruiser, a Large, Deep Space Solar Sail

Solar Cruiser is a small satellite Technology Demonstration Mission (TDM) to mature solar sail propulsion technology using a solar sail larger than 1600 square meters, demonstrating performance as a propulsion system and a stable pointing platform for science observations in an artificial halo orbit sunward of the Sun-Earth Lagrange Point 1 (sub-L1). For the “sailcraft” to meet mission objectives, there are several unique attitude control challenges that the Attitude Determination and Control System (ADCS) must overcome. Large disturbance torques, primarily due to sail deformations coupled with off-sun pointing angles, make it more difficult to maintain adequate controls performance and manage accumulated momentum on the control actuators. The large-amplitude, low-frequency flexible body modes of the sail, in concert with noisy sensors and actuators, make it challenging to mitigate control-structure interactions and maintain fine pointing capabilities. Stability and control performance during sail deployment is complicated by rapidly and widely varying inertias. The Solar Cruiser ADCS effectively addresses these challenges using a simple, traditional control system – including momentum management actuators that use bang-bang control to constrain internal accumulated momentum, a low-pass controls filter and an attitude determination Kalman Filter (KF) that blends multiple star tracker solutions, and a reaction wheel assembly (RWA) with controller gains tuned to specific configurations or inertias. Solar Cruiser’s attitude control approach, and lessons learned from its development, establishes a state of the art of high value to future solar sail missions and other small spacecraft with large deployable structures operating in deep space.

solar sail↗

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.↗