Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized linear model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

Computational and Experimental Study of Energetic Materials in a Counterflow Microgravity Environment

Counterflow diffusion flames are studied for various fuels flowing against decomposition products from solid ammonium perchlorate (AP) pellets in order to obtain fundamental understanding of composite propellant flame structure and chemistry. We illustrate this approach through a combined experimental and numerical study of a fuel mixture consisting of C2H4 CO + H2, and C2H2 + C2H4 flowing against solid AP. For these particular AP-fuel systems, the resulting flame zone simulates the various flame structures that are ex+ to exist between reaction products from Ap crystals and a hydrocarbon binder. As in all our experimental studies, quantitative species and temperature profiles have been measured between the fuel exit and AP surface. Species measured included CN, NH, NO, OH, N2, CO2, CO, H2, CO, HCl, and H2O. Temperature was measured using a thermocouple at the exit, spontaneous Raman scattering measurements throughout the flame, OH rotational population distributions, and NO vibrational population distributions. The burning rate of AP was also measured as a function of strain rate, given by the separation distance between the AP surface and the gaseous hydrocarbon fuel tube exit plane. This distance was nominally set at 5 mm, although studies have been performed for variations in separation distance. The measured 12 scalars are compared with predictions from a detailed gas-phase kinetics model consisting of 86 species and 531 reactions. Model predictions are found to be in good agreement with experiment and illustrate the type of kinetic features that may be expected to occur in propellants when AP particle size distributions are varied. Furthermore, the results constitute the continued development of a necessary database and validation of a comprehensive model for studying more complex AP-solid fuel systems in microgravity. Exploratory studies have also been performed with liquid and solid fuels at normal gravity. Because of melting (and hence dripping) and deep thermal wave penetration into the liquid, these experiments were found feasible, but not used for obtaining quantitative data. Microgravity experiments are needed to eliminate the dripping and boiling phenomena of these systems at normal gravity. Microgravity tests in the NASA Glenn 2.2 second drop tower were performed (1) to demonstrate the feasibility of performing propellant experiments using the NASA Glenn microgravity facilities, (2) to develop the operational procedures for safe handing of the energetic materials and disposal of their toxic combustion by-products and (3) to obtain initial measurements of the AP burning rate and flame structure under microgravity conditions. Experiments were conducted on the CH4/AP system previously studied at normal gravity using a modified design of the counterflow burner and a NASA Glenn Pig Rig, i.e., one of the existing drop rigs for general-purpose usage. In these experiments, the AP burning rate was measured directly with a linear variable differential transducer (LVDT) and video imaging of the flame structure was recorded ignition was achieved by hot wires stretched across the AP surfaces. Initial drop tower combustion data show that with the same burner separation distance and flow conditions of the normal gravity experiments, the AP burning rate is approximately a factor of two lower. This difference is likely a result of radiation effects, but further tests with longer test times need to be conducted to verify that steady state conditions were achieved under microgravity conditions.

Takahashi, Fumiaki↗

Unsteady High Turbulence Effects on Turbine Blade Film Cooling Heat Transfer Performance Using a Transient Liquid Crystal Technique

Unsteady wake effect, with and without trailing edge ejection, on detailed heat transfer coefficient and film cooling effectiveness distributions is presented for a downstream film-cooled gas turbine blade. Tests were performed on a five-blade linear cascade at an exit Reynolds number of 5.3 x 10(exp 5). Upstream unsteady wakes were simulated using a spoke-wheel type wake generator. Coolant blowing ratio was varied from 0.4 to 1.2; air and CO2 were used as coolants to simulate different density ratios. Surface heat transfer and film effectiveness distributions were obtained using a transient liquid crystal technique; coolant temperature profiles were determined with a cold wire technique. Results show that Nusselt numbers for a film cooled blade are much higher compared to a blade without film injection. Unsteady wake slightly enhances Nusselt numbers but significantly reduces film effectiveness versus no wake cases. Nusselt numbers increase only slic,htly but film cooling, effectiveness increases significantly with increasing, blowing ratio. Higher density coolant (CO2) provides higher effectiveness at higher blowing ratios (M = 1.2) whereas lower density coolant (Air) provides higher 0 effectiveness at lower blowing ratios (M = 0.8). Trailing edge ejection generally has more effect on film effectiveness than on the heat transfer, typically reducing film effectiveness and enhancing heat transfer. Similar data is also presented for a film cooled cylindrical leading edge model.

Han, J. C.↗

Global Free-tropospheric NO2 Abundances Derived Using a Cloud Slicing Technique from AURA OMI

We derive free-tropospheric NO2 volume mixing ratios (VMRs) by applying a cloud-slicing technique to data from the Ozone Monitoring Instrument (OMI) on the Aura satellite. In the cloud-slicing approach, the slope of the above-cloud NO2 column versus the cloud scene pressure is proportional to the NO2 VMR. In this work, we use a sample of nearby OMI pixel data from a single orbit for the linear fit. The OMI data include cloud scene pressures from the rotational-Raman algorithm and above-cloud NO2 vertical column density (VCD) (defined as the NO2 column from the cloud scene pressure to the top of the atmosphere) from a differential optical absorption spectroscopy (DOAS) algorithm. We compare OMI-derived NO2 VMRs with in situ aircraft profiles measured during the NASA Intercontinental Chemical Transport Experiment Phase B (INTEX-B) campaign in 2006. The agreement is generally within the estimated uncertainties when appropriate data screening is applied. We then derive a global seasonal climatology of free-tropospheric NO2 VMR in cloudy conditions. Enhanced NO2 in the free troposphere commonly appears near polluted urban locations where NO2 produced in the boundary layer may be transported vertically out of the boundary layer and then horizontally away from the source. Signatures of lightning NO2 are also shown throughout low and middle latitude regions in summer months. A profile analysis of our cloud-slicing data indicates signatures of lightning-generated NO2 in the upper troposphere. Comparison of the climatology with simulations from the global modeling initiative (GMI) for cloudy conditions (cloud optical depth less than10) shows similarities in the spatial patterns of continental pollution outflow. However, there are also some differences in the seasonal variation of free-tropospheric NO2 VMRs near highly populated regions and in areas affected by lightning-generated NOx.

NO2↗

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗

Interannual Variations of Shallow Firn Temperature at Greenland Summit

Firn-temperature profiles are calculated in a thermal model using continuous surface temperatures derived from Automatic Weather Station (AWS) data and passive microwave data in the Greenland Summit region during the period 1987-1999. The results show that significant interannual variations of mean summer (June to August) and annual temperatures occur in the top 15 m, in addition to the normal seasonal cycle of firn temperature. At 5 m depth, the seasonal cycle is damped to 13% of the surface seasonal amplitude, but even at 15 m about 1% or 0.6 C of the seasonal cycle persists. Both summer and mean annual temperatures decrease from 1987 to 1992, followed by a general increasing trend. Interannual variability is 5 C at the surface, but only is only dampened to 3.2 C at 10 m depth and 0.7 C at 15 m depth. Dampening of the interannual variability with depth is slower than dampening of the seasonal cycle, because of the longer time constant of the interannual variation. The warmer spring and summer temperatures experienced in the top 5 m, due to both the seasonal cycle and interannual variations, affect the rate of firn densification, which is non-linearly dependent on temperature. During the 12 year period 1987-1999, the mean annual surface temperature is -29.2 C, and the mean annual 15 m temperature is -30. 1 C, which is more than 1 C warmer than a 15-m borehole temperature representing the period of about 1959 and warmer than the best-fit temperature history by Alley and Koci back to 1500 A.D..

Jun, Li↗

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'propagate-burn-propagate' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons between it and three methods using modified linearized covariance propagation are made. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. In the second method, the maneuver is not modeled but the uncertainty in its magnitude is accounted for by the inclusion of a process noise matrix. In the third method, which is essentially a hybrid of the first two, the nominal portion of the maneuver is included via the state transition matrix while a process noise matrix is used to account for the magnitude uncertainty. Since this method also correctly accounts for the presence of the maneuver in the nominal orbit, it is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, Pc. However, applications for the two other methods exist and are briefly discussed. Despite the fact that the process model ('propagate-burn-propagate') that was studied was very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗

Advanced Modeling of Control-Structure Interaction in Thrust Vector Control Systems

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of 8 mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. Two actuators are used to move each engine in two planes perpendicular to one another (i.e., pitch and yaw). The TVC system design leverages hardware from the Space Shuttle program as well as new hardware designed specifically for the Core Stage. During the development of the SLS TVC system, a family of advanced dynamics models were developed to extend and compliment the simplified quasi-linear “simplex” model historically used for flight control design and stability analysis. The importance of these advanced models became increasingly evident after ambient and hot fire testing of the Core Stage, which revealed a number of findings associated with the dynamic response of the TVC integrated system. Test responses suggested that the TVC did not meet its performance specifications and its step and frequency responses exhibited unexpected departures from prior lab tests and modeled behavior. One driving factor for these results was a higher-than-expected degree of coupling between the TVC system, the engine dynamics, and the Core Stage structure. This paper is the third installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, a new method of modeling rocket vehicle thrust vectoring servoelastic dynamics is presented. In this approach, the load dynamics are replaced by a detailed finite element model containing both the rigid body and elastic modes. A partitioning technique is used to compute the effective compliance from the modal data and obtain accurate simulation results using a reduced number of generalized coordinates. Coupled backup structure and nozzle attach compliance effects on multiple engines are captured in higher fidelity than with a spring approximation, eliciting novel effects due to the complex load paths involved in the Core Stage structure. Validation of the model is demonstrated using a variety of structural/modal, laboratory, and full-scale hot fire test data.

Launch Vehicles↗

Advanced Modeling of Control-Structure Interaction in Thrust Vector Control Systems

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of 8 mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. Two actuators are used to move each engine in two planes perpendicular to one another (i.e., pitch and yaw). The TVC system design leverages hardware from the Space Shuttle program as well as new hardware designed specifically for the Core Stage. During the development of the SLS TVC system, a family of advanced dynamics models were developed to extend and compliment the simplified quasi-linear “simplex” model historically used for flight control design and stability analysis. The importance of these advanced models became increasingly evident after ambient and hot fire testing of the Core Stage, which revealed a number of findings associated with the dynamic response of the TVC integrated system. Test responses suggested that the TVC did not meet its performance specifications and its step and frequency responses exhibited unexpected departures from prior lab tests and modeled behavior. One driving factor for these results was a higher-than-expected degree of coupling between the TVC system, the engine dynamics, and the Core Stage structure. This paper is the third installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, a new method of modeling rocket vehicle thrust vectoring servoelastic dynamics is presented. In this approach, the load dynamics are replaced by a detailed finite element model containing both the rigid body and elastic modes. A partitioning technique is used to compute the effective compliance from the modal data and obtain accurate simulation results using a reduced number of generalized coordinates. Coupled backup structure and nozzle attach compliance effects on multiple engines are captured in higher fidelity than with a spring approximation, eliciting novel effects due to the complex load paths involved in the Core Stage structure. Validation of the model is demonstrated using a variety of structural/modal, laboratory, and full-scale hot fire test data.

Launch Vehicles↗

Evaluating LSM-Based Water Budgets Over a West African Basin Assisted with a River Routing Scheme

Within the framework of the African Monsoon Multidisciplinary Analysis (AMMA) Land Surface Model Intercomparison Project phase 2 (ALMIP-2), this study evaluates the water balance simulated by the Interactions between Soil, Biosphere, and Atmosphere (ISBA) over the upper Oum River basin, in Benin, using a mesoscale river routing scheme (RRS). The RRS is based on the nonlinear Muskingum Cunge method coupled with two linear reservoirs that simulate the time delay of both surface runoff and base flow that are produced by land surface models. On the basis of the evidence of a deep water-table recharge in that region,a reservoir representing the deep-water infiltration (DWI) is introduced. The hydrological processes of the basin are simulated for the 2005-08 AMMA field campaign period during which rainfall and stream flow data were intensively collected over the study area. Optimal RRS parameter sets were determined for three optimization experiments that were performed using daily stream flow at five gauges within the basin. Results demonstrate that the RRS simulates stream flow at all gauges with relative errors varying from -22% to 3% and Nash-Sutcliffe coefficients varying from 0.62 to 0.90. DWI varies from 24% to 67% of the base flow as a function of the sub-basin. The relatively simple reservoir DWI approach is quite robust, and further improvements would likely necessitate more complex solutions (e.g., considering seasonality and soil type in ISBA); thus, such modifications are recommended for future studies. Although the evaluation shows that the simulated stream flows are generally satisfactory, further field investigations are necessary to confirm some of the model assumptions.

Multidisciplinary↗

Second-order closure PBL model with new third-order moments: Comparison with LES data

This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.

Canuto, V. M.↗

Effects of eddy initial conditions on nonlinear forcing of planetary scale waves by amplifying baroclinic eddies

The previous study of Young and Villere concerning growth of planetary scale waves forced by wave-wave interactions of amplifying intermediate scale baroclinic eddies is extended to investigate effects of different eddy initial conditions. A global, spectral, primitive equation model is used for the calculations. For every set of eddy initial conditions considered, growth rates of planetary modes are considerably greater than growth rates computed from linear instability theory for a fixed zonally independent basic state. However, values of growth rates ranged over a factor of 3 depending on the particular set of eddy initial conditions used. Nonlinear forcing of planetary modes via wave-wave coupling becomes more important than baroclinic growth on the basic state at small values of the intermediate-scale modal amplitudes. The relative importance of direct transfer of kinetic energy from intermediate scales of motion to a planetary mode, compared to baroclinic conversion of available potential energy to kinetic energy within that planetary mode, depends on the individual case. In all cases, however, the transfer of either kinetic or available potential energy to the planetary modes was accomplished principally by wave-wave transfer from intermediate scale eddies, rather than from the zonally averaged state. The zonal wavenumber 2 planetary mode was prominent in all solutions, even in those for which eddy initial conditions were such that a different planetary mode was selectively forced at the start. General characteristics of the structural evolution of the planetary wave components of total heat and momentum flux, and modal structures themselves, were relatively insensitive to variations in eddy initial conditions, even though quantitative details varied from case to case.

Young, Richard E.↗

Pluto’s Sputnik Planitia: Composition of geological units from infrared spectroscopy

We have compared spectroscopic data of Sputnik Planitia on Pluto, as acquired by New Horizons’ Linear Etalon Imaging Spectral Array (LEISA) instrument, to the geomorphology as mapped by White et al. (2017) using visible and panchromatic imaging acquired by the LOng-Range Reconnaissance Imager (LORRI) and the Multi-spectral Visible Imaging Camera (MVIC). We have focused on 13 of the geologic units identified by White et al. (2017), which include the plains and mountain units contained within the Sputnik basin. We divided the map of Sputnik Planitia into 15 provinces, each containing one or more geologic units, and we use LEISA to calculate the average spectra of the units inside the 15 provinces. Hapke-based modeling was then applied to the average spectra of the units to infer their surface composition, and to determine if the composition resulting from the modeling of LEISA spectra reflects the geomorphologic analyses of LORRI data, and if areas classified as being the same geologically, but which are geographically separated, share a similar composition. We investigated the spatial distribution of the most abundant ices on Pluto’s surface - CH4, N2, CO, H2O, and a non-ice component presumed to be a macromolecular carbon-rich material, termed a tholin, that imparts a positive spectral slope in the visible spectral region and a negative spectral slope longward of ~1.1 μm. Because the exact nature of the non-ice component is still debated and because the negative spectral slope of the available tholins in the near infrared does not perfectly match the Pluto data, for spectral modeling purposes we reference it generically as the negative spectral slope endmember (NSS endmember). We created maps of variations in the integrated band depth (from LEISA data) and areal mass fraction (from the modeling) of the components. The analysis of correlations between the occurrences of the endmembers in the geologic units led to the observation of an anomalous suppression of the strong CH4 absorption bands in units with compositions that are dominated by H2O ice and the NSS endmember. Exploring the mutual variation of the CH4 and N2 integrated band depths with the abundance of crystalline H2O and NSS endmember revealed that the NSS endmember is primarily responsible for the suppression of CH4 absorptions in mountainous units located along the western edge of Sputnik Planitia. Our spectroscopic analyses have provided additional insight into the geological processes that have shaped Sputnik Planitia. A general increase in volatile abundance from the north to the south of Sputnik Planitia is observed. Such an increase first observed and interpreted by Protopapa et al., 2017 and later confirmed by climate modeling (Bertrand et al., 2018) is expressed geomorphologically in the form of preferential deposition of N2 ice in the upland and mountainous regions bordering the plains of southern Sputnik Planitia. Relatively high amounts of pure CH4 are seen at the southern Tenzing Montes, which are a natural site for CH4 deposition owing to their great elevation and the lower insolation they are presently receiving. The NSS endmember correlates the existence of tholins within certain units, mostly those coating the low-latitude mountain ranges that are co-latitudinal with the tholin-covered Cthulhu Macula. The spectral analysis has also revealed compositional differences between the handful of occurrences of northern non-cellular plains and the surrounding cellular plains, all of which are located within the portion of Sputnik Planitia that is presently experiencing net sublimation of volatiles, and which do not therefore exhibit a surface layer of bright, freshly-deposited N2 ice. The compositional differences between the cellular and non-cellular plains here hint at the effectiveness of convection in entraining and trapping tholins within the body of the cellular plains, while preventing the spread of such tholins to abutting non-cellular plains.

Pluto's Sputnik Planitia↗

Verification of TROPOMI NO2 Product Using OMI NO2 algorithm

We evaluated S5P TROPOspheric Monitoring Instrument (TROPOMI) operational nitrogen dioxide (NO2) product by comparing with the NO2 retrievals from Ozone Monitoring Instrument (OMI) onboard NASA’s Aura satellite. We compared spatially matched NO2 vertical column density (VCD) data from OMI and TROPOMI to identify any discrepancies between the two operational products. We also applied OMI Air-Mass Factor (AMF) algorithm to TROPOMI Slant Column Density (SCD) data (v02.03.01) aiming to create long-term merged NO2 Earth Science Data Record (https://disc.gsfc.nasa.gov/datasets/TROPOMI_MINDS_NO2_1.1/summary ). The algorithm (version 4.0) uses new Geometry-dependent Lambertian Equivalent surface Reflectivity (GLER) product available for each TROPOMI pixel. GLER is pre-calculated using the vector linearized discrete ordinate radiative transfer (VLIDORT) model, which uses as input high-resolution bidirectional reflectance distribution function (BRDF) information from NASA's Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) over land and the wind-dependent Cox–Munk wave-facet slope distribution over water, the latter with a contribution from the water-leaving radiance based on MODIS gap-filled in-water chlorophyll-a data. The GLER data, combined with consistently retrieved cloud parameters, provide improved information for the calculation of the scattering weight profiles. Additional AMF improvements result from using high resolution (0.25o latitude x 0.25o longitude) a priori NO2 profile shapes and other auxiliary information from the Global Modeling Initiative (GMI) Replay simulation sampled using S5P orbital simulator. Our TROPOMI_MINDS_NO2 algorithm employs the stratosphere-troposphere separation scheme, de-striping, and surface snow/ice treatment consistently with the latest OMI NO2 re-processing (version 4). We evaluate the TROPOMI_MINDS_NO2 product by comparing with the S5P operational NO2 product as well as with independent NO2 observations from ground-based Pandora and aircraft observations. Our results indicate that the new TROPOMI_MINDS_NO2 retrievals are generally higher than the S5P operational NO2 product over polluted regions and show improved agreement with independent validation data.

TROPOMI↗

General Tool for Evaluating High-Contrast Coronagraphic Telescope Performance Error Budgets

The Coronagraph Performance Error Budget (CPEB) tool automates many of the key steps required to evaluate the scattered starlight contrast in the dark hole of a space-based coronagraph. The tool uses a Code V prescription of the optical train, and uses MATLAB programs to call ray-trace code that generates linear beam-walk and aberration sensitivity matrices for motions of the optical elements and line-of-sight pointing, with and without controlled fine-steering mirrors (FSMs). The sensitivity matrices are imported by macros into Excel 2007, where the error budget is evaluated. The user specifies the particular optics of interest, and chooses the quality of each optic from a predefined set of PSDs. The spreadsheet creates a nominal set of thermal and jitter motions, and combines that with the sensitivity matrices to generate an error budget for the system. CPEB also contains a combination of form and ActiveX controls with Visual Basic for Applications code to allow for user interaction in which the user can perform trade studies such as changing engineering requirements, and identifying and isolating stringent requirements. It contains summary tables and graphics that can be instantly used for reporting results in view graphs. The entire process to obtain a coronagraphic telescope performance error budget has been automated into three stages: conversion of optical prescription from Zemax or Code V to MACOS (in-house optical modeling and analysis tool), a linear models process, and an error budget tool process. The first process was improved by developing a MATLAB package based on the Class Constructor Method with a number of user-defined functions that allow the user to modify the MACOS optical prescription. The second process was modified by creating a MATLAB package that contains user-defined functions that automate the process. The user interfaces with the process by utilizing an initialization file where the user defines the parameters of the linear model computations. Other than this, the process is fully automated. The third process was developed based on the Terrestrial Planet Finder coronagraph Error Budget Tool, but was fully automated by using VBA code, form, and ActiveX controls.

Marchen, Luis F.↗

Feasibility of Varying Geo-Fence Around an Unmanned Aircraft Operation Based on Vehicle Performance and Wind

Managing trajectory separation is critical to ensuring accessibility, efficiency, and safety in the unmanned airspace. The notion of geo-fences is an emerging concept, where distance buffers enclose individual trajectories and areas of operation in order to manage the airspace. Currently, the Air Traffic Management system for commercial travel defines static distance buffers around the aircraft; however, commercial UASs are envisioned to operate in significantly closer proximity to other UAS requiring a geo-fence for spacing operations. The geo-fence size can be determined based on vehicle performance characteristics, state of the airspace, weather, and other unforeseen events such as emergency or disaster response. Calculation of the geo-fence size could be determined as part of pre-flight planning and during real-time operations. A largely non-homogeneous fleet of UASs will be operating in low altitude and will likely be commercially developed. Due to intellectual property concerns, the operators may not provide detailed specifications of the control system to UTM. In addition, the huge variety of UAS makes modeling each control system prohibitive and flight data for these vehicles may not exist. Therefore, a generalized, simple geo-fence sizing algorithm must be developed such that it does not rely on detailed knowledge of the vehicle control system, accounts for the presence of urban winds, and is sufficiently accurate. In this work, two simple models are investigated to determine its feasibility as an adequate means for calculating the geo-fence size. The vehicle data used in this work are provided by UAS manufactures who have partnered with NASA's UTM project and some publicly available websites. The first model utilizes wind data processed from the NOAA HRRR (Hourly Rapid Refresh) product and Sonar Annemometer data provided by San Jose State. The second model utilizes OpenFOAM which is a CFD code used to generate a wind field for flow around a single building. The key vehicle performance parameters can include UAS response time to disturbances, command to actuation latency, control system rate limits, time to recovery to desired path, and aerodynamics. It was found that the first model provides an initial understanding of geo-fence sizing, but does not provide enough accuracy to provide UTM with an efficient means of scheduling vehicles. The results of the second model reveal that modeling UAS controls systems with a linearized plant and gain scheduled PID controller does not allow capture the UAS flight dynamics within a significant envelope of the wind disturbances.

uncertainty↗

Evaluation of Transport-Equation-Based Transition Models for High-Speed Boundary Layers Using OVERFLOW

Accurate modeling of laminar-turbulent transition is crucial for the design of hypersonic flight systems. However, the current transition models used in production CFD codes are insufficient for high-speed flows. Many extensions to low-speed models have been suggested; however, a thorough verification and validation effort is needed before these models can be used in design settings. Challenges include potentially missing details of the model implementation requirements and/or a complete specification of the input parameters needed to replicate the test findings. A meaningful assessment of the generalization capability of these models is also hindered by a lack of information regarding the specific flow configurations and associated grids employed for model calibration. As a key first step toward model verification, we present an independent assessment of two recently proposed models for high-speed transition, namely, a model within the SST-$\gamma$ framework and a model based on the SST-$\gamma-\nu_L$ equations. These models are implemented in the NASA OVERFLOW 2.3e solver and their performance in predicting first mode, second mode, and crossflow transition has been evaluated for several test cases in the supersonic and hypersonic regimes. Besides the test cases employed by the model developers, which could have also been used for model calibration, the present assessment includes supplementary configurations that contribute to an unbiased assessment of the models. The outcomes presented in this study indicate the potential for the models to be applied to high-speed flight configurations. Key steps toward future improvements to these models are also outlined.

High-speed flow↗

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗

H2, fixed architecture, control design for large scale systems

The H2, fixed architecture, control problem is a classic linear quadratic Gaussian (LQG) problem whose solution is constrained to be a linear time invariant compensator with a decentralized processing structure. The compensator can be made of p independent subcontrollers, each of which has a fixed order and connects selected sensors to selected actuators. The H2, fixed architecture, control problem allows the design of simplified feedback systems needed to control large scale systems. Its solution becomes more complicated, however, as more constraints are introduced. This work derives the necessary conditions for optimality for the problem and studies their properties. It is found that the filter and control problems couple when the architecture constraints are introduced, and that the different subcontrollers must be coordinated in order to achieve global system performance. The problem requires the simultaneous solution of highly coupled matrix equations. The use of homotopy is investigated as a numerical tool, and its convergence properties studied. It is found that the general constrained problem may have multiple stabilizing solutions, and that these solutions may be local minima or saddle points for the quadratic cost. The nature of the solution is not invariant when the parameters of the system are changed. Bifurcations occur, and a solution may continuously transform into a nonstabilizing compensator. Using a modified homotopy procedure, fixed architecture compensators are derived for models of large flexible structures to help understand the properties of the constrained solutions and compare them to the corresponding unconstrained ones.

Mercadal, Mathieu↗