Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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

Alternate Methods of Model Reduction to Avoid Dynamic Modal Truncation Error

Loads analysis is traditionally performed using dynamically reduced models, which provides the benefit to reduce run time. If the reduced model frequency cutoff is not chosen appropriately, the model will lack the dynamic content required to fully represent the response of the non-reduced model. Standard guidance for reduced model frequency content, provided in NASA-STD-5002, is to solve fixed base modes up to a minimum of 1.5x the model frequency content of interest and to employ static modal truncation methods such as residual vectors, the mode acceleration method, and the residual flexibility method to account for the truncated flexibility of the missing modes. Fixed base modes require the selection of a set of degrees of freedom to be constrained which, if not properly selected, may affect the accuracy of the reduced model by excluding some of the dynamic characteristic of the full model. In this case, the standard NASA guidance would be insufficient, but it may not be readily apparent that a portion of the reduced model response is missing. This error was encountered during an independent verification and validation (IV&V) effort, where it was observed that the resulting dynamic response was lower than the inline analysis. In this specific case, despite following the standard NASA model reduction guidelines in the selection of the frequency cutoff, the inline model still did not fully capture the necessary dynamic content. As part of the IV&V, an alternate reduction methodology was employed using an unconstrained mode acceleration method. The original model initially performed a constrained reduction to twice the frequency content of interest before doing a free-free run, employing the mode acceleration method to account for the truncated modes. In contrast in the IV&V, the unconstrained model reduced directly to the needed frequency content of the free-free run, avoiding any interactions between constraints and dynamic content. To verify the model, the original reduction methodology was used to generate a series of Hurty-Craig-Bampton reductions, each with a higher frequency cutoff than the previous. The results were shown to converge once the frequency cutoff increased past eight times the frequency content of interest. At the same frequency cutoff, the results of the Hurty-Craig-Bampton model converged with the results of the unconstrained mode acceleration model. This comparative study provided confidence that the results of the unconstrained modal acceleration reduced model were correct and that the Hurty-Craig-Bampton needed to increase its frequency cutoff to fully capture the dynamic response.

Erin Simmons↗

Comptonization in supercritical winds. I - Spectral evolution. II - Dynamics and observational diagnostics

The physics governing spectral evolution in steady, fully ionized, radiatively driven winds from compact objects is examined. A largely analytic theory is developed to describe the effects of diffusion, adiabatic cooling, and repeated Compton scattering on the emergent radiation spectrum. Analytic forms for the wind temperature profile are derived, and spectral solutions obtained by considering a modified, time-independent Kompaneets equation are presented subject to a variety of initial conditions. A probabilistic formalism which describes the effect of diffusion photons travelling beyond the trapping radius is developed to allow the calculation of spectra observed at infinity. The second paper focuses on the dynamics of supercritical winds and on observational diagnostics which may be used to constrain model parameters in the contexts of specific astrophysical systems, including X-ray burst sources and mass-losing quasars.

Becker, Peter A.↗

The dust coma of Comet Austin (1989c1)

Thermal-infrared (10 and 20 micron) images of Comet Austin were obtained on UT 30.6 Apr., 1.8, 2.8, and 3.6 May 1990. The NASA-Marshall Space Flight Center 20 pixel bolometer array at the NASA 3 meter Infrared Telescope Facility in Hawaii was used. The 10.8 micron (FWHM = 5.3 microns) maps were obtained with maximum dimensions of 113 arcsec (57,500 km) in RA and 45 arcsec (23,000 km) in declination, with a pixel size of 4.2 x 4.2 arcsec. A smaller, 45 x 18 arcsec, map was obtained in the 19.2 micron (FWHM = 5.2 microns) bandpass. At the time of these observations Comet Austin's heliocentric and geocentric distances were 0.7 and 0.5 AU respectively. The peak flux density (within the brightest pixel) was 23 + or - 2 Janskys for the first three dates and only marginally lower the last day; i.e., within the observational uncertainties no evidence was found for day-to-day variability like that observed in Comet Halley. A dynamical analysis of the morphology of the extended dust emission is used to constrain the size distribution and production rate of the dust particles. The results of this analysis are compared with similar studies carried out on comets P/Giacobini-Zinner, P/Brorsen-Metcalf, P/Halley, P/Tempel 2, and Wilson (1987).

Campins, Humberto↗

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗

Autolanding Trajectory Design for the X-34

An Autolanding I-load Program (ALIP) is developed to design unpowered autolanding trajectories for the X-34 Mach 8 vehicle. The trajectory is comprised of geometric flight segments that are based on the shuttle approach and landing design (steep glideslope, circular flare, exponential flare to shallow glideslope). Enforcing physical constraints such as loads, vertical descent rate, continuity and smoothness reduces the design problem to a two point boundary value problem with conditions on the initial and final dynamic pressure. Finding a solution required the development of trajectory simulation techniques that constrained the flight profile to a prescribed geometry. The design methodology can be extended beyond the autolanding flight regime by repeating the series of geometric segments and solving multiple two-point boundary values problems (one for each series). The techniques described in this paper facilitate rapid design of reference trajectories.

Barton, Gregg H.↗

Ionospheric simulation compared with Dynamics Explorer observations for November 22, 1981

The present study uses an extensive DE-2 data base to both constrain inputs to a time-dependent ionospheric model (TDIM) for a simulation of the ionosphere and then check the simulated densities. The investigation was carried out for both a quiet period and a substorm period. The quiet-day study produced very good agreement between modeled and observed electron densities in the topside ionosphere with two significant exceptions: First, across the polar region the DE-2 LANG densities showed fine structure in addition to the overall regional density morphology. Second, a surprising discrepancy arose in the presunrise and midlatitude trough. The TDIM densities were an order of magnitude lower than those observed by DE-2. The substorm study showed remarkably good agreement with the observed densities.

Sojka, J. J.↗

Variety in the Variability of Accreting Supermassive Binary Black Holes

Accreting supermassive binary black holes are key multi-messenger sources for LISA, yet are challenging to simulate realistically as solving the radiation magnetohydrodynamics (MHD) equations over the full dynamic spatio-temporal range of the problem is computationally infeasible at present. We report on our collaboration's progress simulating these systems using general relativistic MHD simulations and dynamic GR. In order to cover a larger temporal range in one set of simulations, we constrain our view to the circumbinary disk region and measure the effects binary mass ratio, accretion disk size, and black hole spin have on the structure and variability of the accretion flow. We particularly emphasize how these parameters influence the overdensity feature, which orbits the binary near the edge of the cavity, since it is responsible for most of the electromagnetic emission's variability and variability is a key signature of a system being a binary. Extending to smaller length scales, we will report on simulations following accretion all the down to the event horizons so that we can begin to investigate how black hole spin affects mini-disk dynamics, accretion rate, and jet power.

Scott C Noble↗

Variety in the Variability of Accreting Supermassive Binary Black Holes

Accreting supermassive binary black holes are key multi-messenger sources for LISA, yet are challenging to simulate realistically as solving the radiation magnetohydrodynamics (MHD) equations over the full dynamic spatio-temporal range of the problem is computationally infeasible at present. We will provide a brief summary of the the progress made in the field to understand these systems theoretically and what new directions groups are pursuing. We will also report on our collaboration's progress to simulate these systems using general relativistic MHD and dynamic GR. In order to cover a larger temporal range in one set of simulations, we constrain our view to the circumbinary disk region and measure how the binary mass ratio, accretion disk size, and black hole spin have on the structure and variability of the accretion flow. We particularly emphasize how these parameters influence the overdensity feature, which orbits the binary near the edge of the cavity, since it is responsible for most of the electromagnetic emission's variability and variability is a key signature of a system being a binary. Extending to smaller length scales, we will report on simulations following accretion all the down to the event horizons so that we can begin to investigate how black hole spin affects mini-disk dynamics, accretion rate, and jet power. The novel computational methods enabling inclusion of the black holes in the domain, including multi-patch methods, will be described.

Scott C Noble↗

Methods of sequential estimation for determining initial data in numerical weather prediction

Numerical weather prediction (NWP) is an initial-value problem for a system of nonlinear differential equations, in which initial values are known incompletely and inaccurately. Observational data available at the initial time must therefore be supplemented by data available prior to the initial time, a problem known as meteorological data assimilation. A further complication in NWP is that solutions of the governing equations evolve on two different time scales, a fast one and a slow one, whereas fast scale motions in the atmosphere are not reliably observed. This leads to the so called initialization problem: initial values must be constrained to result in a slowly evolving forecast. The theory of estimation of stochastic dynamic systems provides a natural approach to such problems. For linear stochastic dynamic models, the Kalman-Bucy (KB) sequential filter is the optimal data assimilation method, for linear models, the optimal combined data assimilation-initialization method is a modified version of the KB filter.

Cohn, S. E.↗

Nonlinear feedback control of multiple robot arms

Multiple coordinated robot arms are modeled by considering the arms: (1) as closed kinematic chains, and (2) as a force constrained mechanical system working on the same object simultaneously. In both formulations a new dynamic control method is discussed. It is based on a feedback linearization and simultaneous output decoupling technique. Applying a nonlinear feedback and a nonlinear coordinate transformation, the complicated model of the multiple robot arms in either formulation is converted into a linear and output decoupled system. The linear system control theory and optimal control theory are used to design robust controllers in the task space. The first formulation has the advantage of automatically handling the coordination and load distribution among the robot arms. In the second formulation, by choosing a general output equation, researchers can superimpose the position and velocity error feedback with the force-torque error feedback in the task space simultaneously.

Tarn, T. J.↗

New nonlinear control algorithms for multiple robot arms

Multiple coordinated robot arms are modeled by considering the arms as closed kinematic chains and as a force-constrained mechanical system working on the same object simultaneously. In both formulations, a novel dynamic control method is discussed. It is based on feedback linearization and simultaneous output decoupling technique. By applying a nonlinear feedback and a nonlinear coordinate transformation, the complicated model of the multiple robot arms in either formulation is converted into a linear and output decoupled system. The linear system control theory and optimal control theory are used to design robust controllers in the task space. The first formulation has the advantage of automatically handling the coordination and load distribution among the robot arms. In the second formulation, it was found that by choosing a general output equation it became possible simultaneously to superimpose the position and velocity error feedback with the force-torque error feedback in the task space.

Tarn, T. J.↗

Control of nonlinear systems with applications to constrained robots and spacecraft attitude stabilization

This thesis is organized in two parts. In Part 1, control systems described by a class of nonlinear differential and algebraic equations are introduced. A procedure for local stabilization based on a local state realization is developed. An alternative approach to local stabilization is developed based on a classical linearization of the nonlinear differential-algebraic equations. A theoretical framework is established for solving a tracking problem associated with the differential-algebraic system. First, a simple procedure is developed for the design of a feedback control law which ensures, at least locally, that the tracking error in the closed loop system lies within any given bound if the reference inputs are sufficiently slowly varying. Next, by imposing additional assumptions, a procedure is developed for the design of a feedback control law which ensures that the tracking error in the closed loop system approaches zero exponentially for reference inputs which are not necessarily slowly varying. The control design methodologies are used for simultaneous force and position control in constrained robot systems. The differential-algebraic equations are shown to characterize the slow dynamics of a certain nonlinear control system in nonstandard singularly perturbed form. In Part 2, the attitude stabilization (reorientation) of a rigid spacecraft using only two control torques is considered. First, the case of momentum wheel actuators is considered. The complete spacecraft dynamics are not controllable. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but a discontinuous feedback control strategy is constructed. Next, the case of gas jet actuators is considered. If the uncontrolled principal axis is not an axis of symmetry, the complete spacecraft dynamics are small time locally controllable. However, the spacecraft attitude cannot be asymptotically stabilized using continuous feedback, but a discontinuous stabilizing feedback control strategy is constructed. If the uncontrolled principal axis is an axis of symmetry, the complete spacecraft dynamics cannot be stabilized. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but again a discontinuous feedback control strategy is constructed.

Krishnan, Hariharan↗

Topography of Impact Structures on the Northern Polar Cap of Mars

The north polar residual ice deposits of mars are thought to be relatively young, based on the reported lack of any fresh impact craters in Viking Orbiter images. A handful of possible impact features were identified, but available data were inconclusive. Determining the number and current topographic characteristics of any craters on the polar residual ice surface is important in constraining the surface age, relative importance of polar cap resurfacing processes and ice flow dynamics, and the role of the cap deposits in the global volatile and climate cycles. Subsequent image processing advances and new Mars Orbiter Laser Altimeter (MOLA) data in the north polar region are now a considerable aid in impact feature identification. This study reconsiders the abundant Viking high-resolution polar images along with the new altimetry data. We examine possible impact features, compare their topography with known mars high-latitude impact features, and use their morphology as a constraint on crater modification processes.

Sakimoto, S. E. H.↗

Colors of Jupiter's large anticyclones and the interaction of a Tropical Red Oval with the Great Red Spot in 2008

The nature and mechanisms producing the chromophore agents that provide color to the upper clouds and hazes of the atmospheres of the giant planets are largely unknown. In recent times, the changes in red coloration that have occurred in large- and medium-scale Jovian anticyclones have been particularly interesting. In late June and early July 2008, a particularly color-intense red tropical oval interacted with the Great Red Spot (GRS) leading to the destruction of the red tropical oval and cloud dispersion. We present a detailed study of the tropical vortices, usually white but sometimes red, and a characterization of their color spectral signatures and dynamics. From the spectral reflectivity in methane bands we study their vertical cloud structure compared to that of the GRS and BA. Using two spectral indices we found a near-correlation between anticyclones cloud top altitudes and red color. We present detailed observations of the interaction of the red oval with the GRS and model simulations of the phenomena that allow us to constrain the relative vertical extent of the vortices. We conclude that the vertical cloud structure, vertical extent and dynamics of Jovian anticyclones are not the causes of their coloration. We propose that the red chromophore forms when background material (a compound or particles) is entrained by the vortex, transforming into red once inside the vortex due to internal conditions, exposure to ultraviolet radiation or to the mixing of two chemical compounds that react inside the vortex, confined by a potential vorticity ring barrier.

red color↗

Constrained Kalman Filtering Via Density Function Truncation 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 an analytic method of incorporating state variable inequality constraints in the Kalman filter. The resultant filter truncates the PDF (probability density function) of the Kalman filter estimate at the known constraints and then computes the constrained filter estimate as the mean of the truncated PDF. The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is demonstrated via simulation results obtained from a turbofan engine model. The turbofan engine model contains 3 state variables, 11 measurements, and 10 component health parameters. It is also shown that the truncated Kalman filter may be a more accurate way of incorporating inequality constraints than other constrained filters (e.g., the projection approach to constrained filtering).

Simon, Dan↗

Extratropical Cyclone Precipitation Life Cycles: A Satellite-Based Analysis

Although extratropical cyclones are the most common midlatitude storms, the relationship between their precipitation life cycle and dynamical strength life cycle has not been thoroughly analyzed. Given that thermodynamic heat exchanges associated with precipitation impact cyclone circulation, there is a need to understand the precipitation/dynamics relationship. Based on Integrated MultisatellitE Retrievals for Global Precipitation MeasurementGPM precipitation and Lagrangian cyclone tracks, the precipitation maximum occurs prior to the dynamical strength maximum 70% of the time. The lag in timing is consistent with the difference in cyclone precipitable water vapor at the two peaks. Conditional subsetting of the cyclone composites shows that if the precipitable water vapor distribution is constrained to be equal throughout the composite life cycle, the precipitation peak occurs very near the time of the peak in cyclone dynamical strength. Thus, the boost in dynamical strength caused by latent heat associated with precipitation manifest itself with little to no time lag.

Extratropical cyclone↗

Dynamics of Numerics & Spurious Behaviors in CFD Computations

The global nonlinear behavior of finite discretizations for constant time steps and fixed or adaptive grid spacings is studied using tools from dynamical systems theory. Detailed analysis of commonly used temporal and spatial discretizations for simple model problems is presented. The role of dynamics in the understanding of long time behavior of numerical integration and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) is explored. The study is complemented with examples of spurious behavior observed in steady and unsteady CFD computations. The CFD examples were chosen to illustrate non-apparent spurious behavior that was difficult to detect without extensive grid and temporal refinement studies and some knowledge from dynamical systems theory. Studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by available computing power. In large scale computations where the physics of the problem under study is not well understood and numerical simulations are the only viable means of solution, extreme care must be taken in both computation and interpretation of the numerical data. The goal of this paper is to explore the important role that dynamical systems theory can play in the understanding of the global nonlinear behavior of numerical algorithms and to aid the identification of the sources of numerical uncertainties in CFD.

Yee, Helen C.↗

Combined Bernstein Polynomial, Optimal Reciprocal Collision Avoidance, Differential Dynamic Programming for Trajectory Replanning and Collision Avoidance for UAM Vehicles

This paper presents an integration of Differential Dynamic Programming (DDP) with the Optimal Reciprocal Collision Avoidance (ORCA) algorithm as the basis for a new algorithm, titled Combined Bernstein Polynomial Optimal Reciprocal Collision Avoidance DDP (COBRA-DDP), for trajectory replanning and collision avoidance for Urban Air Mobility (UAM) vehicles. State-constrained variants of DDP provide the ability to plan trajectories while avoiding obstacles, but these methods require a large increase in computational time per iteration which hinders the overall speed of the algorithm. ORCA utilizes simplified dynamics to recognize potential collisions along a trajectory and provides an optimal velocity for the avoidance of multiple vehicles. These velocity commands, however, may not result in a dynamically feasible trajectory for DDP to plan around. As such, a Bernstein polynomial curve that considers the general dynamic constraints of the vehicle is generated to approximate a trajectory based on the velocity commands. COBRA-DDP optimizes this suggested trajectory via unconstrained DDP to provide a dynamically feasible trajectory that provides collision avoidance. This new trajectory can be applied to the vehicle or used to warm start the state constrained DDP algorithms to decrease computation time. Its benefits and effectiveness of the algorithm are demonstrated on a UAM Vertical Takeoff and Landing (VTOL) vehicle simulation with highly nonlinear dynamics.

Optimal Reciprocal Collision Avoidance↗