Search NASA⌕ Search

SEARCH · Search NASA

Results for “time domain simulation”

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 127 records · Page 7

Direct Numerical Simulation of an Airfoil with Sand Grain Roughness on the Leading Edge

As part of a computational study of acoustic radiation due to the passage of turbulent boundary layer eddies over the trailing edge of an airfoil, the Lattice-Boltzmann method is used to perform direct numerical simulations of compressible, low Mach number flow past an NACA 0012 airfoil at zero degrees angle of attack. The chord Reynolds number of approximately 0.657 million models one of the test conditions from a previous experiment by Brooks, Pope, and Marcolini at NASA Langley Research Center. A unique feature of these simulations involves direct modeling of the sand grain roughness on the leading edge, which was used in the abovementioned experiment to trip the boundary layer to fully turbulent flow. This report documents the findings of preliminary, proof-of-concept simulations based on a narrow spanwise domain and a limited time interval. The inclusion of fully-resolved leading edge roughness in this simulation leads to significantly earlier transition than that in the absence of any roughness. The simulation data is used in conjunction with both the Ffowcs Williams-Hawkings acoustic analogy and a semi-analytical model by Roger and Moreau to predict the farfield noise. The encouraging agreement between the computed noise spectrum and that measured in the experiment indicates the potential payoff from a full-fledged numerical investigation based on the current approach. Analysis of the computed data is used to identify the required improvements to the preliminary simulations described herein.

Ribeiro, Andre F. P.↗

Methods for Combining Payload Parameter Variations with Input Environment

Methods are presented for calculating design limit loads compatible with probabilistic structural design criteria. The approach is based on the concept that the desired limit load, defined as the largest load occuring in a mission, is a random variable having a specific probability distribution which may be determined from extreme-value theory. The design limit load, defined as a particular value of this random limit load, is the value conventionally used in structural design. Methods are presented for determining the limit load probability distributions from both time-domain and frequency-domain dynamic load simulations. Numerical demonstrations of the methods are also presented.

Merchant, D. H.↗

Methods for combining payload parameter variations with input environment

Methods are presented for calculating design limit loads compatible with probabilistic structural design criteria. The approach is based on the concept that the desired limit load, defined as the largest load occurring in a mission, is a random variable having a specific probability distribution which may be determined from extreme-value theory. The design limit load, defined as a particular of this random limit load, is the value conventionally used in structural design. Methods are presented for determining the limit load probability distributions from both time-domain and frequency-domain dynamic load simulations. Numerical demonstrations of the method are also presented.

Merchant, D. H.↗

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

Salamon, N. J.↗

High-Alpha Handling Qualities Flight Research on the NASA F/A-18 High Alpha Research Vehicle

A flight research study of high-angle-of-attack handling qualities has been conducted at the NASA Dryden Flight Research Center using the F/A-18 High Alpha Research Vehicle (HARV). The objectives were to create a high-angle-of-attack handling qualities flight database, develop appropriate research evaluation maneuvers, and evaluate high-angle-of-attack handling qualities guidelines and criteria. Using linear and nonlinear simulations and flight research data, the predictions from each criterion were compared with the pilot ratings and comments. Proposed high-angle-of-attack nonlinear design guidelines and proposed handling qualities criteria and guidelines developed using piloted simulation were considered. Recently formulated time-domain Neal-Smith guidelines were also considered for application to high-angle-of-attack maneuvering. Conventional envelope criteria were evaluated for possible extension to the high-angle-of-attack regime. Additionally, the maneuvers were studied as potential evaluation techniques, including a limited validation of the proposed standard evaluation maneuver set. This paper gives an overview of these research objectives through examples and summarizes result highlights. The maneuver development is described briefly, the criteria evaluation is emphasized with example results given, and a brief discussion of the database form and content is presented.

Keith D Wichman↗

Techniques for measurement of thoracoabdominal asynchrony

Respiratory motion measured by respiratory inductance plethysmography often deviates from the sinusoidal pattern assumed in the traditional Lissajous figure (loop) analysis used to determine thoraco-abdominal asynchrony, or phase angle phi. We investigated six different time-domain methods of measuring phi, using simulated data with sinusoidal and triangular waveforms, phase shifts of 0-135 degrees, and 10% noise. The techniques were then used on data from 11 lightly anesthetized rhesus monkeys (Macaca mulatta; 7.6 +/- 0.8 kg; 5.7 +/- 0.5 years old), instrumented with a respiratory inductive plethysmograph, and subjected to increasing levels of inspiratory resistive loading ranging from 5-1,000 cmH(2)O. L(-1). sec(-1).The best results were obtained from cross-correlation and maximum linear correlation, with errors less than approximately 5 degrees from the actual phase angle in the simulated data. The worst performance was produced by the loop analysis, which in some cases was in error by more than 30 degrees. Compared to correlation, other analysis techniques performed at an intermediate level. Maximum linear correlation and cross-correlation produced similar results on the data collected from monkeys (SD of the difference, 4.1 degrees ) but all other techniques had a high SD of the difference compared to the correlation techniques.We conclude that phase angles are best measured using cross-correlation or maximum linear correlation, techniques that are independent of waveform shape, and robust in the presence of noise. Copyright 2002 Wiley-Liss, Inc.

Non-NASA Center↗

Flight Dynamics Aspects of a Large Civil Tiltrotor Simulation Using Translational Rate Command

An in-depth analysis of a Large Civil Tiltrotor simulation with a Translational Rate Command control law that uses automatic nacelle deflections for longitudinal velocity control and lateral cyclic for lateral velocity control is presented. Results from piloted real-time simulation experiments and offline time and frequency domain analyses are used to investigate the fundamental flight dynamic and control mechanisms of the control law. The baseline Translational Rate Command conferred handling qualities improvements over an attitude command attitude hold control law but in some scenarios there was a tendency to enter PIO. Nacelle actuator rate limiting strongly influenced the PIO tendency and reducing the rate limits degraded the handling qualities further. Counterintuitively, increasing rate limits also led to a worsening of the handling qualities ratings. This led to the identification of a nacelle rate to rotor longitudinal flapping coupling effect that induced undesired pitching motions proportional to the allowable amount of nacelle rate. A modification that applied a counteracting amount of longitudinal cyclic proportional to the nacelle rate significantly improved the handling qualities. The lateral axis of the Translational Rate Command conferred Level 1 handling qualities in a Lateral Reposition maneuver. Analysis of the influence of the modeling fidelity on the lateral flapping angles is presented. It is showed that the linear modeling approximation is likely to have under-predicted the side-force and therefore under-predicted the lateral flapping at velocities above 15 ft/s. However, at lower velocities, and therefore more weakly influenced by the side force modeling, the accelerations that the control law commands also significantly influenced the peak levels of lateral flapping achieved.

Lawrence, Ben↗

Evaluating the Stability of NASA's Space Launch System with Adaptive Augmenting Control

NASA's baseline Space Launch System (SLS) flight control system (FCS) design includes an adaptive augmenting control (AAC) component that modifies the attitude control system response to provide the classical gain-scheduled control architecture with additional performance and robustness. The NASA Engineering and Safety Center (NESC) teamed with the Space Launch System (SLS) Program to perform a comprehensive assessment of the stability and robustness of the FCS with AAC. This paper provides an overview of the approach, specific analysis techniques, and outcomes that were particularly relevant for the SLS Program. Multiple analysis techniques that specifically target the nonlinear AAC were commissioned as part of this assessment, which was completed outside of the Program's standard design analysis cycle. The following analyses were included, with each technique adding its own valuable insights: Lyapunov-based stability analysis, classical stability analysis with static AAC gain variations, circle criterion-based analysis of the FCS with a time-varying element, time-domain stability margin assessment, Monte Carlo simulations with expanded dispersions, and an extensive set of stressing cases. Several of the completed analyses focused on determining whether the inclusion of AAC introduced risk to the FCS, while others quantified the benefits of the adaptive augmentation.

Tannen S Vanzwieten↗

Auralization of Air Vehicle Noise for Community Noise Assessment

This paper serves as an introduction to air vehicle noise auralization and documents the current state-of-the-art. Auralization of flyover noise considers the source, path, and receiver as part of a time marching simulation. Two approaches are offered; a time domain approach performs synthesis followed by propagation, while a frequency domain approach performs propagation followed by synthesis. Source noise description methods are offered for isolated and installed propulsion system and airframe noise sources for a wide range of air vehicles. Methods for synthesis of broadband, discrete tones, steady and unsteady periodic, and a periodic sources are presented, and propagation methods and receiver considerations are discussed. Auralizations applied to vehicles ranging from large transport aircraft to small unmanned aerial systems demonstrate current capabilities.

Rizzi, Stephen A.↗

Simulation-Based Analysis and Prediction of Thrust Vector Servoelastic Coupling

A method of analysis and prediction of servoelastic coupling in launch vehicles is presented, surveying the discovery and subsequent resolution of a predicted servoelastic resonance phenomenon affecting the NASA Space Launch System launch vehicle at specific flight conditions. A physics-based linearized multibody mechanization of the governing equations is combined with first principles analysis to demonstrate that antisymmetric bending of the solid rocket motors leads to a reduction of equivalent viscous modal damping through coupling with the thrust vector control actuators. The sensitivity to parameters and the effects of the resonance phenomenon on flight control performance and stability are confirmed through extensive simulation verification in the time and frequency domain. A novel enhancement in model fidelity that accounts for Coriolis effects of fluid flow on bending within the solid rocket motor case and nozzle is shown to add sufficient damping to reduce the risk of adverse control-structure interaction.

Jeb S Orr↗

Evaluation of a methodology for model identification in the time domain

A model identification methodology for structural dynamics has been applied to simulated vibrational data as a first step in evaluating its accuracy. The evaluation has taken into account a wide variety of factors which affect the accuracy of the procedure. The effects of each of these factors were observed in both the response time histories and the estimates of the parameters of the model by comparing them with the exact values of the system. Each factor was varied independently but combinations of these have also been considered in an effort to simulate real situations. The results of the tests have shown that for the chain model, the procedure yields robust estimates of the stiffness parameters under the conditions studied whenever uniqueness is ensured. When inaccuracies occur in the results, they are intimately related to non-uniqueness conditions inherent in the inverse problem and not to shortcomings in the methodology.

Beck, R. T.↗

Implicit Approaches for Moving Boundaries in a 3-D Cartesian Method

This work considers numerical simulation of three-dimensional flows with time-evolving boundaries. Such problems pose a variety of challenges for numerical schemes, and have received a substantial amount of attention in the recent literature. Since such simulations are unsteady, time-accurate solution of the governing equations is required. In special cases, the body motion can be treated by a uniform rigid motion of the computational domain. For the more general situation of relative-body motion, however, this simplification is unavailable and the simulations require a mechanism for ensuring that the mesh evolves with the moving boundaries. This involves a "remeshing" of the computational domain (either localized or global) at each physical timestep, and places a premium on both the speed and robustness of the remeshing algorithms. This work presents a method which includes unsteady flow simulation, rigid domain motion, and relative body motion using a time-evolving Cartesian grid system in three dimensions.

Murman, Scott M.↗

Verification and Validation of Elastodynamic Simulation Software for Aerospace Research

Physics-based simulation of nondestructive evaluation (NDE) inspection can help to advance the inspectability and reliability of mechanical systems. However, NDE simulations applicable to non-idealized mechanical components often require large compute domains and long run times. This has prompted development of custom NDE simulation software tailored to high performance computing (HPC) hardware. Verification and validation (V&V) is an integral part of developing this software to ensure implementations are robust and applicable to inspection problems, producing tools and simulations suitable for computational NDE research. This presentation addresses factors common to V&V of several elastodynamic simulation codes applicable to ultrasonic NDE. Examples are drawn from in-house simulation software at NASA Langley Research Center, ranging from ensuring reliability in a 1D heterogeneous media wave equation solver to the V&V needs of 3D cluster-parallel elastodynamic software. Factors specific to a research environment are addressed, where individual simulation results can be as relevant as the software product itself. Distinct facets of V&V are discussed including testing to establish software reliability, employing systematic approaches for consistency with fundamental conservation laws, establishing the numerical stability of algorithms, and demonstrating concurrence with empirical data. This talk also addresses V&V practices for small groups of researchers. This includes establishing resources (e.g. time and personnel) for V&V during project planning to mitigate and control the risk of setbacks. Similarly, we identify ways for individual researchers to use V&V during simulation software development itself to both speed up the development process and reduce incurred technical debt.

NDE↗

Computer simulation of multiple pilots flying a modern high performance helicopter

A computer simulation of a human response pilot mechanism within the flight control loop of a high-performance modern helicopter is presented. A human response mechanism, implemented by a low order, linear transfer function, is used in a decoupled single variable configuration that exploits the dominant vehicle characteristics by associating cockpit controls and instrumentation with specific vehicle dynamics. Low order helicopter models obtained from evaluations of the time and frequency domain responses of a nonlinear simulation model, provided by NASA Lewis Research Center, are presented and considered in the discussion of the pilot development. Pilot responses and reactions to test maneuvers are presented and discussed. Higher level implementation, using the pilot mechanisms, are discussed and considered for their use in a comprehensive control structure.

Zipf, Mark E.↗

Eulerian Time-Domain Filtering for Spatial LES

Eulerian time-domain filtering seems to be appropriate for LES (large eddy simulation) of flows whose large coherent structures convect approximately at a common characteristic velocity; e.g., mixing layers, jets, and wakes. For these flows, we develop an approach to LES based on an explicit second-order digital Butterworth filter, which is applied in,the time domain in an Eulerian context. The approach is validated through a priori and a posteriori analyses of the simulated flow of a heated, subsonic, axisymmetric jet.

Pruett, C. David↗

Simulations of Cyclone Sidr in the Bay of Bengal with a High-Resolution Model: Sensitivity to Large-Scale Boundary Forcing

The predictability of Cyclone Sidr in the Bay of Bengal was explored in terms of track and intensity using the Advanced Research Hurricane Weather Research Forecast (AHW) model. This constitutes the first application of the AHW over an area that lies outside the region of the North Atlantic for which this model was developed and tested. Several experiments were conducted to understand the possible contributing factors that affected Sidr s intensity and track simulation by varying the initial start time and domain size. Results show that Sidr s track was strongly controlled by the synoptic flow at the 500-hPa level, seen especially due to the strong mid-latitude westerly over north-central India. A 96-h forecast produced westerly winds over north-central India at the 500-hPa level that were notably weaker; this likely caused the modeled cyclone track to drift from the observed actual track. Reducing the model domain size reduced model error in the synoptic-scale winds at 500 hPa and produced an improved cyclone track. Specifically, the cyclone track appeared to be sensitive to the upstream synoptic flow, and was, therefore, sensitive to the location of the western boundary of the domain. However, cyclone intensity remained largely unaffected by this synoptic wind error at the 500-hPa level. Comparison of the high resolution, moving nested domain with a single coarser resolution domain showed little difference in tracks, but resulted in significantly different intensities. Experiments on the domain size with regard to the total precipitation simulated by the model showed that precipitation patterns and 10-m surface winds were also different. This was mainly due to the mid-latitude westerly flow across the west side of the model domain. The analysis also suggested that the total precipitation pattern and track was unchanged when the domain was extended toward the east, north, and south. Furthermore, this highlights our conclusion that Sidr was influenced from the west side of the domain. The displacement error was significantly reduced after the domain size from the western model boundary was decreased. Study results demonstrate the capability and need of a high-resolution mesoscale modeling framework for simulating the complex interactions that contribute to the formation of tropical cyclones over the Bay of Bengal region

Kumar, Anil↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗

Lightning-related fields in the ionosphere

Lightning-related transient electric fields and conduction current densities due to the redistribution of charge by the return stroke, but not directly associated with the high-intensity current channel, are investigated via a simulation of Maxwell's equations. The investigation presents evidence to suggest that the lightning-related charge redistribution and subsequent relaxation, rather than the high intensity current, is the source of the major component of the transient field and energy coupled to the ionosphere. A time domain finite element analysis is employed in the simulation, using a conductivity profile obtained simultaneously with electric field measurements which agree, in general terms, with the simulations. It is concluded that prior depictions of these phenomena underestimate the fraction of the total energy that couples to and dissipates within the ionosphere and magnetosphere.

Baginski, M. E.↗