Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Aeroacoustic Simulation of Nose Landing Gear on Adaptive Unstructured Grids With FUN3D

Numerical simulations have been performed for a partially-dressed, cavity-closed nose landing gear configuration that was tested in NASA Langley s closed-wall Basic Aerodynamic Research Tunnel (BART) and in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D, developed at NASA Langley Research center, is used to compute the unsteady flow field for this configuration. Starting with a coarse grid, a series of successively finer grids were generated using the adaptive gridding methodology available in the FUN3D code. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these computations. Time-averaged and instantaneous solutions obtained on these grids are compared with the measured data. In general, the correlation with the experimental data improves with grid refinement. A similar trend is observed for sound pressure levels obtained by using these CFD solutions as input to a FfowcsWilliams-Hawkings noise propagation code to compute the farfield noise levels. In general, the numerical solutions obtained on adapted grids compare well with the hand-tuned enriched fine grid solutions and experimental data. In addition, the grid adaption strategy discussed here simplifies the grid generation process, and results in improved computational efficiency of CFD simulations.

Vatsa, Veer N.↗
NASA NTRS2002

Properties of the Residual Stress of the Temporally Filtered Navier-Stokes Equations

The development of a unifying framework among direct numerical simulations, large-eddy simulations, and statistically averaged formulations of the Navier-Stokes equations, is of current interest. Toward that goal, the properties of the residual (subgrid-scale) stress of the temporally filtered Navier-Stokes equations are carefully examined. Causal time-domain filters, parameterized by a temporal filter width 0 less than Delta less than infinity, are considered. For several reasons, the differential forms of such filters are preferred to their corresponding integral forms; among these, storage requirements for differential forms are typically much less than for integral forms and, for some filters, are independent of Delta. The behavior of the residual stress in the limits of both vanishing and in infinite filter widths is examined. It is shown analytically that, in the limit Delta to 0, the residual stress vanishes, in which case the Navier-Stokes equations are recovered from the temporally filtered equations. Alternately, in the limit Delta to infinity, the residual stress is equivalent to the long-time averaged stress, and the Reynolds-averaged Navier-Stokes equations are recovered from the temporally filtered equations. The predicted behavior at the asymptotic limits of filter width is further validated by numerical simulations of the temporally filtered forced, viscous Burger's equation. Finally, finite filter widths are also considered, and a priori analyses of temporal similarity and temporal approximate deconvolution models of the residual stress are conducted.

Pruett, C. D.↗
NASA NTRS2013

Assessments of a Turbulence Model Based on Menter's Modification to Rotta's Two-Equation Model

The main objective of this paper is to construct a turbulence model with a more reliable second equation simulating length scale. In the present paper, we assess the length scale equation based on Menter s modification to Rotta s two-equation model. Rotta shows that a reliable second equation can be formed in an exact transport equation from the turbulent length scale L and kinetic energy. Rotta s equation is well suited for a term-by-term modeling and shows some interesting features compared to other approaches. The most important difference is that the formulation leads to a natural inclusion of higher order velocity derivatives into the source terms of the scale equation, which has the potential to enhance the capability of Reynolds-averaged Navier-Stokes (RANS) to simulate unsteady flows. The model is implemented in the PAB3D solver with complete formulation, usage methodology, and validation examples to demonstrate its capabilities. The detailed studies include grid convergence. Near-wall and shear flows cases are documented and compared with experimental and Large Eddy Simulation (LES) data. The results from this formulation are as good or better than the well-known SST turbulence model and much better than k-epsilon results. Overall, the study provides useful insights into the model capability in predicting attached and separated flows.

Abdol-Hamid, Khaled S.↗
DOE OSTI2025

Inertial collapse of a gas bubble in a shear flow near a rigid wall

Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.

Bhola, Sahil (ORCID:0000000236766302)↗
NASA NTRS1987

Development and application of computational aerothermodynamics flowfield computer codes

Multiple nozzle plume flow field is computed with a 3-D, Navier-Stokes solver. Numerical simulation is performed with a flux-split, two-factor, time asymptotic viscous flow solver of Ying and Steger. The two factor splitting provides a stable 3-D solution procedure under ideal-gas assumptions. An ad-hoc acceleration procedure that shows promise in improving the convergence rate by a factor of three for steady state problems is utilized. Computed solutions to generic problems at various altitude and flight conditions show flow field complexity and three-dimensional effects due to multiple nozzle jet interactions. Viscous, ideal gas solutions for the symmetric nozzle are compared with other numerical solutions.

Venkatapathy, Ethiraj↗
NASA NTRS1988

3-D plume flow computations with an upwind solver

In the present study multi-nozzle plume flows are computed with a three-dimensional, Navier-Stokes solver. Numerical simulations are performed with the flux-split, two-factor, time symptotic, viscous flow solver. The two factor splitting provides a stable three-dimensional solution procedure under ideal-gas assumptions. Viscous, ideal-gas solutions for a thin lip symmetrical nozzle are compared with experimental and numerical solutions. Computed solutions to axisymmetric and three-dimensional, multi-nozzle problems at various altitudes and flight conditions demonstrate flow field complexity and three-dimensional effects.

Venkatapathy, Ethiraj↗
NASA NTRS1991

Methods Of Simulation Of Incompressible Flow

Report discusses methods for numerical simulation of Navier-Stokes equations of viscous, incompressible flow, with emphasis on pseudocompressibility method. Advanced computational techniques reviewed contribute to improved designs by rapidly providing detailed information about conditions at all points in flow fields of proposed designs.

Dwak, D.↗
NASA NTRS1993

Head-on collision of drops: A numerical investigation

The head-on collision of equal sized drops is studied by full numerical simulations. The Navier-Stokes equations are solved for fluid motion both inside and outside the drops using a front tracking/finite difference technique. The drops are accelerated toward each other by a body force that is turned off before the drops collide. When the drops collide, the fluid between them is pushed outward leaving a thin later bounded by the drop surface. This layer gets progressively thinner as the drops continue to deform and in several of the calculations this double layer is artificially removed once it is thin enough, thus modeling rupture. If no rupture takes place, the drops always rebound, but if the film is ruptured the drops may coalesce permanently or coalesce temporarily and then split again.

Nobari, M. R.↗
NASA NTRS1995

Nozzle Thrust Optimization While Reducing Jet Noise

A Bluebell nozzle design concept is proposed for jet noise reduction with minimal thrust loss or even thrust augmentation. A Bluebell nozzle has a sinusoidal lip line edge (chevrons) and a sinusoidal cross section shape with linear amplitude increasing downstream in the divergent nozzle part (corrugations). The experimental tests of several Bluebell nozzle designs have shown nose reduction relative to a convergent-divergent round nozzle with design exhaust number M(e) = 1.5. The best design provides an acoustic benefit near 4dB with about 1 percent thrust augmentation. For subsonic flow ((M(e)= 0.6)), the tests indicated that the present method for design of Bluebell nozzles gives less acoustic benefit and in most cases jet noise increased. The proposed designs incorporate analytical theory and 2D and 3D numerical simulations. Full Navier-Stokes and Euler solvers were utilized. Boundary layer effects were used. Several different designs were accounted for in the Euler applications.

Seiner, J. M.↗
NASA NTRS1996

NPARC Code Upgraded with Two-Equation Turbulence Models

The National PARC (NPARC) Alliance was established by the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center to provide the U.S. aeropropulsion community with a reliable Navier-Stokes code for simulating the nonrotating components of propulsion systems. Recent improvements to the turbulence model capabilities of the NPARC code have significantly improved its capability to simulate turbulent flows. Specifically, the Chien k-epsilon and Wilcox k-omega turbulence models were implemented at Lewis. Lewis researchers installed the Chien k-epsilon model into NPARC to improve the code's ability to calculate turbulent flows with attached wall boundary layers and free shear layers. Calculations with NPARC have demonstrated that the Chien k-epsilon model provides more accurate calculations than those obtained with algebraic models previously available in the code. Grid sensitivity investigations have shown that computational grids must be packed against the solid walls such that the first point off of the wall is placed in the laminar sublayer. In addition, matching the boundary layer and momentum thicknesses entering mixing regions is necessary for an accurate prediction of the free shear-layer growth.

Source record↗
NASA NTRS2005

Partially-Averaged Navier Stokes Model for Turbulence: Implementation and Validation

Partially-averaged Navier Stokes (PANS) is a suite of turbulence closure models of various modeled-to-resolved scale ratios ranging from Reynolds-averaged Navier Stokes (RANS) to Navier-Stokes (direct numerical simulations). The objective of PANS, like hybrid models, is to resolve large scale structures at reasonable computational expense. The modeled-to-resolved scale ratio or the level of physical resolution in PANS is quantified by two parameters: the unresolved-to-total ratios of kinetic energy (f(sub k)) and dissipation (f(sub epsilon)). The unresolved-scale stress is modeled with the Boussinesq approximation and modeled transport equations are solved for the unresolved kinetic energy and dissipation. In this paper, we first present a brief discussion of the PANS philosophy followed by a description of the implementation procedure and finally perform preliminary evaluation in benchmark problems.

Girimaji, Sharath S.↗
NASA NTRS2010

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon similar correlation techniques that have been is use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. Recent entry flight testing performed with the Orbiter Discovery now provides a means to validate this engineering correlation approach to higher confidence. These results only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, and to discuss the characteristics which are desirable in a next generation advancement. The details of the paper will review the experimental facilities and techniques that were utilized to perform the implementation of the Orbiter RTF BLT Vsn 2 prediction capability. Statistically significant results for multiple engineering correlations from a ground testing campaign will be reviewed in order to describe why only certain correlations were selected for complete implementation to support the Shuttle Program. Historical Orbiter flight data on early boundary layer transition due to protruding gap fillers will be described in relation to the selected empirical correlations. In addition, Orbiter entry flight testing results from the BLT Flight Experiment will be discussed in relation to these correlations. Applicability of such correlations to the entry design problem will be reviewed, and finally a perspective on the desirable characteristics for a next generation capability based on high fidelity physical models will be provided.

Campbell, Charles H.↗
NASA NTRS2010

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon the framework of correlation techniques that have been in use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. The basic tenets that guided the formulation and implementation of the Orbiter Return To Flight boundary layer transition prediction capability will be reviewed as a recommended format for future empirical correlation efforts. The validity of this approach has since been demonstrated by very favorable comparison of recent entry flight testing performed with the Orbiter Discovery, which will be graphically summarized. These flight data can provide a means to validate discrete protuberance engineering correlation approaches as well as high fidelity prediction methods to higher confidence. The results of these Orbiter engineering and flight test activities only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information-related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, to illustrate how quantitative statistical evaluations of empirical correlations can and should be used to assess accuracy and to discuss what the authors' perceive as a recent heightened interest in the application of high fidelity numerical modeling of boundary layer transition. Concrete results will also be developed related to empirical boundary layer transition onset correlations. This will include assessment of the discrete protuberance boundary layer transition onset data assembled for the Orbiter configuration during post-Columbia Return To Flight. Assessment of these data will conclude that momentum thickness Reynolds number based correlations have superior coefficients and uncertainty in comparison to roughness height based Reynolds numbers, aka Re(sub k) or Re(sub kk). In addition, linear regression results from roughness height Reynolds number based correlations will be evaluated, leading to a hypothesis that non-continuum effects play a role in the processes associated with incipient boundary layer transition on discrete protuberances.

Campbell, Charles H.↗
NASA NTRS1998

Preliminary Evaluation of Nonlinear Effects on TCA Flutter

The objective of this study is to investigate the effect of nonlinear aerodynamics, especially at high angles-of-attack with leading-edge separation, on the TCA flutter properties at transonic speeds. In order to achieve that objective, flutter simulations with Navier-Stokes CFD must be performed. To this end, time-marching Navier-Stokes solutions are computed for the TCA wing/body configuration at high angles-of-attack in transonic flight regimes. The approach is to perform non-linear flutter calculations on the TCA at two angles-of-attack, the first one being a case with attached flow (a=2.8 degrees) and the second one being a high angle-of-attack case with a wing leading edge vortex (a=12.11 degrees). Comparisons of the resulting histories and frequency damping information for both angles-of-attack will evaluate the impact of high-alpha aerodynamics on flutter.

Arslan, Alan E.↗
NASA NTRSDate not supplied

Dynamic Aeroelastic Flight Dynamic Modeling of Mach 0.745 Transonic Truss-Braced Wing

This paper presents a modeling approach for dynamic aeroelastic flight dynamic analysis of the Mach0.745 Transonic Truss-Braced Wing. The modeling approach is based on a transonic correction method to correct the Theodorsen’s theory for transonic flow. CFD unsteady Reynolds-averaged Navier-Stokes equations (RANS) simulations are conducted using FUN3D solver for a series of wing sections from the Mach 0.745Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen’s theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as functions of the reduced frequency. These derivatives are used to compute the unsteady lift and pitching moment contributions by the angle of attack, pitch rate, and roll rate. They are then approximated using a frequency-domain regression to obtain the dynamic stability derivatives for the Mach 0.745 Transonic Truss-Braced Wing. The structural dynamic mode shapes of the Mach 0.745 TTBW are extracted from a NAS-TRAN finite-element model. These mode shapes are used to compute the generalized unsteady aerodynamic forces. The aerodynamic mass, damping, and stiffness and the aerodynamic lag states are constructed to couple the dynamic aeroelastic contribution to the flight dynamic model of the Mach 0.745 Transonic Truss-Braced Wing. The coupled dynamic aeroelastic flight dynamic equations of motion are formulated. The eigenvalues of the coupled system are computed. All the flight dynamic modes and structural dynamic modes are stable at Match 0.745. The effect of transonic aerodynamics generally causes all the dynamic modes to have lower damping values

Flight Dynamics↗
NASA NTRS2024

Measurements and Computational Analysis of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow Experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This report presents and discusses the results of a recent experiment with the Juncture Flow Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of the test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model angle of incidence. Laser Doppler velocimetry (LDV) measurements were made at three model angles of incidence (0° : fully attached, 1° : incipient separation, and 5° : separated flow) and for each angle, mean-flow and Reynolds-stress data were obtained on the fuselage and at several locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and include model and tunnel wall static pressures, tunnel wall boundary-layer profiles, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. A companion set of Reynolds-averaged Navier-Stokes (RANS) CFD simulations were also performed and those results were used to support the discussion of the experimental test results.

Juncture Flow↗
NASA NTRS2023

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗
NASA NTRSDate not supplied

Classical Flutter Analysis of X-57 Aircraft Mod II, III, and IV Configurations

To support airworthiness assessment of the X-57 electric propulsion demonstrator aircraft, the National Aeronautics and Space Administration (NASA) completed classical flutter analysis of the Mod II, Mod III, and Mod IV aircraft configurations. Aeroelastic analysis was primarily performed using the ZAERO analysis code developed by Zona Technology. Supplemental computational fluid dynamics (CFD) flutter analysis for the Mod III and Mod IV configurations was performed using NASA’s FUN3D software to verify the aeroelastic behavior at deployed flaps conditions. The NASTRAN finite element models (FEMs) used for aeroelastic analysis were developed by NASA. The final Mod II FEM utilized shell elements for the wing and fuselage, and beam elements for the tail and wing control surfaces. For Mod III and Mod IV aircraft models, shell element wing FEMs developed throughout the Mod III/IV wing design cycle were utilized. The fuselage and empennage FEM components were shared between the Mod II, III, and IV aircraft models. ZAERO analyses predicted flutter boundaries were >60% beyond the X-57 flight-test envelope never exceed airspeed (V NE ) for all configurations. Aeroelastic Reynolds-averaged Navier–Stokes (RANS) CFD simulations using FUN3D of the Mod III/IV configurations predicted flap mode instability did not occur until at least speeds of 60% higher than the Mod III/IV maximum flap deployment airspeed (V FE ).

Keerti K. Bhamidipati↗