SEARCH · Engineering Papers

Results for “predictive state estimation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS1996

X-29A Lateral-Directional Stability and Control Derivatives Extracted From High-Angle-of-Attack Flight Data

The lateral-directional stability and control derivatives of the X-29A number 2 are extracted from flight data over an angle-of-attack range of 4 degrees to 53 degrees using a parameter identification algorithm. The algorithm uses the linearized aircraft equations of motion and a maximum likelihood estimator in the presence of state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics over the aircraft at angles of attack above 15 degrees. The results supported the flight-envelope-expansion phase of the X-29A number 2 by helping to update the aerodynamic mathematical model, to improve the real-time simulator, and to revise flight control system laws. Effects of the aircraft high gain flight control system on maneuver quality and the estimated derivatives are also discussed. The derivatives are plotted as functions of angle of attack and compared with the predicted aerodynamic database. Agreement between predicted and flight values is quite good for some derivatives such as the lateral force due to sideslip, the lateral force due to rudder deflection, and the rolling moment due to roll rate. The results also show significant differences in several important derivatives such as the rolling moment due to sideslip, the yawing moment due to sideslip, the yawing moment due to aileron deflection, and the yawing moment due to rudder deflection.

Iliff, Kenneth W.↗

NASA NTRS2003

Validation and Error Characterization for the Global Precipitation Measurement

The Global Precipitation Measurement (GPM) is an international effort to increase scientific knowledge on the global water cycle with specific goals of improving the understanding and the predictions of climate, weather, and hydrology. These goals will be achieved through several satellites specifically dedicated to GPM along with the integration of numerous meteorological satellite data streams from international and domestic partners. The GPM effort is led by the National Aeronautics and Space Administration (NASA) of the United States and the National Space Development Agency (NASDA) of Japan. In addition to the spaceborne assets, international and domestic partners will provide ground-based resources for validating the satellite observations and retrievals. This paper describes the validation effort of Global Precipitation Measurement to provide quantitative estimates on the errors of the GPM satellite retrievals. The GPM validation approach will build upon the research experience of the Tropical Rainfall Measuring Mission (TRMM) retrieval comparisons and its validation program. The GPM ground validation program will employ instrumentation, physical infrastructure, and research capabilities at Supersites located in important meteorological regimes of the globe. NASA will provide two Supersites, one in a tropical oceanic and the other in a mid-latitude continental regime. GPM international partners will provide Supersites for other important regimes. Those objectives or regimes not addressed by Supersites will be covered through focused field experiments. This paper describes the specific errors that GPM ground validation will address, quantify, and relate to the GPM satellite physical retrievals. GPM will attempt to identify the source of errors within retrievals including those of instrument calibration, retrieval physical assumptions, and algorithm applicability. With the identification of error sources, improvements will be made to the respective calibration, assumption, or algorithm. The instrumentation and techniques of the Supersites will be discussed. The GPM core satellite, with its dual-frequency radar and conically scanning radiometer, will provide insight into precipitation drop-size distributions and potentially increased measurement capabilities of light rain and snowfall. The ground validation program will include instrumentation and techniques commensurate with these new measurement capabilities.

Bidwell, Steven W.↗

NASA NTRS2006

Current Sheet Evolution In The Aftermath Of A CME Event

We report on SOHO UVCS observations of the coronal restructuring following a coronal mass ejection (CME) on 2002 November 26, at the time of a SOHO-Ulysses quadrature campaign. Starting about 1.5 hr after a CME in the northwest quadrant, UVCS began taking spectra at 1.7 R, covering emission from both cool and hot plasma. Observations continued, with occasional gaps, for more than 2 days. Emission in the 974.8 A line of [Fe XVIII], indicating temperatures above 6 x 10(exp 6) K, was observed throughout the campaign in a spatially limited location. Comparison with EIT images shows the [Fe XVIII] emission to overlie a growing post-flare loop system formed in the aftermath of the CME. The emission most likely originates in a current sheet overlying the arcade. Analysis of the [Fe XVIII] emission allows us to infer the evolution of physical parameters in the current sheet over the entire span of our observations: in particular, we give the temperature versus time in the current sheet and estimate its density. At the time of the quadrature, Ulysses was directly above the location of the CME and intercepted the ejecta. High ionization state Fe was detected by the Ulysses SWICS throughout the magnetic cloud associated with the CME, although its rapid temporal variation suggests bursty, rather than smooth, reconnection in the coronal current sheet. The SOHO-Ulysses data set provided us with the unique opportunity of analyzing a current sheet structure from its lowest coronal levels out to its in situ properties. Both the remote and in situ observations are compared with predictions of theoretical CME models.

Bemporad, A.↗

NASA NTRSDate not supplied

Comparison of Entry Descent and Landing Aerodynamic Databases with Uncertainty Quantification Developed Using Machine Learning Techniques

When developing the aerodynamic databases for use in trajectory simulations, it is important to develop a system of metrics to qualify which aerodynamic models are best to use. Since aerodynamics are just one input into trajectory simulations, the results of these simulations do not reflect on the quality of the aerodynamic database used. This means that aerodynamic database comparisons must be done offline. While traditional metrics that focus on mean/nominal predictions are a good first step, more robust estimates of the prediction interval become important as more focused uncertainty models are developed. We explore the limitations of evaluating aerodynamic models based purely on nominal-centered response surfaces. Before elaborating and evaluating metrics based on distributed models, the value of evaluating prediction interval and confidence interval are discussed to conclude that prediction intervals are more relevant to the use of trajectory analysis. Several metrics to evaluate the prediction interval are introduced with a focus on the standard calibration metric. Finally, we compare candidate models using both mean and distributed metrics. A finalized candidate model developed using state of the art machine learning methods is compared to a baseline model developed using traditional aerodynamic database modeling techniques.

Aerodynamic Database↗

NASA NTRSDate not supplied

James Webb Space Telescope Navigation Optimization Challenges

The James Webb Space Telescope (JWST) is a NASA flagship mission launched on December 25, 2021. The early orbit phase was highlighted by three midcourse correction burns that were performed to maneuver the vehicle into a libration orbit at the second Sun-Earth-Moon (SEM) libration point (L2). During the coast to L2, the vehicle’s primary observatory and sunshield were deployed. Once the tennis-court-sized sunshield was unfurled and tensioned into place, the large area exposed to solar radiation pressure (SRP) dictated that the SRP force model would need to account for the vehicle’s geometry, reflective properties, and orientation relative to the Sun. Following the insertion of JWST into its L2 libration orbit on January 24, 2022, the commissioning phase for the vehicle’s observatory commenced along with the first cycle of routine station-keeping maneuvers. The NASA Goddard Space Flight Center’s Flight Dynamics Facility (FDF) provides navigation services for the JWST mission. The FDF serves as the prime or backup navigation operations center for more than 30 active missions spanning a wide array of flight regimes performing orbit determination, maneuver planning, trajectory optimization, and tracking data evaluation. Definitive orbit determination for JWST is performed by the FDF using an Extended Kalman Filter (EKF) to estimate the vehicle’s trajectory and reflective properties using Deep Space Network (DSN) TRK-2-34 tracking data and spacecraft attitude telemetry. For the purposes of orbit prediction and station-keeping maneuver targeting, the sensitivity of the SRP force to the vehicle’s orientation requires that predictive attitude information must be incorporated to generate accurate orbit predictions and maneuver plans. Predictive attitude plans are provided to the FDF by the JWST Spacecraft Operations Center (SOC). Short-term (28 day) orbit predictions are propagated using a Short Range Attitude Plan (SRAP) to model the vehicle’s future attitude states for up to a week. SRAP attitude data is highly reliable, as it reflects finalized planned attitude states that the vehicle will be commanded to attain during the ensuing week. A Long Range Attitude Plan (LRAP) can be implemented to model predicted attitude states beyond one week, but these predictions are tentative and subject to revision due to changes in science operations plans. Instead, a conservative approach of applying a Sun-Pointing Neutral (SPN) attitude configuration is utilized for long-term (2 year) orbit predictions where reliable planned attitude data is not available. SPN attitude mode aligns the net SRP force along the JWST-to-Sun vector where it becomes independent of the vehicle’s angle of rotation about this vector, defined as the Sun yaw. Given the dynamic instability of the libration orbit, station-keeping thrust must be applied in either the sunward or anti-sunward direction to balance the resulting orbit. The attitude constraints of the vehicle also impose limits on the available pointing directions for station-keeping thrusters. The thrusters cannot be aligned with the optimal pointing direction for sunward maneuvers, rendering sunward maneuvers to be significantly less fuel efficient than anti-sunward maneuvers. Consequently, station-keeping maneuvers must be targeted to balance the orbit while ensuring that the next maneuver will also be performed in the anti-sunward direction to optimize propellant usage. When targeting a station-keeping maneuver, the attitude and SRP modeling configuration which is applied to the predicted post-maneuver trajectory will dictate the direction of the subsequent station-keeping maneuver, assuming nominal propulsion system performance. If the predicted post-maneuver attitude states result in an SRP model which under-predicts the cumulative SRP impact on the post-maneuver orbit, a targeted anti-sunward maneuver will be larger than necessary and the next maneuver will need to be executed sunward in order to compensate. In contrast, an anti-sunward maneuver which is targeted using an SRP model which over-predicts the cumulative SRP impact will achieve station-keeping while ensuring that the next maneuver will likewise be performed anti-sunward. For this reason, station-keeping maneuvers are targeted while applying the SPN attitude mode, as this mode entails the largest possible SRP area cross-section and therefore a larger modeled cumulative SRP impact. This paper documents the NASA Goddard Space Flight Center’s FDF support for JWST on-orbit operations and the analysis projects undertaken utilizing the experiences and data accumulated throughout the first full year of routine science operations. The results of these analysis efforts have been used to implement improvements to orbit prediction accuracy and maneuver efficiency which have the potential to prolong the lifespan of JWST to continue to conduct ground-breaking infrared astronomy for decades to come.

Flight Dynamics↗

NASA NTRS1975

Attitude determination of a high altitude balloon system. Part 2: Development of the parameter determination process

The attitude of the balloon system is determined as a function of time if: (a) a method for simulating the motion of the system is available, and (b) the initial state is known. The initial state is obtained by fitting the system motion (as measured by sensors) to the corresponding output predicted by the mathematical model. In the case of the LACATE experiment the sensors consisted of three orthogonally oriented rate gyros and a magnetometer all mounted on the research platform. The initial state was obtained by fitting the angular velocity components measured with the gyros to the corresponding values obtained from the solution of the math model. A block diagram illustrating the attitude determination process employed for the LACATE experiment is shown. The process consists of three essential parts; a process for simulating the balloon system, an instrumentation system for measuring the output, and a parameter estimation process for systematically and efficiently solving the initial state. Results are presented and discussed.

Nigro, N. J.↗

NASA NTRS1986

An economical semi-analytical orbit theory for micro-computer applications

An economical algorithm is presented for predicting the position of a satellite perturbed by drag and zonal harmonics J2 through J4. Simplicity being of the essence, drag is modeled as a secular decay rate in the semimajor axis (retarded motion) with the zonal perturbations modeled from a modified version of Brouwers formulas. The algorithm is developed as an alternative on-board orbit predictor; a back up propagator requiring low energy consumption; or a ground based propagator for microcomputer applications (e.g., at the foot of an antenna). An O(J2) secular retarded state partial matrix (matrizant) is also given to employ with state estimation. The theory has been implemented in BASIC on an inexpensive microcomputer, the program occupying under 8K bytes of memory. Simulated trajectory data and real tracking data are employed to illustrate the theory's ability to accurately accommodate oblateness and drag effects.

Gordon, R. A.↗

NASA NTRS2000

Constituent-Based Life Models Being Developed for SiC/SiC Composites

For the successful utilization of ceramic matrix composites (CMC) as hot-section components in advanced aeropropulsion engines, the CMC constituents will need to be tailored and optimized to meet all the critical property demands of each component. Under the High-Speed Research (HSR) and Advanced High-Temperature Engine Materials Technology (HITEMP) Programs, the NASA Glenn Research Center at Lewis Field initiated research to develop mechanistic models for key CMC thermostructural properties. These models would describe the effects of different constituent factors (composition, geometry, and volume fraction) and of potential application conditions (stress, time, temperature, and environment) on these properties. Particular focus was placed on both analytical and numerical modeling of state-of-the-art SiC/SiC composites where the primary load-bearing constituents are stoichiometric SiC fibers in a complex multiphase SiC matrix produced by chemical vapor infiltration and melt infiltration. Recent studies have resulted in computer-based numerical models for the elastic modulus, thermal expansion, and thermal conductivity properties of the SiC/SiC system. Additional studies have generated analytical and empirical models for the time dependence of composite rupture strength at temperatures above 2200 F (1200 C), where CMC's have an important thermostructural advantage over current nickel-based superalloys. These life models utilize thermal activation theory and fiber stress-rupture results measured at Glenn to generate Larson-Miller (L M) plots of fiber rupture strength versus q, a single time- and temperature-dependent parameter. Assuming a worse case in which the SiC matrix is cracked, rupture is then controlled by the time-dependent fracture characteristics of the fiber bundles bridging the matrix cracks. With this as the controlling mechanism, one can then use simple composite theory and the fiber L M plots to predict CMC rupture strength versus the q parameter. The dashed line shows the predicted rupture strength of a SiC/SiC composite that is reinforced by a state-of-the-art stoichiometric SiC fiber. For the q parameter, time is in hours and temperature in degrees kelvin. To generate these predictions, a two-dimensional 0 /90 composite with approximately 16 percent fiber in the applied stress direction and an air test environment were assumed. As such, it is possible to compare the model predictions against limited stress-rupture data for this CMC as shown by the data points. The good agreement confirms the rupture model at least for the selected CMC and test conditions. Thus for this particular SiC/SiC composite, one can estimate a 1000-hr rupture strength of approximately 12 ksi at 2400 F (T = 1588 K and q = 39 700 K). At lower CMC application stresses, the SiC matrix is typically uncracked, so both the fiber and matrix constituents share the composite load. In this case, CMC rupture is controlled by the constituent with the longest rupture time based on the creep rate of the composite. Measured Monkman-Grant plots of rupture time versus creep rate for the two SiC constituents have been used to develop CMC life models for this important application condition. NASA and DOD are currently using this information to establish application and material goals for more advanced CMC's that can be used at even higher temperatures.

DiCarlo, James A.↗

NASA NTRS1980

Oxidation state of the atmosphere and crust of Venus from Pioneer Venus results

Recent spacecraft observations of Venus permit a detailed model of sulfur chemistry in the atmosphere-lithosphere system. Pioneer Venus experiments confirm that, as predicted, COS and H2S are dominant over SO2 in the lower atmosphere, and that the equilibrium concentrations of S2 and S3 are significant. Many criteria serve to bracket the oxidation state of the crust: it is nearly certain that the S2(2-)/SO4(2-) buffer regulates the oxygen fugacity, and that FeO is at least as abundant as F2O3 in crustal silicates. A highly oxidized crust (as, for example, would result from O2 absorption complementary to escape of vast amounts of H2) is incompatible with the gas-phase sulfur chemistry. If the Pioneer Venus mass spectrometer estimates of the abundance of sulfur gases are correct, earth-like models for the bulk composition of Venus are seriously in error, and a far lower FeO content is required for Venus.

Lewis, J. S.↗

NASA NTRS1998

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Elishakoff, Isaac↗

NASA NTRS1995

An Extended Ab Initio and Theoretical Thermodynamics Studies of the Bergman Reaction and the Energy Splitting of the Singlet Ortho-, Meta-, and Para-Benzynes

The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Lindh, Roland↗

NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasizes the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M qi) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in turns of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, Derck↗

NASA NTRS2006

An Intercomparison of Microphysical Retrieval Algorithms for Upper Tropospheric Ice Clouds

The large horizontal extent, location in the cold upper troposphere, and ice composition make cirrus clouds important modulators of the earth's radiation budget and climate. Cirrus cloud microphysical properties are difficult to measure and model because they are inhomogeneous in nature and their ice crystal size distribution and habit are not well characterized. Accurate retrievals of cloud properties are crucial for improving the representation of cloud scale processes in large-scale models and for accurately predicting the earth's future climate. A number of passive and active remote sensing retrieval algorithms exist for estimating the microphysical properties of upper tropospheric clouds. We believe significant progress has been made in the evolution of these retrieval algorithms in the last decade, however, there is room for improvement. Members of the Atmospheric Radiation measurement program (ARM) Cloud properties Working Group are involved in an intercomparison of optical depth(tau), ice water path, and characteristic particle size in clouds retrieved using ground-based instruments. The goals of this intercomparison are to evaluate the accuracy of state-of-the-art algorithms, quantify the uncertainties, and make recommendations for improvement.

Comstock, Jennifer M.↗

NASA NTRSDate not supplied

High-Enthalpy Testing to Validate Simulation of an Aerosol Capture Probe

AERACEPT (Aerosol Rapid Analysis Combined Entry Probe/sonde Technology) is a developing technology enabling in-situ aerosol particle sampling and analysis in a small spacecraft mission envelope. It integrates a passive aerosol sample collection system into a probe’s thermal protection system (TPS) to remove the need for heat shield separation and active descent control (parachutes, gliders, etc.). AERACEPT is being validated against the requirements of the Nephele mission concept, which targets the middle and lower Venus cloud layers. The proposed technology employs 3D Carbon-Carbon (3D-CC) at the probe’s nose to withstand the extreme environments encountered in a Venus entry without producing pyrolysis gases that would contaminate the mission sample collection. Any shape change of this 3D-CC sample inlet will affect airflow through the sampling system, thereby impacting aerosol sample collection bias and efficiency. Accurate prediction of the nose’s material response during entry is there-fore of fundamental importance to the development and practical use of this technology. As a result, the AERACEPT project is planning a test campaign in the newly built PlasmatronX facility with a novel “Open Iso-Q” test article to validate state-of-the-art material response tools for a simulated Venus entry. The PlasmatronX is a 350 kW inductively-coupled plasma facility developed and run by the Center for Hypersonics and Entry Systems Studies at the University of Illinois at Urbana-Champaign. The facility can support ground testing for a variety of planetary destinations by simulating entry conditions in Nitrogen, Air, and Carbon Dioxide. To best simulate the Venus atmosphere and match flight-like recession, the planned AERACEPT test campaign will use Carbon Dioxide as its test gas. AERACEPT has designed two test articles for this campaign: 1) a standard test article whose curved surface approximates a constant applied heat flux (“Closed Iso-Q” model) and 2) a novel inlet test article with a through-hole at the stagnation point (“Open Iso-Q” model). Both article types consist of a 3D-CC sample bonded to a graphite fixture. All articles will be laser scanned before and after testing to estimate shape change and recession, and the “Closed Iso-Q” articles will be instrumented with thermocouples to provide temperature histories. The test article temperatures and inlet shape change will be compared against simulations run with the Porous material Analysis Toolbox based on OpenFOAM (PATO) to validate the AERACEPT 3D-CC material response model. Due to the high thermal conductivity of 3D-CC, multidimensional effects have a large influence on the test article temperature and the standard 1-D material response tools are not sufficient for test planning. The test campaign features (e.g. article design, run duration, heat fluxes) are therefore informed by a series of 3-D simulations using state-of-the-art material response tools such as PATO. The test campaign is scheduled for spring of 2024. Expected results include 1) Key drivers of selected test conditions (e.g. matching mission flight atmosphere & recession), 2) Simulation results that informed test article design, 3) Acquired test data such as recession measurements, temperature histories, and photos, and 4) Comparisons between test results and simulation predictions.

AERACEPT↗

NASA NTRSDate not supplied

High-Enthalpy Testing to Validate Simulation of an Aerosol Capture Probe

AERACEPT↗

NASA NTRS1986

Sequential estimation and satellite data assimilation in meteorology and oceanography

The central theme of this review article is the role that dynamics plays in estimating the state of the atmosphere and of the ocean from incomplete and noisy data. Objective analysis and inverse methods represent an attempt at relying mostly on the data and minimizing the role of dynamics in the estimation. Four-dimensional data assimilation tries to balance properly the roles of dynamical and observational information. Sequential estimation is presented as the proper framework for understanding this balance, and the Kalman filter as the ideal, optimal procedure for data assimilation. The optimal filter computes forecast error covariances of a given atmospheric or oceanic model exactly, and hence data assimilation should be closely connected with predictability studies. This connection is described, and consequences drawn for currently active areas of the atmospheric and oceanic sciences, namely, mesoscale meteorology, medium and long-range forecasting, and upper-ocean dynamics.

Ghil, M.↗

NASA NTRS2000

Controlling Microbial Byproducts using Model-Based Substrate Monitoring and Control Strategies

We have developed a computer-controlled bioreactor system to study various aspects of microbially-mediated nitrogen cycling. The system has been used to investigate methods for controlling microbial denitrification (the dissimilatory reduction of nitrate to N2O and N2) in hydroponic plant growth chambers. Such chambers are key elements of advanced life support systems being designed for use on long duration space missions, but nitrogen use efficiency in them is reduced by denitrification. Control software architecture was designed which permits the heterogeneous control of system hardware using traditional feedback control, and quantitative and qualitative models of various system features. Model-based feed forward control entails prediction of future systems in states and automated regulation of system parameters to achieve desired and avoid undesirable system states. A bacterial growth rate model based on the classic Monod model of saturation kinetics was used to evaluate the response of several individual denitrifying species to varying environmental conditions. The system and models are now being applied to mixed microbial communities harvested from the root zone of a hydroponic growth chamber. The use of a modified Monod organism interaction model was evaluated as a means of achieving more accurate description of the dynamic behavior of the communities. A minimum variance parameter estimation routine was also' used to calibrate the constant parameters in the model by iterative evaluation of substrate (nitrate) uptake and growth kinetics. This representation of processes and interactions aids in the formulation of control laws. The feed forward control strategy being developed will increase system autonomy, reduce crew intervention and limit the accumulation of undesirable waste products (NOx).

Smernoff, David T.↗

NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasises the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters, but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M (dot) q(sub i)) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, D.↗