Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average 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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

The Origin of Titan's External Oxygen: Further Constraints from ALMA Upper Limits on CS and CH2NH

Titan's atmospheric inventory of oxygen compounds (H2O, CO2, CO) are thought to result from photochemistry acting on externally supplied oxygen species (O+, OH, H2O). These species potentially originate from two main sources: (1) cryogenic plumes from the active moon Enceladus and (2) micrometeoroid ablation. Enceladus is already suspected to be the major O+ source, which is required for CO creation. However, photochemical models also require H2O and OH influx to reproduce observed quantities of CO2 and H2O. Here, we exploit sulphur as a tracer to investigate the oxygen source because it has very different relative abundances in micrometeorites (S/O approx. 10(exp -2) and Enceladus' plumes (S/O approx. 10(exp -5). Photochemical models predict most sulphur is converted to CS in the upper atmosphere, so we use Atacama Large Millimeter/submillimeter Array (ALMA) observations at approx. 340 GHz to search for CS emission. We determined stringent CS 3 sigma stratospheric upper limits of 0.0074 ppb (uniform above 100 km) and 0.0256 ppb (uniform above 200 km). These upper limits are not quite stringent enough to distinguish between Enceladus and micrometeorite sources at the 3 sigma level and a contribution from micrometeorites cannot be ruled out, especially if external flux is toward the lower end of current estimates. Only the high flux micrometeorite source model of Hickson et al. can be rejected at 3 sigma. We determined a 3 sigma stratospheric upper limit for CH2NH of 0.35 ppb, which suggests cosmic rays may have a smaller influence in the lower stratosphere than predicted by some photochemical models. Disk-averaged C3H4 and C2H5CN profiles were determined and are consistent with previous ALMA and Cassini/CIRS measurements.

Teanby, N.A.↗
NASA NTRS2018

Numerical Study of Turbulent Separation Bubbles with Varying Pressure Gradient and Reynolds Number

A family of cases each containing a small separation bubble is treated by direct numerical simulation (DNS), varying two parameters: the severity of the pressure gradients, generated by suction and blowing across the opposite boundary, and the Reynolds number. Each flow contains a well-developed entry region with essentially zero pressure gradient, and all are adjusted to have the same value for the momentum thickness, extrapolated from the entry region to the centre of the separation bubble. Combined with fully defined boundary conditions this will make comparisons with other simulations and turbulence models rigorous; we present results for a set of eight Reynolds-averaged Navier–Stokes turbulence models. Even though the largest Reynolds number is approximately 5.5 times higher than in a similar DNS study we presented in 1997, the models have difficulties matching the DNS skin friction very closely even in the zero pressure gradient, which complicates their assessment. In the rest of the domain, the separation location per se is not particularly difficult to predict, and the most definite disagreement between DNS and models is near reattachment. Curiously, the better models tend to cluster together in their predictions of pressure and skin friction even when they deviate from the DNS, although their eddy-viscosity levels are widely different in the outer region near the bubble (or they do not rely on an eddy viscosity). Stratford’s square-root law is satisfied by the velocity profiles, both at separation and reattachment. The Reynolds-number range covers a factor of two, with the Reynolds number based on the extrapolated momentum thickness equal to approximately 1500 and 3000. This allows tentative estimates of the improvements that even higher values will bring to the model comparisons. The solutions are used to assess models through pressure, skin friction and other measures; the flow fields are also used to produce effective eddy-viscosity targets for the models, thus guiding turbulence-modelling work in each region of the flow.

Coleman, G. N.↗
NASA NTRS2021

Spectrophotometric Modeling and Mapping of (101955) Bennu

Using hyperspectral data collected by OVIRS, the visible and infrared spectrometer onboard the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer) spacecraft, we modeled the global average spectrophotometric properties of the carbonaceous asteroid (101955) Bennu and mapped their variations. We restricted our analysis to 0.4–2.5 µm to avoid the wavelengths where thermal emission from the asteroid dominates (>2.5 µm). Bennu has global photometric properties typical of dark asteroids; we found a geometric albedo of 0.046 ± 0.007 and a linear phase slope of 0.024 ± 0.007 mag deg–1 at 0.55 µm. The average spectral slope of Bennu’s normal albedo is –0.0030 µm–1, and the phase reddening parameter is 4.3´10–4 µm–1 deg–1, both over the spectral range of 0.5–2.0 µm. We produced normal albedo maps and phase slope maps at all spectral channels, from which we derived spectral slope and phase reddening maps. Correlation analysis suggests that phase slope variations on Bennu are likely due to photometric roughness variation. A correlation between photometric roughness and thermal roughness is evident, implying that the roughness of Bennu is self-similar on scales from tens of microns to meters. Our analysis reveals latitudinal trends in the spectral color slope and phase reddening on Bennu. The equatorial region appears to be redder than the global average, and the spectral slope decreases towards higher latitudes. Phase reddening on Bennu is relatively weak in the equatorial region and shows an asymmetry between the northern and southern hemispheres. We attributed the latitudinal trend to the geophysical conditions on Bennu that result in a global pattern of mass flow towards the equator.

Jian-Yang Li (李荐扬)↗
NASA NTRS2022

Arizona Water Resources Ii: Utilizing Aerial Imagery and Nasa Earth Observations to Assess Pinyon-Juniper Tree Mortality in Flagstaff, AZ

Pinyon-juniper woodlands (PJW) provide critical and resilient habitat for the local mammal and small bird species of Arizona's northern xeric environment. Drought in Arizona has been persistent for many decades, yet in 2021 PJW experienced a mass tree mortality event at the Wupatki National Monument (WNM) and in other areas across the American Southwest. Previously, researchers at the National Park Service (NPS) and the team from the NASA DEVELOP National Program attempted to quantify the extent of mortality in Northern Arizona between 2015 and 2021 using high resolution National Agricultural Imagery Program (NAIP) aerial photographs. This project aimed to improve the previous term’s methodology and expanded the comparison of the post-mortality event in 2021 to include tree cover assessments for 2017 and 2019. In this iteration, the team utilized NAIP imagery in conjunction with Landscape Fire and Resource Management Planning Tools (LANDFIRE) to calculate the total difference in PJW mortality using an unsupervised classification model trained from multi-date Modified Soil-Adjusted Vegetation Index (MSAVI) and the Visible Atmospherically Resistant Index (VARI) data for the study area. The research also assessed correlations between tree mortality and environmental factors using Western Land Data Assimilation System (WLDAS) modelled climate data. Average PJW mortality from 2015 to 2021 was 21.63% including 19.8% in WNM with the vast majority of dieback occurring between 2019 and 2021. The correlations were weak with the most correlated variables being bare soil evaporation (0.15), rainfall (0.14), groundwater storage (0.13), and wind speed (0.12), perhaps indicating drought as a likely driver of PJW mortality.

pinyon-juniper woodlands↗
NASA NTRSDate not supplied

Role of Arctic Low Clouds in the Arctic Climate System: Isolating the Surface Type Influence on Arctic Low-Clouds

Interactions between sea ice and clouds represent a mechanism through which sea ice influences climate. Understanding the cloud response to the rapidly changing Arctic surface properties is necessary for improving climate simulations and projection of future change inside and outside the Arctic. Towards the goal of understanding the cloud response to sea ice, we composite cloud properties for ice-free, marginal ice zone (MIZ), and ice-covered surfaces during MIZ crossing events to analyze the cloud property differences between surface types. Restricting the analysis to MIZ crossing events enables the isolation of the sea ice effect on clouds from meteorological factors. We find larger cloud fraction and total water concentration below ~1.5 km over ice-free relative to ice-covered surfaces during non-summer months. During summer, the results suggest larger cloud fraction and total water concentration over ice-free surfaces, however differences do not exceed observational uncertainty. Cloud property differences are linked to atmospheric thermodynamic profile differences, namely ice-free surfaces are warmer, moister, less stable, and have more positive surface turbulent fluxes than ice-covered surfaces. Ice-free minus ice-covered cloud property differences scale with surface temperature differences and are only found in the presence of a surface temperature difference. Our results suggest a 0.02 cloud fraction and 0.005 g m 3 total water concentration increase (~5%) at the level of maximum cloud fraction between 2000-2021 due to the observed Arctic sea ice decline in fall, corresponding to ~2 Wm 2 increase in the net surface radiative flux supporting a positive sea ice-cloud radiative feedback in fall and winter and a negative sea ice-cloud radiative feedback in spring. We propose an updated conceptual model where the average surface type influence on cloud properties is mediated by surface temperature differences between ice-free and ice-covered surfaces.

Arctic climate↗
NASA NTRS1996

An efficient approach to ARMA modeling of biological systems with multiple inputs and delays

This paper presents a new approach to AutoRegressive Moving Average (ARMA or ARX) modeling which automatically seeks the best model order to represent investigated linear, time invariant systems using their input/output data. The algorithm seeks the ARMA parameterization which accounts for variability in the output of the system due to input activity and contains the fewest number of parameters required to do so. The unique characteristics of the proposed system identification algorithm are its simplicity and efficiency in handling systems with delays and multiple inputs. We present results of applying the algorithm to simulated data and experimental biological data In addition, a technique for assessing the error associated with the impulse responses calculated from estimated ARMA parameterizations is presented. The mapping from ARMA coefficients to impulse response estimates is nonlinear, which complicates any effort to construct confidence bounds for the obtained impulse responses. Here a method for obtaining a linearization of this mapping is derived, which leads to a simple procedure to approximate the confidence bounds.

Non-NASA Center↗
NASA NTRS2002

A Stochastic Model of Space-Time Variability of Tropical Rainfall: I. Statistics of Spatial Averages

Global maps of rainfall are of great importance in connection with modeling of the earth s climate. Comparison between the maps of rainfall predicted by computer-generated climate models with observation provides a sensitive test for these models. To make such a comparison, one typically needs the total precipitation amount over a large area, which could be hundreds of kilometers in size over extended periods of time of order days or months. This presents a difficult problem since rain varies greatly from place to place as well as in time. Remote sensing methods using ground radar or satellites detect rain over a large area by essentially taking a series of snapshots at infrequent intervals and indirectly deriving the average rain intensity within a collection of pixels , usually several kilometers in size. They measure area average of rain at a particular instant. Rain gauges, on the other hand, record rain accumulation continuously in time but only over a very small area tens of centimeters across, say, the size of a dinner plate. They measure only a time average at a single location. In making use of either method one needs to fill in the gaps in the observation - either the gaps in the area covered or the gaps in time of observation. This involves using statistical models to obtain information about the rain that is missed from what is actually detected. This paper investigates such a statistical model and validates it with rain data collected over the tropical Western Pacific from ship borne radars during TOGA COARE (Tropical Oceans Global Atmosphere Coupled Ocean-Atmosphere Response Experiment). The model incorporates a number of commonly observed features of rain. While rain varies rapidly with location and time, the variability diminishes when averaged over larger areas or longer periods of time. Moreover, rain is patchy in nature - at any instant on the average only a certain fraction of the observed pixels contain rain. The fraction of area covered by rain decreases, as the size of a pixel becomes smaller. This means that within what looks like a patch of rainy area in a coarse resolution view with larger pixel size, one finds clusters of rainy and dry patches when viewed on a finer scale. The model makes definite predictions about how these and other related statistics depend on the pixel size. These predictions were found to agree well with data. In a subsequent second part of the work we plan to test the model with rain gauge data collected during the TRMM (Tropical Rainfall Measuring Mission) ground validation campaign.

Kundu, Prasun K.↗
NASA NTRS1984

Forecasts of the 500 mb height using a dynamically oriented statistical model

The forecast skill of a simple dynamically inspired statistical model of the Northern Hemisphere 500 mb height field is evaluated in spectral and physical space for a variety of forecast lead times (1-32 days) and predictand averaging times (1-32 days). The model includes viscous damping, wave propagation, climatology and implicit stochastic forcing. The largest model skill was found for forecasts of the zonal flow and the largest waves. In general, the largest forecast skills were also associated with the largest forecast error, there being a slight geographic phase shift of the skill with respect to the error. Model skills for climate (time-averaged) forecasts are greater when using instantaneous rather than time averages to forecast time averages. Analysis of model errors suggests areas for improvement in representing forcing terms and model physics. However, the model error fields are largely 'white noise' which suggests that global forecast skills substantially larger than those obtained here are unlikely to be achieved by more sophisticated models.

Roads, J. O.↗
NASA NTRS1985

The zonally averaged transport characteristics of the atmosphere as determined by a general circulation model

Two dimensional modeling has become an established technique for the simulation of the global structure of trace constituents. Such models are simpler to formulate and cheaper to operate than three dimensional general circulation models, while avoiding some of the gross simplifications of one dimensional models. Nevertheless, the parameterization of eddy fluxes required in a 2-D model is not a trivial problem. This fact has apparently led some to interpret the shortcomings of existing 2-D models as indicating that the parameterization procedure is wrong in principle. There are grounds to believe that these shortcomings result primarily from incorrect implementations of the predictions of eddy transport theory and that a properly based parameterization may provide a good basis for atmospheric modeling. The existence of these GCM-derived coefficients affords an unprecedented opportunity to test the validity of the flux-gradient parameterization. To this end, a zonally averaged (2-D) model was developed, using these coefficients in the transport parameterization. Results from this model for a number of contrived tracer experiments were compared with the parent GCM. The generally good agreement substantially validates the flus-gradient parameterization, and thus the basic principle of 2-D modeling.

Plumb, R. A.↗
DOE OSTI2024

Evaluation of FluSight influenza forecasting in the 2021–22 and 2022–23 seasons with a new target laboratory-confirmed influenza hospitalizations

Accurate forecasts can enable more effective public health responses during seasonal influenza epidemics. For the 2021–22 and 2022–23 influenza seasons, 26 forecasting teams provided national and jurisdiction-specific probabilistic predictions of weekly confirmed influenza hospital admissions for one-to-four weeks ahead. Forecast skill is evaluated using the Weighted Interval Score (WIS), relative WIS, and coverage. Six out of 23 models outperform the baseline model across forecast weeks and locations in 2021–22 and 12 out of 18 models in 2022–23. Averaging across all forecast targets, the FluSight ensemble is the 2nd most accurate model measured by WIS in 2021–22 and the 5th most accurate in the 2022–23 season. Forecast skill and 95% coverage for the FluSight ensemble and most component models degrade over longer forecast horizons. In this work we demonstrate that while the FluSight ensemble was a robust predictor, even ensembles face challenges during periods of rapid change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Increasing extreme hourly precipitation risk for New York City after Hurricane Ida

The remnants of Hurricane Ida caused major damage and death in the United States on September 1st, 2021, and 11 people drowned in flooded basement apartments within New York City (NYC). It was catastrophic because the maximum hourly precipitation intensity, recorded as 3.47 inches (88.1 mm) per hour at Central Park, was unprecedentedly high for the NYC region. The stormwater infrastructure in NYC is built for 1.75 inches (44.5 mm) per hour, and so understanding the dynamic risk associated with Ida can inform city planning efforts for climate change’s impact on short duration extreme precipitation events. We contextualize this storm’s record-breaking hourly intensity within the historical record as well as project its risk in the near- to medium-term future using nonstationary stochastic models. These models are conditioned on average temperature (T avg ) and cooling degree day (CDD) projections from three climate models as a covariate, each with a SSP 126 and SSP 370 scenario. The likelihood of such a storm was slowly increasing even before Ida happened, but the projected aggregate reoccurrence risk of an event of Ida’s magnitude over time from the non-stationary models ranges from 4 to 52 times higher than the risk given by the stationary model. Using CDD as a covariate resulted in risks that were more than twice the magnitude than when using T avg . Presenting both covariates provides a broader envelope of uncertainty, which highlights the importance and nuances in the choice of a regionally appropriate covariate for non-stationary risk analysis.

Mossel, Carolien↗
DOE OSTI2023

Reconsidering tympanal-acoustic interactions leads to an improved model of auditory acuity in a parasitoid fly

Although most binaural organisms locate sound sources using neurological structures to amplify the sounds they hear, some animals use mechanically coupled hearing organs instead. One of these animals, the parasitoid fly Ormia ochracea (O. ochracea), has astoundingly accurate sound localization abilities. It can locate objects in the azimuthal plane with a precision of 2°, equal to that of humans, despite an intertympanal distance of only 0.5 mm, which is less than 1/100th of the wavelength of the sound emitted by the crickets that it parasitizes. O. ochracea accomplishes this feat via mechanically coupled tympana that interact with incoming acoustic pressure waves to amplify differences in the signals received at the two ears. In 1995, Miles et al developed a model of hearing mechanics in O. ochracea that represents the tympana as flat, front-facing prosternal membranes, though they lie on a convex surface at an angle from the flies' frontal and transverse planes. The model works well for incoming sound angles less than ±30° but suffers from reduced accuracy (up to 60% error) at higher angles compared to response data acquired from O. ochracea specimens. Despite this limitation, it has been the basis for bio-inspired microphone designs for decades. Here, we present critical improvements to this classic hearing model based on information from three-dimensional reconstructions of O. ochracea's tympanal organ. We identified the orientation of the tympana with respect to a frontal plane and the azimuthal angle segment between the tympana as morphological features essential to the flies' auditory acuity, and hypothesized a differentiated mechanical response to incoming sound on the ipsi- and contralateral sides that depend on these features. We incorporated spatially-varying model coefficients representing this asymmetric response, making a new quasi-two-dimensional (q2D) model. The q2D model has high accuracy (average errors of under 10%) for all incoming sound angles. This improved biomechanical model may inform the design of new microscale directional microphones and other small-scale acoustic sensor systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructural characterization of as-fabricated monolithic plates with boron carbide, aluminum boride, and zirconium boride burnable absorbers

The use of burnable absorbers can be beneficial for nuclear reactors by extending the fuel’s operational cycle, providing additional criticality control, and flattening the power profile. In this work, three burnable absorber materials (boron carbide, aluminum boride, and zirconium boride) embedded in aluminum have been fabricated into foils and clad in AA-6061 for potential use in high performance research reactors. The as-fabricated boron-containing phases were determined using transmission electron microscopy to be AlB2, B4C, and ZrB2. TEM also revealed weak bonding at the B4C-matrix interface. SEM showed a relatively uniform spatial distribution of boron-containing phases for all the candidate materials. Higher porosity was observed in the foil containing ZrB2 in its as-rolled condition. The porosity in the ZrB2 foil was reduced by hot isostatic pressing. The size and shape distributions of the boron-containing phases were analyzed on the criteria of cross-sectional area, perimeter, roundness, circularity, and aspect ratio. A method of converting the 2D burnable absorber dispersoids seen in cross-sectional microscopy images into 3D volumes was derived using both spherical and ellipsoidal geometry models. The difference in calculated burnable absorber dispersoid average volume between the two models ranges from 20% to 100%, which could impact burnable absorber burnout rates due to differences in neutron self-shielding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Structural properties of aqueous grown polydopamine thin films determined by neutron reflectometry

Here, in this work, neutron reflectometry (NR) studies of the bio-mimetic polymer, polydopamine (PDA), deposited from differing initial concentrations of precursor material dopamine hydrochloride for a range of polymerization times, is reported. PDA can form a complex and highly versatile polymer film, with many structure studies having been performed previously, but a comprehensive structural determination of PDA by NR is lacking in the literature. It was found that simple box models were incapable of fully explaining the observed data, necessitating the use of a composite model consisting of the weighted average of both the 1 and 2 box models to fully capture the heterogeneous nature of the PDA film structure. Confocal laser scanning microscopy (CLSM) and atomic force microscopy (AFM) were performed to capture the surface structure and relative mechanical difference between the separate domains of the PDA. The CLSM results provide evidence that the PDA domains are larger than the coherent scattering length of a neutron used in these measurements, 10μm, supporting the need for a composite model. The AFM measurements show complex structure below the coherence length provide physical justification for the two box model. It was determined that film structure and quality are both heavily impacted by the initial concentration of dopamine hydrochloride and the polymerization time, giving confidence to the highly customizable nature of PDA as an adhesion promoting interface treatment in composite systems, such as plastic bonded explosives (PBX).

36 MATERIALS SCIENCE↗
DOE OSTI2024

On Building Predictive Digital Twin Incorporating Wave Predicting Capabilities: Case Study on UMaine Experimental Campaign - FOCAL

The response of floating wind turbines (FWT) are susceptible to stochastic wave variations. For the optimal operation of FWT, a comprehensive understanding of the phaseresolved wave dynamics and the consequential system response is crucial for real-time monitoring and control. A multi-variate, multi-step, long short term memory (MLSTM), a type of recurrent neural network (RNN) is used to capture complex system dynamics for real-time application. Results indicate that the integration of a wave prediction-reconstruction (WRP) model substantially enhances prediction accuracy by 50% on average relative to the baseline model. The improvement is consistent across various wave extremity and prediction horizons, thereby significantly broadening the scope for timely and precise predictive capabilities.

17 WIND ENERGY↗
DOE OSTI2023

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

CFD↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Modeling, Monitoring and Fault Diagnosis of Spacecraft Air Contaminants

Progress and results in the development of an integrated air quality modeling, monitoring, fault detection, and isolation system are presented. The focus was on development of distributed models of the air contaminants transport, the study of air quality monitoring techniques based on the model of transport process and on-line contaminant concentration measurements, and sensor placement. Different approaches to the modeling of spacecraft air contamination are discussed, and a three-dimensional distributed parameter air contaminant dispersion model applicable to both laminar and turbulent transport is proposed. A two-dimensional approximation of a full scale transport model is also proposed based on the spatial averaging of the three dimensional model over the least important space coordinate. A computer implementation of the transport model is considered and a detailed development of two- and three-dimensional models illustrated by contaminant transport simulation results is presented. The use of a well established Kalman filtering approach is suggested as a method for generating on-line contaminant concentration estimates based on both real time measurements and the model of contaminant transport process. It is shown that high computational requirements of the traditional Kalman filter can render difficult its real-time implementation for high-dimensional transport model and a novel implicit Kalman filtering algorithm is proposed which is shown to lead to an order of magnitude faster computer implementation in the case of air quality monitoring.

Ramirez, W. Fred↗