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 451 records · Page 25

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

First Sagittarius A* Event Horizon Telescope Results. VIII. Physical Interpretation of the Polarized Ring

Abstract In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗

Designing open quantum systems with known steady states: Davies generators and beyond

We provide a systematic framework for constructing generic models of nonequilibrium quantum dynamics with a target stationary (mixed) state. Our framework identifies (almost) all combinations of Hamiltonian and dissipative dynamics that relax to a steady state of interest, generalizing the Davies’ generator for dissipative relaxation at finite temperature to nonequilibrium dynamics targeting arbitrary stationary states. We focus on Gibbs states of stabilizer Hamiltonians, identifying local Lindbladians compatible therewith by constraining the rates of dissipative and unitary processes. Moreover, given terms in the Lindbladian not compatible with the target state, our formalism identifies the operations – including syndrome measurements and local feedback – one must apply to correct these errors. Our methods also reveal new models of quantum dynamics: for example, we provide a “measurement-induced phase transition” in which measurable two-point functions exhibit critical (power-law) scaling with distance at a critical ratio of the transverse field and rate of measurement and feedback. Time-reversal symmetry – defined naturally within our formalism – can be broken both in effectively classical and intrinsically quantum ways. Our framework provides a systematic starting point for exploring the landscape of dynamical universality classes in open quantum systems, as well as identifying new protocols for quantum error correction.

Guo, Jinkang [Department of Physics and Center for↗

Exploring the impact of ions on oxygen K-edge X-ray absorption spectroscopy in NaCl solution using the GW-Bethe-Salpeter-equation approach

X-ray absorption spectroscopy (XAS) is a powerful experimental tool to probe the local structure in materials with the core hole excitations. Here, the oxygen K-edge XAS spectra of the NaCl solution and pure water are computed by using a recently developed GW-Bethe-Salpeter equation approach, based on configurations modeled by path-integral molecular dynamics with the deep-learning technique. The neural network is trained on ab initio data obtained with strongly constrained and appropriately normed density functional theory. The observed changes in the XAS features of the NaCl solution, compared to those of pure water, are in good agreement between experimental and theoretical results. We provided detailed explanations for these spectral changes that occur when NaCl is solvated in pure water. Specifically, the presence of solvating ion pairs leads to localization of electron-hole excitons. As a result, our theoretical XAS results support the theory that the effects of the solvating ions on the H-bond network are mainly confined within the first hydration shell of ions, however beyond the shell the arrangement of water molecules remains to be comparable to that observed in pure water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Evolution of radiative thermal instability in a confined medium

Thermally bistable fluid tends to self-organize into clouds of hot and cold material, which are internally uniform and separated by thin conduction fronts. The evolution of these clouds has been studied for isobaric systems, but when pressure is instead treated as a dynamical quantity and allowed to evolve self-consistently, fundamentally different dynamics appear. Such a treatment is necessary in some laboratory plasmas, whose volume is constrained but whose pressure can vary. Here, solutions are derived for the evolution of clouds, accounting for pressure variation and interactions between conduction fronts. Additional stable configurations and secondary instabilities are derived, which may be relevant to fusion plasmas and to the study of photoionized plasma in the laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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