Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent 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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Time-controlled descent guidance in uncertain winds

A procedure has been developed for constructing a statistical model of the altitude-dependent mean wind profile from the historical record of wind measurements at particular locations. The model is constructed by fitting a Markov process, with altitude as the stage variable, to the historical wind data. The wind model, together with the aircraft dynamics and the error characteristics of the navigation system, are incorporated in the design of a state estimator, which gives the minimum variance estimate of the aircraft state and the wind vector. The state and wind estimates are used as inputs to a linear feedback law for guiding the aircraft along the nominal trajectory. An example design of a time-constrained (4D RNAV) descent guidance system is presented, showing tracking accuracy, control activity, and probability of arrival time with and without the wind estimator.

Menga, G.↗
NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗
NASA NTRS2013

Mission Data System Java Edition Version 7

The Mission Data System framework defines closed-loop control system abstractions from State Analysis including interfaces for state variables, goals, estimators, and controllers that can be adapted to implement a goal-oriented control system. The framework further provides an execution environment that includes a goal scheduler, execution engine, and fault monitor that support the expression of goal network activity plans. Using these frameworks, adapters can build a goal-oriented control system where activity coordination is verified before execution begins (plan time), and continually during execution. Plan failures including violations of safety constraints expressed in the plan can be handled through automatic re-planning. This version optimizes a number of key interfaces and features to minimize dependencies, performance overhead, and improve reliability. Fault diagnosis and real-time projection capabilities are incorporated. This version enhances earlier versions primarily through optimizations and quality improvements that raise the technology readiness level. Goals explicitly constrain system states over explicit time intervals to eliminate ambiguity about intent, as compared to command-oriented control that only implies persistent intent until another command is sent. A goal network scheduling and verification process ensures that all goals in the plan are achievable before starting execution. Goal failures at runtime can be detected (including predicted failures) and handled by adapted response logic. Responses can include plan repairs (try an alternate tactic to achieve the same goal), goal shedding, ignoring the fault, cancelling the plan, or safing the system.

Reinholtz, William K.↗
NASA NTRSDate not supplied

Effect of Pressure and Oxidation State on Pnictogen Metal-Silicate Partitioning

The partitioning of trace elements between metal and silicate melts serves as the foundation for understanding the differentiation of a planetary body into a metallic core and silicate mantle. Element activity influences metal-silicate partitioning behavior. Activity coefficients are directly dependent on composition and temperature, and can be indirectly dependent on oxygen fugacity and pressure. Distinguishing the effect of pressure from other variables on the activity coefficients and partitioning is important for understanding the chemical evolution of different planetary bodies during differentiation. In this study, we investigated the influence of pressure on the activity coefficients of Cu, Mo, Pd, Pt, As, Sb, and Bi in Fe-Si metallic liquids. All of these elements exhibit moderate to high activity coefficients in Fe-Si liquids at low pressure, which significantly controls their metal-silicate partitioning behavior. Identifying whether this strong dependence persists at higher pressures is critical to modeling and understanding the chemical consequences of core formation. New experiments at 10 GPa were used to derive activity coefficients for these metals which can be compared to activity coefficients determined at 1 GPa. Experiments were conducted at 10 GPa and 2373K using a 10/5 assembly in the 880-ton multi anvil press at NASA Johnson Space Center. The standard 10/5 COMPRES assembly was slightly modified to accommodate a sample capsule machined from single-crystal MgO that minimizes melt percolation out of the sample volume during the experiment. Experiment starting materials were comprised of 70 wt.% Knippa basalt and 30 wt.% metal. The metal mixture (~85 wt.% Fe) was created by adding the elements of interest (Cu, Mo, Pd, Pt, As, Sb, Bi) to Fe metal powder. Varied amounts of Si metal (0-10 wt.% Si) were added to the metal-silicate mixtures to generate a systematic series of starting materials. For each element, an epsilon interaction parameter in Fe-Si liquid was derived from the results of our 0-10 wt.% Si metal series. To investigate whether pressure influences the trace element activity coefficients in the 1-10 GPa pressure range, we compared our results at 10 GPa to those at 1-4 GPa [1-3]. Our results can also be directly compared to interaction parameters for Au, P, V, Mn, Ga, Zn, Cd, Sn, W, Pb, and Nb previously determined at 10 GPa and 2373K following the same methods [4]. Combined, this suite of interaction parameters will directly inform metal- silicate partitioning between 1 and 10 GPa, and assess whether these values can be extrapolated to modeling differentiation processes at pressures >10 GPa. [1] Righter et al. (2018) GCA 232, 101-123 [2] Righter et al. (2019) MaPS 54, 1379-1394 [3] Steenstra et al. (2020) Icarus, 335, 113391 [4] Righter et al. (2020) Geochem. Persp. Let. 15, 44-49

Kelsey Prissel↗
NASA NTRS1987

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems. II

The results of Goldberg and Tadmor (1985) are extended to achieve improved stability criteria for a large class of approximations to the initial boundary value problem associated with a particular hyperbolic system in a quarter plane. In a stability analysis, it is shown that the entire approximation is stable if and only if the scalar outflow components of its principal part are stable. Thus, the global stability question is reduced to that of a scalar, homogeneous outflow problem. The stability criteria for the reduced problem, which depend both on the basic scheme and the boundary conditions, but very little on the interaction between the two, are stated and used to establish previous examples and new ones, including a host of dissipative and nondissipative examples. There is no difficulty in extending the stability criteria to two-boundary problems and initial-boundary value problems with variable coefficients.

Goldberg, Moshe↗
NASA NTRS1991

A Variational Formalism for the Radiative Transfer Equation: Prelude to Model 3

The MODEL III variational data assimilation model is the third of four general assimilation models designed to blend weather data measured from space based platforms in the meteorological data mainstream in a way that maximizes the information content of the satellite data. Because there are many different observation locations and there are many instruments with different measurement error characteristics, it is also necessary to require that the blending be done to maximize the information content of the data and simultaneously to retain a dynamically consistent and reasonably accurate description of the state of the atmosphere. This is ideally a variational problem for which the data receive relative weights that are inversely proportional to measurement error and are adjusted to satisfy a set of dynamical equations that govern atmospheric processes. The advantage of MODEL III over the previous two models is that radiance, the atmospheric variable measured by satellite, becomes a dependent variable. In the previous versions, mean layer temperatures that had been retrieved from the radiances by some method, were included in the assimilation by substituting them in place of the rawinsonde temperatures. Now both rawinsonde temperatures and satellite radiances are included independently in the assimilation.

Achtemeier, Gary L.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
NASA NTRS2014

Fine Particulate Matter Predictions Using High Resolution Aerosol Optical Depth (AOD) Retrievals

To date, spatial-temporal patterns of particulate matter (PM) within urban areas have primarily been examined using models. On the other hand, satellites extend spatial coverage but their spatial resolution is too coarse. In order to address this issue, here we report on spatial variability in PM levels derived from high 1 km resolution AOD product of Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm developed for MODIS satellite. We apply day-specific calibrations of AOD data to predict PM(sub 2.5) concentrations within the New England area of the United States. To improve the accuracy of our model, land use and meteorological variables were incorporated. We used inverse probability weighting (IPW) to account for nonrandom missingness of AOD and nested regions within days to capture spatial variation. With this approach we can control for the inherent day-to-day variability in the AOD-PM(sub 2.5) relationship, which depends on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and ground surface reflectance among others. Out-of-sample "ten-fold" cross-validation was used to quantify the accuracy of model predictions. Our results show that the model-predicted PM(sub 2.5) mass concentrations are highly correlated with the actual observations, with out-of- sample R(sub 2) of 0.89. Furthermore, our study shows that the model captures the pollution levels along highways and many urban locations thereby extending our ability to investigate the spatial patterns of urban air quality, such as examining exposures in areas with high traffic. Our results also show high accuracy within the cities of Boston and New Haven thereby indicating that MAIAC data can be used to examine intra-urban exposure contrasts in PM(sub 2.5) levels.

aerosol optical depth↗
NASA NTRS2020

Exploring Spatiotemporal Relations Between Soil Moisture, Precipitation, and Streamflow for a Large Set of Watersheds Using Google Earth Engine

An understanding of streamflow variability and its response to changes in climate conditions is essential for water resource planning and management practices that will help to mitigate the impacts of extreme events such as floods and droughts on agriculture and other human activities. This study investigated the relationship between precipitation, soil moisture, and streamflow over a wide range of watersheds across the United States using Google Earth Engine (GEE). The correlation analyses disclosed a strong association between precipitation, soil moisture, and streamflow, however, soil moisture was found to have a higher correlation with the streamflow relative to precipitation. Results indicated different strength of the association depends on the watershed classes and lag times assessments. The perennial watersheds showed higher coherence compared to intermittent watersheds. Previous month precipitation and soil moisture have a stronger influence on the current month streamflow, particularly in the snow-dominated watersheds. Monthly streamflow forecasting models were developed using an autoregressive integrated moving average (ARIMA) and support vector machine (SVM). The results showed that the SVM model generally performed better than the ARIMA model. Overall streamflow forecasting model performance varied considerably among watershed classes, and perennial watersheds tend to exhibit better predictably compared to intermittent watersheds due to lower streamflow variability. The SVM models with precipitation and streamflow inputs performed better than those with streamflow input only. Results indicated that the inclusion of antecedent root-zone soil moisture improved the streamflow forecasting in most of the watersheds, and the largest improvements occurred in the intermittent watersheds. In conclusion, this work demonstrated that knowing the relationship between precipitation, soil moisture, and streamflow in different watershed classes will enhance the understanding of the hydrologic process and can be effectively utilized in improving streamflow forecasting for better satellite-based water resource management strategies.

Nazmus Sazib↗
NASA NTRS1992

Natural versus anthropogenic factors affecting low-level cloud albedo over the North Atlantic

Cloud albedo plays a key role in regulating earth's climate. Cloud albedo depends on column-integrated liquid water content and the density of cloud condensation nuclei, which consists primarily of submicrometer-sized aerosol sulfate particles. A comparison of two independent satellite data sets suggests that, although anthropogenic sulfate emissions may enhance cloud albedo immediately adjacent to the east coast of the United States, over the central North Atlantic Ocean the variability in albedo can be largely accounted for by natural marine and atmospheric processes that probably have remained relatively constant since the beginning of the industrial revolution.

Falkowski, Paul G.↗
NASA NTRS2015

Analytic Thermoelectric Couple Modeling: Variable Material Properties and Transient Operation

To gain a deeper understanding of the operation of a thermoelectric couple a set of analytic solutions have been derived for a variable material property couple and a transient couple. Using an analytic approach, as opposed to commonly used numerical techniques, results in a set of useful design guidelines. These guidelines can serve as useful starting conditions for further numerical studies, or can serve as design rules for lab built couples. The analytic modeling considers two cases and accounts for 1) material properties which vary with temperature and 2) transient operation of a couple. The variable material property case was handled by means of an asymptotic expansion, which allows for insight into the influence of temperature dependence on different material properties. The variable property work demonstrated the important fact that materials with identical average Figure of Merits can lead to different conversion efficiencies due to temperature dependence of the properties. The transient couple was investigated through a Greens function approach; several transient boundary conditions were investigated. The transient work introduces several new design considerations which are not captured by the classic steady state analysis. The work helps to assist in designing couples for optimal performance, and also helps assist in material selection.

Seebeck Effect↗
NASA NTRS1986

Two-phase flows and heat transfer within systems with ambient pressure above the thermodynamic critical pressure

In systems where the design inlet and outlet pressures P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.

Hendricks, R. C.↗
NASA NTRS1986

Two-phase flows and heat transfer within systems with ambient pressure above the thermodynamic critical pressure

In systems where the design inlet and outlet pressure P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.

Hendricks, R. C.↗
NASA NTRS1992

A neural fuzzy controller learning by fuzzy error propagation

In this paper, we describe a procedure to integrate techniques for the adaptation of membership functions in a linguistic variable based fuzzy control environment by using neural network learning principles. This is an extension to our work. We solve this problem by defining a fuzzy error that is propagated back through the architecture of our fuzzy controller. According to this fuzzy error and the strength of its antecedent each fuzzy rule determines its amount of error. Depending on the current state of the controlled system and the control action derived from the conclusion, each rule tunes the membership functions of its antecedent and its conclusion. By this we get an unsupervised learning technique that enables a fuzzy controller to adapt to a control task by knowing just about the global state and the fuzzy error.

Nauck, Detlef↗
NASA NTRS1984

Selection of Slow and Fast Variables in the Three-Dimensional Flight Dynamics

A transformation from the altitude and velocity state variables of three-dimensional flight mechanics to a new set of more desirable variables is found. The new variables provide a greater time-scale separation, decrease system coupling, and give better estimates of the fast-variable values along the reduced solution. One of the new variables is the often-used specific energy, whereas the other variable changes along a given trajectory, depending on the nature of the local reduced solution. Numerical examples are included.

Ardema, M.↗
NASA NTRS2012

Evidence for Quasi-Periodic X-ray Dips from an ULX: Implications for the Binary Motion and the Orbital Inclination

We report results from long-term X-ray (0.3-8.0 keY) monitoring of the ultraluminous X-ray source NGC 5408 X-1 with the Swift/X-Ray Telescope. Our primary results are: (1) the discovery of quasi-periodic dips in the X-ray intensity that recur on average every 243 days, (2) the detection of an energy-dependent (variability amplitude decreases with increasing energy), quasi-sinusoidal X-ray modulation with a period of 112.6 +/- 4 days the amplitude of which decreases during the second half of the light curve and (3) energy spectral evidence for an increase in photoelectric absorption during the last continuous segment of the data, possibly due to a change in the ionization state of the circumbinary material. We interpret the X-ray modulations in the context of binary motion in analogy to that seen in high-inclination low-mass X-ray binaries. If correct, this implies that NGC 5408 X-1 is in a binary with an orbital period of 243 +/- 23 days in contrast to the 115.5 day quasi-sinusoidal period previously reported. In addition, if the X-ray modulation is caused by vertically structured obscuring material in the accretion disk (similar to the phenomenon of dipping LMXBs), this would imply a high value for the inclination of the orbit. A comparison with estimates from accreting X-ray binaries suggests an inclination approx > 60 deg. We note that, in principle, a precessing accretion disk could also produce the observed X-ray modulations.

Pasham, Dheeraj R.↗
NASA NTRS2020

Development of a Non-Iterative Balance Load Prediction Algorithm for the NASA Ames Unitary Plan Wind Tunnel

A non-iterative load prediction algorithm for strain-gage balances was developed for the NASA Ames Unitary Plan Wind Tunnels that computes balance loads from the electrical outputs of the balance bridges and a set of state variables. A state variable could be, for example, a balance temperature difference or the bellows pressure of a flow-through balance. The algorithm directly uses regression models of the balance loads for the load prediction that were obtained by applying global regression analysis to balance calibration data. This choice greatly simplifies both implementation and use of the load prediction process for complex balance configurations as no load iteration needs to be performed. The regression model of a balance load is constructed by using terms from a total of nine term groups. Four term groups are derived from a Taylor Series expansion of the relationship between the load, gage outputs, and state variables. The remaining five term groups are defined by using absolute values of the gage outputs and state variables. Terms from these groups should only be included in the regression model if calibration data from a balance with known bi-directional outputs is analyzed. It is illustrated in detail how global regression analysis may be applied to obtain the coefficients of the chosen regression model of a load component assuming that no linear or massive near-linear dependencies between the regression model terms exist. Data from the machine calibration of a six-component force balance is used to illustrate both application and accuracy of the non-iterative load prediction process.

Ulbrich, Norbert M.↗
NASA NTRS2011

Multi-Wavelength, Multi-Beam, and Polarization-Sensitive Laser Transmitter for Surface Mapping

A multi-beam, multi-color, polarized laser transmitter has been developed for mapping applications. It uses commercial off-the-shelf components for a lowcost approach for a ruggedized laser suitable for field deployment. The laser transmitter design is capable of delivering dual wavelengths, multiple beams on each wavelength with equal (or variable) intensities per beam, and a welldefined state of polarization. This laser transmitter has been flown on several airborne campaigns for the Slope Imaging Multi-Polarization Photon Counting Lidar (SIMPL) instrument, and at the time of this reporting is at a technology readiness level of between 5 and 6. The laser is a 1,064-nm microchip high-repetition-rate laser emitting energy of about 8 microjoules per pulse. The beam was frequency-doubled to 532 nm using a KTP (KTiOPO4) nonlinear crystal [other nonlinear crystals such as LBO (LiB3O5) or periodically poled lithium niobiate can be used as well, depending on the conversion efficiency requirements], and the conversion efficiency was approximately 30 percent. The KTP was under temperature control using a thermoelectric cooler and a feedback monitoring thermistor. The dual-wavelength beams were then spectrally separated and each color went through its own optical path, which consisted of a beam-shaping lens, quarterwave plate (QWP), and a birefringent crystal (in this case, a calcite crystal, but others such as vanadate can be used). The QWP and calcite crystal set was used to convert the laser beams from a linearly polarized state to circularly polarized light, which when injected into a calcite crystal, will spatially separate the circularly polarized light into the two linear polarized components. The spatial separation of the two linearly polarized components is determined by the length of the crystal. A second set of QWP and calcite then further separated the two beams into four. Additional sets of QWP and calcite can be used to further split the beams into multiple orders of two. The spatially separated beams had alternating linearly polarization states; a half-wave plate (HWP) array was then made to rotate the alternating states of A multi-beam, multi-color, polarized laser transmitter has been developed for mapping applications. It uses commercial off-the-shelf components for a lowcost approach for a ruggedized laser suitable for field deployment. The laser transmitter design is capable of delivering dual wavelengths, multiple beams on each wavelength with equal (or variable) intensities per beam, and a welldefined state of polarization. This laser transmitter has been flown on several airborne campaigns for the Slope Imaging Multi-Polarization Photon Counting Lidar (SIMPL) instrument, and at the time of this reporting is at a technology readiness level of between 5 and 6. The laser is a 1,064-nm microchip high-repetition-rate laser emitting energy of about 8 microjoules per pulse. The beam was frequency-doubled to 532 nm using a KTP (KTiOPO4) nonlinear crystal [other nonlinear crystals such as LBO (LiB3O5) or periodically poled lithium niobiate can be used as well, depending on the conversion efficiency requirements], and the conversion efficiency was approximately 30 percent. The KTP was under temperature control using a thermoelectric cooler and a feedback monitoring thermistor. The dual-wavelength beams were then spectrally separated and each color went through its own optical path, which consisted of a beam-shaping lens, quarterwave plate (QWP), and a birefringent crystal (in this case, a calcite crystal, but others such as vanadate can be used). The QWP and calcite crystal set was used to convert the laser beams from a linearly polarized state to circularly polarized light, which when injected into a calcite crystal, will spatially separate the circularly polarized light into the two linear polarized components. The spatial separation of the two linearly polarized components is determined by the length of the crystal. A cond set of QWP and calcite then further separated the two beams into four. Additional sets of QWP and calcite can be used to further split the beams into multiple orders of two. The spatially separated beams had alternating linearly polarization states; a half-wave plate (HWP) array was then made to rotate the alternating states of

Yu, Anthony W.↗