Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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 163 records · Page 9

Wake Vortex Effects Between Urban Air Mobility Aircraft

This study models the wake-induced downwash of an Urban Air Mobility (UAM) aircraft and the wake’s effects on a trailing aircraft’s dynamic stability when intersecting the wake trail. The Proctor wake vortex model is used to compute the wake-induced downwash over a 0.8 × 0.8 nautical mile flight area. The model is complemented by core radius growth, Proctor single phase wake decay, wake propagation, and Lamb wake aging. The wake data is then stored and loaded to the trailing aircraft as an atmospheric condition to affect the dynamics. Preliminary simulations include modeling the downwash from multiple UAM aircraft over the flight area and determining the normalized wake reduction at the ICAO 5 nautical mile separation criteria for a LIGHT aircraft trailing a MEDIUM aircraft.

Wake Vortex↗

Discrete and continuous dynamics modeling of a mass moving on a flexible structure

A general discrete methodology for modeling the dynamics of a mass that moves on the surface of a flexible structure is developed. This problem was motivated by the Space Station/Mobile Transporter system. A model reduction approach is developed to make the methodology applicable to large structural systems. To validate the discrete methodology, continuous formulations are also developed. Three different systems are examined: (1) simply-supported beam, (2) free-free beam, and (3) free-free beam with two points of contact between the mass and the flexible beam. In addition to validating the methodology, parametric studies were performed to examine how the system's physical properties affect its dynamics.

Herman, Deborah Ann↗

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗

Performance of an Optimized Eta Model Code on the Cray T3E and a Network of PCs

In the year 2001, NASA will launch the satellite TRIANA that will be the first Earth observing mission to provide a continuous, full disk view of the sunlit Earth. As a part of the HPCC Program at NASA GSFC, we have started a project whose objectives are to develop and implement a 3D cloud data assimilation system, by combining TRIANA measurements with model simulation, and to produce accurate statistics of global cloud coverage as an important element of the Earth's climate. For simulation of the atmosphere within this project we are using the NCEP/NOAA operational Eta model. In order to compare TRIANA and the Eta model data on approximately the same grid without significant downscaling, the Eta model will be integrated at a resolution of about 15 km. The integration domain (from -70 to +70 deg in latitude and 150 deg in longitude) will cover most of the sunlit Earth disc and will continuously rotate around the globe following TRIANA. The cloud data assimilation is supposed to run and produce 3D clouds on a near real-time basis. Such a numerical setup and integration design is very ambitious and computationally demanding. Thus, though the Eta model code has been very carefully developed and its computational efficiency has been systematically polished during the years of operational implementation at NCEP, the current MPI version may still have problems with memory and efficiency for the TRIANA simulations. Within this work, we optimize a parallel version of the Eta model code on a Cray T3E and a network of PCs (theHIVE) in order to improve its overall efficiency. Our optimization procedure consists of introducing dynamically allocated arrays to reduce the size of static memory, and optimizing on a single processor by splitting loops to limit the number of streams. All the presented results are derived using an integration domain centered at the equator, with a size of 60 x 60 deg, and with horizontal resolutions of 1/2 and 1/3 deg, respectively. In accompanying charts we report the elapsed time, the speedup and the Mflops as a function of the number of processors for the non-optimized version of the code on the T3E and theHIVE. The large amount of communication required for model integration explains its poor performance on theHIVE. Our initial implementation of the dynamic memory allocation has contributed to about 12% reduction of memory but has introduced a 3% overhead in computing time. This overhead was removed by performing loop splitting in some of the high demanding subroutines. When the Eta code is fully optimized in order to meet the memory requirement for TRIANA simulations, a non-negligeable overhead may appear that may seriously affect the efficiency of the code. To alleviate this problem, we are considering implementation of a new algorithm for the horizontal advection that is computationally less expensive, and also a new approach for marching in time.

Kouatchou, Jules↗

CFD Analysis in Advance of the NASA Juncture Flow Experiment

NASA through its Transformational Tools and Technologies Project (TTT) under the Advanced Air Vehicle Program, is supporting a substantial effort to investigate the formation and origin of separation bubbles found on wing-body juncture zones. The flow behavior in these regions is highly complex, difficult to measure experimentally, and challenging to model numerically. Multiple wing configurations were designed and evaluated using Computational Fluid Dynamics (CFD), and a series of wind tunnel risk reduction tests were performed to further down-select the candidates for the final experiment. This paper documents the CFD analysis done in conjunction with the 6 percent scale risk reduction experiment performed in NASA Langley's 14- by 22-Foot Subsonic Tunnel. The combined CFD and wind tunnel results ultimately helped the Juncture Flow committee select the wing configurations for the final experiment.

Flow Experiment↗

CFD Analysis in Advance of the NASA Juncture Flow Experiment

NASA through its Transformational Tools and Technologies Project (TTT) under the Advanced Air Vehicle Program, is supporting a substantial effort to further investigate the formation and origin of separation bubbles found on wing-body juncture zones. The flow behavior in these regions is highly complex, difficult to measure experimentally, and challenging to model numerically. Multiple wing configurations were designed and evaluated using Computational Fluid Dynamics (CFD), and a series of wind tunnel risk reduction tests were performed to further down-select the candidates for the final experiment. This paper documents the CFD analysis done in conjunction with the 6 percent scale risk reduction experiment performed in NASA Langley’s 14- by 22-Foot Subsonic Tunnel. The combined CFD and wind tunnel results ultimately help the Juncture Flow committee select the wing configurations for the final experiment.

Lee, H. C.↗

Arctic Warming in Response to Regional Aerosol Emissions Reductions

This study examines the Arctic surface air temperature response to regional aerosol emissions reductions using three fully coupled chemistry–climate models: National Center for Atmospheric Research-Community Earth System Model version 1, Geophysical Fluid Dynamics Laboratory-Coupled Climate Model version 3 (GFDL-CM3) and Goddard Institute for Space Studies-ModelE version 2. Each of these models was used to perform a series of aerosol perturbation experiments, in which emissions of different aerosol types (sulfate, black carbon (BC), and organic carbon) in different northern mid-latitude source regions, and of biomass burning aerosol over South America and Africa, were substantially reduced or eliminated. We find that the Arctic warms in nearly every experiment, the only exceptions being the U.S. and Europe BC experiments in GFDL-CM3 in which there is a weak and insignificant cooling. The Arctic warming is generally larger than the global mean warming (i.e. Arctic amplification occurs), particularly during non-summer months. The models agree that changes in the poleward atmospheric moisture transport are the most important factor explaining the spread in Arctic warming across experiments: the largest warming tends to coincide with the largest increases in moisture transport into the Arctic. In contrast, there is an inconsistent relationship (correlation) across experiments between the local radiative forcing over the Arctic and the simulated Arctic warming, with this relationship being positive in one model (GFDL-CM3) and negative in the other two. Our results thus highlight the prominent role of poleward energy transport in driving Arctic warming and amplification, and suggest that the relative importance of poleward energy transport and local forcing/feedbacks is likely to be model dependent.

Arctic↗

Modeling and simulations of the DSS 13 antenna control system

A model of the antenna control system for the azimuth and elevation axes of the Deep Space Station (DSS) 13 antenna is developed. This model is used for simulation of elevation and azimuth dynamics, cross-coupled dynamics, and radio-frequency pointing error due to both input commands and wind disturbances. This model also serves as a tool for the antenna controller design. A modal state-space model of the antenna structure was obtained from its finite-element model with a free-rotating tipping structure and alidade. Model reduction techniques applied separately for the antenna structure, elevation, and azimuth drives, and rate-loop model reduce the system order to one-third of that of the original, while preserving its significant dynamic properties. Extensive simulation results illustrate properties of the model.

Gawronski, W.↗

Application of Quasi-Static Modal Analysis to an Orion Multi-Purpose Crew Vehicle Test

Bolted structural joints often exhibit load-dependent stiffness and energy dissipation that leads to nonlinear, amplitude dependent frequency and damping in the structure. As an alternative to direct integration of the nonlinear equations of motion, quasi-static modal analysis (QSMA) determines the dependence of frequency and damping on response amplitude using loading behavior from nonlinear static analyses. QSMA has previously been demonstrated to substantially reduce computational cost and maintain accuracy relative to full nonlinear dynamic simulation. This work explores the applicability of QSMA to a complex, large-scale aerospace structure. QSMA is employed to analyze a nonlinear model of test hardware developed to support the Orion Multi-Purpose Crew Vehicle program, which exhibited nonlinear behavior during dynamic testing at flight-like load levels. In addition to the extraction of amplitude-dependent frequency and damping curves, a Bouc-Wen hysteresis model was used in conjunction with the quasi-static results to develop nonlinear, uncoupled, time-domain modal equations of motion for the structure. Excellent agreement was observed between the reduced and full-order nonlinear models, encouraging future employment of QSMA to support accurate and efficient model reduction of structures with bolted joint nonlinearities.

QSMA↗

Development of procedures for calculating stiffness and damping properties of elastomers in engineering applications. Part 1: Verification of basic methods

The primary aim of this investigation was verification of basic methods which are to be used in cataloging elastomer dynamic properties (stiffness and damping) in terms of viscoelastic model constants. These constants may then be used to predict dynamic properties for general elastomer shapes and operating conditions, thereby permitting optimum application of elastomers as energy absorption and/or energy storage devices in the control of vibrations in a broad variety of applications. The efforts reported involved: (1) literature search; (2) the design, fabrication and use of a test rig for obtaining elastomer dynamic test data over a wide range of frequencies, amplitudes, and preloads; and (3) the reduction of the test data, by means of a selected three-element elastomer model and specialized curve fitting techniques, to material properties. Material constants thus obtained have been used to calculate stiffness and damping for comparison with measured test data. These comparisons are excellent for a number of test conditions and only fair to poor for others. The results confirm the validity of the basic approach of the overall program and the mechanics of the cataloging procedure, and at the same time suggest areas in which refinements should be made.

Chiang, T.↗

Gravity Recovery and Interior Laboratory Mission (GRAIL) Orbit Determination

Launched on 10 September 2011 from the Cape Canaveral Air Force Station, Florida, the twin-spacecraft Gravity Recovery and Interior Laboratory (GRAIL) has the primary mission objective of generating a lunar gravity map with an unprecedented resolution via the Ka-band Lunar Gravity Ranging System (LGRS). After successfully executing nearly 30 maneuvers on their six-month journey, Ebb and Flow (aka GRAIL-A and GRAIL-B) established the most stringent planetary formation orbit on 1 March 2012 of approximately 30 km x 90 km in orbit size. This paper describes the orbit determination (OD) filter configurations, analyses, and results during the Trans-Lunar Cruise, Orbit Period Reduction, and Transition to Science Formation phases. The maneuver reconstruction strategies and their performance will also be discussed, as well as the navigation requirements, major dynamic models, and navigation challenges. GRAIL is the first mission to generate a full high-resolution gravity field of the only natural satellite of the Earth. It not only enables scientists to understand the detailed structure of the Moon but also further extends their knowledge of the evolutionary histories of the rocky inner planets. Robust and successful navigation was the key to making this a reality.

Orbit Period Reduction (OPR)↗

Computational fluid dynamic control

A general technique is presented for modeling fluid, or gas, dynamic systems specifically for the development of control systems. The numerical methods which are generally used in computational fluid dynamics are borrowed to create either continuous-time or discrete-time models of the particular fluid system. The resulting equations can be either left in a nonlinear form, or easily linearized about an operating point. As there are typically very many states in these systems, the usual linear model reduction methods can be used on them to allow a low-order controller to be designed. A simple example is given which typifies many internal flow control problems. The resulting control is termed computational fluid dynamic control.

Hartley, Tom T.↗

State reduction for semi-Markov reliability models

Semi-Markov processes have proved to be an effective and convenient tool to construct models of systems that achieve reliability by redundancy and reconfiguration. These models are able to depict complex system architectures and to capture the dynamics of fault arrival and system recovery. A disadvantage of this approach is that the models can be extremely large, which poses both a model construction and a computational problem. Techniques are needed to reduce the model size. Because these systems are used in critical applications where failure can be expensive, there must be an analytically derived bound for the error produced by the model reduction technique. Automatic model generation programs have been written to help the reliability analyst produce models of complex systems. Because of the importance of these programs, the model reduction technique needs to be precise and easily implemented. This paper presents a model reduction technique called trimming that can be applied to a popular class of systems. An error bound for the trimming procedure is derived that uses readily available system parameters. The trimming procedure is precisely described and appears easy to implement in a model generation program.

Reduced order systems↗

A Reduced Dimension Static, Linearized Kalman Filter and Smoother

An approximate Kalman filter and smoother, based on approximations of the state estimation error covariance matrix, is described. Approximations include a reduction of the effective state dimension, use of a static asymptotic error limit, and a time-invariant linearization of the dynamic model for error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. Examples of use come from TOPEX/POSEIDON.

Kalman filter smoother modeling ocean simulation↗

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

An Approximate Kalman Filter for Ocean Data Assimilation; An Example with an Idealized Gulf Stream Model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

['Kalman↗

Dynamical and photometric investigation of cometary type 2 tails

Photographic observations of the antitail of Comet Kohoutek (1973f) are photometrically studied and the results compared with the working model of an antitail. The applied technique of photometric reduction is described and the radial and transverse profiles of the antitail, corrected for the effects of the variable intensity of the sky background, are derived. The most important result reached so far is a quantitative confirmation of the previously suggested hypothesis, arguing that dust particles in the antitail suffered a significant loss in radius due to evaporation near the perihelion passage. Only particles initially larger than 0.1 to 0.15 mm in diameter survived. Numerically, however, this result is still tentative, because the dynamical effect exerted by particle evaporation remains to be accounted for.

Sekanina, Z.↗

A continuum model for dynamic analysis of the Space Station

Dynamic analysis of the International Space Station using MSC/NASTRAN had 1312 rod elements, 62 beam elements, 489 nodes and 1473 dynamic degrees of freedom. A realtime, man-in-the-loop simulation of such a model is impractical. This paper discusses the mathematical model for realtime dynamic simulation of the Space Station. Several key questions in structures and structural dynamics are addressed. First, to achieve a significant reduction in the number of dynamic degrees of freedom, a continuum equivalent representation of the Space Station truss structure which accounted for the unsymmetry of the basic configuration and resulted in the coupling of extensional and transverse deformation, is developed. Next, dynamic equations for the continuum equivalent of the Space Station truss structure are formulated using a matrix version of Kane's dynamical equations. Flexibility is accounted for by using a theory that accommodates extension, bending in two principal planes and shear displacement. Finally, constraint equations suitable for dynamic analysis of flexible bodies with closed loop configuration are developed and solution of the resulting system of equations is based on the zero eigenvalue theorem.

Thomas, Segun↗