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.
Three new options were incorporated into an existing computer program for the design and analysis of low speed airfoils. These options permit the analysis of airfoils having variable chord (variable geometry), a boundary layer displacement iteration, and the analysis of the effect of single roughness elements. All three options are described in detail and are included in the FORTRAN IV computer program.
Experimental measurements in the inlet of a transonic turbine blade cascade showed unacceptable pitchwise flow non-uniformity. A three-dimensional, Navier-Stokes computational fluid dynamics (CFD) analysis of the imbedded bellmouth inlet in the facility was performed to identify and eliminate the source of the flow non-uniformity. The blockage and acceleration effects of the blades were accounted for by specifying a periodic static pressure exit condition interpolated from a separate three-dimensional Navier-Stokes CFD solution of flow around a single blade in an infinite cascade. Calculations of the original inlet geometry showed total pressure loss regions consistent in strength and location to experimental measurements. The results indicate that the distortions were caused by a pair of streamwise vortices that originated as a result of the interaction of the flow with the imbedded bellmouth. Computations were performed for an inlet geometry which eliminated the imbedded bellmouth by bridging the region between it and the upstream wall. This analysis indicated that eliminating the imbedded bellmouth nozzle also eliminates the pair of vortices, resulting in a flow with much greater pitchwise uniformity. Measurements taken with an installed redesigned inlet verify that the flow non-uniformity has indeed been eliminated.
An axisymmetric full Navier-Stokes computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy, grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex. Although the one-dimensional isentropic nozzle plume results look reasonable, they fail to capture the sonic boom shock train in the highly underexpanded nozzle flow.
A computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the Sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy. grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex.
Magnetic nozzles are used in many laboratory experiments in which plasma flows are to be confined, cooled. accelerated, or directed. At present, however, there is no generally accepted theoretical description that explains the phenomena of plasma detachment from an externally-imposed magnetic field. This is an important problem in the field of plasma propulsion, where the ionized gas must detach from the applied, solenoidal magnetic field to realize thrust production. In this paper we simulate a plasma flowing in the presence of an applied magnetic field using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The new tool employs the same formulation of the governing equation set. but retains the axial and radial components of magnetic field and the azimuthal component of velocity that were neglected in other works. We aim to compare the computational results with the various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kurtki and Okada. We model this geometry because there ts a substantial amount of experimental data that can be compared to our computations, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how the 3 scale with different input parameters.
Magnetic nozzles are used in many laboratory experiments in which plasma flows are to be confined, cooled, accelerated, or directed. At present, however, there is no generally accepted theoretical description that explains the phenomena of plasma expansion in and detachment from an externally-imposed magnetic field. The latter is an especially important problem in the field of plasma propulsion, where the ionized gas must detach from the applied, solenoidal magnetic field to realize thrust production. In this paper we simulate a plasma flowing in the presence of an applied magnetic field using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The new tool employs the same formulation of the governing equation set, but retains the axial and radial components of magnetic field and the azimuthal component of velocity that were neglected. We aim to compare the computational results with the various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver, and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. 4 A schematic showing the setup used in those experiments is shown. We model this geometry because there is a substantial amount of experimental data that can be compared to our computations, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters.
Accurate heats of formation are computed for PHn and PHn+, for n=1-3, and for PHn-, for n= 1-2. The geometries and vibrational frequencies are determined at the B3LYP level of theory. The energetics are determined at the CCSD(T) level of theory. Basis set limit values are obtained by extrapolation. Spin-orbit effects are taken from experiment. The temperature dependence of the heat of formation, heat capacity, and entropy are computed for the temperature range 300 to 4000 K and fit to a polynomial.
Based on multilevel partitioning, this paper develops a structural parallelizable solution methodology that enables a significant reduction in computational effort and memory requirements for very large scale linear and nonlinear steady and transient thermal (heat conduction) models. Due to the generality of the formulation of the scheme, both finite element and finite difference simulations can be treated. Diverse model topologies can thus be handled, including both simply and multiply connected (branched/perforated) geometries. To verify the methodology, analytical and numerical benchmark trends are verified in both sequential and parallel computer environments.
Results of detailed structural studies, using elastic neutron scattering, of the freezing of liquid O2 and D2 in porous vycor glass, are presented. The experimental studies have been complemented by computer simulations of the dynamics of freezing of a Lennard-Jones liquid in narrow channels bounded by molecular walls. Results point to a new simple physical interpretation of freezing in confined geometries.
A gradient-based shape optimization based on quasi-analytical sensitivities has been extended for practical three-dimensional aerodynamic applications. The flow analysis has been rendered by a fully implicit, finite-volume formulation of the Euler and Thin-Layer Navier-Stokes (TLNS) equations. Initially, the viscous laminar flow analysis for a wing has been compared with an independent computational fluid dynamics (CFD) code which has been extensively validated. The new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4 with coarse- and fine-grid based computations performed with Euler and TLNS equations. The influence of the initial constraints on the geometry and aerodynamics of the optimized shape has been explored. Various final shapes generated for an identical initial problem formulation but with different optimization path options (coarse or fine grid, Euler or TLNS), have been aerodynamically evaluated via a common fine-grid TLNS-based analysis. The initial constraint conditions show significant bearing on the optimization results. Also, the results demonstrate that to produce an aerodynamically efficient design, it is imperative to include the viscous physics in the optimization procedure with the proper resolution. Based upon the present results, to better utilize the scarce computational resources, it is recommended that, a number of viscous coarse grid cases using either a preconditioned bi-conjugate gradient (PbCG) or an alternating-direction-implicit (ADI) method, should initially be employed to improve the optimization problem definition, the design space and initial shape. Optimized shapes should subsequently be analyzed using a high fidelity (viscous with fine-grid resolution) flow analysis to evaluate their true performance potential. Finally, a viscous fine-grid-based shape optimization should be conducted, using an ADI method, to accurately obtain the final optimized shape.
Computational predictions of ice accretion on flying aircraft most commonly rely on modeling in 2D. These 2D methods treat an aircraft geometry either as wing-like with infinite span, or as an axisymmetric body. Recently, fully 3D methods have been introduced that model an aircraft's true 3D shape. Because 3D methods are more computationally expensive than 2D methods, 2D methods continue to be widely used. However, a 3D method allows investigation of whether it is valid to continue applying 2D methods to a finite wing. The extent of disagreement between LEWICE, a 2D method, and LEWICE3D, a 3D method, in calculating local collection efficiencies at the leading edge of finite wings is investigated.
The observations of the Jovian aurora made by the Voyager UV spectrometer (UVS) during 1979 are analyzed, with special consideration given to the model fitting process. Several different estimates of the function specifying auroral arc geometry were tried, with the Broadfoot et al. (1981) estimate, found by probing with the tip of the UVS slit, giving the best agreement; this agreement is slightly better than the agreement with the Io torus footprint computed by Roederer et al. (1977). The results suggest that the UVS observations are more sensitive to the surface field geometry than are the Voyager flyby in situ observations. The results of a study of intensity maximum positions indicate that the particles exciting the spatially variable portion of the aurora are drifting west, implying either that these particles are electrons or that they are positive ions drifting east more slowly than the corotation lag of the Io torus region carries them to the west. The latter case is most consistent with the high-energy charged particle measurements.
The sublaminate or ply-level analysis of composite structures is presently undertaken by a computational procedure yielding the stresses in regions affected by delaminations, transverse cracks, and discontinuities that are related to material properties, geometries, and loads. Attention is given to layers or groups of layers that are immediately affected by flaws; these are analyzed as if they were homogeneous bodies in equilibrium, in isolation from the rest of the laminate. Computed stresses agree with those from a three-dimensional FEM analysis.
A computer code called NCOREL (for Nonconical Relaxation) has been developed to solve for supersonic full potential flows over complex geometries. The method first solves for the conical at the apex and then marches downstream in a spherical coordinate system. Implicit relaxation techniques are used to numerically solve the full potential equation at each subsequent crossflow plane. Many improvements have been made to the original code including more reliable numerics for computing wing-body flows with multiple embedded shocks, inlet flow through simulation, wake model and entropy corrections. Line relaxation or approximate factorization schemes are optionally available. Improved internal grid generation using analytic conformal mappings, supported by a simple geometric Harris wave drag input that was originally developed for panel methods and internal geometry package are some of the new features.
Several finite difference schemes are applied to the stress and free vibration analysis of homogeneous isotropic and layered orthotropic shells of revolution. The study is based on a form of the Sanders-Budiansky first-approximation linear shell theory modified such that the effects of shear deformation and rotary inertia are included. A Fourier approach is used in which all the shell stress resultants and displacements are expanded in a Fourier series in the circumferential direction, and the governing equations reduce to ordinary differential equations in the meridional direction. While primary attention is given to finite difference schemes used in conjunction with first order differential equation formulation, comparison is made with finite difference schemes used with other formulations. These finite difference discretization models are compared with respect to simplicity of application, convergence characteristics, and computational efficiency. Numerical studies are presented for the effects of variations in shell geometry and lamination parameters on the accuracy and convergence of the solutions obtained by the different finite difference schemes. On the basis of the present study it is shown that the mixed finite difference scheme based on the first order differential equation formulation and two interlacing grids for the different fundamental unknowns combines a number of advantages over other finite difference schemes previously reported in the literature.
A numerical analysis was developed to determine the airloads on helicopter rotors operating under near-hovering flight conditions capable of producing impulsive noise. A computer program was written in which the solutions for the rotor tip vortex geometry, inflow, aeroelastic response, and airloads are solved in a coupled manner at sequential time steps, with or without the influence of an imposed steady ambient wind or transient gust. The program was developed for future applications in which predicted airloads would be incorporated in an acoustics analysis to attempt to predict and analyze impulsive noise (blade slap). The analysis was applied to a hovering full-scale rotor for which impulsive noise was recorded in the presence of ambient wind. The predicted tip vortex coordinates are in reasonable agreement with the test data, and the blade airload solutions converged to a periodic behavior for an imposed steady ambient wind conditions.
The cooling of post-flare plasmas is discussed and the formation of loop prominences is explained as due to a thermal instability. A one-dimensional model was developed for active loop prominences. Only the motion and heat fluxes parallel to the existing magnetic fields are considered. The relevant size scales and time scales are such that single-fluid MHD equations are valid. The effects of gravity, the geometry of the field and conduction losses to the chromosphere are included. A computer code was constructed to solve the model equations. Basically, the system is treated as an initial value problem (with certain boundary conditions at the chromosphere-corona transition region), and a two-step time differencing scheme is used.
An analysis of the geometry of the burning in both star-cylindrical port interface regions and regions of partially inhibited slots is presented. Some characteristics parameters are defined and illustrated. Methods are proposed for calculating burning areas which functionally depend only on the total distance burned. According to this method, several points are defined where abrupt changes in geometry occur, and these are tracked throughout the burn. Equations are developed for computing port perimeter and port area at pre-established longitudinal positions. Some common formulas and some newly developed formulas are then used to compute burning surface area and port volume. Some specific results are presented for the solid rocket motor committed to the space shuttle project.