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 271 records · Page 15

Natural Rolling Responses of a Delta Wing in Transonic and Subsonic Flows

The unsteady, three-dimensional, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the natural rolling response of slender delta wings of zero thickness at moderate to high angles of attack, to transonic and subsonic flows. The governing equations of fluid flow and dynamics of the present multi-disciplinary problem are solved using the time-accurate solution of the NS equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The main focus is to analyze the effect of Mach number and angle of attack on the leading edge vortices and their breakdown, the resultant rolling motion, and overall aerodynamic response of the wing. Three cases demonstrate the natural response of a 65 deg swept, cropped delta wing in a transonic flow with breakdown of the leading edge vortices and an 80 deg swept delta wing in a subsonic flow undergoing either damped or self-excited limit-cycle rolling oscillations as a function of angle of attack. Comparisons with an experimental investigation completes this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret A.↗

Turbulent fluid motion 2: Scalars, vectors, and tensors

The author shows that the sum or difference of two vectors is a vector. Similarly the sum of any two tensors of the same order is a tensor of that order. No meaning is attached to the sum of tensors of different orders, say u(sub i) + u(sub ij); that is not a tensor. In general, an equation containing tensors has meaning only if all the terms in the equation are tensors of the same order, and if the same unrepeated subscripts appear in all the terms. These facts will be used in obtaining appropriate equations for fluid turbulence. With the foregoing background, the derivation of appropriate continuum equations for turbulence should be straightforward.

Deissler, Robert G.↗

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗

An equation of state for oxygen and nitrogen

Preliminary equations of state are presented for oxygen and nitrogen which provide accurate representations of the available P-density-T data for both fluids. The equation for nitrogen is applicable for temperatures from 70 K to 1300 K at pressures to 10,000 atmospheres, and the equation for oxygen for temperatures from 70 K to 323 K at pressures to 350 atmospheres. Deviations of calculated densities from representative experimental data are included. A volume-explicit equation of state for oxygen to be used in estimating density values in the range of applicability of the equation of state is also presented.

Stewart, R. B.↗

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unsteady, Transonic Flow Around Delta Wings Undergoing Coupled and Natural Modes Response: A Multidisciplinary Problem

The unsteady, three-dimensional Navier-Stokes equations coupled with the Euler equations of rigid-body dynamics are sequentially solved to simulate and analyze the aerodynamic response of a high angle of attack delta wing undergoing oscillatory motion. The governing equations of fluid flow and dynamics of the multidisciplinary problem are solved using a time-accurate solution of the laminar, unsteady, compressible, full Navier- Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The primary model under consideration consists of a 65 deg swept, sharp-edged, cropped delta wing of zero thickness at 20 deg angle of attack. In a freestream of Mach 0.85 and Reynolds number of 3.23 x 10(exp 6), the flow over the upper surface of the wing develops a complex shock system which interacts with the leading-edge primary vortices producing vortex breakdown. The effect of the oscillatory motion of the wing on the vortex breakdown and overall aerodynamic response is detailed to provide insight to the complicated physics associated with unsteady flows and the phenomenon of wing rock. Forced sinusoidal single and coupled mode rolling and pitching motion is presented for the wing in a transonic freestream. The Reynolds number, frequency of oscillation, and the phase angle are varied. Comparison between the single and coupled mode forced rolling and pitching oscillation cases illustrate the effects of coupling the motion. This investigation shows that even when coupled, forced rolling oscillation at a reduced frequency of 2(pi) eliminates the vortex breakdown which results in an increase in lift. The coupling effect for in phase forced oscillations show that the lift coefficient of the pitching-alone case and the rolling-moment coefficient of the rolling-alone case dominate the resulting response. However, with a phase lead in the pitching motion, the coupled motion results in a non-periodic response of the rolling moment. The second class of problems involve releasing the wing in roll to respond to the flowfield. Two models of sharp-edged delta wings, the previous 65 deg swept model and an 80 deg swept, sharp-edged delta wing, are used to observe the aerodynamic response of a wing free to roll in a transonic and subsonic freestream, respectively. These cases demonstrate damped oscillations, self-sustained limit cycle oscillations, and divergent rolling oscillations. Ultimately, an active control model using a mass injection system was applied on the surface of the wing to suppress the self-sustained limit cycle oscillation known as wing rock. Comparisons with experimental investigations complete this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret Anne↗

New Approaches to Inferences for Steep-Sided Domes on Venus

New mathematical approaches for the relaxation and emplacement of viscous lava domes are presented and applied to steep-sided domes on Venus. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry for two distinct scenarios. In the first scenario, dome relaxation is explored assuming a constant volume of fluid (i.e. lava) has been rapidly emplaced onto the surface. Cooling of lava is represented by a time-variable viscosity and singularities inherent in previous models for dome relaxation have been eliminated. At the onset of relaxation, bulk dynamic viscosities lie in the range between 10(exp 10) - 10(exp 16) Pa s, consistent with basaltic-andesite to rhyolitic compositions. Plausible relaxation times range from 5 to 5000 years, depending on initial lava viscosity. The first scenario, however, is only valid during the final stages of dome relaxation and does not consider the time taken for lava to be extruded onto the surface. In the second scenario, emplacement and growth of a steep-sided dome is considered when the volume of lava on the surface increases over time (i.e. time-variable volume approach). The volumetric flow rate may depend on an arbitrary power of the dome thickness, thus embracing Newtonian as well as other rheologies for describing terrestrial and planetary mass flows. The approach can be used to distinguish between basic flow rate models for fluid emplacement. The formalism results in radial expansion of a dome proportional to t(sup 1/2), consistent with the diffusive nature of the governing equation. The flow at the front is shown to thicken as the front advances for a constant rate of lava supply. Emplacement times are intimately correlated with the bulk rheology. Comparison of the theoretical profiles with the shape of a typical dome on Venus indicates that a Newtonian bulk rheology is most appropriate, consistent with prior studies. However, results here suggest a bulk dynamic viscosity of 10(exp 12) - 10(exp 13) Pa s and emplacement times of approximately 2-16 years. Both scenarios investigated give emplacement times significantly less than prior estimates and compositions consistent with basaltic andesite.

steep-sided domes↗

Prediction of hydrodynamics and chemistry of confined turbulent methane-air frames in a two concentric tube combustor

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a two-concentric-tube combustor is presented. A numerical procedure for the solution of the governing differential equations is described and models for chemical-equilibrium and chemical-kinetics calculations are presented. The chemical-equilibrium model is used to characterize the hydrocarbon reactions. The chemical-kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor considered consists of two coaxial ducts. Concentric streams of gaseous fuel and air enter the inlet duct at one end; the flow then reverses and flows out through the outer duct. Two sample cases with specified inlet and boundary conditions are considered and the results are discussed.

Markatos, N. C.↗

Normal modes of a compound drop

The modes are characterized by their frequency, the attendant displacement of fluid boundaries, and the flow pressure fields within the fluids. The drops consist of three fluids; a core fluid, a fluid shell surrounding the core, and a host fluid surrounding the shell. These fluids are assumed to be inviscid and incompressible, and the core and the shell to be concentric. The theory is obtained by linearization of the equations of fluid motion to the lowest order of nonlinearity that yields the normal modes. Numerical values of mode frequencies and the associated relative displacements of the fluid boundaries are presented for several specific systems, and the results compared with observations.

Saffren, M.↗

A pressure correction method for the calculation of compressible chemical reacting flows

A recently developed noniterative method for the solution of the transient fluid flow equations at all speed is extended to handle chemical reacting flows. The species conservation equations are loosely coupled into the predictor/multicorrector sequence of the solution procedure. A split-operator method separates the chemical kinetics terms from the fluid-dynamical terms, as well as an implicit differencing method enhance the numerical stability. The method was applied for turbulent diffusion flame calculations and for the analyses of high pressure, axisymmetric turbulent hypersonic nozzle flows. The diffusion flame results were compared with a similar pressure method for fast chemistry integration scheme without operator-splitting. Simulations of the nozzle flow indicated that the nonideal intermolecular effects must be included in the analysis and design of high pressure hypersonic nozzle.

Chen, Z. J.↗

Three-dimensional unstructured method for flows past bodies in 6-DOF relative motion

A three dimensional, unstructured-mesh methodology was developed to simulate unsteady flows past bodies in relative motion, where the trajectory was determined from the instantaneous aerodynamics. The method coupled the equations of fluid flow and those of rigid-body dynamics, and captured the time-dependent interference between stationary and moving boundaries. The unsteady, compressible Euler equations were solved on dynamic, unstructured meshes by an explicit, finite-volume, upwind method. The grid adaptation was performed within a window placed around the moving body. The Euler equations of dynamics were solved by a Runge-Kutta integration scheme. The flow solver and the adaptation scheme were validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic pitching motion, then comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom (DOF) motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, K. P.↗

Development of an integrated BEM (Boundary Element Mesh) for hot fluid-structure interaction

A boundary integral representation for a coupled approach to fluid flow and solid deformation problems associated with the design of hot-section components such as those in the Space Shuttle Main Engine is discussed. The formulation is based on the fundamental analytical solution of the Navier-Stokes equation for fluid velocity in an infinite domain. This fundamental solution was obtained by decomposing a Navier-Stokes equation into vorticity and dilation transport equations. A boundary integral involving convolutions in time was then constructed in which the convective terms appear in the volume integral.

Dargush, G. F.↗

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗

Molecular hydrodynamic theory of the velocity autocorrelation function

The velocity autocorrelation function (VACF) encapsulates extensive information about a fluid’s molecular-structural and hydrodynamic properties. We address the following fundamental question: How well can a purely hydrodynamic description recover the molecular features of a fluid as exhibited by the VACF? To this end, we formulate a bona fide hydrodynamic theory of the tagged-particle VACF for simple fluids. Our approach is distinguished from previous efforts in two key ways: collective hydrodynamic modes and tagged-particle self-motion are modeled by linear hydrodynamic equations; the fluid’s spatial velocity power spectrum is identified as a necessary initial condition for the momentum current correlation. This formulation leads to a natural physical interpretation of the VACF as a superposition of products of quasinormal hydrodynamic modes weighted commensurately with the spatial velocity power spectrum, the latter of which appears to physically bridge continuum hydrodynamical behavior and discrete-particle kinetics. The methodology yields VACF calculations quantitatively on par with existing approaches for liquid noble gases and alkali metals. Furthermore, we obtain a new, hydrodynamic form of the self-intermediate scattering function whose description has been extended to low densities where the Schmidt number is of order unity; various calculations are performed for gaseous and supercritical argon to support the general validity of the theory. Excellent quantitative agreement is obtained with recent MD calculations for a dense supercritical Lennard-Jones fluid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

Circularly polarized inertial wave vectors in rotating fluids

The Navier-Stokes equations for a rotating fluid are harmonically analyzed for planar motion in an infinite half-space. All solutions are shown to be a sum of two inertial wave vectors, one circularly polarized to the left (CPL) and the other circularly polarized to the right (CPR). These basic solutions are therefore presented in the same nomenclature and form as that found useful by experimentalists in analyzing flow data (called 'rotary spectra'). The CPL wave acts counter to the Coriolis force and consequently has a slower phase speed and larger damping than the CPR wave. At resonance (forcing frequency = Coriolis frequency) the CPR wave has an infinite phase speed and no damping and is the important component leading to the singular nature of the solutions for certain boundary conditions. All possible resonant singularities are explicitly shown. The unsteady development of these unbounded (limited space structure), cyclic (no time origin or structure) flows is presented to show that with time structure the resonant singularities evolve in a self-similar manner.

Cogley, A. C.↗

Progress in unstructured-grid methods development for unsteady aerodynamic applications

The development of unstructured-grid methods for the solution of the equations of fluid flow and what was learned over the course of the research are summarized. The focus of the discussion is on the solution of the time-dependent Euler equations including spatial discretizations, temporal discretizations, and boundary conditions. An example calculation with an implicit upwind method using a CFL number of infinity is presented for the Boeing 747 aircraft. The results were obtained in less than one hour CPU time on a Cray-2 computer, thus, demonstrating the speed and robustness of the capability. Additional calculations for the ONERA M6 wing demonstrate the accuracy of the method through the good agreement between calculated results and experimental data for a standard transonic flow case.

Batina, John T.↗