SEARCH · Engineering Papers
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.
Prony coefficients found by computer program, without trial and error. Computational method and computer program developed to exploit full potential of Prony's interpolation method in analysis of experimental data on relaxation modules of viscoelastic material. Prony interpolation curve chosen to give least-squares best fit to "B-spline" interpolation of experimental data.
An approximate numerical method has been developed for computing the entire transonic flow field over an airfoil, including the boundary layer, the shock-wave/boundary-layer interaction, and the shock-induced-separation. The inviscid flow, the turbulent boundary layer in the weak interaction zone, and the strong interaction elliptic region are solved separately. The inviscid flow is represented by the solution of the transonic small-disturbance equation, while the boundary layer is solved by the generalized Galerkin method, and the strong interaction region (including the wake) is computed from the time-averaged Navier-Stokes equations. A second-order interaction is computed iteratively using the inviscid pressure distribution as a boundary condition for the boundary layer and Navier-Stokes codes, and the viscous displacement thickness to modify the airfoil geometry for the inviscid computation. The computed pressure distributions and shock-wave profiles are compared with experimental data. Agreement is very good except that the pressure recovery on the trailing edge is greater than the experimental data.
The oxidation of stoichiometric mixtures of benzene and toluene behind incident shock waves was studied for a temperature range from 1700 to 2800 K and a pressure range from 1.1 to 1.7 atm. The concentration of CO and CO2 produced were measured as well as the product of the oxygen atom and carbon monoxide concentrations. Comparisons between the benzene experimental data and results calculated by use of a reaction mechanism published in the open literature were carried out. With some additional reactions and changes in rate constants to reflect the pressure-temperature range of the experimental data, a good agreement was achieved between computed and experimental results. A reaction mechanism was developed for toluene oxidation based on analogous rate steps from the benzene mechanism. Measurements of NOx levels in an actual flame device, a jet-stirred combustor, were reproduced successfully by use of the reaction mechanism developed from the shock-tube experiments on toluene. These experimental measurements of NOx levels were reproduced from a computer simulation of a jet-stirred combustor.
A computational technique for computing the three-dimensional inviscid flow over blunt bodies having large regions of embedded subsonic flow is detailed. Results, which were obtained using the CDC Cyber 203 vector processing computer, are presented for several analytic shapes with some comparison to experimental data. Finally, windward surface pressure computations over the first third of the Space Shuttle vehicle are compared with experimental data for angles of attack between 25 and 45 degrees.
A combined experimental and numerical study has been conducted on the hypersonic shock-wave turbulent-boundary layer interaction at Mach 8.2 generated by a single fin of angles alpha = 10 and 15 deg. Three models are considered: (1) the 3D compressible Reynolds-averaged Navier-Stokes (RANS) equations using the k-epsilon turbulence model, (2) the 3D RANS using the Rodi turbulence model, and (3) the conical RANS using the Baldwin-Lomax algebraic turbulence model. The computations are compared with various experimental data. The computations using models (1) and (2) show quantitatively very similar results and very good agreement with experimental data for surface pressure and skin friction. Comparison with boundary layer profiles of pitot pressure and yaw angle are also generally good, but the peak surface heat transfer is overestimated by up to 48 percent. The effect of the laminar boundary layer on the fin is restricted to the immediate vicinity of the fin surface. Conical calculations using model (3) show substantially poorer agreement with experiment.
The influence of fiber waviness and matrix nonlinearity on the compressive behavior of continuous fiber composites is studied. A micromechanics model based on the kinematics of the fibers has been developed to predict the behavior of unidirectional composites with initially wavy fibers under compressive loads. The initial waviness has been idealized as sinusoidal. Nonlinear shear behavior of the matrix has been included. Because the shear strain in the matrix is a function of position, the deformed shape of the fibers differs from the initial shape. Therefore, the deformed fiber shape is represented by a sine series. Fiber waviness and matrix shear properties have been determined experimentally for T650-42/Radel C under several conditions. These experimental data have been used as input into the analytical model. Predicted compressive response is in reasonable agreement with experimental data.
Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.
Land mobile satellite communications at Ka-band (30/20 GHz) are attracting an increasing interest among researchers because of the frequency band availability and the possibility of small earth station designs. However, communications at the Ka-band pose significant challenges in the system designs due to severe channel impairments. Because only very limited experimental data for mobile applications at Ka-band is available, this paper studies the channel characteristics based on experimental data at L-band (1.6/1.5 GHz) and the use of frequency scaling. The land mobile satellite communication channel at Ka-band is modelled as log-normal Rayleigh fading channel. The first and second-order statistics of the fading channel are studied. The performance of a coherent BPSK system over the fading channel at L-band and K-band is evaluated theoretically and validated by computer simulations. Conclusions on the communication channel characteristics and system performance at L-band and Ka-band are presented.
The dissociation of WF6 and the related singly-charged cations and anions into the lower fluorides and fluorine atoms has been investigated theoretically using density functional theory (B3LYP) and relativistic effective core potentials, with estimates of spin-orbit effects included using a simple model. The inclusion of spin-orbit is essential for a correct description of the thermochemistry. The total atomization energy of the neutral and anionic WF6 is reproduced to within 25 kcal/mol, but comparison of individual bond dissociation energies with available experimental data shows discrepancies of up to 10 kcal/mol. The results are nevertheless useful to help resolve discrepancies in experimental data and provide estimates of missing data.
This paper describes the process used for estimating flow-through balance momentum tares. The interaction of jet engine exhausts on the BOEINGERA Hybrid Wing Body (HWB) was simulated in the NFAC 40x80 wind tunnel at NASA Ames using a pair of turbine powered simulators (TPS). High-pressure air was passed through a flow-through balance and manifold before being delivered to the TPS units. The force and moment tares that result from the internal shear and pressure distribution were estimated using CFD. Validation of the CFD simulations for these complex internal flows is a challenge, given limited experimental data due to the complications of the internal geometry. Two CFD validation efforts are documented, and comparisons with experimental data from the final model installation are provided.
This white paper focuses on the scientific rationale and motivation for deploying phase change processes-based spacecraft thermal management systems for space missions [1,2]. Phase change processes need to be implemented in the design of various life support systems, sub-critical cryogenic and thermal control systems, taking advantage of the latent heat of phase change so as to significantly increase the heat transfer rate per unit fluid mass. Experimental data are critical to understanding gravitational effects on complex phase change processes using novel materials and working fluids. Experimental data can be acquired by either designing new testing platforms or performing experiments with simulant fluids on existing microgravity platforms. With the acquired data, theoretical flow physics models and numerical simulations can then be validated for use under long-duration microgravity conditions. The acquired data and validated models will enable the deployment of game-changing thermal management systems with remarkably improved efficiency and stability in heat acquisition, transport, and rejection for long-duration space missions.
The sloshing of propellants can affect the stability of a spacecraft and the integrity of the tank structure. Undesirable sloshing can be controlled by the addition of anti-slosh baffles, and the spacing and configuration of baffles are driven by damping requirements. The structural design of the baffle is determined after consideration of many factors, such as the strength and rigidity needed to support the baffle for its lifetime. Therefore, knowledge of distributed pressure loading is important for detailed structural design. In addition, the resultant force and moment produced by the distributed pressure are of direct importance to the design of a vehicle’s control system. Previous experimental investigations have been conducted to determine the liquid pressure loads and slosh damping associated with a rigid ring baffle. The results suggested that when the nondimensional velocity parameter is larger than 1.0, the theories agree with the test. However, when the velocity parameter is less than 1.0, all theories are nonconservative and under-predict the pressure loads. The present study has derived a maximum pressure load on the slosh baffle based on the energy conservation principle. It is verified from CFD that pressure in the slosh flow field can be decomposed into static and transient components. The CFD results confirm that there is a phase shift in pressure across the baffle, which depends on the fluid damping. Higher damping leads to a higher phase shift. The CFD investigation further verifies the proposed theory: the maximum pressure load occurs when the phase shift is 90 degrees. A comparison of the present computational results to the previous comprehensive experimental data validates the maximum pressure theory. When the baffle is submerged, the maximum pressure theory envelopes all the experimental data points.
A new approach is presented for solving the problem of predicting dynamic stall characteristics. In this approach, the unsteady aerodynamic characteristics are related theoretically to static aerodynamic characteristics which are readily available for a great number of airfoil shapes. Using these static experimental data as an input, the developed analytical method predicts dynamic stall characteristics that are in good agreement with available experimental data. The analysis is also extended to include frequency and amplitude modulation effects.
Equations based on the mixtures rule are presented for predicting the physical, thermal, hygral, and mechanical properties of unidirectional intraply hybrid composites (UIHC) from the corresponding properties of their constituent composites. Bounds were derived for uniaxial longitudinal strengths, tension, compression, and flexure of UIHC. The equations predict shear and flexural properties which agree with experimental data from UIHC. Use of these equations in a composites mechanics computer code predicted flexural moduli which agree with experimental data from various intraply hybrid angleplied laminates (IHAL). It is indicated, briefly, how these equations can be used in conjunction with composite mechanics and structural analysis during the analysis/design process.
An Euler flow solver and a thin layer Navier-Stokes flow solver were used to numerically simulate the supersonic leeside flow fields over delta wings which were observed experimentally. Three delta wings with 75, 67.5, and 60 deg leading edge sweeps were computed over an angle-of-attack range of 4 to 20 deg at a Mach number 2.8. The Euler code and Navier-Stokes code predict equally well the primary flow structure where the flow is expected to be separated or attached at the leading edge based on the Stanbrook-Squire boundary. The Navier-Stokes code is capable of predicting both the primary and the secondary flow features for the parameter range investigated. For those flow conditions where the Euler code did not predict the correct type of primary flow structure, the Navier-Stokes code illustrated that the flow structure is sensitive to boundary layer model. In general, the laminar Navier-Stokes solutions agreed better with the experimental data, especially for the lower sweep delta wings. The computational results and a detailed re-examination of the experimental data resulted in a refinement of the flow classifications. This refinement in the flow classification results in the separation bubble with the shock flow type as the intermediate flow pattern between separated and attached flows.
An approximate solution technique has been developed for three-dimensional, inviscid, hypersonic flows. The method uses Maslen's explicit pressure equation and the assumption of approximate stream surfaces in the shock layer. This approximation represents a simplification of Maslen's asymmetric method. The solution procedure involves iteratively changing the shock shape in the subsonic-transonic region until the correct body shape is obtained. Beyond this region, the shock surface is determined by using a marching procedure. Results are presented herein for a paraboloid and elliptic cone at angle of attack. Calculated surface pressure distributions, shock shapes, and property profiles are compared with experimental data and finite-difference solutions of the Euler equations. Comparisons of the results of the present method with experimental data and detailed predictions are very good. Since the present method provides a very rapid computational procedure, it can be used for parametric or preliminary design applications. One useful application would be to incorporate a heating procedure for aerothermal studies.
A computational model is presented for the prediction of solid/liquid phase change energy transport including the influence of free convection fluid flow in the liquid phase region. The computational model considers the velocity components of all non-liquid phase change material control volumes to be zero but fully solves the coupled mass-momentum problem within the liquid region. The thermal energy model includes the entire domain and uses an enthalpy like model and a recently developed method for handling the phase change interface nonlinearity. Convergence studies are performed and comparisons made with experimental data for two different problem specifications. The convergence studies indicate that grid independence was achieved and the comparison with experimental data indicates excellent quantitative prediction of the melt fraction evolution. Qualitative data is also provided in the form of velocity vector diagrams and isotherm plots for selected times in the evolution of both problems. The computational costs incurred are quite low by comparison with previous efforts on solving these problems.
A numerical study based on the 3D Reynolds-averaged Navier-Stokes equation has been conducted to investigate the detailed flow physics inside a transonic compressor. 3D shock structure, shock-boundary layer interaction, flow separation, radial mixing, and wake development are all investigated at design and off-design conditions. Experimental data based on laser anemometer measurements are used to assess the overall quality of the numerical solution. An additional experimental study to investigate end-wall flow with a hot-film was conducted, and these results are compared with the numerical results. Detailed comparison with experimental data indicates that the overall features of the 3D shock structure, the shock-boundary layer interaction, and the wake development are all calculated very well in the numerical solution. The numerical results are further analyzed to examine the radial mixing phenomena in the transonic compressor. A thin sheet of particles is injected in the numerical solution upstream of the compressor. The movement of particles is traced with a 3D plotting package. This numerical survey of tracer concentration reveals the fundamental mechanisms of radial transport in this transonic compressor.