Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled problems”

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 577 records · Page 32

Analytical aeroelastic stability considerations and conversion loads for an XV-15 tilt-rotor in a wind tunnel simulation

A rotorcraft analysis is conducted to assess tilt-rotor stability and conversion loads for the XV-15 rotor with metal blades within its specified test envelope. A 38-DOF flutter analysis based on the code by Johnson (1988) is developed to simulate a wind-tunnel test in which the rotor torque is constant and thereby study stability. The same analytical model provides the simulated loads including hub loads, blade loads, and oscillatory pitch-link loads with attention given to the nonuniform inflow through the proprotor in the presence of the wing. Tilt-rotor stability during the cruise mode is found to be sensitive to coupling effects in the control system stiffness, and a stability problem is identified in the XV-15 Advanced Technology Blades. The present analysis demonstrates that the tilt-rotor is stable within the specified test envelope of the NASA 40 x 80-ft wind tunnel.

Kottapalli, Sesi↗

Multidisciplinary optimization of a controlled space structure using 150 design variables

A general optimization-based method for the design of large space platforms through integration of the disciplines of structural dynamics and control is presented. The method uses the global sensitivity equations approach and is especially appropriate for preliminary design problems in which the structural and control analyses are tightly coupled. The method is capable of coordinating general purpose structural analysis, multivariable control, and optimization codes, and thus, can be adapted to a variety of controls-structures integrated design projects. The method is used to minimize the total weight of a space platform while maintaining a specified vibration decay rate after slewing maneuvers.

James, Benjamin B.↗

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(C0)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule 04(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in 04(+) has also been studied using the Brueckner coupled cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗

Application of CFD Analysis to Design Support and Problem Resolution for ASRM and RSRM

The use of Navier-Stokes CFD codes to predict the internal flow field environment in a solid rocket motor is a very important analysis element during the design phase of a motor development program. These computational flow field solutions uncover a variety of potential problems associated with motor performance as well as suggesting solutions to these problems. CFD codes have also proven to be of great benefit in explaining problems associated with operational motors such as in the case of the pressure spike problem with the STS-54B flight motor. This paper presents results from analyses involving both motor design support and problem resolution. The issues discussed include the fluid dynamic/mechanical stress coupling at field joints relative to significant propellant deformations, the prediction of axial and radial pressure gradients in the motor associated with motor performance and propellant mechanical loading, the prediction of transition of the internal flow in the motor associated with erosive burning, the accumulation of slag at the field joints and in the submerged nozzle region, impingement of flow on the nozzle nose, and pressure gradients in the nozzle region of the motor. The analyses presented in this paper have been performed using a two-dimensional axisymmetric model. Fluent/BFC, a three dimensional Navier-Stokes flow field code, has been used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical pressure-velocity coupling algorithm. Wall functions are used to represent the character of the viscous sub-layer flow, and an adjusted k-epsilon turbulence model especially configured for mass injection internal flows, is used to model the growth of turbulence in the motor port. Conclusions discussed in this paper consider flow field effects on the forward, center, and aft propellant grains except for the head end star grain region of the forward propellant segment. The field joints and the submerged nozzle are discussed as well. Conclusions relative to both the design evaluation of the ASRM and the RSRM scenarios explaining the pressure spikes were based on the flow field solutions presented in this paper.

Dill, Richard A.↗

Parallelization of a Parabolized Navier-Stokes Solver with a Design Optimizer

The design of future supersonic aircraft, such as the High Speed Civil Transport (HSCT), will rely heavily on computational methods for aircraft design and the prediction of the complex aerodynamics encountered in flight. Parabolized Navier-Stokes (PNS) equation flow solvers are recognized as efficient and accurate computational tools for the solution of supersonic and hypersonic flow-fields, while design optimizers have the potential to be valuable tools within the overall design process. Presently, however, the execution of the flow solver in conjunction with a design optimizer presents a computationally intensive and formidable problem. To meet the challenges and increasing demand for multidisciplinary numerical tools which are faster, more robust and provide greater functionality, alternative strategies are explored to increase computational throughput by coupling a design optimizer and flow solver in a parallel processing environment. To address this problem the parallel processing of a PNS flow solver with a nonlinear constraint design optimizer is investigated as an alternative computational method.

Pallis, J. M.↗

Improvements in Block-Krylov Ritz Vectors and the Boundary Flexibility Method of Component Synthesis

A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence relationship, proposed by Wilson, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is based upon the boundary flexibility vectors of the component. Improvements have been made in the formulation of the initial seed to the Krylov sequence, through the use of block-filtering. A method to shift the Krylov sequence to create Ritz vectors that will represent the dynamic behavior of the component at target frequencies, the target frequency being determined by the applied forcing functions, has been developed. A method to terminate the Krylov sequence has also been developed. Various orthonormalization schemes have been developed and evaluated, including the Cholesky/QR method. Several auxiliary theorems and proofs which illustrate issues in component mode synthesis and loss of orthogonality in the Krylov sequence have also been presented. The resulting methodology is applicable to both fixed and free- interface boundary components, and results in a general component model appropriate for any type of dynamic analysis. The accuracy is found to be comparable to that of component synthesis based upon normal modes, using fewer generalized coordinates. In addition, the block-Krylov recurrence algorithm is a series of static solutions and so requires significantly less computation than solving the normal eigenspace problem. The requirement for less vectors to form the component, coupled with the lower computational expense of calculating these Ritz vectors, combine to create a method more efficient than traditional component mode synthesis.

Carney, Kelly Scott↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

An Indirect Mixed-Sensitivity Approach to Microgravity Vibration Isolation: The Exploitation of Kinematic Coupling In Frequency-weighting Design-Filter Selections

Many space-science experiments need an active isolation system to provide them with the requisite microgravity environment. The isolation systems planned for use with the International Space Station have been appropriately modeled using relative position relative velocity, and acceleration states. In theory, frequency-weighting design filters can be applied to these state-space models, in order to develop optimal H2 or mixed-norm controllers with desired stability and performance characteristics. In practice, however, the kinematic coupling among the various states can lead, through the associated frequency-weighting-filters, to conflicting demands on the Riccati design "machinery." The results can be numerically ill-conditioned regulator and estimator Riccati equations and/or reduced intuition in the design process. In addition, kinematic coupling can result in a redundancy in the demands imposed by the frequency weights. Failure properly to account for this type of coupling can lead to an unnecessary increase in controller dimensionality and, in turn, controller complexity. This paper suggests a rational approach to the assignment of frequency-weighting design filters, in the presence of the kinematic coupling among states that exists in the microgravity vibration isolation problem.

Hampton, R. David↗

Modeling the Diurnal Tides in the MLT Region with the Doppler Spread Parameterization of Gravity Waves

Special Session: SA03 The mesosphere/lower thermosphere region: Structure, dynamics, composition, and emission. Ground based and satellite observations in the upper mesosphere and lower thermosphere (MLT) reveal large seasonal variations in the horizontal wind fields of the diurnal and semidiurnal tides. To provide an understanding of the observations, we discuss results obtained with our Numerical Spectral Model (NMS) that incorporates the gravity wave Doppler Spread Parameterization (DSP) of Hines. Our model reproduces many of the salient features observed, and we discuss numerical experiments that delineate the important processes involved. Gravity wave momentum deposition and the seasonal variations in the tidal excitation contribute primarily to produce the large equinoctial amplitude maxima in the diurnal tide. Gravity wave induced variations in eddy viscosity, not accounted for in the model, have been shown by Akmaev to be important too. For the semidiurnal tide, with amplitude maximum observed during winter solstice, these processes also contribute, but filtering by the mean zonal circulation is more important. A deficiency of our model is that it cannot reproduce the observed seasonal variations in the phase of the semidiurnal tide, and numerical experiments are being carried out to diagnose the cause and to alleviate this problem. The dynamical components of the upper mesosphere are tightly coupled through non-linear processes and wave filtering, and this may constrain the model and require it to reproduce in detail the observed phenomenology.

Mayr, H. G.↗

The Transient Dendritic Solidification Experiment (TDSE)

Dendritic solidification is a common mode of solidification. It is also an important model problem in nonequilibrium physics and pattern formation physics. Current theories couple the transfer of latent heat with selection mechanisms at the interface. Measurements of succinonitrile (SCN) dendrites in microgravity show reasonable agreement between heat transfer predictions and experiment. However, data and analysis for assessing interfacial physics theories are less definitive. We are studying, and will present data on, transient effects in dendritic growth of SCN. We employ the Clapeyron pressure/melting temperature effect to make a rapid change in a sample's hydrostatic pressure, and thereby rapidly change the specimen's melting temperature, forcing the dendrite to select a new steady-state. These initial measurements show some surprising and non-intuitive effects.

Koss, M. B.↗

An Indirect Mixed-Sensitivity Approach to Microgravity Vibration Isolation: The Exploitation of Kinematic Coupling In Frequency-Weighting Design-Filter Selections

Many space science experiments need an active isolation system to provide them with the requisite microgravity environment. The isolation systems planned for use with the International Space Station have been appropriately modeled using relative position, relative velocity, and acceleration states. In theory, frequency design filters can be applied to these state-space models, in order to develop optimal H, or mixed-norm controllers with desired stability- and performance characteristics. In practice. however, the kinematic coupling among the various states can lead, through the associated frequency-weighting-filters, to conflicting demands on the Riccati design "machinery." The results can be numerically ill-conditioned regulator and estimator Riccati equations and/or reduced intuition in the design process. In addition, kinematic coupling can result in a redundancy in the demands imposed by the frequency weights. Failure properly to account for this type of coupling can lead to an unnecessary increase in controller dimensionality and, in turn, controller complexity. This paper suggests a rational approach to the assignment of frequency weighting design filters, in the presence of the kinematic coupling among states that exists in the microgravity vibration isolation problem.

Hampton, R. David↗

Atmospheric Correction for Satellite Ocean Color Radiometry

This tutorial is an introduction to atmospheric correction in general and also documentation of the atmospheric correction algorithms currently implemented by the NASA Ocean Biology Processing Group (OBPG) for processing ocean color data from satellite-borne sensors such as MODIS and VIIRS. The intended audience is graduate students or others who are encountering this topic for the first time. The tutorial is in two parts. Part I discusses the generic atmospheric correction problem. The magnitude and nature of the problem are first illustrated with numerical results generated by a coupled ocean-atmosphere radiative transfer model. That code allow the various contributions (Rayleigh and aerosol path radiance, surface reflectance, water-leaving radiance, etc.) to the topof- the-atmosphere (TOA) radiance to be separated out. Particular attention is then paid to the definition, calculation, and interpretation of the so-called "exact normalized water-leaving radiance" and its equivalent reflectance. Part I ends with chapters on the calculation of direct and diffuse atmospheric transmittances, and on how vicarious calibration is performed. Part II then describes one by one the particular algorithms currently used by the OBPG to effect the various steps of the atmospheric correction process, viz. the corrections for absorption and scattering by gases and aerosols, Sun and sky reflectance by the sea surface and whitecaps, and finally corrections for sensor out-of-band response and polarization effects. One goal of the tutorial-guided by teaching needs- is to distill the results of dozens of papers published over several decades of research in atmospheric correction for ocean color remote sensing.

MODIS↗

Incremental States for Precise On-Orbit Relative Knowledge in Formation Flight

High precision and close proximity formation flight is an enabling technology for future space missions and requires an on-board relative navigation capability that is accurate to the mm-level and robust to formation parameters. Common estimation techniques linearize the entire formation about one spacecraft’s position, resulting in degraded accuracy due to linearization errors when separation distances become large. Additionally, methods which decouple absolute and relative state estimates usually require ad-hoc methods to incorporate the estimates together. This work discusses an alternative ”incremental” formulation of the relative navigation problem which is in variant to formation size, robust to coupling between absolute and relative dynamics, and can undergo similarity transformations to smoothly incorporate either absolute or relative information without numerical issues. A specific example of this architecture is presented in the context of formation navigation using Carrier-Differential Global Positioning System measurements, and is compared to a traditional leader-linearized filter. In the presence of accurate measurements, the incremental architecture is shown to reduce linearization errors and mean estimate errors by up to three orders of magnitude without the incorporation of new sensor information.

Seubert, Carl↗

Sensitivity Study of Dust-Induced Surface Erosion during Martian Planetary Entry

Large, global dust storms that periodically form in the Martian atmosphere pose a threat to NASA missions due to a lack of understanding of how they will affect the heat loading and mass loss of the entry vehicle’s thermal protection system (TPS). Previous studies have laid out viable methodologies for estimating the amount of surface recession from a significant dust storm, but rely on underlying modeling assumptions that are difficult to validate. The current work utilizes the DUST library, which is a Lagrangian particle solver under development at NASA Ames Research Center. DUST includes variable-fidelity fluid-particle coupling and has been tested for a wide range of problem types. The DUST library is used in conjunction with the US3D Naiver-Stokes solver. Using these tools, a sensitivity study is performed to determine the influence of the underlying dust modeling (drag coefficient, Nusselt number, and size distribution) on the surface erosion through a hypersonic Martian entry trajectory. The surface erosion was found to be insensitive to the drag and heating models, but was affected significantly by the size distributions parameters. The total recession was comparable to the thermochemical ablation predicted by material response codes, but the lack of a particle erosion model for PICA limits the predictive capability of the current work.

Thermal Protection Systems↗

An historical survey of computational methods in optimal control.

Review of some of the salient theoretical developments in the specific area of optimal control algorithms. The first algorithms for optimal control were aimed at unconstrained problems and were derived by using first- and second-variation methods of the calculus of variations. These methods have subsequently been recognized as gradient, Newton-Raphson, or Gauss-Newton methods in function space. A much more recent addition to the arsenal of unconstrained optimal control algorithms are several variations of conjugate-gradient methods. At first, constrained optimal control problems could only be solved by exterior penalty function methods. Later algorithms specifically designed for constrained problems have appeared. Among these are methods for solving the unconstrained linear quadratic regulator problem, as well as certain constrained minimum-time and minimum-energy problems. Differential-dynamic programming was developed from dynamic programming considerations. The conditional-gradient method, the gradient-projection method, and a couple of feasible directions methods were obtained as extensions or adaptations of related algorithms for finite-dimensional problems. Finally, the so-called epsilon-methods combine the Ritz method with penalty function techniques.

Polak, E.↗

An investigation of helicopter dynamic coupling using an analytical model

Many attempts have been made in recent years to predict the off-axis response of a helicopter to control inputs, and most have had little success. Since physical insight is limited by the complexity of numerical simulation models, this paper examines the off-axis response problem using an analytical model, with the goal of understanding the mechanics of the coupling. A new induced velocity model is extended to include the effects of wake distortion from pitch rate. It is shown that the inclusion of these results in a significant change in the lateral flap response to a steady pitch rate. The proposed inflow model is coupled with the full rotor/body dynamics, and comparisons are made between the model and flight test data for a UH-60 in hover. Results show that inclusion of induced velocity variations due to shaft rate improves correlation in the pitch response to lateral cycle inputs.

Keller, Jeffrey D.↗

The semi-discrete Galerkin finite element modelling of compressible viscous flow past an airfoil

A method is developed to solve the two-dimensional, steady, compressible, turbulent boundary-layer equations and is coupled to an existing Euler solver for attached transonic airfoil analysis problems. The boundary-layer formulation utilizes the semi-discrete Galerkin (SDG) method to model the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby permitting the use of a uniform finite element grid which provides high resolution near the wall and automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes, through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past NACA 0012 and RAE 2822 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack. All results show good agreement with experiment, and the coupled code proved to be a computationally-efficient and accurate airfoil analysis tool.

Meade, Andrew J., Jr.↗

Intricacies of frustrated magnetism in the Kondo metal YbAgGe

The combination of localized magnetic moments, their frustration and interaction with itinerant electrons is a key challenge of condensed matter physics. Frustrated magnetic interactions promote degenerate ground states with enhanced fluctuations, a topic that is predominantly studied in magnetic insulators. The coupling between itinerant and localized electrons in metals add complexity to the problem, and is presently formulated only for extreme cases in which the itinerant electrons mediate exchange between localized spins (RKKY interaction) or suppress the formation of magnetic moments (Kondo screening). Here, we report an in-depth experimental study of the distorted Kagome metal YbAgGe, unravelling the open questions of how frustration, localized magnetism and itinerant electrons are intertwined in frustrated Kondo metals. We find that coupled itinerant and localized electrons give rise to dynamic magnetic correlations below T* ≈ 20 K. At lower temperature, frustrated magnetic interactions establish anisotropic magnetic short-range correlations that culminate into antiferromagnetic long-range order below T N = 0.68 K with a significantly reduced modulated magnetic moment. We show that local moment Hamiltonians can yield limited understanding of the microscopic behaviour in frustrated metals, and prompt the extension of more sophisticated model Hamiltonians incorporating itinerant effects.

Mazzone, Daniel G. [PSI Center for Neutron and Muo↗