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 289 records · Page 16

New Observations of the Crab Nebula and Pulsar

We present a phase-resolved study of the X-ray spectrum of the Crab Pulsar, using data obtained in a special mode with the Chandra X-ray Observatory. The superb angular resolution easily enables discerning the Pulsar from the surrounding nebulosity, even at pulse minimum. We find that the Pulsar's X-ray spectral index varies sinusoidally with phase---except over the same phase range for which rather abrupt changes in optical polarization magnitude and position angle have been reported. In addition, we use the X-ray data to constrain the surface temperature for various neutron-star equations of state and atmospheres. Finally, we present new data on dynamical variations of structure within the Nebula.

Weisskopf, Martin C.↗

VLTI-GRAVITY Measurements of Cool Evolved Stars I. Variable Photosphere and Extended Atmosphere of the Mira Star R Peg

Context. Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Aim. Here we measure the variability of the oxygen-rich Mira star R Peg in near-continuum and molecular bands. Methods. We used near-infrared K-band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the Very Large Telescope Interferometer (VLTI). Results. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.24 m are steadily increasing with values of 8.7 0.1 mas, 9.4 0.1 mas, 9.8 0.1 mas, and 9.9 0.1 mas at visual phases of 0.15, 0.36, 0.45, 0.53, respectively. UD diameters at a bandpass around 2.05 m, dominated by water vapour, follow the near-continuum variability at larger UD diameters between 10.7 mas and 11.7 mas. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters, with values between 12.3 mas and 11.7 mas. Conclusions. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes in the extended atmosphere along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.

photospheric radius↗

On spurious behavior of CFD simulations

Spurious behavior in underresolved grids and/or semi-implicit temporal discretizations for four computational fluid dynamics (CFD) simulations are studied. The numerical simulations consist of (a) a 1-D chemically relaxed nonequilibrium model, (b) the direct numerical simulation (DNS) of 2-D incompressible flow over a backward facing step, (c) a loosely-coupled approach for a 2-D fluid-structure interaction, and (d) a 3-D compressible unsteady flow simulation of vortex breakdown in delta wings. Using knowledge from dynamical systems theory, various types of spurious behaviors that are numerical artifacts were systematically identified. These studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by the available computing power. In large scale computations underresolved grids, semi-implicit procedures, loosely-coupled implicit procedures, and insufficiently long time integration in DNS are most often unavoidable. Consequently, care must be taken in both computation and in interpretation of the numerical data. The results presented confirm the important role that dynamical systems theory can play in the understanding of the nonlinear behavior of numerical algorithms and in aiding the identification of the sources of numerical uncertainties in CFD.

Numerical Solution↗

Quantifying Emergent Fluid Dynamics Using Reynolds-Interpolated Fluid Reduced-order Models

Fluid reduced-order models (ROMs) which capture the flow physics within the problem's physical domain are usually constrained in accuracy to only the parameter points, e.g. Reynolds and Mach numbers, at which reference data was provided. Interpolation-focused quantity-of-interest ROMs are often structured differently and fail to provide flow volume data with the same quality - if at all. In this paper, techniques which reside at the intersection of these two ROM schools - flow physics ROMs which can be interpolated within a parameter space of interest - are explored. Using a combination of existing and novel techniques, emergent physics are identified using a fluid ROM at parameter points which are not provided in the ROM's training data.

uncertainty quantification↗

Quantifying Emergent Fluid Dynamics Using Reynolds-Interpolated Fluid Reduced-order Models

Fluid reduced-order models (ROMs) which capture the flow physics within the problem's physical domain are usually constrained in accuracy to only the parameter points, e.g. Reynolds and Mach numbers, at which reference data was provided. Interpolation-focused quantity-of-interest ROMs are often structured differently and fail to provide flow volume data with the same quality - if at all. In this paper, techniques which reside at the intersection of these two ROM schools - flow physics ROMs which can be interpolated within a parameter space of interest - are explored. Using a combination of existing and novel techniques, emergent physics are identified using a fluid ROM at parameter points which are not provided in the ROM's training data.

uncertainty quantification↗

Experimental Validation of a Thermoelastic Model for SMA Hybrid Composites

This study presents results from experimental validation of a recently developed model for predicting the thermomechanical behavior of shape memory alloy hybrid composite (SMAHC) structures, composite structures with an embedded SMA constituent. The model captures the material nonlinearity of the material system with temperature and is capable of modeling constrained, restrained, or free recovery behavior from experimental measurement of fundamental engineering properties. A brief description of the model and analysis procedures is given, followed by an overview of a parallel effort to fabricate and characterize the material system of SMAHC specimens. Static and dynamic experimental configurations for the SMAHC specimens are described and experimental results for thermal post-buckling and random response are presented. Excellent agreement is achieved between the measured and predicted results, fully validating the theoretical model for constrained recovery behavior of SMAHC structures.

Turner, Travis L.↗

Digital robust control law synthesis using constrained optimization

Development of digital robust control laws for active control of high performance flexible aircraft and large space structures is a research area of significant practical importance. The flexible system is typically modeled by a large order state space system of equations in order to accurately represent the dynamics. The active control law must satisy multiple conflicting design requirements and maintain certain stability margins, yet should be simple enough to be implementable on an onboard digital computer. Described here is an application of a generic digital control law synthesis procedure for such a system, using optimal control theory and constrained optimization technique. A linear quadratic Gaussian type cost function is minimized by updating the free parameters of the digital control law, while trying to satisfy a set of constraints on the design loads, responses and stability margins. Analytical expressions for the gradients of the cost function and the constraints with respect to the control law design variables are used to facilitate rapid numerical convergence. These gradients can be used for sensitivity study and may be integrated into a simultaneous structure and control optimization scheme.

Mukhopadhyay, Vivekananda↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Earth system dynamics↗

Prolate spheroidal slosh model for fluid motion

Mathematical model, designed for zero gravity conditions, analyzes dynamic effects of large amplitude fluid motion interior to a rigid body. It has two advantages over other mathematical models: (1) constrains slosh motion to given region in natural manner, and (2) allows equilibrium position of slosh mass to be anywhere on slosh surface.

Flanders, H. A.↗

Imaging spectrophotometry of the nuclear outflow of NGC 1068

This observational program (in conjunction with R. B. Tully (IfA, Honolulu), and J. Bland (Rice U., Houston)) aims to constrain the kinematic organization and dominant excitation mechanisms of ionized gas in active galaxies. More generally, researchers are interested in the dynamics of radiative, supersonic flows in the Interstellar Medium (ISM). Imaging Fabry-Perot interferometers and low-noise Charge Coupled Devices (CCDs) are used for complete spatial coverage of the complex gas distribution in circumnuclear narrow-line regions (NLRs). Extranuclear emission line widths in NLRs can exceed 3000 km s(-1), so to avoid inter-order confusion researchers use an etalon of 4000 km s(-1) free spectral range to map (N II) lambda lambda 6548, 6583 and H alpha. To maximize spatial resolution, researchers select nearby active systems independent of luminosity but known to possess interesting morphologies and/or high-velocity extranuclear ionized gas. Monochromatic images Full Width Half Maximum (FWHM) approx. 65 km s(-1) have thus far been obtained in 1 second or better seeing at the U. Hawaii 2.2m, CFH 3.6m, and CTIO 4.0m telescopes. These are stacked into grids of line profiles, of spectrophotometric quality, at sub-arcsecond increments across a 3 second field. To handle the approx. 20,000 to 300,000 useful spectra that arise from a typical night's work, researchers have developed a complete analysis and reduction package for VAX and Sun image workstations. Reduction involves parametrization of approx. 10 to the 8th raw data points to a few maps (e.g., velocities of each kinematic subsystem, continuum-free line fluxes) containing approx. 10 to the 5th pixels. Researchers identify kinematic and structural symmetries by examining these maps and the point to point variations of the synthesized line profiles. The combination of monochromatic images and full spatial sampling of line profiles has allowed them to isolate such symmetries and has led to reliable kinematic deprojections.

Cecil, Gerald↗

Historical Carbon Dioxide Emissions Caused by Land-Use Changes are Possibly Larger than Assumed

The terrestrial biosphere absorbs about 20% of fossil-fuel CO2 emissions. The overall magnitude of this sink is constrained by the difference between emissions, the rate of increase in atmospheric CO2 concentrations, and the ocean sink. However, the land sink is actually composed of two largely counteracting fluxes that are poorly quantified: fluxes from land-use change andCO2 uptake by terrestrial ecosystems. Dynamic global vegetation model simulations suggest that CO2 emissions from land-use change have been substantially underestimated because processes such as tree harvesting and land clearing from shifting cultivation have not been considered. As the overall terrestrial sink is constrained, a larger net flux as a result of land-use change implies that terrestrial uptake of CO2 is also larger, and that terrestrial ecosystems might have greater potential to sequester carbon in the future. Consequently, reforestation projects and efforts to avoid further deforestation could represent important mitigation pathways, with co-benefits for biodiversity. It is unclear whether a larger land carbon sink can be reconciled with our current understanding of terrestrial carbon cycling. Our possible underestimation of the historical residual terrestrial carbon sink adds further uncertainty to our capacity to predict the future of terrestrial carbon uptake and losses.

Arneth, A.↗

A Generalized Framework for Constrained Design Optimization of General Supersonic Configurations Using Adjoint Based Sensitivity Derivatives

The Aeronautics Research Mission Directorate (ARMD) sent out an NASA Research Announcement (NRA) for proposals soliciting research and technical development. The proposed research program was aimed at addressing the desired milestones and outcomes of ROA (ROA-2006) Subtopic A.4.1.1 Advanced Computational Methods. The second milestone, SUP.1.06.02 Robust, validated mesh adaptation and error quantification for near field Computational Fluid Dynamics (CFD), was addressed by the proposed research. Additional research utilizing the direct links to geometry through a CAD interface enabled by this work will allow for geometric constraints to be applied and address the final milestone, SUP2.07.06 Constrained low-drag supersonic aerodynamic design capability. The original product of the proposed research program was an integrated system of tools that can be used for the mesh mechanics required for rapid high fidelity analysis and for design of supersonic cruise vehicles. These Euler and Navier-Stokes volume grid manipulation tools were proposed to efficiently use parallel processing. The mesh adaptation provides a systematic approach for achieving demonstrated levels of accuracy in the solutions. NASA chose to fund only the mesh generation/adaptation portion of the proposal. So this report describes the completion of the proposed tasks for mesh creation, manipulation and adaptation as it pertains to sonic boom prediction of supersonic configurations.

Karman, Steve L., Jr.↗

Constraining the Hydration of X-Ray Amorphous and Clay Mineral Phases in Gale Crater, Mars

Clay minerals and salts in the Gale crater stratigraphy investigated by the Mars Science Laboratory (MSL) Curiosity rover preserve information about sediment sources, depositional environment, and post-depositional alteration history. The compositions of amorphous components in samples have bearing on the redox potential, pH, and chemistry of alteration fluids, and on the heterogeneity of these conditions. Interlayer clay hydration in present-day Gale crater could be an indication of high water/rock ratios of late-stage fluid alteration events, which increases the potential habitability of these post-depositional environments. The MSL Dynamic Albedo of Neutrons instrument (DAN) measures bulk hydration and here we attempt to determine the hydration states of both clay minerals and amorphous components throughout the Curiosity traverse to constrain the abundance of individual amorphous phases and determine if clay interlayer water is present today. These implications are also more broadly applicable to similar sedimentary successions outside Gale crater. In this study, we compare hydration results from DAN to mineralogical results from the MSL Chemistry and Mineralogy instrument (CheMin) for 16 different fluvio-lacustrine sample locations. These samples typically contain one or more hydrated crystalline phases, a hydrated X-ray amorphous fraction composed of multiple amorphous phases, and one or more clay mineral phases which CheMin has consistently measured as collapsed (i.e., no interlayer water). However, Gale samples are known to dehydrate during sample acquisition and handling, and because DAN measures hydration to a depth of tens of cm, we can test for the presence of clay interlayer water in situ by comparing DAN bulk hydration to CheMin mineralogy. We can also constrain the hydration of the samples’ amorphous fractions, and with geochemical data from the ChemCam and APXS instruments, this will allow us to better define the phases which compose this fraction.

S Czarnecki↗

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗

Input Constrained M-MRAC for Multirotors Operating in an Urban Environment

The paper presents a modified model reference adaptive control (M-MRAC) for multi-input multi-output nonlinear dynamical systems with time varying parametric uncertainties and bounded external disturbances. It uses a prediction model to rapidly generate adaptive estimates of the system's uncertainties with adjustable errors that converge to a small neighborhood of the origin. A sufficient condition is derived to specify the region of attraction in the space of initialization errors, design parameters and external commends. The approach is applied to thrust controlled multi-rotor air vehicles operating in an urban environment. It is shown that the designed controller can provide a good tracking of a given trajectory in the unknown urban wind field, assuming that the maximum thrust generated by the rotors is known. The performance of the algorithms are demonstrated in simulations.

Stepanyan, Vahram↗

Flutter-Constrained Optimization with the Linearized Frequency-Domain Approach

Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.

Aeroelasticity↗

Flutter-Constrained Optimization with the Linearized Frequency-Domain Approach

Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.

Aeroelasticity↗

Interplanetary trajectory design for the Mariner Mark II Comet Rendezvous/Asteroid Flyby mission

The primary purpose of the first mission utilizing a Mariner Mark II spacecraft will be to rendezvous with the short-period comet Wild 2 and conduct a proximal examination during the comet's quiescent and active phases. The purpose of the present study is to elucidate, using an accurate dynamic model, some of the characteristics of the interplanetary trajectory that are relevant to the mission design. Optimal transfer trajectories are generated, constrained by the Shuttle/Centaur launch vehicle capability and the total spacecraft mass and performance. Asteroid flyby opportunities are examined, and the effects of selected asteroid encounters on mission performance are discussed.

Myers, M. R.↗