Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compressible flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
NASA NTRS1975

Interaction of strong shock wave with turbulent boundary layer

The reported investigation represents an extension of the time-dependent solution of separated laminar flows based on the complete Navier-Stokes equations reported by MacCormack (1971) and Carter (1973). The current study includes turbulence models in conjunction with the compressible flow equations. The calculations start with a uniform flow except for values imposed along the upstream and outer boundaries. The basic numerical method is discussed along with questions concerning the exponential accuracy and the resolution of the viscous sublayer in a compressed region.

Baldwin, B. S.↗
NASA NTRS1943

The flow of a compressible fluid past a curved surface

An iteration method is employed to obtain the flow of a compressible fluid past a curved surface. The first approximation which leads to the Prandtl-Glauert rule, is based on the assumption that the flow differs but little from a pure translation. The iteration process then consists in improving this first approximation in order that it will apply to a flow differing from pure translatory motion to a greater degree. The method fails when the Mach number of the undisturbed stream reaches unity but permits a transition from subsonic to supersonic conditions without the appearance of a compression shock. The limiting value at which potential flow no longer exits is indicated by the apparent divergence of the power series representing the velocity of the fluid at the surface of the solid boundary.

Kaplan, Carl↗
NASA NTRS1974

Conical flow near singular rays

The steady flow of an ideal gas past a conical body is investigated by the method of matched asymptotic expansions, with particular emphasis on the flow near the singular ray occurring in linearized theory. The first-order problem governing the flow in this region is formulated, leading to the equation of Kuo, and an approximate solution is obtained in the case of compressive flow behind the main front. This solution is compared with the results of previous investigations with a view to assessing the applicability of the Lighthill-Whitham theories.

Zahalak, G. I.↗
NASA NTRS1992

Comments on the present state and future directions of PDF methods

The one point probability density function (PDF) method is examined in light of its use in actual engineering problems. The PDF method, although relatively complicated, appears to be the only format available to handle the nonlinear stochastic difficulties caused by typical reaction kinetics. Turbulence modeling, if it is to play a central role in combustion modeling, has to be integrated with the chemistry in a way which produces accurate numerical solutions to combustion problems. It is questionable whether the development of turbulent models in isolation from the peculiar statistics of reactant concentrations is a fruitful line of development as far as propulsion is concerned. There are three issues for which additional viewgraphs are prepared: the one point pdf method; the amplitude mapping closure; and a hybrid strategy for replacing a full two point pdf treatment of reacting flows by a single point pdf and correlation functions. An appeal is made for the establishment of an adequate data base for compressible flow with reactions for Mach numbers of unity or higher.

Obrien, E. E.↗
NASA NTRS1993

Adaptive-mesh algorithms for computational fluid dynamics

The basic goal of adaptive-mesh algorithms is to distribute computational resources wisely by increasing the resolution of 'important' regions of the flow and decreasing the resolution of regions that are less important. While this goal is one that is worthwhile, implementing schemes that have this degree of sophistication remains more of an art than a science. In this paper, the basic pieces of adaptive-mesh algorithms are described and some of the possible ways to implement them are discussed and compared. These basic pieces are the data structure to be used, the generation of an initial mesh, the criterion to be used to adapt the mesh to the solution, and the flow-solver algorithm on the resulting mesh. Each of these is discussed, with particular emphasis on methods suitable for the computation of compressible flows.

Powell, Kenneth G.↗
NASA NTRS1990

Embedded function methods for supersonic turbulent boundary layers

The development of embedded functions to represent the mean velocity and total enthalpy distributions in the wall layer of a supersonic turbulent boundary layer is considered. The asymptotic scaling laws (in the limit of large Reynolds number) for high speed compressible flows are obtained to facilitate eventual implementation of the embedded functions in a general prediction method. A self-consistent asymptotic structure is derived, as well as a compressible law of the wall in which the velocity and total enthalpy are logarithmic within the overlap zone, but in the Howarth-Dorodnitsyn variable. Simple outer region turbulence models are proposed (some of which are modifications of existing incompressible models) to reflect the effects of compressibility. As a test of the methodology and the new turbulence models, a set of self-similar outer region profiles is obtained for constant pressure flow; these are then coupled with embedded functions in the wall layer. The composite profiles thus obtained are compared directly with experimental data and good agreement is obtained for flows with Mach numbers up to 10.

He, J.↗
NASA NTRS1992

Conservative Grid-Interface Algorithm For Computing Flows

Best features of structured- and unstructured-grid methods combined. Gaps and overlaps between zonal grids eliminated by grid-interface algorithm, which generates single interfacial grid and corrects fluxes of flow quantities accordingly. Incorporated into two three-dimensional Navier-Stokes finite-volume codes and tested in computations of incompressible and compressible flows about simple bodies. Good numerical results obtained. General enough to be incorporated into other finite-volume codes without restrictions on complexities of shapes of bodies and zonal interfaces.

Klopfer, G. H.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1975

Finite element computation of a viscous compressible free shear flow governed by the time dependent Navier-Stokes equations

A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.

Cooke, C. H.↗
NASA NTRS1985

Numerical simulation of boundary layer excitation by surface heating/cooling

The concept of active control of growing disturbances in an unstable compressible flow by using time periodic, localized surface heating is studied numerically. The simulations are calculated by a fourth-order accurate solution of the compressible, laminar Navier-Stokes equations. Fourth-order accuracy is particularly important for this problem because the solution must be computed over many wavelengths. The numerical results demonstrate the growth of an initially small fluctuation into the nonlinear regime where a local breakdown into smaller scale disturbances can be observed. It is shown that periodic surface heating over a small strip can reduce the level of the fluctuation provided that the phase of the heating current is properly chosen.

Bayliss, A.↗
NASA NTRS1986

Numerical simulation of boundary-layer excitation by surface heating/cooling

The concept of active control of growing disturbances in an unstable compressible flow by using time periodic, localized surface heating is studied numerically. The simulations are calculated by a fourth-order accurate solution of the compressible, laminar Navier-Stokes equations. Fourth-order accuracy is particularly important for this problem because the solution must be computed over many wavelengths. The numerical results demonstrate the growth of an initially small fluctuation into the nonlinear regime where a local breakdown into smaller scale disturbances can be observed. It is shown that periodic surface heating over a small strip can reduce the level of the fluctuation provided that the phase of the heating current is properly chosen.

Bayliss, A.↗
NASA NTRS1993

Experimental investigations of the time and flow-direction responses of shear-stress-sensitive liquid crystal coatings

Time and flow-direction responses of shearstress-sensitive liquid crystal coatings were explored experimentally. For the time-response experiments, coatings were exposed to transient, compressible flows created during the startup and off-design operation of an injector-driven supersonic wind tunnel. Flow transients were visualized with a focusing Schlieren system and recorded with a 1000 frame/sec color video camera. Liquid crystal responses to these changing-shear environments were then recorded with the same video system, documenting color-play response times equal to, or faster than, the time interval between sequential frames (i.e., 1 millisecond). For the flow-direction experiments, a planar test surface was exposed to equal-magnitude and known-direction surface shear stresses generated by both normal and tangential subsonic jet-impingement flows. Under shear, the sense of the angular displacement of the liquid crystal dispersed (reflected) spectrum was found to be a function of the instantaneous direction of the applied shear. This technique thus renders dynamic flow reversals or flow divergences visible over entire test surfaces at image recording rates up to 1 KHz. Extensions of the technique to visualize relatively small changes in surface shear stress direction appear feasible.

Reda, Daniel C.↗
NASA NTRS1994

Time and flow-direction responses of shear-styress-sensitive liquid crystal coatings

Time and flow-direction responses of shear-stress liquid crystal coatings were exploresd experimentally. For the time-response experiments, coatings were exposed to transient, compressible flows created during the startup and off-design operation of an injector-driven supersonic wind tunnel. Flow transients were visualized with a focusing schlieren system and recorded with a 100 frame/s color video camera.

Reda, Daniel C.↗
NASA NTRS1993

A numerical analysis of supersonic flow over an axisymmetric afterbody

The performance of a pressure-based algorithm and the k-epsilon type turbulence models, including the original formulation proposed by Jones and Launder (1972), model extensions for refined treatment of time scales, and model corrections developed for treating compressible flows is assessed by analyzing the near-wake region of a Mach 2.5 flow over an axisymmetric cylindrical afterbody. The mechanisms responsible for the different levels of performance of the various turbulence models for the flow are investigated. It is shown that the k-epsilon models extended to allow for better response to the mean strain rate and corrected for compressibility result in better flowfield predictions than the original k-epsilon model.

Tucker, P. K.↗
NASA NTRS1992

Laminar hypersonic flow over a compression using the HANA code

A hypersonic blunt-body code, HANA (Hypersonic Analysis for Aerobrakes), has been used to model the flow with compression corner angles of 15 and 24 degrees. A comparison is made between experimental flow data and numerical results of previous investigations which include surface pressures, skin friction coefficients, heat transfer rates, and flow field geometries.

Fay, John F.↗
NASA NTRS1988

Flow-thermal-structural study of aerodynamically heated leading edges

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated leading edges is presented. The Navier-Stokes equations for high speed compressible flow, the energy equation, and the quasi-static equilibrium equations for the leading edge are solved using a single finite element approach in one integrated, vectorized computer program called LIFTS. The fluid-thermal-structural coupling is studied for Mach 6.47 flow over a 3-inch diameter cylinder for which the flow behavior and the aerothermal loads are calibrated by experimental data. Issues of the thermal-structural response are studied for hydrogen cooled, super thermal conducting leading edges subjected to intense aerodynamic heating.

Dechaumphai, Pramote↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗