Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical fluid dynamics”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

A three-dimensional model of co-rotating streams in the solar wind. 2: Hydrodynamic streams

Theoretical aspects of corotating solar wind dynamics on a global scale are explored by means of numerical simulations executed with a nonlinear, inviscid, adiabatic, single-fluid, three-dimensional (3-D) hydrodynamic formulation. A simple, hypothetical 3-D stream structure is defined on a source surface located at 35 solar radius and carefully documents its evolution to 1 AU under the influence of solar rotation. By manipulating the structure of this prototype configuration at the source surface, it is possible to elucidate the factors most strongly affecting stream evolution: (1) the intrinsic correlations among density, temperature, and velocity existing near the source; (2) the amplitude of the stream; (3) the longitudinal breadth of the stream; (4) the latitudinal breadth of the stream; and (5) the heliographic latitude of the centroid of the stream.

Pizzo, V. J.↗
NASA NTRS1981

Comparisons of some large scientific computers

In 1975, the National Aeronautics and Space Administration (NASA) began studies to assess the technical and economic feasibility of developing a computer having sustained computational speed of one billion floating point operations per second and a working memory of at least 240 million words. Such a powerful computer would allow computational aerodynamics to play a major role in aeronautical design and advanced fluid dynamics research. Based on favorable results from these studies, NASA proceeded with developmental plans. The computer was named the Numerical Aerodynamic Simulator (NAS). To help insure that the estimated cost, schedule, and technical scope were realistic, a brief study was made of past large scientific computers. Large discrepancies between inception and operation in scope, cost, or schedule were studied so that they could be minimized with NASA's proposed new compter. The main computers studied were the ILLIAC IV, STAR 100, Parallel Element Processor Ensemble (PEPE), and Shuttle Mission Simulator (SMS) computer. Comparison data on memory and speed were also obtained on the IBM 650, 704, 7090, 360-50, 360-67, 360-91, and 370-195; the CDC 6400, 6600, 7600, CYBER 203, and CYBER 205; CRAY 1; and the Advanced Scientific Computer (ASC). A few lessons learned conclude the report.

Credeur, K. R.↗
NASA NTRS1991

CFD Activity at Aerojet Related to Seals and Fluid Film Bearing

Computational Fluid Dynamics (CFD) activities related to seals and fluid film bearings are presented. Among the topics addressed are the following: Aerovisc Numeric and its capabilities; Recent Seal Applications; and Future Code Developments.

Bache, George E.↗
NASA NTRS1992

Development of a CFD code for internal flows in liquid fueled engines

To support the design efforts of engines in which liquid propellants are used, one is often required to analyze incompressible, two-dimensional or axisymmetric flows within ducts and cavities with rotating walls and complicated geometries. The steady-state solution is of interest in most cases. This code is intended to provide a tool for efficient computational fluid dynamics (CFD) analysis of such flow problems, taking advantage of the artificial compressibility concept and the Beam-Warming numerical scheme modified for second order accurate, implicit boundary conditions. These concepts ensure a stable, robust, and accurate algorithm due to the reduced speed of sound and the accuracy of the boundary conditions. The code is dedicated only to two-dimensional or axisymmetric flows with or without swirl. Three-dimensional computation is excluded to increase efficiency and speed. This paper briefly presents the theory of the code, as well as several benchmark applications with comparison to well known analytical solutions. In all these test cases, the code produced accurate results.

Dakhoul, Youssef↗
NASA NTRS1991

h-p adaptive finite element methods in computational fluid dynamics

The principal ideas of h-p adaptive finite element methods for fluid dynamics problems are discussed. Applications include acoustics, compressible Euler and both compressible and incompressible Navier-Stokes equations. Several numerical examples illustrate the presented concepts.

Oden, J. T.↗
NASA NTRS1993

Basic studies of baroclinic flows

A fully nonlinear 3-dimensional numerical model (GEOSIM), previously developed and validated for several cases of geophysical fluid flow, has been used to investigate the dynamical behavior of laboratory experiments of fluid flows similar to those of the Earth's atmosphere. The phenomena investigated are amplitude vacillation, and the response of the fluid system to uneven heating and cooling. The previous year's work included hysteresis in the transition between axisymmetric and wave flow. Investigation is also continuing of the flows in the Geophysical Fluid Flow Cell (GFFC), a low-gravity Spacelab experiment. Much of the effort in the past year has been spent in validation of the model under a wide range of external parameters including nonlinear flow regimes. With the implementation of a 3-dimensional upwind differencing scheme, higher spectral resolution, and a shorter time step, the model has been found capable of predicting the majority of flow regimes observed in one complete series of baroclinic annulus experiments of Pfeffer and co-workers. Detailed analysis of amplitude vacillation has revealed that the phase splitting described in the laboratory experiments occurs in some but not all cases. Through the use of animation of the models output, a vivid 3-dimensional view of the phase splitting was shown to the audience of the Southeastern Geophysical Fluid Dynamics Conference in March of this year. A study on interannual variability was made using GEOSIM with periodic variations in the thermal forcing. Thus far, the model has not predicted a chaotic behavior as observed in the experiments, although there is a sensitivity in the wavenumber selection to the initial conditions. Work on this subject, and on annulus experiments with non-axisymmetric thermal heating, will continue. The comparison of GEOSIM's predictions will result from the Spacelab 3 GFFC experiments continued over the past year, on a 'back-burner' basis. At this point, the study (in the form of a draft of a journal article) is nearly completed. The results from GEOSIM compared very well with the experiments, and the use of the model allows the demonstration of flow mechanics that were not possible with the experimental data. For example, animation of the model output shows that the forking of the spiral bands is a transient phenomenon, due to the differential east-west propagation of convection bands from different latitudes.

Miller, Timothy L.↗
NASA NTRS1995

Comparison of tropical pacific temperature and current simulations with two vertical mixing schemes embedded in an ocean general circulation model and reference to observations

The Pacanowski-Philander (PP) and Mellor-Yamada (MY) parameterization models of vertical mixing by turbulent processes were embedded in the Geophysical Fluid Dynamics Laboratory high-resolution ocean general circulation model of the tropical Pacific Ocean. All other facets of the numerical simulations were the same. Simulations were made for the 1987-1988 period. At the equator the MY simulation produced near-surface temperatures more uniform with depth, a deeper thermocline, a deeper core speed of the Equatorial Undercurrent, and a South Equatorial Current with greater vertical thickness compared with that computed with the PP method. Along 140 deg W, between 5 deg N and 10 deg N, both simulations were the same. Moored buoy current and temperature observations had been recorded by the Pacific Marine Environmental Laboratory at three sites (165 deg E, 140 deg W, 110 deg W) along the equator and at three sites (5 deg N, 7 deg N, 9 deg N) along 140 deg W. Simulated temperatures were lower than those observed in the near-surface layer and higher than those observed in the thermocline. Temperature simulations were in better agreement with observations compared to current simulations. At the equator, PP current and temperature simulations were more representative of the observations than MY simulations.

Halpern, David↗
NASA NTRS2003

Assessing Impacts of Global Warming on Tropical Cyclone Tracks

A new approach is proposed to assess the possible impacts of the global climate change on tropical cyclone (TC) tracks in the western North Pacific (WNP) basin. The idea is based on the premise that the future change of TC track characteristics is primarily determined by changes in large-scale environmental steering flows. It is demonstrated that the main characteristics of the current climatology of TC tracks can be derived from the climatological mean velocity field of TC motion by using a trajectory model. The climatological mean velocity of TC motion, which is composed of the large-scale steering and beta drift, is determined on each grid of the basin. The mean beta drift is estimated from the best track data, and the mean large-scale steering flow is computed from the NCEP/NCAR reanalysis for the current climate state. The derived mean beta drift agrees well with the results of previous observational and numerical studies in terms of its direction and magnitude. The outputs of experiments A2 and B2 of the Geophysical Fluid Dynamics Laboratory (GFDL) R30 climate model suggest that the subtropical high will be persistently weak over the western part of the WNP or shift eastward during July-September in response to the future climate change. By assuming that the mean beta drift in the future climate state is unchanged, the change in the general circulation by 2059 will decrease the TC activities in the WNP, but favor a northward shift of typical TC tracks. As a result, the storm activities in the South China Sea will decrease by about 12%, while the Japan region will experience an increase of TCs by 12-15%. During the period of 2000-2029, the tropical storms that affect the China region will increase by 5-6%, but return to the current level during 2030-2059. It is also suggested that, during the period of 2030-2059 tropical storms will more frequently affect Japan and the middle latitude region of China given that the formation locations remain the same as in the current climate state.

Wu, Li-Guang↗
NASA NTRS2007

An Introduction to Thermodynamic Performance Analysis of Aircraft Gas Turbine Engine Cycles Using the Numerical Propulsion System Simulation Code

This document is intended as an introduction to the analysis of gas turbine engine cycles using the Numerical Propulsion System Simulation (NPSS) code. It is assumed that the analyst has a firm understanding of fluid flow, gas dynamics, thermodynamics, and turbomachinery theory. The purpose of this paper is to provide for the novice the information necessary to begin cycle analysis using NPSS. This paper and the annotated example serve as a starting point and by no means cover the entire range of information and experience necessary for engine performance simulation. NPSS syntax is presented but for a more detailed explanation of the code the user is referred to the NPSS User Guide and Reference document (ref. 1).

Jones, Scott M.↗
NASA NTRS2018

Generalised Primal-Dual Grids for Unstructured Co-Volume Schemes

The generation of high-quality staggered unstructured grids is considered, leading to the development of a new optimisation-based strategy designed to construct weighted 'Regular-Power' tessellations appropriate for co-volume type numerical discretisation techniques. This new framework aims to extend the conventional Delaunay-Voronoi primal-dual structure; seeking to assemble generalised orthogonal tessellations with enhanced geometric quality. The construction of these grids is motivated by the desire to improve the performance and accuracy of numerical methods based on unstructured co-volume type schemes, including various staggered grid techniques for the simulation of fluid dynamics and hyperbolic transport. In this study, a new hybrid optimisation strategy is proposed; seeking to optimise the geometry, topology and weights associated with general, two-dimensional Regular-Power tessellations using a combination of gradient-ascent and energy-based techniques. The performance of this new method is tested experimentally, with a range of complex, multi-resolution primal-dual grids generated for various coastal and regional ocean modelling applications.

power diagrams↗
DOE OSTI2020

Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes

Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.

42 ENGINEERING↗
NASA NTRS1996

Numerical Simulation of Flow in a Whirling Annular Seal and Comparison With Experiments

The turbulent flow field inside a whirling annular seal was simulated by using SCISEAL, a three-dimensional computational fluid dynamics code. The rotor center described a circular synchronous whirl. A rotating frame transformation was used to make the problem quasi-steady. The flow field at an axial Reynolds number of 24,000 and a Taylor number of 6600 was simulated. The standard kappa-epsilon model with wall functions and the low-Reynolds-number model were used to treat turbulence. An experimentally measured velocity field was used at the inlet boundary. Numerical predictions of the velocities and the stator wall pressures compared well with experimental data. Both turbulence models yielded nearly the same results. The capability of the SCISEAL code to analyze this complex flow field was demonstrated; the isotropic turbulence models performed adequately on this nonisotropic turbulence flow.

Athavale, M. M.↗
NASA NTRS1991

Nonlinear effects in a plain journal bearing. I - Analytical study. II - Results

In the first part of this work, a numerical model is presented which couples the variable-property Reynolds equation with a rotor-dynamics model for the calculation of a plain journal bearing's nonlinear characteristics when working with a cryogenic fluid, LOX. The effects of load on the linear/nonlinear plain journal bearing characteristics are analyzed and presented in a parametric form. The second part of this work presents numerical results obtained for specific parametric-study input variables (lubricant inlet temperature, external load, angular rotational speed, and axial misalignment). Attention is given to the interrelations between pressure profiles and bearing linear and nonlinear characteristics.

Choy, F. K.↗
DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗
DOE OSTI2022

Editorial: Particulate Matter Emissions From Conventional and Reformulated Fuel Combustion: Advances in Experiments and Simulations

Particulate matter (PM or soot) emissions generated by combustion of conventional and reformulated fuels adversely impact the environment and human health (Bond et al., 2013; Landrigan et al., 2018; World Health Organization (WHO), 2016), generating public awareness and motivating efforts towards the mitigation of their harmful effects. Soot formation is one of the most complex phenomena in combustion, involving interactions between combustion chemistry, fluid mechanics, mass/heat transport, and particle dynamics, spanning different spatial and temporal scales. A fundamental understanding of soot formation process is thus necessary to achieve a strong reduction of PM emissions and design cleaner and more efficient combustion systems. This has motivated a longstanding and ongoing research activity aimed at improving our understanding of the physical and chemical processes involved in soot formation, well reviewed in these recent works (Wang and Chung, 2019; Michelsen et al., 2020; Martin et al., 2022). Despite the wide interest, the transition from gas-phase molecules to incipient soot particles is still elusive and the successive particle growth and oxidation processes are far from being fully understood, especially in conditions relevant to real-world applications. Recent advances in combustion PM emission diagnostic and computational capabilities helped in improving the predictability of fundamental chemical and aerosol models for practical applications, thus tackling some of the above-mentioned challenges. The aim of this research topic is to display the ongoing research efforts in addressing the existing gaps on particulate formation from various fuel sources and in conditions typical of practical combustion applications (e.g. flames, engines, pool fires), through both experimental and numerical approaches.

42 ENGINEERING↗
DOE OSTI2020

A space-time parallel algorithm with adaptive mesh refinement for computational fluid dynamics

This work describes a space-time parallel algorithm with space-time adaptive mesh refinement (AMR). AMR with subcycling is added to multigrid reduction-in-time (MGRIT) in order to provide solution efficient adaptive grids with a reduction in work performed on coarser grids. This algorithm is achieved by integrating two software libraries: XBraid (Parallel time integration with multigrid. https://computation.llnl.gov/projects/parallel-timeintegration-multigrid) and Chombo (Chombo software package for AMR applications—design document, 2014). The former is a parallel time integration library using multigrid and the latter is a massively parallel structured AMR library. Employing this adaptive space-time parallel algorithm is Chord (Comput Fluids 123:202–217, 2015), a computational fluid dynamics (CFD) application code for solving compressible fluid dynamics problems. For the same solution accuracy, speedups are demonstrated from the use of space-time parallelization over the time-sequential integration on Couette flow and Stokes’ second problem. On a transient Couette flow case, at least a 1.5× speedup is achieved, and with a time periodic problem, a speedup of up to 13.7× over the time-sequential case is obtained. In both cases, the speedup is achieved by adding processors and exploring additional parallelization in time. The numerical experiments show the algorithm is promising for CFD applications that can take advantage of the time parallelism. Future work will focus on improving the parallel performance and providing more tests with complex fluid dynamics to demonstrate the full potential of the algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2013

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗
NASA NTRS1991

Three-dimensional flows in a transonic compressor rotor

This study involves an experimental and numerical investigation of the three-dimensional flows in a transonic compressor rotor. A variety of data which could be used, in a complementary fashion, to validate/calibrate the computational fluid dynamics turbomachinery code and improve understanding of the flow physics, were acquired. Detailed radial survey data which consisted of total pressure, total temperature, static pressure and flow angle were obtained at stations upstream and downstream of the rotor blade. Detailed velocity and turbulence profiles were obtained upstream of the rotor and used as the upstream boundary conditions for the numerical analysis. Calibrated flush-mounted hot film probes were used to measure wall shear stress on the hub and casing walls upstream of the rotor. The blade-to-blade shear-stress angle distributions were obtained at two axial locations on the rotor casing, using flush-mounted hot film probes. A numerical analysis conducted using a three-dimensional Navier-Stokes code was compared with the experimental results.

Reid, Lonnie↗