Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

The development of optimal control laws for orbiting tethered platform systems

A mathematical model of the open and closed loop in orbit plane dynamics of a space platform-tethered-subsatellite system is developed. The system consists of a rigid platform from which an (assumed massless) tether is deploying (retrieving) a subsatellite from an attachment point which is, in general, offset from the platform's mass center. A Langrangian formulation yields equations describing platform pitch, subsatellite tetherline swing, and varying tether length motions. These equations are linearized about the nominal station keeping motion. Control can be provided by both modulation of the tether tension level and by a momentum type platform-mounted device; system controllability depends on the presence of both control inputs. Stability criteria are developed in terms of the control law gains, the platform inertia ratio, and tether offset parameter. Control law gains are obtained based on linear quadratic regulator techniques. Typical transient responses of both the state and required control effort are presented.

Bainum, P. M.↗
NASA NTRS1986

The development of optimal control laws for orbiting tethered platform systems

A mathematical model of the open and closed loop in-orbit plane dynamics of a space platform-tethered-subsatellite system is developed. The system consists of a rigid platform from which an (assumed massless) tether is deploying (retrieving) a subsatellite from an attachment point which is, in general, offset from the platform's mass center. A Lagrangian formulation yields equations describing platform pitch, subsatellite tether-line swing, and varying tether length motions. These equations are linearized about the nominal station keeping motion. Control can be provided by both modulation of the tether tension level and by a momentum type platform-mounted device; system controllability depends on the presence of both control inputs. Stability criteria are developed in terms of the control law gains, the platform inertia ratio, and tether offset parameter. Control law gains are obtained based on linear quadratic regulator techniques. Typical transient responses of both the state and required control effort are presented.

Bainum, P. M.↗
NASA NTRS1996

Effects of modeling errors on trajectory predictions in air traffic control automation

Air traffic control automation synthesizes aircraft trajectories for the generation of advisories. Trajectory computation employs models of aircraft performances and weather conditions. In contrast, actual trajectories are flown in real aircraft under actual conditions. Since synthetic trajectories are used in landing scheduling and conflict probing, it is very important to understand the differences between computed trajectories and actual trajectories. This paper examines the effects of aircraft modeling errors on the accuracy of trajectory predictions in air traffic control automation. Three-dimensional point-mass aircraft equations of motion are assumed to be able to generate actual aircraft flight paths. Modeling errors are described as uncertain parameters or uncertain input functions. Pilot or autopilot feedback actions are expressed as equality constraints to satisfy control objectives. A typical trajectory is defined by a series of flight segments with different control objectives for each flight segment and conditions that define segment transitions. A constrained linearization approach is used to analyze trajectory differences caused by various modeling errors by developing a linear time varying system that describes the trajectory errors, with expressions to transfer the trajectory errors across moving segment transitions. A numerical example is presented for a complete commercial aircraft descent trajectory consisting of several flight segments.

Jackson, Michael R. C.↗
NASA NTRS1999

Nonlinear Local Bending Response and Bulging Factors for Longitudinal Cracks in Pressurized Cylindrical Shells

Results of a geometrically nonlinear finite element parametric study to determine curvature correction factors or "bulging factors" that account for increased stresses due to curvature for longitudinal cracks in unstiffened pressurized cylindrical shells are presented. Geometric parameters varied in the study include the shell radius, the shell wall thickness, and the crack length. The major results are presented in graphs of the bulging factor as a function of the applied load and as a function of geometric parameters that include the shell radius, the shell thickness and the crack length. The computed bulging factors are compared with solutions based on linear shallow shell theory, and with semi-empirical solutions that approximately account for the nonlinear deformation in the vicinity of the crack. The effect of biaxial loads on the computed bulging factors is also discussed.

Rose, Cheryl A.↗
NASA NTRS2011

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗
DOE OSTI2023

Gas–Liquid Chemical Reactions with Nanosecond Pulses: Role of Frequency and Pulse Delivery Modes

The effects of changes in the mode of delivery of nanosecond pulses in a gas-liquid plasma reactor on the formation of hydrogen peroxide, H 2 O 2 , and, as an indicator for ·OH radicals, the decoloration of methylene blue, MB, were determined for pulse delivery by a) increasing frequency with uniform pulses (5-50 kHz), b) variation of the time between bursts of pulses (burst period), c) changing the inner burst frequency (1 over the time between the pulses in the burst), and d) variation of number of pulses in a burst (N-cycles). H 2 O 2 peroxide formation was not affected by the method of pulse delivery in the range of parameters studied here and all data followed an approximately linear increase in H 2 O 2 production rate with discharge power. In contrast, the MB decoloration rate was affected by the burst modes. In terms of discharge power, the MB decoloration rate was highest for the uniform pulse mode; however, the linear trend in increase of MB decoloration with power when the burst period was varied, suggest that at higher power the burst mode may be more effective than the uniform pulsing. Consideration of the per pulse decoloration with energy per pulse and with number of pulses suggest that the burst mode can affect reactions differently from applying a uniform pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

3D pattern formation of a protein–membrane suspension

Many essential cellular processes, including cell division and the establishment of cell polarity during embryogenesis, are regulated by pattern-forming proteins. These proteins often need to bind to a substrate, such as the cell membrane, onto which they interact and form two-dimensional (2D) patterns. It is unclear how the membrane’s continuity and dimensionality impact pattern formation. Here, we address this gap using the MinDE system, a prototypical example of pattern-forming membrane proteins. We show that when the lipid substrate is fragmented into submicrometer-sized diffusive liposomes, adenosine triphosphate-driven protein–protein interactions generate three-dimensional (3D) spatially extended patterns, despite the complete loss of membrane continuity. Remarkably, these 3D patterns emerge at scales four orders of magnitude larger than the individual liposomes. By systematically varying protein concentration, liposome size, and density, we observed and characterized a variety of 3D dynamical patterns not seen on continuous 2D membranes, including traveling waves, dynamical spirals, and a coexistence phase. Simulations and linear stability analysis of a coarse-grained model revealed that the physical properties of the dispersed membrane effectively rescale both the protein–membrane binding rates and diffusion, two key parameters governing pattern formation and wavelength selection. These findings highlight the robustness of Min’s pattern-forming ability, suggesting that protein–membrane suspensions could serve as an adaptable template for studying out-of-equilibrium self-organization in 3D, beyond in vivo contexts.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Ionic Polymer-Metal Composites (IPMCs) as Biomimetic Sensors, Actuators and Artificial Muscles: A Review

This paper presents an introduction to ionic polymer-metal composites and some mathematical modeling pertaining to them. It further discusses a number of recent findings in connection with ion-exchange polymer-metal composites (IPMCS) as biomimetic sensors and actuators. Strips of these composites can undergo large bending and flapping displacement if an electric field is imposed across their thickness. Thus, in this sense they are large motion actuators. Conversely by bending the composite strip, either quasi-statically or dynamically, a voltage is produced across the thickness of the strip. Thus, they are also large motion sensors. The output voltage can be calibrated for a standard size sensor and correlated to the applied loads or stresses. They can be manufactured and cut in any size and shape. In this paper first the sensing capability of these materials is reported. The preliminary results show the existence of a linear relationship between the output voltage and the imposed displacement for almost all cases. Furthermore, the ability of these IPMCs as large motion actuators and robotic manipulators is presented. Several muscle configurations are constructed to demonstrate the capabilities of these IPMC actuators. This paper further identifies key parameters involving the vibrational and resonance characteristics of sensors and actuators made with IPMCS. When the applied signal frequency varies, so does the displacement up to a critical frequency called the resonant frequency where maximum deformation is observed, beyond which the actuator response is diminished. A data acquisition system was used to measure the parameters involved and record the results in real time basis. Also the load characterizations of the IPMCs were measured and it was shown that these actuators exhibit good force to weight characteristics in the presence of low applied voltages. Finally reported are the cryogenic properties of these muscles for potential utilization in an outer space environment of a few Torrs and temperatures of the order of - 140 degrees Celsius. These muscles are shown to work quite well in such harsh cryogenic environments and thus present a great potential as sensors and actuators that can operate at cryogenic temperatures.

Shahinpoor, M.↗
NASA NTRS2003

Kalman Filtering with Inequality Constraints for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops two analytic methods of incorporating state variable inequality constraints in the Kalman filter. The first method is a general technique of using hard constraints to enforce inequalities on the state variable estimates. The resultant filter is a combination of a standard Kalman filter and a quadratic programming problem. The second method uses soft constraints to estimate state variables that are known to vary slowly with time. (Soft constraints are constraints that are required to be approximately satisfied rather than exactly satisfied.) The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is proven theoretically and shown via simulation results. The use of the algorithm is demonstrated on a linearized simulation of a turbofan engine to estimate health parameters. The turbofan engine model contains 16 state variables, 12 measurements, and 8 component health parameters. It is shown that the new algorithms provide improved performance in this example over unconstrained Kalman filtering.

Simon, Dan↗
NASA NTRS1979

Further comments on the effects of vacuum birefringence on the polarization of X-rays emitted from magnetic neutron stars

The birefringence of the vacuum in the presence of strong (of the order of 1 teragauss) magnetic fields will in general affect the polarization of X-rays propagating through these fields. Two of the four Stokes parameters will vary so rapidly with wavelength as to be 'washed out' and unobservable, but the remaining two parameters will be unaffected. These results show that one conclusion of an earlier work is incorrect: Polarized X-ray emission from the surface of a magnetic neutron star will not in general be completely depolarized by the effects of vacuum birefringence. In particular, this birefringence has no effect on the linear polarization of cyclotron emission from the poles of magnetic neutron stars, and a similar result holds for synchrotron emission. More general cases of the propagation of polarized X-rays in magnetic fields are also discussed.

Chanan, G. A.↗
NASA NTRS2016

A Quasi-Global Approach to Improve Day-Time Satellite Surface Soil Moisture Anomalies through the Land Surface Temperature Input

Passive microwave observations from various spaceborne sensors have been linked to the soil moisture of the Earth's surface layer. A new generation of passive microwave sensors are dedicated to retrieving this variable and make observations in the single theoretically optimal L-band frequency (1-2 GHz). Previous generations of passive microwave sensors made observations in a range of higher frequencies, allowing for simultaneous estimation of additional variables required for solving the radiative transfer equation. One of these additional variables is land surface temperature, which plays a unique role in the radiative transfer equation and has an influence on the final quality of retrieved soil moisture anomalies. This study presents an optimization procedure for soil moisture retrievals through a quasi-global precipitation-based verification technique, the so-called Rvalue metric. Various land surface temperature scenarios were evaluated in which biases were added to an existing linear regression, specifically focusing on improving the skills to capture the temporal variability of soil moisture. We focus on the relative quality of the day-time (01:30 pm) observations from the Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E), as these are theoretically most challenging due to the thermal equilibrium theory, and existing studies indicate that larger improvements are possible for these observations compared to their night-time (01:30 am) equivalent. Soil moisture data used in this study were retrieved through the Land Parameter Retrieval Model (LPRM), and in line with theory, both satellite paths show a unique and distinct degradation as a function of vegetation density. Both the ascending (01:30 pm) and descending (01:30 am) paths of the publicly available and widely used AMSR-E LPRM soil moisture products were used for benchmarking purposes. Several scenarios were employed in which the land surface temperature input for the radiative transfer was varied by imposing a bias on an existing regression. These scenarios were evaluated through the Rvalue technique, resulting in optimal bias values on top of this regression. In a next step, these optimal bias values were incorporated in order to re-calibrate the existing linear regression, resulting in a quasi-global uniform LST relation for day-time observations. In a final step, day-time soil moisture retrievals using the re-calibrated land surface temperature relation were again validated through the Rvalue technique. Results indicate an average increasing Rvalue of 16.5%, which indicates a better performance obtained through the re-calibration. This number was confirmed through an independent Triple Collocation verification over the same domain, demonstrating an average root mean square error reduction of 15.3%. Furthermore, a comparison against an extensive in situ database (679 stations) also indicates a generally higher quality for the re-calibrated dataset. Besides the improved day-time dataset, this study furthermore provides insights on the relative quality of soil moisture retrieved from AMSR-E's day- and night-time observations.

soil moisture↗
NASA NTRS1992

Moment method analysis of linearly tapered slot antennas: Low loss components for switched beam radiometers

A Moment Method Model for the radiation pattern characterization of single Linearly Tapered Slot Antennas (LTSA) in air or on a dielectric substrate is developed. This characterization consists of: (1) finding the radiated far-fields of the antenna; (2) determining the E-Plane and H-Plane beamwidths and sidelobe levels; and (3) determining the D-Plane beamwidth and cross polarization levels, as antenna parameters length, height, taper angle, substrate thickness, and the relative substrate permittivity vary. The LTSA geometry does not lend itself to analytical solution with the given parameter ranges. Therefore, a computer modeling scheme and a code are necessary to analyze the problem. This necessity imposes some further objectives or requirements on the solution method (modeling) and tool (computer code). These may be listed as follows: (1) a good approximation to the real antenna geometry; and (2) feasible computer storage and time requirements. According to these requirements, the work is concentrated on the development of efficient modeling schemes for these type of problems and on reducing the central processing unit (CPU) time required from the computer code. A Method of Moments (MoM) code is developed for the analysis of LTSA's within the parameter ranges given.

Koeksal, Adnan↗
NASA NTRS1978

Study of the effect of static/dynamic Coulomb friction variation at the tape-head interface of a spacecraft tape recorder by non-linear time response simulation

A description is presented of six simulation cases investigating the effect of the variation of static-dynamic Coulomb friction on servo system stability/performance. The upper and lower levels of dynamic Coulomb friction which allowed operation within requirements were determined roughly to be three times and 50% respectively of nominal values considered in a table. A useful application for the nonlinear time response simulation is the sensitivity analysis of final hardware design with respect to such system parameters as cannot be varied realistically or easily in the actual hardware. Parameters of the static/dynamic Coulomb friction fall in this category.

Mukhopadhyay, A. K.↗
NASA NTRS1989

Scattering of light by stochastically rough particles

The single particle phase function and the linear polarization for large stochastically deformed spheres have been calculated by Monte Carlo simulation using the geometrical optics approximation. The radius vector of a particle is assumed to obey a bivariate lognormal distribution with three free parameters: mean radius, its standard deviation and the coherence length of the autocorrelation function. All reflections/refractions which include sufficient energy have been included. Real and imaginary parts of the refractive index can be varied without any restrictions. Results and comparisons with some earlier less general theories are presented. Applications of this theory to the photometric properties of atmosphereless bodies and interplanetary dust are discussed.

Peltoniemi, Jouni I.↗
NASA NTRS1993

The onset of electrohydrodynamic instability in isoelectric focusing

The onset of electrohydrodynamic motion associated with the imposition of an electric field across a thin layer of liquid has been investigated for the case in which the electrical conductivity varies linearly over the depth of the layer. The variation of the conductivity is due to concentration gradients in the charge-carrying solutes and its spatio-temporal evolution is represented by a convective-diffusion equation. When the viscous relaxation time is long compared to the time for charge relaxation, the analysis reveals that the neutral stability curves for the layer can be characterized by three dimensionless parameters: Ra(sub e) is equivalent to d(epsilon)E(sup 2, sub 0) Delta sigma/mu K(sub eff) sigma(sub 0), an electrical Rayleigh number; delta sigma/sigma(sub 0), the relative conductivity increment; and alpha, the transverse wave number of the disturbance. Here d is the thickness, epsilon is the dielectric constant, and mu is the viscosity of layer, E(sub 0) is the applied field strength at the lower conductivity boundary, and K(sub eff) is an effective diffusivity associated with the Brownian motion of the charge-carrying solutes. With viscous-stress-free boundaries, at which the electrical conductivity and the normal component of the electric field are prescribed, the critical Ra(sub e) is 1.504 x 10(exp 4) at a critical transverse wave number of 1.97 when Delta sigma/sigma(sub 0) is 10. As Delta sigma/sigma(sub 0) increases, the critical Ra(sub e) increases and shifts to shorter wavelength disturbances; the critical imposed field strength, however, passes through a minimum because the lower-conductivity boundary exerts a considerable stabilizing influence in the presence of steep conductivity gradients. Similar trends were obtained for liquid layers with rigid boundaries.

Baygents, James C.↗
NASA NTRS2004

A 3D Convective Model for the Jovian Wind Bands

In an earlier paper, we proposed that Jupiter's alternating wind bands are a manifestation of the global interaction between rotation and convection in a shallow layer. The model, however, was obtained from linearization of the 2D equations of motions. At HKUST/Hong Kong, we are now trying to study this problem by rigorous numerical simulation. Using a three-dimensional spectral numerical code, we compute models for the outermost layer of Jupiter's convective envelope. Two cases have been studied. In one the atmospheric pressure varies from 1 to 23 bar, and in the other from 1 to 115 bar. The physical parameters (internal energy flux, rotation rate) are chosen to be close to those expected, but solar heating, chemistry, as well as dynamical influences from deeper layers are ignored. The models generate wind field patterns that contain alternating jet streams with resemblance to the Jovian bands. Instantaneous values of the mean zonal flow at the equator reach 80 m/sec. Yet the mean meridional flows are less than 1% of such value. The meridional temperature profile at the cloud top level also shows a double hump structure of a few degrees (as observed) in the subtropics. Though there is not complete quantitative agreement (caused perhaps by neglected effects like solar radiation), these models demonstrate, in principle, the feasibility of generating a Jovian type wind pattern through the interaction of fast rotation and convection in a thin shell.

Mayr, H. G.↗
DOE OSTI2022

Thermodynamics of the Inner Heliosheath

We derive annual sky maps of the proton temperature in the inner heliosheath (IHS), and track their temporal evolution over the years 2009–2016 of Interstellar Boundary Explorer observations. Other associated thermodynamic parameters also determined are the density, kappa (the parameter that characterizes kappa distributions), temperature rate, polytropic index, and entropy. We exploit the theory of kappa distributions and their connection with polytropes, to (i) express a new polytropic quantity Π that remains invariant along streamlines where temperature and density may vary, (ii) parameterize the proton flux in terms of the Π invariant and kappa, and (iii) derive the temperature and density, respectively, from the slope and intercept of the linear relationship between kappa and logarithm of Π. We find the following thermodynamic characteristics: (1) temperature sky maps and histograms shifted to their lowest values in 2012 and their highest in 2015; (2) temperature negatively correlated with density, reflecting the subisothermal polytropic behavior; (3) temperature positively correlated with kappa, revealing characteristics of the mechanism responsible for generating kappa distributions; (4) processes in IHS are subisothermal tending toward isobaric, consistent with previously published results; (5) linear relationship between kappa and polytropic indices, revealing characteristics of the particle potential energy; and (6) entropy positively correlated with polytropic index, aligned with the underlying theory that entropy increases toward the isothermal state where the kappa distribution reduces to the Maxwell–Boltzmann description.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Error assessments of widely-used orbit error approximations in satellite altimetry

From simulations, the orbit error can be assumed to be a slowly varying sine wave with a predominant wavelength comparable to the Earth's circumference. Thus, one can derive analytically the error committed in representing the orbit error along a segment of the satellite ground track by a bias; by a bias and tilt (linear approximation); or by a bias, tilt, and curvature (quadratic approximation). The result clearly agrees with what is obvious intuitively, i.e., (1) the fit is better with more parameters, and (2) as the length of the segment increases, the approximation gets worse. But more importantly, it provides a quantitative basis to evaluate the accuracy of past results and, in the future, to select the best approximation according to the required precision and the efficiency of various approximations.

Tai, Chang-Kou↗