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 973 records · Page 54

Starvation effects on the hydrodynamic lubrication of rigid nonconformal contacts in combined rolling and normal motion

The effect of inlet starvation on the hydrodynamic lubrication of lightly loaded rigid nonconformal contacts in combined rolling and normal motion is determined through a numerical solution of the Reynolds' equation for an isoviscous, incompressible lubricant. Starvation is effected by systematically reducing the fluid inlet level. The pressures are taken to be ambient at the inlet meniscus boundary and Reynolds' boundary condition is applied for film rupture in the exit region. Results are presented for the dynamic performance of the starved contacts in combined rolling and normal motion for both normal approach and separation. During normal approach the dynamic load ratio (i.e. ratio of dynamic to steady state load capacity) increases considerably with increase in the inlet starvation. The effect of starvation on the dynamic peak pressure ratio is relatively small. Further, it has been observed that with increasing starvation, film thickness effects become significant in the dynamic behavior of the nonconformal contacts. For significantly starved contacts the dynamic load ratio increases with increase in film thickness during normal approach and a similar reduction is observed during separation. A similar effect is noted for the dynamic peak pressure ratio.

Ghosh, M. K.↗

Starvation effects on the hydrodynamic lubrication of rigid nonconformal contacts in combined rolling and normal motion

The effect of inlet starvation on the hydrodynamic lubrication of lightly loaded rigid nonconformal contacts in combined rolling and normal motion is determined through a numerical solution of the Reynolds' equation for an isoviscous, incompressible lubricant. Starvation is effected by systematically reducing the fluid inlet level. The pressures are taken to be ambient at the inlet meniscus boundary and Reynolds' boundary condition is applied for film rupture in the exit region. Results are presented for the dynamic performance of the starved contacts in combined rolling and normal motion for both normal approach and separation. During normal approach the dynamic load ratio (i.e. ratio of dynamic to steady state load capacity) increases considerably with increase in the inlet starvation. The effect of starvation on the dynamic peak pressure ratio is relatively small. Further, it has been observed that with increasing starvation, film thickness effects become significant in the dynamic behavior of the nonconformal contacts. For significantly starved contacts the dynamic load ratio increases with increase in film thickness during normal approach and a similar reduction is observed during separation. A similar effect is noted for the dynamic peak pressure ratio.

Ghosh, M. K.↗

Numerical study of chemically reacting flows using an LU scheme

A new computational fluid dynamic code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite volume, lower-upper symmetric successive overrelaxation scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully-coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step finite rate chemistry model whereas turbulence is simulated by a Baldwin-Lomax algebraic model. The validity of the code is demonstrated by comparing the numerical calculations with both experimental data and previous calculations of a cold flow helium injection into a straight channel and premixed hydrogen-air reacting flows in a ramped duct. The code is then used to calculate the mixing and chemical reactions of a hydrogen jet transversely injected into a supersonic airstream. Results are presented describing the flow field, the recirculation regions in front and behind the injector, and the chemical reactions.

Shuen, Jian Shun↗

Numerical study of chemically reacting flows using an LU scheme

A new computational fluid dynamic code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite volume, lower-upper symmetric successive overrelaxation scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully-coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step finite rate chemistry model whereas turbulence is simulated by a Baldwin-Lomax algebraic model. The validity of the code is demonstrated by comparing the numerical calculations with both experimental data and previous calculations of a cold flow helium injection into a straight channel and premixed hydrogen-air reacting flows in a ramped duct. The code is then used to calculate the mixing and chemical reactions of a hydrogen jet transversely injected into a supersonic airstream. Results are presented describing the flow field, the recirculation regions in front and behind the injector, and the chemical reactions.

Shuen, Jian Shun↗

Progress in the development of PDF turbulence models for combustion

A combined Monte Carlo-computational fluid dynamic (CFD) algorithm was developed recently at Lewis Research Center (LeRC) for turbulent reacting flows. In this algorithm, conventional CFD schemes are employed to obtain the velocity field and other velocity related turbulent quantities, and a Monte Carlo scheme is used to solve the evolution equation for the probability density function (pdf) of species mass fraction and temperature. In combustion computations, the predictions of chemical reaction rates (the source terms in the species conservation equation) are poor if conventional turbulence modles are used. The main difficulty lies in the fact that the reaction rate is highly nonlinear, and the use of averaged temperature produces excessively large errors. Moment closure models for the source terms have attained only limited success. The probability density function (pdf) method seems to be the only alternative at the present time that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus may be the only viable approach for more accurate turbulent combustion calculations. Assumed pdf's are useful in simple problems; however, for more general combustion problems, the solution of an evolution equation for the pdf is necessary.

Hsu, Andrew T.↗

Boost-phase discrimination research

The final report describes the combined work of the Computational Chemistry and Aerothermodynamics branches within the Thermosciences Division at NASA Ames Research Center directed at understanding the signatures of shock-heated air. Considerable progress was made in determining accurate transition probabilities for the important band systems of NO that account for much of the emission in the ultraviolet region. Research carried out under this project showed that in order to reproduce the observed radiation from the bow shock region of missiles in their boost phase it is necessary to include the Burnett terms in the constituent equation, account for the non-Boltzmann energy distribution, correctly model the NO formation and rotational excitation process, and use accurate transition probabilities for the NO band systems. This work resulted in significant improvements in the computer code NEQAIR that models both the radiation and fluid dynamics in the shock region.

Langhoff, Stephen R.↗

Lox droplet vaporization in a supercritical forced convective environment

A systematic investigation has been conducted to study the effects of ambient flow conditions (i.e. pressure and velocity) on supercritical droplet gasification in a forced-convective environment. The model is based on the time-dependent conservation equations in axisymmetric coordinates, and accommodates thermodynamic nonidealities and transport anomalies. In addition, an efficient scheme for evaluating thermophysical properties over the entire range of fluid thermodynamic states is established. The analysis allows a thorough examination of droplet behavior during its entire lifetime, including transient gasification, dynamic deformation, and shattering. A parametric study of droplet vaporization rate in terms of ambient pressure and Reynolds number is also conducted.

Hsiao, Chia-Chun↗

Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. These algorithms are now used in the practical engineering design of devices such as cars and airplanes as well as theoretical studies of complex phenomena such as fluid turbulence. In past years, limited computer resources have motivated the development of efficient numerical methods in computational fluid dynamics (CFD) utilizing structured meshes. These meshes are comprised of systematic arrays of quadrilateral or hexahedral cells. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Structured meshes also permit the use of highly efficient solution techniques such as alternating direction implicit (ADI) iteration schemes or multigrid. Following the dramatic improvement in computing speed in recent years, emphasis has shifted towards the design of algorithms capable of treating complex geometries. The automatic generation of structured grids about complex geometries is problematic. Unstructured grids offer one promising alternative technique for treating these general geometries. Unstructured meshes have irregular connectivity and usually contain triangles and/or quadrilaterals in two dimensions and tetrahedra and/or hexahedra in three dimensions. The generation and use of unstructured grids poses new challenges in computational fluid dynamics. This is true for both grid generation as well as for the design of algorithms for flow solution. The purpose of these notes is to present recent developments in the unstructured grid generation and flow solution technology.

Barth, T. J.↗

Numerical Calculation of the Drag Force Acting on a Solid Particle Pushed by a Solid/Liquid Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: Abstract The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: F(sub D, sup num) =sqoare root of 3(pi)(eta)V(sub p)(R(sub p)/d)(sup 10(gamma)/3). This is to compared to the solution provided by the lubrication theory: F(sub D), sup theor) = 3 square root of 2(pi)(eta)V(sub p)(R(sub p)/d)(sub 3/2), where F(sub D) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma approx. = 0.577), and d is the width of the gap between the particle and the SLI. The domain on which the equations proposed by the lubrication theory are relevant is clearly identified in the paper. The numerical model was then validated against the classical lubrication theory within the domain of its validity.

Catalina, Adrian V.↗

Factorizable Upwind Schemes: the Triangular Unstructured Grid Formulation

The upwind factorizable schemes for the equations of fluid was introduced recently. They facilitate achieving the Textbook Multigrid Efficiency (TME) and are expected also to result in the solvers of unparalleled robustness. The approach itself is very general. Therefore, it may well become a general framework for the large-scale Computational Fluid Dynamics. In this paper we outline the triangular grid formulation of the factorizable schemes. The derivation is based on the fact that the factorizable schemes can be expressed entirely using vector notation, without explicitly mentioning a particular coordinate frame. We describe the resulting discrete scheme in detail and present some computational results verifying the basic properties of the scheme/solver.

Nielsen, Eric J.↗

Factorizable Upwind Schemes: The Triangular Unstructured Grid Formulation

The upwind factorizable schemes for the equations of fluid were introduced recently. They facilitate achieving the Textbook Multigrid Efficiency (TME) and are expected also to result in the solvers of unparalleled robustness. The approach itself is very general. Therefore, it may well become a general framework for the large-scale, Computational Fluid Dynamics. In this paper we outline the triangular grid formulation of the factorizable schemes. The derivation is based on the fact that the factorizable schemes can be expressed entirely using vector notation. without explicitly mentioning a particular coordinate frame. We, describe the resulting discrete scheme in detail and present some computational results verifying the basic properties of the scheme/solver.

Sidilkover, David↗

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) is a Space Shuttle Orbiter support piece of the Ground Servicing Equipment (GSE). The purpose of the UCTS is to provide two support services to the Space Shuttle Orbiter immediately after landing at the Shuttle Landing Facility. The Simulation uses GSE Models to stand in for the actual systems to support testing of SCCS systems s:luring their development. As an intern at KSC, my assignment was to develop a model component for the UCTS. I was given a fluid component (drier) to model in Matlab. The drier was a Catch All replaceable core type filter-drier. The filter-drier provides maximum protection for the thermostatic expansion valve and solenoid valve from dirt that may be in the system. The filter-drier also protects the valves from freezing up. I researched fluid dynamics to understand the function of my component. I completed training for UNIX and Simulink to help aid in my assignment. The filter-drier was modeled by determining affects it has on the pressure, velocity and temperature of the system. I used Bernoulli's Equation to calculate the pressure and velocity differential through the dryer. I created my model filter-drier 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.

Fields, Christina M.↗

Reference Solutions for Benchmark Turbulent Flows in Three Dimensions

A grid convergence study is performed to establish benchmark solutions for turbulent flows in three dimensions (3D) in support of turbulence-model verification campaign at the Turbulence Modeling Resource (TMR) website. The three benchmark cases are subsonic flows around a 3D bump and a hemisphere-cylinder configuration and a supersonic internal flow through a square duct. Reference solutions are computed for Reynolds Averaged Navier Stokes equations with the Spalart-Allmaras turbulence model using a linear eddy-viscosity model for the external flows and a nonlinear eddy-viscosity model based on a quadratic constitutive relation for the internal flow. The study involves three widely-used practical computational fluid dynamics codes developed and supported at NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions computed with these three codes on families of consistently refined grids are presented. Grid-to-grid and code-to-code variations are described in detail.

Diskin, Boris↗

Propellant Sloshing Parameter Extraction from CFD Analysis

Propellant slosh is a potential source of disturbance critical to the stability of space vehicle. The sloshing dynamics is typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters required by the mechanical model include natural frequency of the sloshing, sloshing mass, sloshing mass center coordinates, and critical damping coefficient. During the 1960 s US space program, these parameters were either computed from analytical solution for simple geometry or by experimental testing for the sub-scaled configurations. The purpose of this work is to demonstrate the soundness of a CFD approach in modeling the detailed fluid dynamics of tank sloshing and the excellent accuracy in extracting mechanical properties for different tank configurations and at different fill levels. The validation studies included straight cylinder against analytical solution, and sub-scaled Centaur LOX and LH2 tanks with and without baffles against experimental results. This effort shows that CFD technology can provide accurate mechanical parameters for any tank configuration, and is especially valuable to the future design of propellant tanks, as there is no previous experimental data available for the same size and configuration.

Yang, H. Q.↗

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗

Progress and supercomputing in computational fluid dynamics; Proceedings of U.S.-Israel Workshop, Jerusalem, Israel, December 1984

Current developments and future trends in the application of supercomputers to computational fluid dynamics are discussed in reviews and reports. Topics examined include algorithm development for personal-size supercomputers, a multiblock three-dimensional Euler code for out-of-core and multiprocessor calculations, simulation of compressible inviscid and viscous flow, high-resolution solutions of the Euler equations for vortex flows, algorithms for the Navier-Stokes equations, and viscous-flow simulation by FEM and related techniques. Consideration is given to marching iterative methods for the parabolized and thin-layer Navier-Stokes equations, multigrid solutions to quasi-elliptic schemes, secondary instability of free shear flows, simulation of turbulent flow, and problems connected with weather prediction.

Murman, E. M.↗

The 3DGRAPE book: Theory, users' manual, examples

A users' manual for a new three-dimensional grid generator called 3DGRAPE is presented. The program, written in FORTRAN, is capable of making zonal (blocked) computational grids in or about almost any shape. Grids are generated by the solution of Poisson's differential equations in three dimensions. The program automatically finds its own values for inhomogeneous terms which give near-orthogonality and controlled grid cell height at boundaries. Grids generated by 3DGRAPE have been applied to both viscous and inviscid aerodynamic problems, and to problems in other fluid-dynamic areas. The smoothness for which elliptic methods are known is seen here, including smoothness across zonal boundaries. An introduction giving the history, motivation, capabilities, and philosophy of 3DGRAPE is presented first. Then follows a chapter on the program itself. The input is then described in detail. A chapter on reading the output and debugging follows. Three examples are then described, including sample input data and plots of output. Last is a chapter on the theoretical development of the method.

Sorenson, Reese L.↗

Euler and Potential Experiment/CFD Correlations for a Transport and Two Delta-Wing Configurations

A selection of successes and failures of Computational Fluid Dynamics (CFD) is discussed. Experiment/CFD correlations involving full potential and Euler computations of the aerodynamic characteristics of four commercial transport wings and two low aspect ratio, delta wing configurations are shown. The examples consist of experiment/CFD comparisons for aerodynamic forces, moments, and pressures. Navier-Stokes equations are not considered.

Hicks, R. M.↗