Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluid equations”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Development of a linearized unsteady Euler analysis for turbomachinery blade rows

A linearized unsteady aerodynamic analysis for axial-flow turbomachinery blading is described in this report. The linearization is based on the Euler equations of fluid motion and is motivated by the need for an efficient aerodynamic analysis that can be used in predicting the aeroelastic and aeroacoustic responses of blade rows. The field equations and surface conditions required for inviscid, nonlinear and linearized, unsteady aerodynamic analyses of three-dimensional flow through a single, blade row operating within a cylindrical duct, are derived. An existing numerical algorithm for determining time-accurate solutions of the nonlinear unsteady flow problem is described, and a numerical model, based upon this nonlinear flow solver, is formulated for the first-harmonic linear unsteady problem. The linearized aerodynamic and numerical models have been implemented into a first-harmonic unsteady flow code, called LINFLUX. At present this code applies only to two-dimensional flows, but an extension to three-dimensions is planned as future work. The three-dimensional aerodynamic and numerical formulations are described in this report. Numerical results for two-dimensional unsteady cascade flows, excited by prescribed blade motions and prescribed aerodynamic disturbances at inlet and exit, are also provided to illustrate the present capabilities of the LINFLUX analysis.

Verdon, Joseph M.↗
DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗
NASA NTRS1998

Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft

This report compiles the various research activities conducted under the auspices of the NASA Grant NAG3-1026, "Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft" during the period of April 1989 to April 1994. The effort involved the development of multigrid based algorithms and computer programs for the calculation of the flow and temperature fields generated by Short Take-off and Vertical Landing (STOVL) aircraft, while hovering in ground proximity. Of particular importance has been the interaction of the exhaust jets with the head wind which gives rise to the hot gas ingestion process. The objective of new STOVL designs to reduce the temperature of the gases ingested into the engine. The present work describes a solution algorithm for the multi-dimensional elliptic partial-differential equations governing fluid flow and heat transfer in general curvilinear coordinates. The solution algorithm is based on the multigrid technique which obtains rapid convergence of the iterative numerical procedure for the discrete equations. Initial efforts were concerned with the solution of the Cartesian form of the equations. This algorithm was applied to a simulated STOVL configuration in rectangular coordinates. In the next phase of the work, a computer code for general curvilinear coordinates was constructed. This was applied to model STOVL geometries on curvilinear grids. The code was also validated in model problems. In all these efforts, the standard k-Epsilon model was used.

Vanka, S. P.↗
NASA NTRS1987

Comparison of generalized Reynolds and Navier Stokes equations for flow of a power law fluid

This paper compares a finite element solution of a modified Reynolds equation with a finite difference solution of the Navier-Stokes equation for a power law fluid. Both the finite element and finite difference formulation are reviewed. Solutions to spiral flow in parallel and conical geometries are compared. Comparison with experimental results are also given. The effects of the assumptions used in the Reynolds equation are discussed.

Mullen, R. L.↗
NASA NTRS1986

Generalization and extension of the law of acoustic energy conservation in a nonuniform flow

An exact conservation equation is derived which generalizes the familiar acoustic energy equations. The new relation is valid for arbitrary disturbances to a viscous, compressible flow. It is suggested by a development of the acoustic energy equation by means of a regular perturbation expansion of the general energy equation of fluid mechanics. A perturbation energy density and flux are defined and identified as the exact physical quantities whose leading order perturbation representations are the usual acoustic energy density and flux. The conservation equation governing the perturbation energy quantities is shown to yield previously known results for several special cases.

Myers, M. K.↗
NASA NTRS2008

An Initial Investigation of the Effects of Turbulence Models on the Convergence of the RK/Implicit Scheme

A three-stage Runge-Kutta (RK) scheme with multigrid and an implicit preconditioner has been shown to be an effective solver for the fluid dynamic equations. This scheme has been applied to both the compressible and essentially incompressible Reynolds-averaged Navier-Stokes (RANS) equations using the algebraic turbulence model of Baldwin and Lomax (BL). In this paper we focus on the convergence of the RK/implicit scheme when the effects of turbulence are represented by either the Spalart-Allmaras model or the Wilcox k-! model, which are frequently used models in practical fluid dynamic applications. Convergence behavior of the scheme with these turbulence models and the BL model are directly compared. For this initial investigation we solve the flow equations and the partial differential equations of the turbulence models indirectly coupled. With this approach we examine the convergence behavior of each system. Both point and line symmetric Gauss-Seidel are considered for approximating the inverse of the implicit operator of the flow solver. To solve the turbulence equations we use a diagonally dominant alternating direction implicit (DDADI) scheme. Computational results are presented for three airfoil flow cases and comparisons are made with experimental data. We demonstrate that the two-dimensional RANS equations and transport-type equations for turbulence modeling can be efficiently solved with an indirectly coupled algorithm that uses the RK/implicit scheme for the flow equations.

Swanson, R. C.↗
NASA NTRS1994

An introduction to generalized functions with some applications in aerodynamics and aeroacoustics

In this paper, we start with the definition of generalized functions as continuous linear functionals on the space of infinitely differentiable functions with compact support. The concept of generalization differentiation is introduced next. This is the most important concept in generalized function theory and the applications we present utilize mainly this concept. First, some of the results of classical analysis, such as Leibniz rule of differentiation under the integral sign and the divergence theorem, are derived using the generalized function theory. It is shown that the divergence theorem remains valid for discontinuous vector fields provided that the derivatives are all viewed as generalized derivatives. This implies that all conservation laws of fluid mechanics are valid as they stand for discontinuous fields with all derivatives treated as generalized deriatives. Once these derivatives are written as ordinary derivatives and jumps in the field parameters across discontinuities, the jump conditions can be easily found. For example, the unsteady shock jump conditions can be derived from mass and momentum conservation laws. By using a generalized function theory, this derivative becomes trivial. Other applications of the generalized function theory in aerodynamics discussed in this paper are derivation of general transport theorems for deriving governing equations of fluid mechanics, the interpretation of finite part of divergent integrals, derivation of Oswatiitsch integral equation of transonic flow, and analysis of velocity field discontinuities as sources of vorticity. Applications in aeroacoustics presented here include the derivation of the Kirchoff formula for moving surfaces,the noise from moving surfaces, and shock noise source strength based on the Ffowcs Williams-Hawkings equation.

Farassat, F.↗
NASA NTRS1983

Nonlinear aerodynamic modeling of flap oscillations in transonic flow - A numerical validation

The regime of validity of a nonlinear aerodynamic force and moment formulation, based on concepts from nonlinear functional analysis and applicable to a transonic airfoil with a deflecting flap, is investigated. A time-dependent finite difference technique is used to evaluate the aerodynamic data of the formulation in terms of specified, characteristic motions. Flap-motion histories are generated from the flap inertial equations of motion, with aerodynamic reactions specified by the moment formulation. The motion histories depicting the cases of decaying and growing flap oscillations are compared with histories generated through simultaneous, coupled solution of the fluid-dynamic equations and flap inertial equations of motion. The range of applicability of the formulation is discussed.

Chyu, W. J.↗
NASA NTRS1987

The structure of the time-dependent interstellar shocks and grain destruction in the interstellar medium

A new theoretical analysis of the structure of interstellar shocks and the grain dynamics of these shocks is presented. The basic hydrodynamic equations for J-shocks in interstellar gas are given in which shock-driving pressure is allowed to be weakly time-dependent and the grains are treated as a separate two-dimensional fluid. Specific equations for the grain dynamics in the cooling postshock gas are derived. An analytic theory for the propagation of a shock driven into an interstellar cloud by the blast wave of a supernova remnant which sweeps over the cloud is developed. An improved calculation of the grain charge in postshock gas is described, giving a simple analytic approximation for the results. The consequences of including these processes in the numerical code of Seab and Shull (1983) are addressed, including the effect of realistic magnetic fields in low-density gas.

Mckee, Christopher F.↗
DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗
NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗
NASA NTRS1992

A new method for transonic static aeroelasticity problems

A new method has been developed to calculate the steady flow and structural deformations for fluid/structure interaction problems. The discretized fluid dynamic and structural equations are regarded as a single set of coupled, nonlinear, algebraic equations. The equilibrium solution is directly obtained using Newton's method. The governing equations used for the fluid flow are the two-dimensional Navier-Stokes equations, and a finite-element model is used to represent the structure. This paper describes the analytical method and presents sample calculations demonstrating the technique. The results show rapid convergence and good agreement with experimental data.

Felker, Fort F.↗
NASA NTRS2010

Multi-Fluid Block-Adaptive-Tree Solar Wind Roe-Type Upwind Scheme: Magnetospheric Composition and Dynamics During Geomagnetic Storms, Initial Results

The magnetosphere contains a significant amount of ionospheric O{+}, particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the BATS-R-US model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multi-fluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single- fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multi-fluid MHD model (with outflow) gives comparable results to the multi-species MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multi-species and multi-fluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H{+} and O(+}, which is not possible when utilizing the other techniques considered.

Gkocer, A.↗
NASA NTRS2003

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗
NASA NTRS1987

Force and moment rotordynamic coefficients for pump-impeller shroud surfaces

Governing equations of motion are derived for a bulk-flow model of the leakage path between an impeller shroud and a pump housing. The governing equations consist of a path-momentum, a circumferential - momentum, and a continuity equation. The fluid annulus between the impeller shroud and pump housing is assumed to be circumferentially symmetric when the impeller is centered; i.e., the clearance can vary along the pump axis but does not vary in the circumferential direction. A perturbation expansion of the governing equations in the eccentricity ratio yields a set of zeroth and first-order governing equations. The zeroth-order equations define the leaking rate and the circumferential and path velocity distributions and pressure distributions for a centered impeller position. The first-order equations define the perturbations in the velocity and pressure distributions due to either a radial-displacement perturbation or a tilt perturbation of the impeller. Integration of the perturbed pressure and shear-stress distribution acting on the rotor yields the reaction forces and moments acting on the impeller face.

Childs, Dara W.↗
NASA NTRS1983

Role of thermal conduction in the acceleration of the solar wind

The thermal and suprathermal processes involving solar wind electrons are discussed from a theoretical point of view. A model for the electron distribution function, f(e), based on the solutions of the Boltzmann equation in Krook's approximation is outlined: the angular and energy dependences of for various distances from the Sun between the coronal base and the Earth are presented with the express purpose of obtaining the radial profile of the heat flux generated by the suprathermals. The basic fluid-dynamical equations for the solar wind plasma as a whole along any given magnetic field line are solved. Results are in good agreement with observations. The predicted density and temperature profiles for positive ions exhibit the characteristics of recent measurements, both in the corona (above the coronal holes) and between 0.3 and 1 A.U. (the Helios region). It is concluded that the suprathermal electrons play an important role in the formation and the dynamics of the solar wind.

Olbert, S.↗
DOE OSTI2022

How slowly can the early Universe expand?

When the expansion of the Universe is dominated by a perfect fluid with equation of state parameter w and a sound speed c s satisfying w = $c^{2}_{s}$ ≤ 1 , the Hubble parameter H and time t satisfy the bound Ht ≥ 1/3 . There has been recent interest in “ultraslow” expansion laws with Ht <1/3 (sometimes described as “fast expanding” models). Here, we examine various models that can produce ultraslow expansion: scalar fields with negative potentials, barotropic fluids, braneworld models, or a loitering phase in the early Universe. Scalar field models and barotropic models for ultraslow expansion are unstable to evolution toward w = 1 or w → ∞ in the former case and w → ∞ in the latter case. Braneworld models can yield ultraslow expansion but require an expansion law beyond the standard Friedman equation. Loitering early universe models can produce a quasistatic expansion phase in the early Universe but require an exotic negative-density component. These results suggest that appeals to an ultraslow expansion phase in the early Universe should be approached with some caution, although the loitering early universe may be worthy of further investigation. These results do not apply to ultraslow contracting models.

79 ASTRONOMY AND ASTROPHYSICS↗