Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Modeling of transient two-component flow using a four-point implicit method

The four-point, centered implicit scheme that is extensively used in open channel flow simulation is shown to be applicable to rapid and slow pressure transient problems in conduits with nearly single phase and two-phase flows. It is only necessary to choose the proper weighting factor value, theta, of the Courant number. For rapid pressure transients such as waterhammer, the implicit method can yield reasonable results with limited numerical dispersion and attenuation if theta is only slightly greater than the critical value of 0.5. For slower pressure gradients in single and two-phase flows, reasonable numerical solutions may be achieved for Courant number values as high as 20.

Wiggert, D. C.↗
NASA NTRS1993

Navier-Stokes, dynamics and aeroelastic computations for vortical flows, buffet and flutter applications

Research on Navier-Stokes, dynamics, and aeroelastic computations for vortical flows, buffet, and flutter applications was performed. Progress during the period from 1 Oct. 1992 to 30 Sep. 1993 is included. Papers on the following topics are included: vertical tail buffet in vortex breakdown flows; simulation of tail buffet using delta wing-vertical tail configuration; shock-vortex interaction over a 65-degree delta wing in transonic flow; supersonic vortex breakdown over a delta wing in transonic flow; and prediction and control of slender wing rock.

Kandil, Osama A.↗
NASA NTRS1995

Comparison of Implicit Schemes for the Incompressible Navier-Stokes Equations

For a computational flow simulation tool to be useful in a design environment, it must be very robust and efficient. To develop such a tool for incompressible flow applications, a number of different implicit schemes are compared for several two-dimensional flow problems in the current study. The schemes include Point-Jacobi relaxation, Gauss-Seidel line relaxation, incomplete lower-upper decomposition, and the generalized minimum residual method preconditioned with each of the three other schemes. The efficiency of the schemes is measured in terms of the computing time required to obtain a steady-state solution for the laminar flow over a backward-facing step, the flow over a NACA 4412 airfoil, and the flow over a three-element airfoil using overset grids. The flow solver used in the study is the INS2D code that solves the incompressible Navier-Stokes equations using the method of artificial compressibility and upwind differencing of the convective terms. The results show that the generalized minimum residual method preconditioned with the incomplete lower-upper factorization outperforms all other methods by at least a factor of 2.

Rogers, Stuart E.↗
NASA NTRS1994

Merging Intelligent Systems Technologies with CFD Analysis Strategies: Prototype Development at NASA Ames

A collaborative team of researchers from fields of Computational Fluid Dynamics (CFD), fluid physics, computer architectures, and computer science and knowledge engineering have begun work on a prototype system that addresses several of industry's concerns in using NASA-developed CFD codes as part of the design cycle. A major problem exists in the application of CFD technologies within the aeronautics design cycle due primarily to misunderstandings in the ranges of applicability of the various solver codes or turbulence models. Features that arise during the CFD solution process need to be discriminated and recognized as actual flow features with physical support in the geometry and flow conditions of the problem being solved, or as numerical or non-physical errors arising from mis-application of solver code and its parameters, gridding strategies, or discretization. interpolations. The fundamental concept is to develop an intelligent computational system that can accept the engineer's definition of the problem and construct an optimal CFD solution. To do this requires capturing both the knowledge of how to apply the various CFD tools and how to adapt the application of those tools to flow structures as they evolve during the flow simulation. Embedded within this adaptive system approach is the additional desire to automatically identify and quantify the quality of resolution of the pertinent flow structures, be they genuine or error-induced, and then to adjust the solution strategy accordingly. This paper discusses the status of that prototyping effort.

Thompson, David E.↗
NASA NTRS2004

Effect of Detonation through a Turbine Stage

Pulse detonation engines (PDE) have been investigated as a more efficient means of propulsion due to its constant volume combustion rather than the more often used constant pressure combustion of other propulsion systems. It has been proposed that a hybrid PDE-gas turbine engine would be a feasible means of improving the efficiency of the typical constant pressure combustion gas turbine cycle. In this proposed system, multiple pulse detonation tubes would replace the conventional combustor. Also, some of the compressor stages may be removed due to the pressure rise gained across the detonation wave. The benefits of higher thermal efficiency and reduced compressor size may come at a cost. The first question that arises is the unsteadiness in the flow created by the pulse detonation tubes. A constant pressure combustor has the advantage of supplying a steady and large mass flow rate. The use of the pulse detonation tubes will create an unsteady mass flow which will have currently unknown effects on the turbine located downstream of the combustor. Using multiple pulse detonation tubes will hopefully improve the unsteadiness. The interaction between the turbine and the shock waves exiting the tubes will also have an unknown effect. Noise levels are also a concern with this hybrid system. These unknown effects are being investigated using TURBO, an unsteady turbomachinery flow simulation code developed at Mississippi State University. A baseline case corresponding to a system using a constant pressure combustor with the same mass flow rate achieved with the pulse detonation hybrid system will be investigated first.

Ellis, Matthew T.↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2015

Magnetohydrodynamics using path or stream functions

Magnetization in highly conductive plasmas is ubiquitous to astronomical systems. Flows in such media can be described by three path functions Λ{sub α}, or, for a steady flow, by two stream functions λ{sub κ} and an additional field such as mass density ρ, velocity v, or travel time Δt. While typical analyses of a frozen magnetic field B are problem-specific and involve nonlocal gradients of the fluid element position x(t), we derive the general, local (in Λ or λ space) solution B=(∂x/∂Λ{sub α}){sub t} B-tilde {sub α}ρ/ ρ-tilde , where Lagrangian constants denoted by a tilde are directly fixed at a boundary hypersurface H-tilde on which B is known. For a steady flow, ρ-tilde B/ρ=(∂x/∂λ{sub κ}){sub Δt} B-tilde {sub κ}+v B-tilde {sub 3}/ v-tilde ; here the electric field E∼( B-tilde {sub 2}∇λ{sub 1}− B-tilde {sub 1}∇λ{sub 2})/ ρ-tilde depends only on λ{sub κ} and the boundary conditions. Illustrative special cases include compressible axisymmetric flows and incompressible flows around a sphere, showing that viscosity and compressibility enhance the magnetization and lead to thicker boundary layers. Our method is especially useful for directly computing electric fields, and for addressing upstream magnetic fields that vary in spacetime. We thus estimate the electric fields above heliospheres and magnetospheres, compute the draping of magnetic substructure around a planetary body, and demonstrate the resulting inverse polarity reversal layer. Our analysis can be immediately incorporated into existing hydrodynamic codes that are based on stream or path functions, to passively evolve the electromagnetic fields in a simulated flow. Furthermore, in such a prescription the electromagnetic fields are frozen onto the grid, so it may be developed into a fully magnetohydrodynamic, efficient simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Permeability simulation of kidney-bean shaped carbon fibers

Cost saving efforts in carbon fiber manufacturing may result in irregularly shaped cross-sections, often as a kidney-bean shape. Here, in this research, flow simulations were performed to determine transverse permeability of these kidney-bean shaped fibers. Three different unit-cell packing arrangements for the kidney-bean shaped fibers and unit-cell orientations of 0°, + 90°, and – 90° were considered over the fiber volume fraction range 0.30 – 0.75. The kidney-bean shaped fibers had two distinct hexagonally packed unit-cell geometries because of the fiber’s asymmetry with different packing efficiency. The maximum fiber packing fraction for kidney-bean shaped fibers was 0.923, which was larger than the 0.907 for circular fibers. When compared to circular fibers with hexagonal packing, the kidney-bean shaped fibers had up to a 74% reduction in permeability. The ±90° orientations for KB-fibers with hexagonal packing had an infiltration time up to 3.86 longer than the hexagonal packed circular fibers for the same volume fraction.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Modeling of heat transfer, fluid flow and solidification microstructure of nickel-base superalloy fabricated by laser powder bed fusion

Laser-Powder Bed Fusion (L-PBF), an additive manufacturing process, produces a distinctive microstructure that closely resembles the weld metal microstructure but at a much finer scale. The solidification parameters, particularly temperature gradient and solidification rate, are important to study the as-built microstructure. In the present paper, a computational framework with meso-scale resolution is developed for L-PBF of Inconel ® 718 (IN718), a Ni-base superalloy. The framework combines a powder packing model based on Discrete Element Method and a 3-D transient heat and fluid flow simulation. The latter, i.e., the molten pool model, captures the interaction between laser beam and individual powder particles including free surface evolution, surface tension and evaporation. The solidification parameters, calculated from the temperature fields, are used to assess the solidification morphology and grain size using existing theoretical models. The IN718 coupon built by L-PBF are characterized using optical and scanning electron microscopies. Finally, the experimental data of molten pool size and solidification microstructure are compared to the corresponding simulation results.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Transonic aerodynamic design experience

Advancements have occurred in transonic numerical simulation that place aerodynamic performance design into a relatively well developed status. Efficient broad band operating characteristics can be reliably developed at the conceptual design level. Recent aeroelastic and separated flow simulation results indicate that systematic consideration of an increased range of design problems appears promising. This emerging capability addresses static and dynamic structural/aerodynamic coupling and nonlinearities associated with viscous dominated flows.

Bonner, E.↗
NASA NTRS2015

CFD Simulations of the Space Launch System Ascent Aerodynamics and Booster Separation

This paper presents details of Computational Fluid Dynamic modeling of the Space Launch System during ascent. The primary focus of the paper is the flow simulation of the vehicle during ascent using the Overflow Navier-Stokes code. Computations of 739 first-stage flight conditions covering a range of Mach numbers, angles of attack, and roll angles were computed. The overset grid system contained 375 million grid points, and over 28 mil- lion CPU hours were used in the simulations. The simulations were run on the Pleiades supercomputer at the NASA Advanced Supercomputer Center at Ames Research Center. The data products from this work include integrated line-loads, surface pressure coefficients, venting pressures, and protuberance air-loads. Detailed comparisons were made of the aerodynamic performance predicted by Overflow and the wind-tunnel derived aero- dynamic database. A small number of the cases were run with two different turbulence models and with two differencing schemes. These results were used to quantify the sensitivity to the choice of the turbulence model and to the differencing scheme. The paper also introduces an effort to use the inviscid, unstructured Cartesian solver Cart3D to compute the aerodynamics during booster separation. Adaptive mesh refinement is being used to enable accurate simulations of sixteen booster-separation-motor plumes. The use of this tool is explored in preparation for building a booster-separation aerodynamic database.

Ascent Aerodynamics↗
NASA NTRSDate not supplied

Simulations of the Modular Axisymmetric Scramjet Test Rig Under Reacting Flow Conditions

Simulations of the Air Force Research Laboratory (AFRL) Modular Axisymmetric Scramjet Test Rig (MASTeR) are presented. MASTeR is a parametric test article capable of investigating various scramjet cavity flameholder designs and fueling strategies with the goal to characterize and optimize flameholding capability. In the current work, three cavity aspect ratios, and two fueling strategies (upstream and in-cavity) with ethylene at a nominal facility pressure and temperature condition are evaluated. The reacting simulations are performed, and the resulting flow characteristics are discussed. For each configuration and fueling strategy, cavity residence time, entrainment rate, and fuel-air equivalence ratio are computed. The MASTeR geometry is defined in Engineering Sketch Pad (ESP) and the simulations use a sketch-to-solution (S2S) automated unstructured grid adaptation tool in VULCAN-CFD. This tool automatically generates a simulation grid from the ESP geometry and systematically adapts it to the numerical solution based on the Hessian error estimate of a specified flow field parameter. Reynolds averaged simulations (RAS) are used with a two-equation linear eddy viscosity and diffusivity model. The resulting database can be compared with the experimental data as those become available, and explored to develop models for cavity performance for scramjet propulsion design applications.

hypersonics↗
DOE OSTI2017

Analysis of operator splitting errors for near-limit flame simulations

High-fidelity simulations of ignition, extinction and oscillatory combustion processes are of practical interest in a broad range of combustion applications. Splitting schemes, widely employed in reactive flow simulations, could fail for stiff reaction–diffusion systems exhibiting near-limit flame phenomena. The present work first employs a model perfectly stirred reactor (PSR) problem with an Arrhenius reaction term and a linear mixing term to study the effects of splitting errors on the near-limit combustion phenomena. Analysis shows that the errors induced by decoupling of the fractional steps may result in unphysical extinction or ignition. The analysis is then extended to the prediction of ignition, extinction and oscillatory combustion in unsteady PSRs of various fuel/air mixtures with a 9-species detailed mechanism for hydrogen oxidation and an 88-species skeletal mechanism for n-heptane oxidation, together with a Jacobian-based analysis for the time scales. The tested schemes include the Strang splitting, the balanced splitting, and a newly developed semi-implicit midpoint method. Results show that the semi-implicit midpoint method can accurately reproduce the dynamics of the near-limit flame phenomena and it is second-order accurate over a wide range of time step size. For the extinction and ignition processes, both the balanced splitting and midpoint method can yield accurate predictions, whereas the Strang splitting can lead to significant shifts on the ignition/extinction processes or even unphysical results. With an enriched H radical source in the inflow stream, a delay of the ignition process and the deviation on the equilibrium temperature are observed for the Strang splitting. On the contrary, the midpoint method that solves reaction and diffusion together matches the fully implicit accurate solution. The balanced splitting predicts the temperature rise correctly but with an over-predicted peak. For the sustainable and decaying oscillatory combustion from cool flames, both the Strang splitting and the midpoint method can successfully capture the dynamic behavior, whereas the balanced splitting scheme results in significant errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

Investigation A: Numerical solutions of driven vortices of binary fluid in a cylindrical geometry

The steady state laminar motion of a viscous, incompressible and binary fluid is studied for a rotating flow in a cylindrical geometry. The mathematical model employed is a cold flow simulation of the fluid mechanics of the light-bulb concept of the gaseous core nuclear engine. A numerical treatment is developed for the rotating flow which includes a description of the nuclear fuel addition. The problem is formulated with the complete Navier-Stokes equations in order to show the interaction between the fuel addition, the main flow, and the boundary layer flow in an accurate manner. The results presented show holdup of the nuclear fuel for the case of steady fuel addition.

Shin, Y. W.↗
NASA NTRS1982

A simple finite difference procedure for the vortex controlled diffuser

A simple prediction procedure for sudden expansion incompressible flows is developed and applied to the vortex controlled diffuser. Transient Navier-Stokes equations of an incompressible fluid are solved by means of their associated finite difference equations in terms of the primitive pressure velocity variables. A computer code is developed using a laminar flow simulation with free slip or no slip wall boundary conditions. In addition, predicted results confirm that effectiveness increases with increases in duct length and bleed flow rate

Busnaina, A. A.↗
NASA NTRS2004

Development of a Flexible Framework for Hypersonic Navier-Stoke Space Shuttle Orbiter Meshes

A flexible framework constructing block structured volume grids for hypersonic Navier-Strokes flow simulations was developed for the analysis of the Shuttle Orbiter Columbia. The development of the framework, which was partially basedon the requirements of the primary flow solvers used resulted in an ability to directly correlate solutions contributed by participating groups on a common surface mesh. A foundation was built through the assessment of differences between differnt solvers, which provided confidence for independent assessment of other damage scenarios by team members. The framework draws on the experience of NASA Langley and NASA Ames Research Centers in structured grid generation, and consists of a grid generation, and consist of a grid generation process implemented through a division of responsibilities. The nominal division of labor consisted of NASA Johnson Space Center coordinating the damage scenarios to be analyzed by the Aerothermodynamics Columbia Accident Investigation (ACAI) team, Ames developing the surface grids that described the computational volume about the Orbiter, and Langley improving grid quality of Ames generated data and constructing the final computational volume grids. Distributing the work among the participant in th ACAI team resulted in significantl less time required to construct complete meshes than possible by any individual participant. The approach demonstrated that the One-NASA grid generation team could sustain the demand of for five new meshes to explore new damage scenarios within an aggressive time-line.

Alter, Stephen J.↗
DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2002

The Fractional Step Method Applied to Simulations of Natural Convective Flows

This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The development of the momentum and continuity equations for natural convection in a fluid, a permeable medium, and in a binary alloy undergoing directional solidification will be presented. Finally, results for natural convection in a pure liquid, natural convection in a medium with a constant permeability, and for directional solidification will be presented.

Westra, Douglas G.↗