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

Land Surface Data Assimilation and the Gulf Coast Sea-Breeze

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The focus of this paper is to examine how the satellite assimilation technique impacts simulations of near-surface meteorology on the 4- to 48-hour time scale. The technique is implemented within the PSU/NCAR MM5 V34 and applied on a 36-km CONUS domain and a 12-km nest centered over Oklahoma. Two Friday periods during 2000 were selected for study, one during the warm season (May-June) and the other in the cold season (Nov-Dec). Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint show that the assimilation technique can improve numerical simulations on both regional and continental scales. Comparison of the simulated surface energy and water fluxes with observations from Energy Balance Bowen Ratio measurements taken at CART/ARM sites in the Southern Great Plains showed the assimilation technique produces a realistic argument of the Bowen ratio in both seasons. Additional runs with the Oregon State University land surface scheme as implemented within MM5 will be made and comparisons with assimilation results will be presented at the meeting.

Lapenta, William M.↗

An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder

Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.

Mankbadi, M. R.↗

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

Psychoacoustic Analysis of Synthesized Jet Noise

An aircraft noise synthesis capability is being developed so the annoyance caused by proposed aircraft can be assessed during the design stage. To make synthesized signals as realistic as possible, high fidelity simulation is required for source (e.g., engine noise, airframe noise), propagation and receiver effects. This psychoacoustic study tests whether the jet noise component of synthesized aircraft engine noise can be made more realistic using a low frequency oscillator (LFO) technique to simulate fluctuations in level observed in recordings. Jet noise predictions are commonly made in the frequency domain based on models of time-averaged empirical data. The synthesis process involves conversion of the frequency domain prediction into an audible pressure time history. However, because the predictions are time-invariant, the synthesized sound lacks fluctuations observed in recordings. Such fluctuations are hypothesized to be perceptually important. To introduce time-varying characteristics into jet noise synthesis, a method has been developed that modulates measured or predicted 1/3-octave band levels with a (<20Hz) LFO. The LFO characteristics are determined through analysis of laboratory jet noise recordings. For the aft emission angle, results indicate that signals synthesized using a generic LFO are perceived as more similar to recordings than those using no LFO, and signals synthesized with an angle-specific LFO are more similar to recordings than those synthesized with a generic LFO.

Okcu, Selen↗

Tangle-Free Finite Element Mesh Motion for Ablation Problems

In numerical simulations involving boundaries that evolve in time, the primary challenge is updating the computational mesh to reflect the physical changes in the domain. In particular, the fundamental objective for any such \mesh motion" scheme is to maintain mesh quality and suppress unphysical geometric anamolies and artifacts. External to a physical process of interest, mesh motion is an added component that determines the specifics of how to move the mesh given certain limited information from the main system. This paper develops a set of boundary conditions designed to eliminate tangling and internal collision within the context of PDE-based mesh motion (linear elasticity). These boundary conditions are developed for two- and three-dimensional meshes. The paper presents detailed algorithms for commonly occuring topological scenarios and explains how to apply them appropriately. Notably, the techniques discussed herein make use of none of the specifics of any particular formulation of mesh motion and thus are more broadly applicable. The two-dimensional algorithms are validated by an extensive verification procedure. Finally, many examples of diverse geometries in both two- and three-dimensions are shown to showcase the capabilities of the tangle-free boundary conditions.

Droba, Justin↗

Wall-Resolved Large Eddy Simulations of Transonic Shock-Induced Flow Separation

This paper reports the wall-resolved large eddy simulations of shock-induced boundary layer separation over an axisymmetric bump for a flow Mach number of 0.875 and a chord-based Reynolds number of 2.763 million. The incoming boundary layer has a momentum-thickness Reynolds number of 6600 at one and a half chord lengths upstream of the leading edge. The calculations simulate the experiment by Bachalo and Johnson (AIAA Journal, Vol. 24, No. 3, 1986), except that the tunnel walls are ignored and the simulations are performed assuming free air with as many as 24 billion grid points. The effects of domain span, grid resolution and time step on the predictions are examined. The results are found to show some sensitivity to the studied parameters. Owing to the outer boundary conditions, the predicted surface pressure distribution as well as the flow separation and reattachment locations tend to agree better with the experimental results from the larger (6 × 6 ft) tunnel than those from the smaller (2 × 2 ft) tunnel. The predicted Reynolds shear stress profiles in the separated region differ by as much as 31%from the experimental results that were only obtained in the smaller tunnel. The most accurate surface pressure distribution obtained in this study lies within the scatter of the measurements taken in the two facilities.

Uzun, Ali↗

The Attenuation of a Detonation Wave by an Aircraft Engine Axial Turbine Stage

A Constant Volume Combustion Cycle Engine concept consisting of a Pulse Detonation Combustor (PDC) followed by a conventional axial turbine was simulated numerically to determine the attenuation and reflection of a notional PDC pulse by the turbine. The multi-stage, time-accurate, turbomachinery solver TURBO was used to perform the calculation. The solution domain consisted of one notional detonation tube coupled to 5 vane passages and 8 rotor passages representing 1/8th of the annulus. The detonation tube was implemented as an initial value problem with the thermodynamic state of the tube contents, when the detonation wave is about to exit, provided by a 1D code. Pressure time history data from the numerical simulation was compared to experimental data from a similar configuration to verify that the simulation is giving reasonable results. Analysis of the pressure data showed a spectrally averaged attenuation of about 15 dB across the turbine stage. An evaluation of turbine performance is also presented.

VanZante, Dale↗

High-power graphic computers for visual simulation: a real-time--rendering revolution

Advances in high-end graphics computers in the past decade have made it possible to render visual scenes of incredible complexity and realism in real time. These new capabilities make it possible to manipulate and investigate the interactions of observers with their visual world in ways once only dreamed of. This paper reviews how these developments have affected two preexisting domains of behavioral research (flight simulation and motion perception) and have created a new domain (virtual environment research) which provides tools and challenges for the perceptual psychologist. Finally, the current limitations of these technologies are considered, with an eye toward how perceptual psychologist might shape future developments.

NASA Discipline Space Human Factors↗

Real-Time Dynamic Modeling - Data Information Requirements and Flight Test Results

Practical aspects of identifying dynamic models for aircraft in real time were studied. Topics include formulation of an equation-error method in the frequency domain to estimate non-dimensional stability and control derivatives in real time, data information content for accurate modeling results, and data information management techniques such as data forgetting, incorporating prior information, and optimized excitation. Real-time dynamic modeling was applied to simulation data and flight test data from a modified F-15B fighter aircraft, and to operational flight data from a subscale jet transport aircraft. Estimated parameter standard errors, prediction cases, and comparisons with results from a batch output-error method in the time domain were used to demonstrate the accuracy of the identified real-time models.

Morelli, Eugene A.↗

Development of a High-Fidelity Simulation Environment for Shadow-Mode Assessments of Air Traffic Concepts

This paper will describe the purpose, architecture, and implementation of a gate-to-gate, high-fidelity air traffic simulation environment called the Shadow Mode Assessment using Realistic Technologies for the National Airspace System (SMART-NAS) Test Bed.The overarching purpose of the SMART-NAS Test Bed (SNTB) is to conduct high-fidelity, real-time, human-in-the-loop and automation-in-the-loop simulations of current and proposed future air traffic concepts for the Next Generation Air Transportation System of the United States, called NextGen. SNTB is intended to enable simulations that are currently impractical or impossible for three major areas of NextGen research and development: Concepts across multiple operational domains such as the gate-to-gate trajectory-based operations concept; Concepts related to revolutionary operations such as the seamless and widespread integration of large and small Unmanned Aerial System (UAS) vehicles throughout U.S. airspace; Real-time system-wide safety assurance technologies to allow safe, increasingly autonomous aviation operations. SNTB is primarily accessed through a web browser. A set of secure support services are provided to simplify all aspects of real-time, human-in-the-loop and automation-in-the-loop simulations from design (i.e., prior to execution) through analysis (i.e., after execution). These services include simulation architecture and asset configuration; scenario generation; command, control and monitoring; and analysis support.

human-in-the-loop↗

Formation Flying Control of Multiple Spacecraft

The problem of coordination and control of multiple spacecraft (MS) moving in formation is considered. Here, each MS is modeled by a rigid body with fixed center of mass. First, various schemes for generating the desired formation patterns are discussed, Then, explicit control laws for formation-keeping and relative attitude alignment based on nearest neighbor-tracking are derived. The necessary data which must be communicated between the MS to achieve effective control are examined. The time-domain behavior of the feedback-controlled MS formation for typical low-Earth orbits is studied both analytically and via computer simulation. The paper concludes with a discussion of the implementation of the derived control laws, and the integration of the MS formation coordination and control system with a proposed inter-spacecraft communication/computing network.

Hadaegh, F. Y.↗

Performance Dependences of Multiplication Layer Thickness for InP/InGaAs Avalanche Photodiodes Based on Time Domain Modeling

InP/InGaAs avalanche photodiodes (APDs) are being widely utilized in optical receivers for modern long haul and high bit-rate optical fiber communication systems. The separate absorption, grading, charge, and multiplication (SAGCM) structure is an important design consideration for APDs with high performance characteristics. Time domain modeling techniques have been previously developed to provide better understanding and optimize design issues by saving time and cost for the APD research and development. In this work, performance dependences on multiplication layer thickness have been investigated by time domain modeling. These performance characteristics include breakdown field and breakdown voltage, multiplication gain, excess noise factor, frequency response and bandwidth etc. The simulations are performed versus various multiplication layer thicknesses with certain fixed values for the areal charge sheet density whereas the values for the other structure and material parameters are kept unchanged. The frequency response is obtained from the impulse response by fast Fourier transformation. The modeling results are presented and discussed, and design considerations, especially for high speed operation at 10 Gbit/s, are further analyzed.

Xiao, Yegao↗

Receptivity of Hypersonic Boundary Layers Due to Acoustic Disturbances over Blunt Cone

The transition process induced by the interaction of acoustic disturbances in the free-stream with boundary layers over a 5-degree straight cone and a wedge with blunt tips is numerically investigated at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves the Navier-Stokes equations are solved in axisymmetric coordinates. The governing equations are solved using the 5th -order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, acoustic disturbances are introduced at the outer boundary of the computational domain and unsteady simulations are performed. Generation and evolution of instability waves and the receptivity of boundary layer to slow and fast acoustic waves are investigated. The mean flow data are compared with the experimental results. The results show that the instability waves are generated near the leading edge and the non-parallel effects are stronger near the nose region for the flow over the cone than that over a wedge. It is also found that the boundary layer is much more receptive to slow acoustic wave (by almost a factor of 67) as compared to the fast wave.

Kara, K.↗

The Fast Scattering Code (FSC): Validation Studies and Program Guidelines

The Fast Scattering Code (FSC) is a frequency domain noise prediction program developed at the NASA Langley Research Center (LaRC) to simulate the acoustic field produced by the interaction of known, time harmonic incident sound with bodies of arbitrary shape and surface impedance immersed in a potential flow. The code uses the equivalent source method (ESM) to solve an exterior 3-D Helmholtz boundary value problem (BVP) by expanding the scattered acoustic pressure field into a series of point sources distributed on a fictitious surface placed inside the actual scatterer. This work provides additional code validation studies and illustrates the range of code parameters that produce accurate results with minimal computational costs. Systematic noise prediction studies are presented in which monopole generated incident sound is scattered by simple geometric shapes - spheres (acoustically hard and soft surfaces), oblate spheroids, flat disk, and flat plates with various edge topologies. Comparisons between FSC simulations and analytical results and experimental data are presented.

Tinetti, Ana F.↗

A Microfabricated 8-40 GHz Dual-Polarized Reflector Feed

Planar antennas based on tightly coupled dipole arrays (also known as a current sheet antenna or CSA) are amenable for use as electronically scanned phased arrays. They are capable of performance nearing a decade of bandwidth. These antennas have been demonstrated in many implementations at frequencies below 18 GHz. This paper describes the implementation using a relatively new multi-layer microfabrication process resulting in a small, 6x6 element, dual-linear polarized array with beamformer that operates from 8 to 40 GHz. The beamformer includes baluns that feed the dual-polarized differential antenna elements and reactive splitter networks that also cover the full frequency range of operation. This antenna array serves as a reflector feed for a multi-band instrument designed to measure snow water equivalent (SWE) from airborne platforms. The instrument has both radar and radiome try capability at multiple frequencies. Scattering-parameter and time-domain measurements have been used to characterize the array feed. Radiation patterns of the antenna have been measured and are compared to simulation. To the best of the authors' knowledge, this work represents the most integrated multi-octave millimeter-wave antenna feed fabricated to date.

current sheet antenna↗

Dynamic Modeling from Flight Data with Unknown Time Skews

A method for estimating dynamic model parameters from flight data with unknown time skews is described and demonstrated. The method combines data reconstruction, nonlinear optimization, and equation-error parameter estimation in the frequency domain to accurately estimate both dynamic model parameters and the relative time skews in the data. Data from a nonlinear F-16 aircraft simulation with realistic noise, instrumentation errors, and arbitrary time skews were used to demonstrate the approach. The approach was further evaluated using flight data from a subscale jet transport aircraft, where the measured data were known to have relative time skews. Comparison of modeling results obtained from time-skewed and time-synchronized data showed that the method accurately estimates both dynamic model parameters and relative time skew parameters from flight data with unknown time skews.

Morelli, Eugene A.↗

Application of radiative boundary conditions to nonhydrostatic primitive equation models

A procedure is presented for applying a radiative boundary condition to the lateral boundary of a limited domain model. The model considered, the Terminal Area Simulation System (TASS), has nonhydrostatic compressible and unsteady governing equations. The primary applications to date have been studies of downbursts, wind shear, thunderstorms and mass fires. A time-splitting integration is used with TASS, with the acoustic terms being integrated separately with a smaller time step. A first-order closure approximation is employed to ensure subgrid turbulence closure. TASS applications for two-dimensional axisymmetric and three-dimensional simulations are described. The radiative boundary conditions were applied in the two-dimensional case to the temperature, pressure and non-normal velocity components. Solutions in both cases avoided runaway circulations and experienced minimal domain-wide mass loss. Also, distorsion was minimized at the boundaries.

Proctor, F. H.↗

3D Multistage Simulation of Each Component of the GE90 Turbofan Engine

A 3D multistage simulation of each component of the GE90 Turbofan engine has been made. This includes 49 blade rows. A coupled simulation of all blade rows will be made very soon. The simulation is running using two levels of parallelism. The first level is on a blade row basis with information shared using files. The second level is using a grid domain decomposition with information shared using MPI. Timings will be shown for running on the SP2, an SGI Origin and a distributed system of HP workstations. On the HP workstations, the CHIMP version of MPI is used, with queuing supplied by LSF (Load Sharing Facility). A script-based control system is used to ensure reliability. An MPEG movie illustrating the flow simulation of the engine has been created using PV3, a parallel visualization library created by Bob Haimes of MIT. PVM is used to create a virtual machine from 10 HP workstations and display on an SGI workstation. A representative component simulation will be compared to rig data to demonstrate its usefulness in turbomachinery design and analysis.

Turner, Mark↗