Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equation of 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 811 records · Page 45

Design of an Object-Oriented Turbomachinery Analysis Code: Initial Results

Performance prediction of turbomachines is a significant part of aircraft propulsion design. In the conceptual design stage, there is an important need to quantify compressor and turbine aerodynamic performance and develop initial geometry parameters at the 2-D level prior to more extensive Computational Fluid Dynamics (CFD) analyses. The Object-oriented Turbomachinery Analysis Code (OTAC) is being developed to perform 2-D meridional flowthrough analysis of turbomachines using an implicit formulation of the governing equations to solve for the conditions at the exit of each blade row. OTAC is designed to perform meanline or streamline calculations; for streamline analyses simple radial equilibrium is used as a governing equation to solve for spanwise property variations. While the goal for OTAC is to allow simulation of physical effects and architectural features unavailable in other existing codes, it must first prove capable of performing calculations for conventional turbomachines.OTAC is being developed using the interpreted language features available in the Numerical Propulsion System Simulation (NPSS) code described by Claus et al (1991). Using the NPSS framework came with several distinct advantages, including access to the pre-existing NPSS thermodynamic property packages and the NPSS Newton-Raphson solver. The remaining objects necessary for OTAC were written in the NPSS framework interpreted language. These new objects form the core of OTAC and are the BladeRow, BladeSegment, TransitionSection, Expander, Reducer, and OTACstart Elements. The BladeRow and BladeSegment consumed the initial bulk of the development effort and required determining the equations applicable to flow through turbomachinery blade rows given specific assumptions about the nature of that flow. Once these objects were completed, OTAC was tested and found to agree with existing solutions from other codes; these tests included various meanline and streamline comparisons of axial compressors and turbines at design and off-design conditions.

compressors↗

Design of an Object-Oriented Turbomachinery Analysis Code: Initial Results

Performance prediction of turbomachines is a significant part of aircraft propulsion design. In the conceptual design stage, there is an important need to quantify compressor and turbine aerodynamic performance and develop initial geometry parameters at the 2-D level prior to more extensive Computational Fluid Dynamics (CFD) analyses. The Object-oriented Turbomachinery Analysis Code (OTAC) is being developed to perform 2-D meridional flowthrough analysis of turbomachines using an implicit formulation of the governing equations to solve for the conditions at the exit of each blade row. OTAC is designed to perform meanline or streamline calculations; for streamline analyses simple radial equilibrium is used as a governing equation to solve for spanwise property variations. While the goal for OTAC is to allow simulation of physical effects and architectural features unavailable in other existing codes, it must first prove capable of performing calculations for conventional turbomachines. OTAC is being developed using the interpreted language features available in the Numerical Propulsion System Simulation (NPSS) code described by Claus et al (1991). Using the NPSS framework came with several distinct advantages, including access to the pre-existing NPSS thermodynamic property packages and the NPSS Newton-Raphson solver. The remaining objects necessary for OTAC were written in the NPSS framework interpreted language. These new objects form the core of OTAC and are the BladeRow, BladeSegment, TransitionSection, Expander, Reducer, and OTACstart Elements. The BladeRow and BladeSegment consumed the initial bulk of the development effort and required determining the equations applicable to flow through turbomachinery blade rows given specific assumptions about the nature of that flow. Once these objects were completed, OTAC was tested and found to agree with existing solutions from other codes; these tests included various meanline and streamline comparisons of axial compressors and turbines at design and off-design conditions.

object-oriented programming↗

Grid Convergence for Three Dimensional Benchmark Turbulent Flows

Grid convergence studies are performed to establish reference solutions for benchmark three dimensional turbulent flows in support of the ongoing turbulence model verification and validation e ort at the Turbulence Modeling Resource website curated by NASA. The bench- mark cases are a subsonic flow around a hemisphere cylinder and a transonic flow around the ONERA M6 wing with a sharp trailing edge. The study applies widely-used computational fluid dynamics codes developed and supported at the NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions are computed for the Reynolds-Averaged Navier-Stokes equations with the Spalart-Allmaras turbulence model on families of consistently-refined grids composed of different types of cells. Coarse-to- ne and code-to-code solution variation is described in detail.

Diskin, Boris↗

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗

Invariant two-dimensional structures of vorticity field in nonviscous fluid

It is shown that a fluid flow which are invariant under translations along their vorticity field directions can be described in terms of a formally 2D system. Its nonviscous behavior of the scalar vorticity is governed by Hamiltonian dynamics, similar to the 2D line vortex case. Applying a statistical mechanical method, the equilibrium distribution of the vorticity field for such a general 2D system is shown to obey a nonlinear partial differential equation that is a generalization of the sinh-Poisson equation derived by Montgomery and Joyce (1973, 1974). Thus, if the domain containing the fluid is finite, similar negative temperature states may occur and the distribution of the vorticity intensity may become localized in space. The results are believed to have applications in studying large-scale partial and temporal properties of flows having small but nonzero viscosity. The analysis can probably be extended to other fluid systems with curved magnetic field lines.

Chen, Hudong↗

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗

Preliminary experimental results for a cryogenic brush seal configuration

Preliminary fluid nitrogen flow data are reported for a five-brush, ceramic-coated-rub-runner brush seal system, where the brushes and the rub runner were placed at each end of a centrally pressurized multifunction tester ('back-to-back' set of brushes) and tested at rotor speeds of 0, 10, 18, and 22.5 krpm. After testing, both the brushes and the ceramic-coated rub runner appeared pristine. The coating withstood both the thermomechanical and dynamic loadings with minor wear track scarring. The bristle tips showed some indication of material shearing (smearing) wear. The Ergun (1952) porous flow equation was applied to the brush seal data. The Ergun relation fit the gaseous data but was in poor agreement with the fluid results. The brush seal exit conditions were two phase. Two-phase, choked-flow design charts were applied but required one data point at each rotor speed to define the flow and area coefficients. Reasonable agreement between prediction and data was found, but such methods are not to be construed as two-phase-flow brush seal analyses.

Carlile, J. A.↗

Modeling of InP metalorganic chemical vapor deposition

The growth of InP by metalorganic chemical vapor deposition (MOCVD) in a horizontal reactor is being modeled with a commercially available computational fluid dynamics modeling code. The mathematical treatment of the MOCVD process has four primary areas of concern: 1) transport phenomena, 2) chemistry, 3) boundary conditions, and 4) numerical solution methods. The transport processes involved in CVD are described by conservation of total mass, momentum, energy, and atomic species. Momentum conservation is described by a generalized form of the Navier-Stokes equation for a Newtonian fluid and laminar flow. The effect of Soret diffusion on the transport of particular chemical species and on the predicted deposition rate is examined. Both gas-phase and surface chemical reactions are employed in the model. Boundary conditions are specified at the inlet and walls of the reactor for temperature, fluid flow and chemical species. The coupled set of equations described above is solved by a finite difference method over a nonuniform rectilinear grid in both two and three dimensions. The results of the 2-D computational model is presented for gravity levels of zero- and one-g. The predicted growth rates at one-g are compared to measured growth rates on fused silica substrates.

Black, Linda R.↗

PAGOSA Theory Manual

PAGOSA is a computational fluid dynamics program developed at LANL for the study of high-speed compressible flow and high-rate material deformation. PAGOSA is a three-dimensional Eulerian finite-difference code, solving problems with a wide variety of EOS, material strength, and explosive modeling options. This document presents the finite difference equations that are used in the PAGOSA continuum mechanics computer code. This program is especially intended to be used for the numerical simulation of the interactions of gases, fluids, and solids.

36 MATERIALS SCIENCE↗

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model verification studies using accurate measurements of spin up

The reliability and accuracy of the numerical code for spin up flows in a cylinder by comparing the numerical results against high resolution laser Doppler velocimeter (LDV) measurements of the azimuthal flows were checked. A computer code to generate numerical solution for axisymmetric rotating fluid in a cylinder was obtained and amended for routine use at MSFC. The numerical simulations used the Navier-Stokes equations in axisymmetric form and employed finite difference techniques on both constant and variable grids. The numerical solutions are analyzed to gain further insight into the fundamental questions analyzed in rotating fluid dynamics.

Hyun, J. M.↗

A finite element solver for 3-D compressible viscous flows

Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.

Reddy, K. C.↗

Aerothermodynamic environment predictions in support of the Aeroassist Flight Experiment

A computational fluid dynamic (CFD) simulation has been made of the flow past an Aeroassist Flight Experiment (AFE) aerobrake model with a payload and sting. The simulation involves solving the complete Navier-Stokes equations in three dimensions. Comparisons with data taken in a Mach 10 wind tunnel test are made and agreement is shown to be very good. Points of comparison include the aerobrake surface pressure and heat transfer rate, the sting surface streamlines, and the sting heat transfer rate. The dependence of the solution on the grid is also explored. Finally, predictions are made for conditions which correspond to a hypersonic shock tunnel force and moment test to be carried out later this year.

Fay, John F.↗

CFD Simulations Of Sonic Booms In Near And Mid Fields

Report discusses computational fluid dynamics (CFD) to simulate generation and propagation of sonic booms in near- and mid-field regions of supersonic flows about simplified bodies representative of advanced airplanes. Parabolized Navier-Stokes equations integrated by implicit, approximate-factorization, finite-volume algorithm in which crossflow inviscid fluxes evaluated by Roe's flux-difference-splitting scheme. Near-field solutions obtained by applying algorithm to flows immediately surrounding bodies. Solutions transferred to computer codes based on Whitham"s F-function theory for extrapolation to far-field.

Cheung, Samson H.↗

Computational fluid dynamics application to hypersonic flow over a Martian entry vehicle (MEV) - A correlation with experiment

A class of vehicles for a mission to Mars are analyzed for aerodynamic characteristics using advanced Computational Fluid Dynamics (CFD). The general configuration is a modified cone-conical-frustum geometry where the nose radius has a large influence on the flowfield. Inviscid-compressible flow using the Euler equations and viscous-compressible flow using the thin-layer Navier-Stokes equations is applied to the configuration. The surface modeling, grid generation and application of state-of-the-art CFD software are described. The effects of nose radius, angle of attack, and hypersonic velocity on the flight characteristics of the vehicle are discussed. The numerical simulations demonstrate the merits of the inviscid and viscous software. Results are compared with wind tunnel experiments.

Stewart, John E.↗

Developing a Model Component

The Spaceport Command and Control System (SCCS) Simulation Computer Software Configuration Item (CSCI) is responsible for providing simulations to support test and verification of SCCS hardware and software. The Universal Coolant Transporter System (UCTS) was a Space Shuttle Orbiter support piece of the Ground Servicing Equipment (GSE). The initial purpose of the UCTS was to provide two support services to the Space Shuttle Orbiter immediately after landing at the Shuttle Landing Facility. The UCTS is designed with the capability of servicing future space vehicles; including all Space Station Requirements necessary for the MPLM Modules. The Simulation uses GSE Models to stand in for the actual systems to support testing of SCCS systems during their development. As an intern at Kennedy Space Center (KSC), my assignment was to develop a model component for the UCTS. I was given a fluid component (dryer) to model in Simulink. I completed training for UNIX and Simulink. The dryer is a Catch All replaceable core type filter-dryer. The filter-dryer provides maximum protection for the thermostatic expansion valve and solenoid valve from dirt that may be in the system. The filter-dryer also protects the valves from freezing up. I researched fluid dynamics to understand the function of my component. The filter-dryer was modeled by determining affects it has on the pressure and velocity of the system. I used Bernoulli's Equation to calculate the pressure and velocity differential through the dryer. I created my filter-dryer model in Simulink and wrote the test script to test the component. I completed component testing and captured test data. The finalized model was sent for peer review for any improvements. I participated in Simulation meetings and was involved in the subsystem design process and team collaborations. I gained valuable work experience and insight into a career path as an engineer.

Fields, Christina M.↗

Full Stall Simulations of a Redesigned Ventilation Fan for the ISS

The concept of a stall is studied rigorously in the aerospace industry. From a design standpoint, instabilities such as stall are undesirable– operation in the stall regime has a tremendous impact on aerodynamic performance as well as structural integrity. In extreme cases, operating in stall conditions can cause failure. Further, stall can cause a loss of lift on aircraft wings or a loss of thrust in aircraft engines. In any case, stall continues to be a topic of interest in the aerospace industry. The present work aims to analyze the stall characteristics of a ventilation fan that was recently designed for the International Space Station (ISS). Although the ventilation fan has a rotor-stator design, this paper considers a rotor-only configuration. The FUN3D Computational Fluid Dynamic (CFD) solver developed by NASA Langley Research Center was used to simulate the operational characteristics of the ventilation fan. FUN3D solves the Unsteady Reynolds-Averaged Naiver-Stokes (URANS) equations using implicit time marching and a dynamic overset grid. The FUN3D solver was originally written for exterior flow fields; however, this work represents an extension of the FUN3D solver to turbomachinery or interior flow fields. The FUN3D results for the rotor-only ventilation fan accurately captured the operational characteristics inherent to compressors– the results shared similar performance trends when compared to the experimental results for the rotor-stator case. A peak adiabatic efficiency of 96% occurred at a MFR of 105.2 CFM, the minimum aerodynamically stable point. A computationally stable stall occurred at a mass flow rate of 43.8 CFM where the adiabatic efficiency dropped to 69%.

CFD↗