SEARCH · Search NASA

Results for “State Variable”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Minimum-time acceleration of aircraft turbofan engines

Minimum-time accelerations of the F100 turbofan engine are presented. A piecewise-linear engine model, having three state variables and four control variables, is used to obtain the minimum-time solutions. The linear model which applies at a given time in the trajectory is determined by calculating a normalized 'distance' from the current state to the equilibrium state associated with each linear model. The linear model associated with the closest equilibrium point is then used. The control histories for the minimum-time solutions are used as input to a nonlinear simulation of the F100 engine to verify the accuracy of the piecewise-linear solutions.

Teren, F.↗

NASA NTRS1997

Visualization of Unsteady Computational Fluid Dynamics

The current compute environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array and the J90 cluster) provide the required computation bandwidth for CFD calculations of transient problems. If we follow the traditional computational analysis steps for CFD (and we wish to construct an interactive visualizer) we need to be aware of the following: (1) Disk space requirements. A single snap-shot must contain at least the values (primitive variables) stored at the appropriate locations within the mesh. For most simple 3D Euler solvers that means 5 floating point words. Navier-Stokes solutions with turbulence models may contain 7 state-variables. (2) Disk speed vs. Computational speeds. The time required to read the complete solution of a saved time frame from disk is now longer than the compute time for a set number of iterations from an explicit solver. Depending, on the hardware and solver an iteration of an implicit code may also take less time than reading the solution from disk. If one examines the performance improvements in the last decade or two, it is easy to see that depending on disk performance (vs. CPU improvement) may not be the best method for enhancing interactivity. (3) Cluster and Parallel Machine I/O problems. Disk access time is much worse within current parallel machines and cluster of workstations that are acting in concert to solve a single problem. In this case we are not trying to read the volume of data, but are running the solver and the solver outputs the solution. These traditional network interfaces must be used for the file system. (4) Numerics of particle traces. Most visualization tools can work upon a single snap shot of the data but some visualization tools for transient problems require dealing with time.

Haimes, Robert↗

DOE OSTI2016

Materials science in the time domain using Bragg coherent diffraction imaging

Materials are generally classified by a phase diagram which displays their properties as a function of external state variables, typically temperature and pressure. A new dimension that is relatively unexplored is time: a rich variety of new materials can become accessible in the transient period following laser excitation from the ground state. The timescale of nanoseconds to femtoseconds, is ripe for investigation using x-ray free-electron laser (XFEL) methods. There is no shortage of materials suitable for time-resolved materials-science exploration. Oxides alone represent most of the minerals making up the Earth's crust, catalysts, ferroelectrics, corrosion products and electronically ordered materials such as superconductors, to name a few. Some of the elements have metastable phase diagrams with predicted new phases. There are some examples known already: an oxide 'hidden phase' living only nanoseconds and an electronically ordered excited phase of fullerene C 60 , lasting only femtoseconds. In a completely general way, optically excited states of materials can be probed with Bragg coherent diffraction imaging, both below the damage threshold and in the destructive regime. Lastly, prospective methods for carrying out such XFEL experiments are discussed.

36 MATERIALS SCIENCE↗

NASA NTRS1992

Volume 2: Explicit, multistage upwind schemes for Euler and Navier-Stokes equations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using monotonic upstream schemes for conservation laws (MUSCL)-type differencing to obtain state variables at cell interface. Variable interpolations were written in the k-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator. Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-state predictor-corrector schemes, and multistage time-stepping schemes. The coefficients of the multistage time-stepping schemes have been modified successfully to achieve better performance with upwind differencing. A technique was developed to optimize the coefficients for good high-frequency damping at relatively high CFL numbers. Local time-stepping, implicit residual smoothing, and multigrid procedure were added to the explicit time stepping scheme to accelerate convergence to steady-state. The developed algorithm was implemented successfully in a multi-block code, which provides complete topological and geometric flexibility. The only requirement is C degree continuity of the grid across the block interface. The algorithm has been validated on a diverse set of three-dimensional test cases of increasing complexity. The cases studied were: (1) supersonic corner flow; (2) supersonic plume flow; (3) laminar and turbulent flow over a flat plate; (4) transonic flow over an ONERA M6 wing; and (5) unsteady flow of a compressible jet impinging on a ground plane (with and without cross flow). The emphasis of the test cases was validation of code, and assessment of performance, as well as demonstration of flexibility.

Elmiligui, Alaa↗

DOE OSTI2023

Risk-Constrained Reinforcement Learning for Inverter-Dominated Power System Controls

Here, this paper develops a risk-aware controller for grid-forming inverters (GFMs) to minimize large frequency oscillations in GFM inverter-dominated power systems. To tackle the high variability from loads/renewables, we incorporate a mean-variance risk constraint into the classical linear quadratic regulator (LQR) formulation for this problem. The risk constraint aims to bound the time-averaged cost of state variability and thus can improve the worst-case performance for large disturbances. The resulting risk-constrained LQR problem is solved through the dual reformulation to a minimax problem, by using a reinforcement learning (RL) method termed as stochastic gradient-descent with max-oracle (SGDmax). In particular, the zero-order policy gradient (ZOPG) approach is used to simplify the gradient estimation using simulated system trajectories. Numerical tests conducted on the IEEE 68-bus system have validated the convergence of our proposed SGDmax for GFM model and corroborate the effectiveness of the risk constraint in improving the worst-case performance while reducing the variability of the overall control cost.

Frequency control↗

DOE OSTI2021

The Atmospheric Carbon and Transport (ACT)-America Mission

The Atmospheric Carbon and Transport (ACT)-America NASA Earth Venture Suborbital Mission set out to improve regional atmospheric greenhouse gas (GHG) inversions by exploring the intersection of the strong GHG fluxes and vigorous atmospheric transport that occurs within the midlatitudes. In this study, two research aircraft instrumented with remote and in situ sensors to measure GHG mole fractions, associated trace gases, and atmospheric state variables collected 1,140.7 flight hours of research data, distributed across 305 individual aircraft sorties, coordinated within 121 research flight days, and spanning five 6-week seasonal flight campaigns in the central and eastern United States. Flights sampled 31 synoptic sequences, including fair-weather and frontal conditions, at altitudes ranging from the atmospheric boundary layer to the upper free troposphere. The observations were complemented with global and regional GHG flux and transport model ensembles. We found that midlatitude weather systems contain large spatial gradients in GHG mole fractions, in patterns that were consistent as a function of season and altitude. We attribute these patterns to a combination of regional terrestrial fluxes and inflow from the continental boundaries. These observations, when segregated according to altitude and air mass, provide a variety of quantitative insights into the realism of regional CO 2 and CH 4 fluxes and atmospheric GHG transport realizations. The ACT-America dataset and ensemble modeling methods provide benchmarks for the development of atmospheric inversion systems. As global and regional atmospheric inversions incorporate ACT-America’s findings and methods, we anticipate these systems will produce increasingly accurate and precise subcontinental GHG flux estimates.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1971

Ride comfort control in large flexible aircraft

The problem of ameliorating the discomfort of passengers on a large air transport subject to flight disturbances is examined. The longitudinal dynamics of the aircraft, including effects of body flexing, are developed in terms of linear, constant coefficient differential equations in state variables. A cost functional, penalizing the rigid body displacements and flexure accelerations over the surface of the aircraft is formulated as a quadratic form. The resulting control problem, to minimize the cost subject to the state equation constraints, is of a class whose solutions are well known. The feedback gains for the optimal controller are calculated digitally, and the resulting autopilot is simulated on an analog computer and its performance evaluated.

Warren, M. E.↗

DOE OSTI2025

Comparison Study of Conventional and Adaptive Mesh Refinement in Organic Material Decomposition Models

This study compares conventional mesh refinement techniques, specifically Uniform Mesh Refinement (UMR), with a new Adaptive Mesh Refinement (AMR) method, applied to Organic Material Decomposition (OMD) models. The proposed benefit of AMR is that only areas that require refinement, based on minimizing a specific field gradient, are refined thus decreasing model wall time compared to conventional UMR methods. This work specifically focuses on comparing UMR and AMR methods on decomposing (both No-Flow and Porous-Flow material models) Polymeric Methylene Diisocyanate (PMDI) polyurethane foam. Throughout the work, the geometry increased in complexity to assess the refinement methods performance at varying levels geometric intricacy. While AMR has been shown to work well in a variety of applications, the UMR approach proved to be computationally faster, for many of the geometries and foam decomposition models, than AMR. However, it was observed that at higher levels of refinement, greater than 3 UMR, AMR begins to be computationally better. Additionally, the settings used to perform AMR greatly impact its performance, and lessons learned, in terms of OMD models, are shared. Due to physics involved in material decomposition, specifically the evolution of state variables, these problems don’t fully benefit from the advantages of AMR.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Laser-induced quenching of metastability at the Mott insulator to metal transition

There is growing interest in strongly correlated insulator thin films because the intricate interplay of their intrinsic and extrinsic state variables causes memristive behavior that might be used for biomimetic devices in the emerging field of neuromorphic computing. Here, in this study, we find that laser irradiation tends to drive V 2 ⁢O 3 from supercooled/superheated metastable states toward thermodynamic equilibrium, most likely in a nonthermal way. We study thin films of the prototypical Mott-insulator V2⁢O3, which show spontaneous phase separation into metal-insulator herringbone domains during the Mott transition. Here, we use low-temperature microscopy to investigate how these metal-insulator domains can be modified by scanning a focused laser beam across the thin film surface. We find that the response depends on the thermal history: When the thin film is heated from below the Mott transition temperature, the laser beam predominantly induces metallic domains. On the contrary, when the thin film is cooled from a temperature above the transition, the laser beam predominantly induces insulating domains. Very likely, the V 2 ⁢O 3 thin film is in a superheated or supercooled state, respectively, during the first-order phase transition, and the perturbation by a laser beam drives these metastable states into stable ones. This way, the thermal history is locally erased. Our findings are supported by a phenomenological model with a laser-induced lowering of the energy barrier between the metastable and equilibrium states.

Materials Science↗

NASA NTRS1989

Model reduction in the physical coordinate system

In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.

Yae, K. Harold↗

NASA NTRS2000

Timing analysis by model checking

The safety of modern avionics relies on high integrity software that can be verified to meet hard real-time requirements. The limits of verification technology therefore determine acceptable engineering practice. To simplify verification problems, safety-critical systems are commonly implemented under the severe constraints of a cyclic executive, which make design an expensive trial-and-error process highly intolerant of change. Important advances in analysis techniques, such as rate monotonic analysis (RMA), have provided a theoretical and practical basis for easing these onerous restrictions. But RMA and its kindred have two limitations: they apply only to verifying the requirement of schedulability (that tasks meet their deadlines) and they cannot be applied to many common programming paradigms. We address both these limitations by applying model checking, a technique with successful industrial applications in hardware design. Model checking algorithms analyze finite state machines, either by explicit state enumeration or by symbolic manipulation. Since quantitative timing properties involve a potentially unbounded state variable (a clock), our first problem is to construct a finite approximation that is conservative for the properties being analyzed-if the approximation satisfies the properties of interest, so does the infinite model. To reduce the potential for state space explosion we must further optimize this finite model. Experiments with some simple optimizations have yielded a hundred-fold efficiency improvement over published techniques.

Naydich, Dimitri↗

NASA NTRS2010

Land Surface Emission Modeling to Support Physical Precipitation Retrievals

Peters-Lidard, Christina D.↗

NASA NTRS2011

Global Precipitation Measurement, Validation, and Applications Integrated Hydrologic Validation to Improve Physical Precipitation Retrievals for GPM

Land surface modeling and data assimilation can provide dynamic land surface state variables necessary to support physical precipitation retrieval algorithms over land. It is well-known that surface emission, particularly over the range of frequencies to be included in the Global Precipitation Measurement Mission (GPM), is sensitive to land surface states, including soil properties, vegetation type and greenness, soil moisture, surface temperature, and snow cover, density, and grain size. In order to investigate the robustness of both the land surface model states and the microwave emissivity and forward radiative transfer models, we have undertaken a multi-site investigation as part of the NASA Precipitation Measurement Missions (PMM) Land Surface Characterization Working Group. Specifically, we will demonstrate the performance of the Land Information System (LIS; http://lis.gsfc.nasa.gov; Peters-Lidard et aI., 2007; Kumar et al., 2006) coupled to the Joint Center for Satellite Data Assimilation (JCSDA's) Community Radiative Transfer Model (CRTM; Weng, 2007; van Deist, 2009). The land surface is characterized by complex physical/chemical constituents and creates temporally and spatially heterogeneous surface properties in response to microwave radiation scattering. The uncertainties in surface microwave emission (both surface radiative temperature and emissivity) and very low polarization ratio are linked to difficulties in rainfall detection using low-frequency passive microwave sensors (e.g.,Kummerow et al. 2001). Therefore, addressing these issues is of utmost importance for the GPM mission. There are many approaches to parameterizing land surface emission and radiative transfer, some of which have been customized for snow (e.g., the Helsinki University of Technology or HUT radiative transfer model;) and soil moisture (e.g., the Land Surface Microwave Emission Model or LSMEM).

Peters-Lidar, Christa D.↗

NASA NTRS2019

Analyzing Machine Learning Predictions of Passive Microwave Brightness Temperature Spectral Difference Over Snow-Covered Terrain in High Mountain Asia

Snow is an important component of the terrestrial freshwater budget in high mountainAsia (HMA) and contributes to the runoff in Himalayan rivers through snowmelt. Despitethe importance of snow in HMA, considerable spatiotemporal uncertainty exists across the different estimates of snow water equivalent for this region. In order to better estimate snow water equivalent, radiative transfer models are often used in conjunction with microwave brightness temperature measurements. In this study, the efficacy of support vector machines (SVMs), a machine learning technique, to predict passive microwave brightness temperature spectral difference (1Tb) as a function of geophysical variables (snow water equivalent, snow depth, snow temperature, and snow density) is explored through a sensitivity analysis. The use of machine learning (as opposed to radiative transfer models) is a relatively new and novel approach for improving snow water equivalent estimates. The Noah-MP land surface model within the NASALand Information System framework is used to simulate the hydrologic cycle over HMA and model geophysical variables that are then used for SVM training. The SVMsserve as a nonlinear map between the geophysical space (modeled in Noah-MP) andthe observation space (1Tb as measured by the radiometer). Advanced MicrowaveScanning Radiometer-Earth Observing System measured passive microwave brightness temperatures over snow-covered locations in the HMA region are used as training data during the SVM training phase. Sensitivity of well-trained SVMs to each Noah-MP modeled state variable is assessed by computing normalized sensitivity coefficients. Sensitivity analysis results generally conform with the known first-order physics. Input states that increase volume scattering of microwave radiation, such as snow density and snow water equivalent, exhibit a plurality of positive normalized sensitivity coefficients. In general, snow temperature was the most sensitive input to the SVM predictions. The sensitivity of each state is location and time dependent. The signs of normalized sensitivity coefficients that indicate physical irrationality are ascribed to significant cross-correlation between Noah-MP simulated states and decreased SVM prediction capability at specific locations due to insufficient training data. SVM prediction pitfalls do exist that serve to highlight the limitations of this particular machine learning algorithm.

high mountain Asia↗

DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗

DOE OSTI2022

Physics-based compact modeling of electro-thermal memristors: Negative differential resistance, local activity, and non-local dynamical bifurcations

Leon Chua's Local Activity theory quantitatively relates the compact model of an isolated nonlinear circuit element, such as a memristor, to its potential for desired dynamical behaviors when externally coupled to passive elements in a circuit. However, the theory's use has often been limited to potentially unphysical toy models and analyses of small-signal linear circuits containing pseudo-elements (resistors, capacitors, and inductors), which provide little insight into required physical, material, and device properties. Furthermore, the Local Activity concept relies on a local analysis and must be complemented by examining dynamical behavior far away from the steady-states of a circuit. In this work, we review and study a class of generic and extended one-dimensional electro-thermal memristors (i.e., temperature is the sole state variable), re-framing the analysis in terms of physically motivated definitions and visualizations to derive intuitive compact models and simulate their dynamical behavior in terms of experimentally measurable properties, such as electrical and thermal conductance and capacitance and their derivatives with respect to voltage and temperature. Within this unified framework, we connect steady-state phenomena, such as negative differential resistance, and dynamical behaviors, such as instability, oscillations, and bifurcations, through a set of dimensionless nonlinearity parameters. In particular, we reveal that the reactance associated with electro-thermal memristors is the result of a phase shift between oscillating current and voltage induced by the dynamical delay and coupling between the electrical and thermal variables. We thus, demonstrate both the utility and limitations of local analyses to understand non-local dynamical behavior. Critically for future experimentation, the analyses show that external coupling of a memristor to impedances within modern sourcing and measurement instruments can dominate the response of the total circuit, making it impossible to characterize the response of an uncoupled circuit element for which a compact model is desired. However, these effects can be minimized by proper understanding of the Local Activity theory to design and utilize purpose-built instruments.

Brown, Timothy D. (ORCID:0000000347068355)↗

NASA NTRS1978

Human operator identification model and related computer programs

Four computer programs which provide computational assistance in the analysis of man/machine systems are reported. The programs are: (1) Modified Transfer Function Program (TF); (2) Time Varying Response Program (TVSR); (3) Optimal Simulation Program (TVOPT); and (4) Linear Identification Program (SCIDNT). The TV program converts the time domain state variable system representative to frequency domain transfer function system representation. The TVSR program computes time histories of the input/output responses of the human operator model. The TVOPT program is an optimal simulation program and is similar to TVSR in that it produces time histories of system states associated with an operator in the loop system. The differences between the two programs are presented. The SCIDNT program is an open loop identification code which operates on the simulated data from TVOPT (or TVSR) or real operator data from motion simulators.

Kessler, K. M.↗

NASA NTRS2022

Perspective on Satellite-Based Land Data Assimilation to Estimate Water Cycle Components in an Era of Advanced Data Availability and Model Sophistication

The beginning of the 21st century is marked by a rapid growth of land surface satellite data and model sophistication. This offers new opportunities to estimate multiple components of the water cycle via satellite-based land data assimilation (DA) across multiple scales. By resolving more processes in land surface models and by coupling the land, the atmosphere, and other Earth system compartments, the observed information can be propagated to constrain additional unobserved variables. Furthermore, access to more satellite observations enables the direct constraint of more and more components of the water cycle that are of interest to end users. However, the finer level of detail in models and data is also often accompanied by an increase in dimensions, with more state variables, parameters, or boundary conditions to estimate, and more observations to assimilate. This requires advanced DA methods and efficient solutions. One solution is to target specific observations for assimilation based on a sensitivity study or coupling strength analysis, because not all observations are equally effective in improving subsequent forecasts of hydrological variables, weather, agricultural production, or hazards through DA. This paper offers a perspective on current and future land DA development, and suggestions to optimally exploit advances in observing and modeling systems.

Data Assimilation↗