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 271 records · Page 15

Turbulent flow in a channel with a wall with progressive waves

The effects of prescribed wall motion on turbulent channel flow were examined with the objective of understanding the drag mechanisms and obtaining a possible means for drag reduction. A computer program for large-eddy simulation (LES) of turbulent flow in a channel was adapted to treat the case of time-varying wall motion using linearized boundary conditions. The code was applied to flow in a channel with prescribed wall motion on one wall. It was found that the results can be explained by a model based on a Stokes-like layer near the moving wall. This provides a theory which can be used to guide the selection of parameters. Considering both viscous and pressure drag components, it was found that net drag reduction is possible under certain conditions of wall motion. However, the work required to drive the wall motion can be greater than the reduced through-flow energy requirement. The major effects of the wall motion are restricted to the viscous sublayer of the turbulent flow.

Kuhn, G. D.↗

The role of shock induced trailing-edge separation in limit cycle oscillations

The potential role of shock induced trailing edge separation (SITES) in limit cycle oscillations (LCO) was established. It was shown that the flip-flop characteristics of transition to and from SITES as well as its hysteresis could couple with wing modes with torsional motion and low damping. This connection led to the formulation of a very simple nonlinear math model using the linear equations of motion with a nonlinear step forcing function with hysteresis. A finite difference solution with time was developed and calculations were made for the F-111 TACT were used to determine the step forcing function due to SITES transition. Since no data were available for the hysteresis, a parameter study was conducted allowing the hysteresis effect to vary. Very small hysteresis effects, which were within expected bounds, were required to obtain reasonable response levels that essentially agreed with flight test results. Also in agreement with wind tunnel tests, LCO calculations for the 1/6 scale F-111 model showed that the model should have not experienced LCO.

Cunningham, Atlee M., Jr.↗

Flashback Studies of High-Hydrogen Flames using High-Speed OH Planar Laser-Induced Fluorescence

Flame flashback is a crucial challenge during the application of high-hydrogen fuels for the development of next-generation lean premixed combustors. The fundamentals of stable flame configuration, flashback phenomena, and spatiotemporal information of flame dynamics in a premixed, swirl-stabilized burner are investigated experimentally at atmospheric pressure conditions using high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) measurements. The inlet conditions are varied systematically with respect to equivalence ratio (ϕ), hydrogen percentage in a hydrogen-methane mixture (50%–90%) and pre-mixer velocities and the experiments are performed on the burner modified to provide optical access to the premixing section. Signal interferences such as saturation effects, background chemiluminescence, and laser non-uniformity issues were minimized or corrected during data acquisition and subsequent data processing steps. A stable flame characterization was performed at 20 kHz-repetition-rate for varied pre-mixer velocities, hydrogen content and ϕ in the reactant mixture and lift-off length (L) investigation showed dependence on each of these parameters. Flame dynamics such as flame flashback events, local ignition, local extinction, flame curvature as well as the main reaction zone of the combustion mechanism are well characterized. Flashback ϕ showed an expected linearly increasing trend with increasing pre-mixer velocity and decreasing hydrogen content and the conclusions drawn aligned well with detailed L investigation. Such spatially and temporally resolved kHz-rate OH-PLIF is a promising technique to observe rapidly occurring flashback dynamics and can further validate turbulence-chemistry interaction models of swirling flames.

Parajuli, Pradeep↗

Study of Quadrupole Fringe Fields in the Interaction Region of the Hadron Storage Ring of the Electron Ion Collider

Fringe fields in quadrupole magnets are usually neglected in studies of beam dynamics at accelerators. However, the extreme optical parameters present in the final focus of a collider such as the Electron–Ion Collider (EIC) may give rise to effects that should not be overlooked. The calculation of quadrupole fringe fields presented in this study follows the procedure outlined in Ref. [1], specialized to the case of a straight reference orbit (i.e., with no dipole field component). A right-handed Cartesian coordinate system is employed, with the $z$-axis aligned with the quadrupole axis and $x$ and $y$ denoting the horizontal and vertical transverse coordinates, respectively. The magnetic quadrupole field gradient, $G = \partial B_y / \partial x$, transitions from its peak value inside the quadrupole—where it is nearly independent of the longitudinal coordinate $z$—to zero at some distance beyond the magnet edge. Consequently, $G$ is treated as a function of $z$. The region over which this variation occurs is defined as the quadrupole fringe field region. The study begins with the development of a description of the magnetic field in the fringe region using a power-series expansion in the transverse coordinates $x$and $y$, consistent with the longitudinally varying gradient. A model for the $z$-dependence of the gradient is then proposed and adjusted to reproduce magnetic field data obtained from three-dimensional field calculations. To evaluate the impact of the fringe fields on beam dynamics, the corresponding vector potential is derived and incorporated into a Hamiltonian formulation of particle motion. The significance of the fringe fields is quantified by calculating the amplitude-dependent tune shift from the Hamiltonian. Using linear beam optics parameters of the Hadron Storage Ring (HSR) of the EIC, the tune shift due to the fringe fields of all quadrupole magnets in the IR-6 interaction region is evaluated. Finally, the resulting tune shifts are compared with those arising from other nonlinear field components present in the HSR.

43 PARTICLE ACCELERATORS↗

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New Directions in Asymptotically Stable Finite-dimensional Adaptive Control of Linear Distributed Parameter Systems

Distributed Parameter Systems (DPS), such as systems described by partial differential equations, require infinite-dimensional state space descriptions to correctly model their dynamical behavior. However, any adaptive control algorithm must be finite-dimensional in order to be implemented via on-line digital computers. Finite-dimensional adaptive control of linear DPS requires stability analysis of nonlinear, time-varying, infinite-dimensional systems. The structure of nonadaptive finite-dimensional control of linear DPS is summarized as it relates to the existence of limiting systems for adaptive control. Two candidate schemes for finite-dimensional adaptive control of DPS are described and critical issues in infinite-dimensional stability analysis are discussed, in particular, the invariance principle, center manifold theory, and relationships between input-output and internal stability.

Balas, M. J.↗

Effects of Elastic Edge Restraints and Initial Prestress on the Buckling Response of Compression-Loaded Composite Panels

A parametric study of the effects of test-fixture-induced initial prestress and elastic edge restraints on the prebuckling and buckling responses of a compression-loaded, quasi-isotropic curved panel is presented. The numerical results were obtained by using a geometrically nonlinear finite element analysis code with high-fidelity models. The results presented show that a wide range of prebuckling and buckling behavior can be obtained by varying parameters that represent circumferential loaded-edge restraint and rotational unloaded-edge restraint provided by a test fixture and that represent the mismatch in specimen and test-fixture radii of curvature. For a certain range of parameters, the panels exhibit substantial nonlinear prebuckling deformations that yield buckling loads nearly twice the corresponding buckling load predicted by a traditional linear bifurcation buckling analysis for shallow curved panels. In contrast, the results show another range of parameters exist for which the nonlinear prebuckling deformations either do not exist or are relatively benign, and the panels exhibit buckling loads that are nearly equal to the corresponding linear bifurcation buckling load. Overall, the results should also be of particular interest to scientists, engineers, and designers involved in simulating flight-hardware boundary conditions in structural verification and certification tests, involved in validating structural analysis tools, and interested in tailoring buckling performance.

Hilburger, Mark W.↗

Compression-Loaded Composite Panels With Elastic Edge Restraints and Initial Prestress

A parametric study of the effects of test-fixture-induced initial prestress and elastic edge restraints on the prebuckling and buckling responses of a compression-loaded, quasi-isotropic curved panel is presented. The numerical results were obtained by using a geometrically nonlinear finite element analysis code with high-fidelity models. The results presented show that a wide range of prebuckling and buckling behavior can be obtained by varying parameters that represent circumferential loaded-edge restraint and rotational unloaded-edge restraint provided by a test fixture and that represent the mismatch in specimen and test-fixture radii of curvature. For a certain range of parameters, the panels exhibit substantial nonlinear prebuckling deformations that yield buckling loads nearly twice the corresponding buckling load predicted by a traditional linear bifurcation buckling analysis for shallow curved panels. In contrast, the results show another range of parameters exist for which the nonlinear prebuckling deformations either do not exist or are relatively benign, and the panels exhibit buckling loads that are nearly equal to the corresponding linear bifurcation buckling load. Overall, the results should be of particular interest to scientists, engineers, and designers involved in simulating flight-hardware boundary conditions in structural verification and certification tests, involved in validating structural analysis tools, and interested in tailoring buckling performance.

Hilburger, Mark W.↗

Canopy reflectance modelling of semiarid vegetation

Three different types of remote sensing algorithms for estimating vegetation amount and other land surface biophysical parameters were tested for semiarid environments. These included statistical linear models, the Li-Strahler geometric-optical canopy model, and linear spectral mixture analysis. The two study areas were the National Science Foundation's Jornada Long Term Ecological Research site near Las Cruces, NM, in the northern Chihuahuan desert, and the HAPEX-Sahel site near Niamey, Niger, in West Africa, comprising semiarid rangeland and subtropical crop land. The statistical approach (simple and multiple regression) resulted in high correlations between SPOT satellite spectral reflectance and shrub and grass cover, although these correlations varied with the spatial scale of aggregation of the measurements. The Li-Strahler model produced estimated of shrub size and density for both study sites with large standard errors. In the Jornada, the estimates were accurate enough to be useful for characterizing structural differences among three shrub strata. In Niger, the range of shrub cover and size in short-fallow shrublands is so low that the necessity of spatially distributed estimation of shrub size and density is questionable. Spectral mixture analysis of multiscale, multitemporal, multispectral radiometer data and imagery for Niger showed a positive relationship between fractions of spectral endmembers and surface parameters of interest including soil cover, vegetation cover, and leaf area index.

Franklin, Janet↗

Empirical Model of Puncture Energy for Metals

The purpose of this work is to fit a previously developed empirical equation for puncture energy to simulation data. The conservative puncture energy equation could be used to expedite the process of performing calculations in the development of safety measures, avoiding the need to create complex finite element models for specific puncture scenarios. A total of 108 simulations are developed by varying coupon thickness, coupon material, probe shape, and probe diameter. The simulations are comprised of a low-velocity probe puncturing the coupons, from which the probe kinetic energy change is calculated. The empirical equation is fit to the dimensions, material properties, and energy results using a non-linear least-squares regression method within Python, which determines the two constant parameters for each fit. More statistically significant fit results are achieved by separating the data by probe shape and coupon material.

42 ENGINEERING↗

Lidar observations of banded convection during BLX83

Lidar observations of clear-air convection during the 1983 Boundary Layer Experiment (BLX83) reveal the presence of elongated, parallel regions of updrafts marked by enhanced aerosol backscattering. These linear (banded) aerosol structures were observed over a 2-hr period during a cloud-free morning. During this period, the depth of the convective boundary layer (CBL) increased from 100 to 1300 m. Wind speeds averaged over the depth of the CBL varied between 0 and 2 m/s, while the wind direction varied over a range of 160 deg. The CBL instability parameter, -Zi/L, increased from approximately 25 (weakly unstable) to 250 (strongly unstable). The spacings of the elongated, parallel plumes scaled with the CBL height. These findings suggest that secondary circulations in the form of horizontal roll vortices were present under conditions not normally associated with roll vortices. The lines of aerosol structures aligned much more closely (within 15 deg) with the direction of the vertical shear to the horizontal wind through the depth of the CBL than with either the surface wind, mean CBL wind, or the wind at an altitude of 1.1 Zi.

Ferrare, R. A.↗

Rapid time variations and radiative instabilities of astrophysical masers

The time-varying intensities are obtained for astrophysical masers that are radiatively unstable. Numerical integrations of the time-dependent, nonlinear equations of radiative transfer are performed with the usual approximation of a linear maser. At long times after changes in the physical conditions, the intensity of maser radiation reaches an asymptotic behavior and oscillates permanently in these idealized calculations with a period that is related to the length of the maser divided by the speed of light. The intensity varies by more than a factor of 10. These intensities depend upon the same four parameters as we originally found to determine the regime for radiative instabilities based on a stability analysis of the steady state. A detailed comparison is made between the predictions of the stability analysis and the time variations of the intensities. Calculations are performed for interacting pairs as well as isolated, individual masers.

Scappaticci, Gerardo A.↗

Radar Scattering and Block Size Properties of Lunar Crater Ejecta From Mini-RF and LROC NAC Data

A major objective of the Mini-RF experiment is to distinguish lunar surfaces that may contain water/ice deposits [1,2]. Better understanding of the backscattering properties of craters of varying age and size is crucial for interpreting data received from the Mini-RF. The Mini-RF transmits a circularly polarized RF electromagnetic energy and coherently receives orthogonal linear polarization echoes [1]. The Mini- RF maps in two separate bands ( =12.6 and 4.5 cm) at a high resolution mode of 30 m/pixel [1]. Given the variables mentioned, the four stokes parameters are reconstructed. The Circular Polarization Ratio (CPR) is calculated for the purposes of understanding subsurface and surface roughness. The CPR is determined from reflections acquired from the ratio of power of the transmitted radio wave in same sense to the reflected radio wave in the opposite sense [1]. Ice in the permanently shadowed regions (PSRs) would be transparent to radar, but the inclusions of materials and imperfections would cause the radio wave to reflect multiple times [3], enhancing the number of same sense reflections and increasing the CPR. In addition, ice also displays the coherent backscatter opposition effect (CBOE), an interferrometric addition of same sense backscatter that further increases the CPR of ice targets [7]. High CPR values also correlate to multiple reflections and are typically associated with very rough surfaces [3]. The average dry lunar surface has a CPR in the range of 0.2-0.4 at 48deg incidence [3]. The purpose of this study is to begin to quantify degrees of surface wavelength-scale roughness with CPR and to understand how such surface roughness is created and gradually destroyed by erosion on the lunar surface. Another goal is to identify and isolate the possible causes of high CPR within the shadowed areas of anomalous polar craters [3]. All the studied craters are non-polar, so that we can see into their interiors in NAC images. The idea is to understand what controls blockiness in these craters so that we can rule out rocks (and rule in ice) for the anomalous polar dark ones [3].

Spudis, P. D.↗

Development and Integration of an Advanced Stirling Convertor Linear Alternator Model for a Tool Simulating Convertor Performance and Creating Phasor Diagrams

A simple model of the Advanced Stirling Convertors (ASC) linear alternator and an AC bus controller has been developed and combined with a previously developed thermodynamic model of the convertor for a more complete simulation and analysis of the system performance. The model was developed using Sage, a 1-D thermodynamic modeling program that now includes electro-magnetic components. The convertor, consisting of a free-piston Stirling engine combined with a linear alternator, has sufficiently sinusoidal steady-state behavior to allow for phasor analysis of the forces and voltages acting in the system. A MATLAB graphical user interface (GUI) has been developed to interface with the Sage software for simplified use of the ASC model, calculation of forces, and automated creation of phasor diagrams. The GUI allows the user to vary convertor parameters while fixing different input or output parameters and observe the effect on the phasor diagrams or system performance. The new ASC model and GUI help create a better understanding of the relationship between the electrical component voltages and mechanical forces. This allows better insight into the overall convertor dynamics and performance.

Metscher, Jonathan F.↗

Development and Integration of an Advanced Stirling Convertor Linear Alternator Model for a Tool Simulating Convertor Performance and Creating Phasor Diagrams

A simple model of the Advanced Stirling Convertors (ASC) linear alternator and an AC bus controller has been developed and combined with a previously developed thermodynamic model of the convertor for a more complete simulation and analysis of the system performance. The model was developed using Sage, a 1-D thermodynamic modeling program that now includes electro-magnetic components. The convertor, consisting of a free-piston Stirling engine combined with a linear alternator, has sufficiently sinusoidal steady-state behavior to allow for phasor analysis of the forces and voltages acting in the system. A MATLAB graphical user interface (GUI) has been developed to interface with the Sage software for simplified use of the ASC model, calculation of forces, and automated creation of phasor diagrams. The GUI allows the user to vary convertor parameters while fixing different input or output parameters and observe the effect on the phasor diagrams or system performance. The new ASC model and GUI help create a better understanding of the relationship between the electrical component voltages and mechanical forces. This allows better insight into the overall convertor dynamics and performance.

Metscher, Jonathan F.↗

A Bayesian approach to tracking patients having changing pharmacokinetic parameters

This paper considers the updating of Bayesian posterior densities for pharmacokinetic models associated with patients having changing parameter values. For estimation purposes it is proposed to use the Interacting Multiple Model (IMM) estimation algorithm, which is currently a popular algorithm in the aerospace community for tracking maneuvering targets. The IMM algorithm is described, and compared to the multiple model (MM) and Maximum A-Posteriori (MAP) Bayesian estimation methods, which are presently used for posterior updating when pharmacokinetic parameters do not change. Both the MM and MAP Bayesian estimation methods are used in their sequential forms, to facilitate tracking of changing parameters. Results indicate that the IMM algorithm is well suited for tracking time-varying pharmacokinetic parameters in acutely ill and unstable patients, incurring only about half of the integrated error compared to the sequential MM and MAP methods on the same example.

Bayes Theorem↗

Failure of the method of slowly varying amplitude and phase for non-linear, singular oscillators

It is shown that the method of slowly varying amplitude and phase yields erroneous results in the study of the mathematical properties of nonlinear singular oscillator systems. The analytical solution is described in which the phase function is constant and for which a special limiting behavior exists when the wavelength is zero. The previous method based on the condition of boundedness cannot be satisfied for nonlinear singular characteristics, and the erroneous designation of the expansion parameter is identified.

Mickens, R. E.↗

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.↗