Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling and filtering”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A stochastic covariance shrinkage approach to particle rejuvenation in the ensemble transform particle filter

Rejuvenation in particle filters is necessary to prevent the collapse of the weights when the number of particles is insufficient to properly sample the high-probability regions of the state space. Rejuvenation is often implemented in a heuristic manner by the addition of random noise that widens the support of the ensemble. This work aims at improving canonical rejuvenation methodology by the introduction of additional prior information obtained from climatological samples; the dynamical particles used for importance sampling are augmented with samples obtained from stochastic covariance shrinkage. A localized variant of the proposed method is developed. Numerical experiments with the Lorenz '63 model show that modified filters significantly improve the analyses for low dynamical ensemble sizes. Furthermore, localization experiments with the Lorenz '96 model show that the proposed methodology is extendable to larger systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Estimation With Range Dependent Sensor Model

This paper focuses on the improvement of target tracking accuracy taking into account the sensor’s reading degradation as the relative range increases. The approach is based on the assumption that range measurement error depends on the actual range. Specifically, we model the measurement error as a proportional to the actual range term plus a zero-mean, Gaussian distributed, and uncorrelated process. We present two approaches to tracking a target with constant velocity model. The first approach uses Extended Kalman Filter (EKF) framework to estimate target's states and the unknown proportionality coefficient with linearization on the observation model on the predicted states at each time step. The second one uses a measurement conversion method to estimate the target's scaled by the unknown constant term states, which is shown to be unbiased. This conversion results in linear state and observation models, hence the standard Kalman filter can be applied. The actual target state is computed by application of back scaling with the estimate of unknown scale factor. We evaluate the approaches in desktop simulations.

Target tracking↗
NASA NTRS1995

Assimilation of TOPEX Sea Level Measurements with a Reduced-Gravity, Shallow Water Model of the Tropical Pacific Ocean

Sea surface height variability measured by TOPEX is analyzed in the tropical Pacific Ocean by way of assimilation into a wind-driven, reduced-gravity, shallow water model using an approximate Kalman filter and smoother. The analysis results in an optimal fit of the dynamic model to the observations, providing it dynamically consistent interpolation of sea level and estimation of the circulation. Nearly 80% of the expected signal variance is accounted for by the model within 20 deg of the equator, and estimation uncertainty is substantially reduced by the voluminous observation. Notable features resolved by the analysis include seasonal changes associated with the North Equatorial Countercurrent and equatorial Kelvin and Rossby waves. Significant discrepancies are also found between the estimate and TOPEX measurements, especially near the eastern boundary. Improvements in the estimate made by the assimilation are validated by comparisons with independent tide gauge and current meter observations. The employed filter and smoother are based on approximately computed estimation error covariance matrices, utilizing a spatial transformation and an symptotic approximation. The analysis demonstrates the practical utility of a quasi-optimal filter and smoother.

Fukumori, Ichiro↗
NASA NTRS1991

Integrated adaptive filtering and design for control experiments of flexible structures

A novel method is presented of identifying a state space model and a state estimator for linear stochastic systems from input and output data. The method is primarily based on the relations between the state space model and the finite difference model for linear stochastic systems derived through projection filters. It is proven that least squares identification of a finite difference model converges to the model derived from the projection filters. System pulse response samples are computed from the coefficients of the finite difference model. In estimating the corresponding state estimator gain, a z-domain method is used. First the deterministic component of the output is subtracted out, and then the state estimator gain is obtained by whitening the remaining signal. Experimental example is used to illustrate the feasibility of the method.

Huang, Jen-Kuang↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1997

System Dynamic Analysis of a Wind Tunnel Model with Applications to Improve Aerodynamic Data Quality

The research investigates the effect of wind tunnel model system dynamics on measured aerodynamic data. During wind tunnel tests designed to obtain lift and drag data, the required aerodynamic measurements are the steady-state balance forces and moments, pressures, and model attitude. However, the wind tunnel model system can be subjected to unsteady aerodynamic and inertial loads which result in oscillatory translations and angular rotations. The steady-state force balance and inertial model attitude measurements are obtained by filtering and averaging data taken during conditions of high model vibrations. The main goals of this research are to characterize the effects of model system dynamics on the measured steady-state aerodynamic data and develop a correction technique to compensate for dynamically induced errors. Equations of motion are formulated for the dynamic response of the model system subjected to arbitrary aerodynamic and inertial inputs. The resulting modal model is examined to study the effects of the model system dynamic response on the aerodynamic data. In particular, the equations of motion are used to describe the effect of dynamics on the inertial model attitude, or angle of attack, measurement system that is used routinely at the NASA Langley Research Center and other wind tunnel facilities throughout the world. This activity was prompted by the inertial model attitude sensor response observed during high levels of model vibration while testing in the National Transonic Facility at the NASA Langley Research Center. The inertial attitude sensor cannot distinguish between the gravitational acceleration and centrifugal accelerations associated with wind tunnel model system vibration, which results in a model attitude measurement bias error. Bias errors over an order of magnitude greater than the required device accuracy were found in the inertial model attitude measurements during dynamic testing of two model systems. Based on a theoretical modal approach, a method using measured vibration amplitudes and measured or calculated modal characteristics of the model system is developed to correct for dynamic bias errors in the model attitude measurements. The correction method is verified through dynamic response tests on two model systems and actual wind tunnel test data.

Buehrle, Ralph David↗
NASA NTRS2019

Global Evolution of Solar Magnetic Fields and Prediction of Solar Activity Cycles

Prediction of solar activity cycles is challenging because the physical processes inside the Sun involve a broad range of multiscale dynamics that no model can reproduce, and the available observations are highly limited and cover mostly surface layers. Helioseismology makes it possible to probe solar dynamics in the convective zone, but variations in the differential rotation and meridional circulation are currently available for only two solar activity cycles. It has been demonstrated that sunspot observations, which cover over 400 years, can be used to calibrate the Parker-Kleeorin-Ruzmaikin model and that the Ensemble Kalman Filter (EnKF) method can be used to link the model magnetic fields to sunspot observations to make reliable predictions of a following cycle. However, for more accurate predictions, it is necessary to use actual observations of the solar magnetic fields, which are available for only four solar cycles. This raises the question of how limitations in observational data and model uncertainties affect predictive capabilities and implies the need for the development of new forecast methodologies and validation criteria. In this presentation, I will discuss the influence of the limited number of available observations on the accuracy of EnKF estimates of solar cycle parameters.

Kitiashvili, Irina N.↗
NASA NTRS1999

Fidelity Analysis of Sampled Imaging Systems

Many modeling, simulation and performance analysis studies of sampled imaging systems are inherently incomplete because they are conditioned on a discrete-input, discrete-output model that only accounts for blurring during image acquisition and additive noise. For those sampled imaging systems where the effects of digital image acquisition, digital filtering and reconstruction are significant, the modeling, simulation and performance analysis should be based on a more comprehensive continuous-input, discrete-processing, continuous-output end-to-end model. This more comprehensive model should properly account for the low-pass filtering effects of image acquisition prior to sampling, the potentially important noiselike effects of the aliasing caused by sampling, additive noise due to device electronics and quantization, the generally high-boost filtering effects of digital processing, and the low-pass filtering effects of image reconstruction. This model should not, however, be so complex as to preclude significant mathematical analysis, particularly the mean-square (fidelity) type of analysis so common in linear system theory. We demonstrate that, although the mathematics of such a model is more complex, the increase in complexity is not so great as to prevent a complete fidelity-metric analysis at both the component level and at the end-to-end system level: that is, computable mean-square-based fidelity metrics are developed by which both component-level and system-level performance can be quantified. In addition, we demonstrate that system performance can be assessed qualitatively by visualizing the output image as the sum of three component images, each of which relates to a corresponding fidelity metric. The cascaded, or filtered, component accounts for the end-to-end system filtering of image acquisition, digital processing, and image reconstruction; the random noise component accounts for additive random noise, modulated by digital processing and image reconstruction filtering; and the aliased noise component accounts for the frequency folding effect of sampling, modulated by digital processing and image reconstruction filtering.

Park, Stephen K.↗
DOE OSTI2022

Design of a supervisory control system for autonomous operation of advanced reactors

Advanced reactors to be deployed in the coming decades will face deregulated energy markets, and may adopt flexible operation to boost profitability. To aid in the transition from baseload to flexible operation paradigm, autonomous operation is sought. This work focuses on the control aspect of autonomous operation. Specifically, a hierarchical control system is designed to support constraint enforcement during routine operational transients. Within the system, data-driven modeling, physics-based state observation, and classical control algorithms are integrated to provide an adaptable and robust solution. A 320 MW Fluoride-cooled High-temperature Pebble-bed Reactor is the design basis for demonstrating the proposed control system. The hierarchical control system consists of a supervisory layer and low-level layer. The supervisory layer receives requests to change the system's operating conditions (e.g., the current reactor power to meet a load -follow), and accepts or rejects them based on constraints that have been assigned. Constraints are issued to keep the plant within an optimal operating region. The low-level layer interfaces with the actuators of the system to fulfill requested changes, while maintaining tracking and regulation duties. Further, to accept requests at the supervisory layer, the Reference Governor algorithm was adopted. To model the dynamics of the reactor, a system identification algorithm, Dynamic Mode Decomposition, was utilized. To estimate the evolution of process variables that cannot be directly measured (e.g., the propagation of delayed neutron precursors), the Unscented Kalman Filter, incorporating a nonlinear model of nuclear dynamics, was adopted. The composition of these algorithms led to a numerical demonstration of constraint enforcement during a 40% power drop transient (at a rate of 5 %/min). Uncontrolled secondary-side temperatures were successfully constrained. Adaptability of the proposed system was demonstrated by modifying the constraint values, and enforcing them during the transient. Robustness was also demonstrated by enforcing constraints under noisy environments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1985

Geophysical Interpretation of Venus Gravity Data

The subsurface distribution of Venus was investigated through the analysis of the data from Pioneer Venus Orbiter (PVO). In particular, the Doppler tracking data were used to map the gravitational potential. These were compared to the topographic data from the PVO radar (ORAD). In order to obtain an unbiased comparison, the topography data obtained from the PVO-ORAD were filtered to introduce distortions which are the same as those of the gravity models. Both the gravity and filtered topography maps are derived by two stage processes with a common second stage. In the first stage, the topography was used to calculate a corresponding spacecraft acceleration under the assumptions that the topography has a uniform given density and no compensation. In the second stage, the acceleration measures found in the first stage were passed through a linear inverter to yield maps of gravity and topography. Because these maps are the result of the same inversion process, they contain the same distortion; a comparison between them is unbiased to first order.

Reasenberg, R. D.↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS1971

The effects of random accelerations on estimation accuracy with applications to the Mariner 1969 relativity experiment

A simple point of view is presented for establishing the relationship between a general modeling error and estimation accuracy. Tentative models are given for the random, nongravitational forces thought likely to affect the Mariner 6 and 7 spacecraft. The navigation data from these spacecraft were used to determine the relativity parameter Gamma. The effect that random accelerations have on the classical least squares filter, which assumes the accelerations do not exist, is calculated, and it is demonstrated that their presence seriously affects the estimates, particularly when an extended data arc is employed. Finally, some preliminary exploration using sequential-type filters, which attempt to model and thus ameliorate the effects of these same accelerations, is reported.

Curkendall, D. W.↗
NASA NTRS2021

In-Flight Turbulence Emulation for Fixed-Wing Aircraft using Equivalent Multisine Excitations

A procedure is discussed for designing control-equivalent turbulence input (CETI) excitations to emulate a fixed-wing aircraft response to turbulence while flying through calm air. These excitations are added to the actuator commands to produce the same motion of the aircraft that would have resulted from flying in actual turbulence. Although CETI models are usually shaping filters driven by white noise inputs, the equivalent control inputs in this report were constructed as multisine excitations, which enables mimicking the response to turbulence with arbitrary power spectra, within hardware limitations. Practical aspects of the approach are discussed, such as actuator dynamics, model uncertainty, actual in-flight turbulence, and others. The procedure is demonstrated using a flight dynamics simulation model of a subscale generic jet transport aircraft.

control-equivalent turbulence inputs (CETI)↗
NASA NTRS1979

Velocity structure of the shallow lunar crust

Data from the thumper shots of the Apollo 14 and Apollo 16 active seismic experiments, testing whether the velocity variation in the shallow lunar crust (depths less than or equal to 10 m) can be represented by a self-compacting-power-layer or by a constant-velocity-layer model, are analyzed. Although filtering and stacking improved the S/N ratios, it was found that measuring the arrival times or amplitudes of arrivals beyond 32 m was not possible. The data quality precluded a definitive distinction between the power-law velocity variation and the layered-velocity model. Furthermore, it was found that the shallow lunar regolith is made up of fine particles, which supports the idea of a 1/6 power-velocity model. Analysis of the amplitudes of first arrivals revealed large errors in the data due to variations in the geophone sensitivities and shot strengths; a least-squares method, that uses data redundancy was employed to eliminate them.

Gangi, A. F.↗
DOE OSTI2023

Multi-regime mixing modeling for local extinction and re-ignition in turbulent non-premixed flame by using LES/FDF method

Local extinction and re-ignition occur in turbulent non-premixed combustion when the Damköhler number is not large enough and the combustion is not fully mixing controlled. The occurrence of local extinction introduces locally extinguished flame holes followed by a premixed flame propagation toward the hole center to potentially reignite the flame. The co-existence of the non-premixed and premixed combustion regimes complicates the modeling since traditional combustion models are mostly for a single regime. In this work, we examine the effect of multi-regime mixing modeling in the transported filtered density function (FDF) method on the predictions of local extinction and reignition. Predictions of local extinction and reignition remain a challenge for the FDF method despite the progress made in the past. To account for the multi-regime combustion, two different mixing timescale models for non-premixed and premixed combustion are combined. A flame index based on the gradients of fuel and oxidizer is used to define a weighting factor to blend the two mixing timescale models. A turbulent jet non-premixed flame with substantial local extinction, the Sydney piloted jet flame L, is adopted as a test case to examine the performance of the multi-regime model in large-eddy simulation/FDF modeling. It is found that the traditional non-premixed mixing timescale model when combined with the modified Curl mixing leads to global extinction for the Sydney flame L without the presence of the premixed combustion regime. After accounting for the multi-regime combustion with proper detection of the different combustion regimes, the predictions for the flame statistics and the amount of local extinction are significantly improved. Furthermore, it suggests the need of including multi-regime combustion for the predictions of local extinction and re-ignition in turbulent non-premixed combustion configurations.

42 ENGINEERING↗
NASA NTRS1991

A system identification model for adaptive nonlinear control

A system identification model that combines generalized-spline function approximation with a nonlinear control system is described. The complete control system contains three main elements: a nonlinear-inverse-dynamic control law that depends on a comprehensive model of the plant, a state estimator whose outputs drive the control law, and a function approximation scheme that models the system dynamics. The system-identification task, which combines an extended Kalman filter with a function approximator modeled as an artificial neural network, is considered. The results of an application of the identification techniques to a nonlinear transport aircraft model are presented.

Linse, Dennis J.↗