Search NASA⌕ Search

SEARCH · Search NASA

Results for “unsteady simulations”

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 199 records · Page 11

Improved Depiction Of Computed Unsteady Flows

PLOT4D and STREAKER programs increase ease of production and utility of animated depictions of unsteady flows simulated numerically by advanced computational fluid dynamics. Programs developed to maximize information displayed while minimizing amount stored. PLOT4D, an extended version of PLOT3D, used to depict steady three-dimensional flows, incorporates features making it more computationally efficient than to visualize unsteady flow as succession of steady flows, animating frame by frame. PLOT4D postprocesses numerical output of fluid-dynamical computation into maps of flow quantities on two-coordinate surfaces depicting cut through flow field or conform to body surface. STREAKER generates streaklines, close approximations of patterns traced by streams of dye or smoke released into simulated flows. Streaklines defined via spatial succession of endpoints of pathlines of massless particles released into flow. STREAKER computes endpoints by fourth-order Runge-Kutta integration of velocity vectors.

Smith, Merritt H.↗

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Newman, James C., III↗

Large-Amplitude, High-Rate Roll Oscillations of a 65 deg Delta Wing at High Incidence

The IAR/WL 65 deg delta wing experimental results provide both detail pressure measurements and a wide range of flow conditions covering from simple attached flow, through fully developed vortex and vortex burst flow, up to fully-stalled flow at very high incidence. Thus, the Computational Unsteady Aerodynamics researchers can use it at different level of validating the corresponding code. In this section a range of CFD results are provided for the 65 deg delta wing at selected flow conditions. The time-dependent, three-dimensional, Reynolds-averaged, Navier-Stokes (RANS) equations are used to numerically simulate the unsteady vertical flow. Two sting angles and two large- amplitude, high-rate, forced-roll motions and a damped free-to-roll motion are presented. The free-to-roll motion is computed by coupling the time-dependent RANS equations to the flight dynamic equation of motion. The computed results are compared with experimental pressures, forces, moments and roll angle time history. In addition, surface and off-surface flow particle streaks are also presented.

Chaderjian, Neal M.↗

Comparison of Computed and Measured Vortex Evolution for a UH-60A Rotor in Forward Flight

A Computational Fluid Dynamics (CFD) simulation using the Navier-Stokes equations was performed to determine the evolutionary and dynamical characteristics of the vortex flowfield for a highly flexible aeroelastic UH-60A rotor in forward flight. The experimental wake data were acquired using Particle Image Velocimetry (PIV) during a test of the fullscale UH-60A rotor in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel. The PIV measurements were made in a stationary cross-flow plane at 90 deg rotor azimuth. The CFD simulation was performed using the OVERFLOW CFD solver loosely coupled with the rotorcraft comprehensive code CAMRAD II. Characteristics of vortices captured in the PIV plane from different blades are compared with CFD calculations. The blade airloads were calculated using two different turbulence models. A limited spatial, temporal, and CFD/comprehensive-code coupling sensitivity analysis was performed in order to verify the unsteady helicopter simulations with a moving rotor grid system.

Evolution↗

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗

Numerical Investigation of Multi-Rotor Configurations

The development of a conservative, parallel blade-element rotor model and its implementation into an adaptive Cartesian method is briefly reviewed. The model is first validated and then applied in edgewise forward flight for a variety of configurations. Validation studies comparing torque required for an isolated rotor match XV-15 wind tunnel trends well. Simulations of a simple quadrotor in “×” and “+” configurations are performed, demonstrating how the trailing rotor(s) lose thrust because of unfavorable aerodynamic interactions with the forward rotor(s), while power required remains near constant. A modified “+1” configuration is also examined, where the lateral rotors are brought inboard by 25%. This configuration showed a further 2% performance improvement over the “+” arrangement. The final paper will also include validation of a tandem rotor configuration with wind tunnel data as well as results and analysis of unsteady quadcopter simulations. In addition, a 4th configuration (× will be analyzed.

Multirotor↗

RCA: Modeling Tools for CLmax Prediction

TACP06 is aimed at further research in maturing eddy-resolving modeling tools, in their accuracy and efficiency, and demonstrating the tools for application to the prediction of aircraft maximum lift (CLmax). The goal is to achieve CLmaxprediction accuracy of the same level as in aircraft certification flight tests. This will require comparison of various modeling approaches against experimental results from the planned wind tunnel tests, down selecting an approach, further maturing the technology and validating against flight test data. Advances in computational fluid dynamics (CFD) over the last several decadeshas fundamentally changed the aerospace design process. Advanced simulation capabilities not only enable reductions in ground-based and in-flight testing requirements, but also provide added physical insight, enable superior designs at reduced cost and risk, and open upnew frontiers in aerospace vehicle design and performance. The NASA sponsored CFD Vision 2030 Study, while highlighting these accomplishments, brought out several challenges and deficiencies in the computational technology and developed a research roadmap for advancing the state-of-the-art required for enabling NASA missions in aeronautics and space applications. Based on TC TACP01 research completed in May 2018, it became clear that unsteady flow simulation capability is needed for expanding the role of CFD across the entire flight envelope to enable design of future advanced aircraft and space vehicles, and that capability will be the primary outcome of this technical challenge.

Mujeeb R Malik↗

A Comparative Study of Simulated and Measured Main Landing Gear Noise for Large Civil Transports

Computational results for the NASA 26%-scale model of a six-wheel main landing gear with and without a toboggan-shaped noise reduction fairing are presented. The model is a high-fidelity representation of a Boeing 777-200 aircraft main landing gear. A lattice Boltzmann method was used to simulate the unsteady flow around the model in isolation. The computations were conducted in free-air at a Mach number of 0.17, matching a recent acoustic test of the same gear model in the Virginia Tech Stability Wind Tunnel in its anechoic configuration. Results obtained on a set of grids with successively finer spatial resolution demonstrate the challenge in resolving/capturing the flow field for the smaller components of the gear and their associated interactions, and the resulting effects on the high-frequency segment of the farfield noise spectrum. Farfield noise spectra were computed based on an FWH integral approach, with simulated pressures on the model solid surfaces or flow-field data extracted on a set of permeable surfaces enclosing the model as input. Comparison of these spectra with microphone array measurements obtained in the tunnel indicated that, for the present complex gear model, the permeable surfaces provide a more accurate representation of farfield noise, suggesting that volumetric effects are not negligible. The present study also demonstrates that good agreement between simulated and measured farfield noise can be achieved if consistent post-processing is applied to both physical and synthetic pressure records at array microphone locations.

Konig, Benedikt↗

Computational strategies for three-dimensional flow simulations on distributed computer systems

An increasing amount of research activity in computational fluid dynamics has been devoted to the development of efficient algorithms for parallel computing systems. The increasing performance to price ratio of engineering workstations has led to research to development procedures for implementing a parallel computing system composed of distributed workstations. This thesis proposal outlines an ongoing research program to develop efficient strategies for performing three-dimensional flow analysis on distributed computing systems. The PVM parallel programming interface was used to modify an existing three-dimensional flow solver, the TEAM code developed by Lockheed for the Air Force, to function as a parallel flow solver on clusters of workstations. Steady flow solutions were generated for three different wing and body geometries to validate the code and evaluate code performance. The proposed research will extend the parallel code development to determine the most efficient strategies for unsteady flow simulations.

Weed, Richard Allen↗

Unsteady Phenomena During Operation of the SSME Fuel Flowmeter

This report describes a part of the analysis carried in support of the SSME (Space Shuttle Main Engine) Fuel Flowmeter redesign, addressing an intensely researched phenomenon known as "shifting" of the flowmeter constant value. It consists of a sudden change in the flowmeter indication, which occurs simultaneously with the onset of an oscillatory variation of the rotor speed. The change in the flowmeter indications does not correspond to a real change in the volumetric flow through the device. Several causes have been investigated in detail, in the past, without conclusive evidence towards a cause of this phenomenon. The present analysis addresses the flow physics through the flowmeter by assembling results from 3-D CFD (computational fluid dynamics) calculations, airfoil C(sub D)/C(sub L) performance curves and mass moment of inertia characteristics of the rotor into a synergistic calculation which simulates the unsteady regime of the flowmeter operation. The results show that the 4-bladed rotor interacts with the periodic flow pattern created behind the flow straightener upstream in a manner that generates a steady, periodic fluctuation in the rotor's speed. The amplitude of this fluctuation is significantly smaller than the 0.5% of mean speed threshold which constitutes a flight operational limit. When manufacturing errors occur, however, the fluctuations are amplified and can generate a significant apparent change in the flowmeter indication. Two types of possible fabrication errors-which can occur even for parts fabricated within the accepted tolerances for the blade airfoil-are presented, together with their effect on the flowmeter operation.

Marcu, Bogdan↗

CFD Simulations for the Effect of Unsteady Wakes on the Boundary Layer of a Highly Loaded Low-Pressure Turbine Airfoil (L1A)

This report is the third part of a three-part final report of research performed under an NRA cooperative Agreement contract. The first part was published as NASA/CR-2012-217415. The second part was published as NASA/CR-2012-217416. The study of the very high lift low-pressure turbine airfoil L1A in the presence of unsteady wakes was performed computationally and compared against experimental results. The experiments were conducted in a low speed wind tunnel under high (4.9%) and then low (0.6%) freestream turbulence intensity for Reynolds number equal to 25,000 and 50,000. The experimental and computational data have shown that in cases without wakes, the boundary layer separated without reattachment. The CFD was done with LES and URANS utilizing the finite-volume code ANSYS Fluent (ANSYS, Inc.) under the same freestream turbulence and Reynolds number conditions as the experiment but only at a rod to blade spacing of 1. With wakes, separation was largely suppressed, particularly if the wake passing frequency was sufficiently high. This was validated in the 3D CFD efforts by comparing the experimental results for the pressure coefficients and velocity profiles, which were reasonable for all cases examined. The 2D CFD efforts failed to capture the three dimensionality effects of the wake and thus were less consistent with the experimental data. The effect of the freestream turbulence intensity levels also showed a little more consistency with the experimental data at higher intensities when compared with the low intensity cases. Additional cases with higher wake passing frequencies which were not run experimentally were simulated. The results showed that an initial 25% increase from the experimental wake passing greatly reduced the size of the separation bubble, nearly completely suppressing it.

Vinci, Samuel, J.↗

Simulating a Three-Dimensional Flow in Pipes

Numerical method simulates in unsteady, incompressible flow in axisymmetric pipe in three dimensions. Method relies on velocity-vector expansion technique in which each vector function in expansion set is divergence-free and satisfies boundary conditions for viscous flow.

Leonard, A.↗

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 2nd, California State University, Long Beach, CA, January 17-20, 1983, Proceedings

The present conference covers topics concerning the measurement and calculation of interactive flows, together with problems posed by subsonic and transonic wings, missiles, and ships. Discussions are presented on the time-dependent finite difference simulation of unsteady interactive flows, Navier-Stokes equation methods, numerical solutions for spatially periodic boundary layers, the application of unsteady laminar tripple deck theory to viscous-inviscid interaction, the coupling of boundary layer and Euler equation solutions, and viscous-inviscid flow interactions. Also discussed are leading and trailing edge flows, three-dimensional wing flows, small disturbance calculations including entropy corrections, an inviscid computational method for tactical missiles, and boundary layer and flow separation characteristics of bodies of revolution at incidence.

Source record↗

Dynamic response characteristics of a circulation control rotor model pneumatic system

Numerical and experimental simulation of unsteady airflow through the control valve and slotted air duct of a circulation control rotor is described. The numerical analysis involves the solution of the quasi-one-dimensional compressible fluid-dynamic equations in the blade air duct together with the coupled isentropic flow equations for flow into the blade through the valve and out of the blade through the Coanda slot. Numerical solutions are compared with basic experimental results obtained for a mockup of a circulation control rotor and its pneumatic valving system. The pneumodynamic phenomena that were observed are discussed with particular emphasis on the characteristic system time lags associated with the response of the flow variables to transient and periodic control valve inputs.

Watkins, C. B.↗

Simulation of self-induced unsteady motion in the near wake of a Joukowski airfoil

The unsteady Navier-Stokes analysis is shown to be capable of analyzing the massively separated, persistently unsteady flow in the post-stall regime of a Joukowski airfoil for an angle of attack as high as 53 degrees. The analysis has provided the detailed flow structure, showing the complex vortex interaction for this configuration. The aerodynamic coefficients for lift, drag, and moment were calculated. So far only the spatial structure of the vortex interaction was computed. It is now important to potentially use the large-scale vortex interactions, an additional energy source, to improve the aerodynamic performance.

Ghia, K. N.↗

Numerical simulation of excited jet mixing layers

A numerical simulation of unsteady flow in jet mixing layers, both with and without external excitation, has been performed by solving the time-dependent compressible Navier-Stokes equations. Computations were performed on a CRAY X-MP computer using MacCormick's explicit finite difference algorithm. Different excitation methods were investigated and were shown to be very effective in controlling the well organized periodic production, shedding and pairing of large scale vortex structures. It is found that pressure excitation was generally more effective than temperature excitation, and that grid refinement results in substantial improvement in the resolution of unsteady features. The location and orientation, in addition to the frequency, of the excitation source are shown to have a significant influence on the production and interaction of large scale vortex structures in the jet mixing layer.

Scott, J. N.↗