SEARCH · Search NASA
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.
The breakdown of the conservative potential approximation occurs for all the airfoils tested. It develops as soon as the shock waves appear in the flow field. Since shock waves are not properly represented by the potential approximation, it is conjectured that the breakdown is due to the isentropic shock jump condition of the potential approximation.
Patched-grid algorithms for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver were investigated. Conservative and non-conservative approaches for performing zonal interpolations were implemented. The latter approach yields the most flexible technique in that it can handle both patched and overlaid grids. Results for a two-dimensional blunt body problem show that either approach yield accurate steady-state shock locations and jump conditions. In addition, calculations of the turbulent flow through a hypersonic inlet on a three-zone grid show that the numerical prediction is in good agreement with the experimental results. Through the use of a generalized coordinate transformation at the zonal interface between two or more blocks, the algorithm can be applied to highly stretched viscous grids and to arbitrarily-shaped zonal boundaries. Applications were made to the F-18 aircraft at subsonic, high-alpha conditions, in support of the NASA High-Alpha Research Program. The calculations were compared to ground-based and flight test experiments and were used as a guide to understanding the ground-based tests, which are laminar and transitional, and their relationship to flight. Calculations about a complete reconnaissance aircraft were also performed in order to further demonstrate the capability of the patched-grid algorithm.
The influence of electron thermal conduction on radiative shock structure is studied for both one- and two-temperature plasmas. The dimensionless ratio of the conductive length to the cooling length determines whether or not conduction is important, and shock jump conditions with conduction are established for a collisionless shock front. Approximate solutions are obtained, with the assumptions that the ionization state of the gas is constant and the cooling rate is a function of temperature alone. In the absence of magnetic fields, these solutions indicate that conduction noticeably influences normal-abundance interstellar shocks with velocities 50-100 km/s and dramatically affects metal-dominated shocks over a wide range of shock velocities.
A theory of evaporation, ignition, and burning of moderately dense spherical drop clusters has been developed. The theory takes into account burning of premixed air and fuel internal to the cluster at ignition and subsequent burning of fuel emitted from the cluster by a flame sheet surrounding it. The model considers interdrop interaction, momentum exchange between drops and gas, and turbulent exchange processes between the cluster and its surroundings. Calcualtions are performed for varying initial air to fuel mass ratios, initial cluster radii, ambient gas temperatures, and initial drop temperatures. Results are presented for fuel burn fractions at ignition and at the moment of drop disappearance, as well as jump conditions at ignition.
The paper presents recent work on evolution of recurrent solar wind structures in the outer heliosphere. Corotating shocks, corotating interaction regions, and merged interaction regions are studied, and an MHD simulation model in which the jump conditions at all shock crossings satisfy the Rankine-Hugoniot solution is examined. Simulation results which describe the evolution of idealized recurrent solar wind structures between 14 AU and the termination shock are reported which show that merged interaction regions belonging to two to four consecutive solar rotations coalesce with each other between 30 AU and the termination shock to reshape the heliospheric structures. In the outer heliosphere, forward and reverse interplanetary shocks which continuously heat the solar wind plasma are investigated.
The coronal slow shock was predicted to exist embedded in large coronal holes at 4 to 10 solar radii. A three-fluid model was used to study the jumps in minor ions propertes across the coronal slow shock. The jump conditions were formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. A simularity solution was also obtained for heavy ions. The results show that on the downstream side of the coronal slow shock the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.
The model used by the present unsteady cascade flow calculation method is based on the linearized unsteady Euler equations and accounts for blade loading and geometry, shock motion, and wake motion. The steady flow that must be determined prior to the unsteady flow is ascertained by means of a Newton iteration procedure. A noteworthy feature of the procedure is the use of shock-fitting to determine steady and unsteady shock positions; the use of the Euler equations in conjunction with the Rankine-Hugoniot shock-jump conditions is found to directly model the generation of entropy and vorticity at shocks. Results are presented for both channel and cascade flows.
A theory of evaporation, ignition, and burning of moderately dense spherical drop clusters has been developed. The theory takes into account burning of premixed air and fuel internal to the cluster at ignition and subsequent burning of fuel emitted from the cluster by a flame sheet surrounding it. The model considers interdrop interaction, momentum exchange between drops and gas, and turbulent exchange processes between the cluster and its surroundings. Calculations are performed for varying initial air-to-fuel-mass ratios, initial cluster radii, ambient gas temperatures and initial drop temperatures. Results are presented for ratios of fuel mass burned to fuel mass lost from the cluster between drop ignition and drop disappearance, fuel burned fractions at ignition and at the moment of drop disappearance, and jump conditions at ignition.
The objective is to better characterize the macroscopic properties of the interhemisphere plasma flow by solving a more complete set of hydrodynamic equations than that solved previously. Specifically, the ion continuity, momentum and energy equations were solved for the plasma flow along the closed magnetic field lines. During the initial stage of the supersonic outflow in the equatorial region, the ions cool substantially. Using the hydrodynamic model for the large-scale plasma flow, the dynamics of shocks was examined which form in the geomagnetic flux tubes during the early stages of refilling. These shocks are more like those forming in neutral gases than the electrostatic shocks driven by microinstabilities involving ion-ion interaction. Therefore, the shocks seen in the hydrodynamic model are termed as hydrodynamic shocks. Such shocks are generally unsteady and therefore the usual shock jump conditions given by Rankine-Hugoniot relations are not strictly applicable to them. The density, flow velocity and temperature structures associated with the shocks are examined for both asymmetrical and symmetrical flows. In the asymmetrical flow the outflow from one of two conjugate ionospheres is dominant. On the other hand, in the symmetrical case outflows from the two ionospheric sources are identical. Both cases are treated by a two-stream model. In the late type of flow, the early-time refilling shows a relaxation type of oscillation, which is driven by the large-scale interactions between the two identical streams. After this early stage, the resulting temperature structure shows some interesting features. In the equatorial region the streams are isothermal, but in the off-equatorial regions the streams have quite different temperatures, and also densities and flow velocities. The dense and slow stream is found to be warmer than the low-density fast stream. In the late stage of refilling, the temperature is found to steadily increase from the conjugate ionospheres towards the equator; the equatorial temperature is found to be as high as about 8000 K compared to the ionospheric temperature of 3600 K.
An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.
Panel methods are numerical schemes for solving (the Prandtl-Glauert equation) for linear, inviscid, irrotational flow about aircraft flying at subsonic or supersonic speeds. The tools at the panel-method user's disposal are (1) surface panels of source-doublet-vorticity distributions that can represent nearly arbitrary geometry, and (2) extremely versatile boundary condition capabilities that can frequently be used for creative modeling. Panel-method capabilities and limitations, basic concepts common to all panel-method codes, different choices that were made in the implementation of these concepts into working computer programs, and various modeling techniques involving boundary conditions, jump properties, and trailing wakes are discussed. An approach for extending the method to nonlinear transonic flow is also presented. Three appendices supplement the main test. In appendix 1, additional detail is provided on how the basic concepts are implemented into a specific computer program (PANAIR). In appendix 2, it is shown how to evaluate analytically the fundamental surface integral that arises in the expressions for influence-coefficients, and evaluate its jump property. In appendix 3, a simple example is used to illustrate the so-called finite part of the improper integrals.
The coronal slow shock has been predicted to exist embedded in large coronal holes at 4-10 solar radii. In this paper, a three-fluid model is used to study the jumps in minor ion properties across a slow shock such as the coronal slow shock. The jump conditions are formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. Also obtained is a similarity solution for heavy ions. The results show that, on the downstream side of the slow shock, the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe, in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.
An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.
One-dimensional particle-in-cell plasma simulations are used to examine the mechanical structure and thermalization properties of collisionless relativistic shock waves in electron-positron plasmas. Shocks propagating perpendicularly to the magnetic field direction are considered. It is shown that these shock waves exist, and that they are completely parameterized by the ratio of the upstream Poynting flux to the upstream kinetic energy flux. The way in which the Rankine-Hugoniot shock jump conditions are modified by the presence of wave fluctuations is shown, and they are used to provide a macroscopic description of these collisionless shock flows. The results of a 2D simulation that demonstrates the generality of these results beyond the assumption of the 1D case are discussed. It is suggested that the thermalization mechanism is the formation of a synchrotron maser by the coherently reflected particles in the shock front. Because the downstream medium is thermalized, it is argued that perpendicular shocks in pure electron-positron plasmas are not candidates as nonthermal particle accelerators.
The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in this work. A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models.
The Hollweg et al. (1990) analysis of MHD surface waves in a stationary equilibrium is extended. The conservation laws and jump conditions at Alfven and slow resonance points obtained by Sakurai et al. (1990) are generalized to include an equilibrium flow, and the assumption that the Eulerian perturbation of total pressure is constant is recovered as the special case of the conservation law for an equilibrium with straight magnetic field lines and flow along the magnetic field lines. It is shown that the conclusions formulated by Hollweg et al. are still valid for the straight cylindrical case. The effect of curvature is examined.
The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.
A two-dimensional model is developed for simulating heat pipes transient performance following changes in the input/rejection power or in the evaporator/condenser temperatures. The model employs the complete form of governing equations and momentum and energy jump conditions at the liquid-vapor interface. Although the model is capable of handling both cylindrical and rectangular geometries, the results reported are for a circular heat pipe with liquid lithium as the working fluid. The model incorporates a variety of other working fluids, such as water, ammonia, potassium, sodium and mercury, and offers combinations of isothermal, isoflux, convective and radiative heating/cooling conditions in the evaporator and condenser regions of the heat pipe. Results presented are for lithium heat pipes with exponential heating of the evaporator and isothermal cooling of the condenser.