Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface codes”

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 1,729 records · Page 96

Thermal study of the Missouri River in North Dakota using infrared imagery

Studies of infrared imagery obtained from aircraft at 305- to 1,524-meter altitudes indicate the feasibility of monitoring thermal changes attributable to the operation of thermal electric plants and storage reservoirs, as well as natural phenomena such as tributary inflow and ground water seeps in large rivers. No identifiable sources of ground water inflow below the surface of the river could be found in the imagery. The thermal patterns from the generating plants and the major tributary inflow are readily apparent in imagery obtained from an altitude of 305 meters. Portions of the tape-recorded imagery were processed in a color-coded quantization to enhance the displays and to attach quantitative significance to the data. The study indicates a marked decrease in water temperature in the Missouri River prior to early fall and a moderate increase in temperature in late fall because of the Lake Sakakawea impoundment.

Crosby, O. A.↗

Space environmental interactions with biased spacecraft surfaces

Large, high-voltage space power systems are being proposed for future space missions. These systems must operate in the charged-particle environment of space, and interactions between this environment and the high-voltage surfaces are possible. Ground simulation testing has indicated that dielectric surfaces that usually surround biased conductors can influence these interactions. For positive voltages greater than 100 V, it has been found that the dielectrics contribute to discharges. Using these experimental results a large, high-voltage power system operating in geosynchronous orbit was analyzed with the NASCAP code. Results of this analysis indicated that very strong electric fields exist in these power systems. A technology investigation is required to understand the interactions and develop techniques to alleviate any impact on power system performance.

Stevens, N. J.↗

Results of a flight experiment on the catalytic efficiency of the Space Shuttle heat shield

An experiment was performed on the second Space Shuttle flight to verify data from arc-heated wind-tunnel tests indicating that the thermal protection tiles used on the Shuttle's surface are noncatalytic to the recombination of dissociated air. The experiment indicates that the tiles are indeed noncatalytic, and explains the lower than expected temperatures observed in the first two Shuttle flights. The data are reasonably well predicted with computations, using a chemically reactive boundary-layer code. The experimental setup and the computational analysis are described.

Rakich, J. V.↗

A fast viscous correction method for full-potential transonic wing analysis

An analysis of the transonic flowfield around a three-dimensional wing is carried out using a strip method. Attention is given to the boundary layer growth in the streamwise direction. A viscous correction technique is defined for the TWING code for solving the full potential equations. A viscous ramp at the base of a shock is superimposed on the boundary layer displacement thickness generated by an integral boundary layer method. A relationship is then obtained between the effective displacement thickness and a vertical component of the surface velocity, a transpirational boundary condition. The viscous correction is found to be unnecessary in weak shock conditions but gives a better shock position and pressure distribution in a strong shock condition when compared with data from an ONERA M6 airfoil and the Hinson and Burdges (1980) Wing A.

Lee, S. C.↗

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for Mach number 3.5 and Reynolds number of 2,000,000. The simulation is done by integrating the viscous Navier-Stokes equations. These equations govern the unsteady, viscous, compressible and heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. The O-type and H-type grids have been used for this study, and results are compared. The governing equations are solved by the MacCormack time-split method, and the results are compared with other theoretical and experimental results. The codes are written in FORTRAN, vectorized and currently run on the CDC Vector Processing System (VPS-32) computer.

Abolhassani, Jamshid S.↗

Analytical study of the liquid phase transient behavior of a high temperature heat pipe

The transient operation of the liquid phase of a high temperature heat pipe is studied. The study was conducted in support of advanced heat pipe applications that require reliable transport of high temperature drops and significant distances under a broad spectrum of operating conditions. The heat pipe configuration studied consists of a sealed cylindrical enclosure containing a capillary wick structure and sodium working fluid. The wick is an annular flow channel configuration formed between the enclosure interior wall and a concentric cylindrical tube of fine pore screen. The study approach is analytical through the solution of the governing equations. The energy equation is solved over the pipe wall and liquid region using the finite difference Peaceman-Rachford alternating direction implicit numerical method. The continuity and momentum equations are solved over the liquid region by the integral method. The energy equation and liquid dynamics equation are tightly coupled due to the phase change process at the liquid-vapor interface. A kinetic theory model is used to define the phase change process in terms of the temperature jump between the liquid-vapor surface and the bulk vapor. Extensive auxiliary relations, including sodium properties as functions of temperature, are used to close the analytical system. The solution procedure is implemented in a FORTRAN algorithm with some optimization features to take advantage of the IBM System/370 Model 3090 vectorization facility. The code was intended for coupling to a vapor phase algorithm so that the entire heat pipe problem could be solved. As a test of code capabilities, the vapor phase was approximated in a simple manner.

Roche, Gregory Lawrence↗

Computations of unsteady multistage compressor flows in a workstation environment

High-end graphics workstations are becoming a necessary tool in the computational fluid dynamics environment. In addition to their graphic capabilities, workstations of the latest generation have powerful floating-point-operation capabilities. As workstations become common, they could provide valuable computing time for such applications as turbomachinery flow calculations. This report discusses the issues involved in implementing an unsteady, viscous multistage-turbomachinery code (STAGE-2) on workstations. It then describes work in which the workstation version of STAGE-2 was used to study the effects of axial-gap spacing on the time-averaged and unsteady flow within a 2 1/2-stage compressor. The results included time-averaged surface pressures, time-averaged pressure contours, standard deviation of pressure contours, pressure amplitudes, and force polar plots.

Gundy-Burlet, Karen L.↗

A determination of the external forces required to move the benchmark active controls testing model in pure plunge and pure pitch

In view of the strong need for a well-documented set of experimental data which is suitable for the validation and/or calibration of modern Computational Fluid Dynamics codes, the Benchmark Models Program was initiated by the Structural Dynamics Division of the NASA Langley Research Center. One of the models in the program, the Benchmark Active Controls Testing Model, consists of a rigid wing of rectangular planform with a NACA 0012 profile and three control surfaces (a trailing-edge control surface, a lower-surface spoiler, and an upper-surface spoiler). The model is affixed to a flexible mount system which allows only plunging and/or pitching motion. An approximate analytical determination of the forces required to move this model, with its control surfaces fixed, in pure plunge and pure pitch at a number of test conditions is included. This provides a good indication of the type of actuator system required to generate the aerodynamic data resulting from pure plunging and pure pitching motion, in which much interest was expressed. The analysis makes use of previously obtained numerical results.

Dcruz, Jonathan↗

Relativistic tidal interaction of a white dwarf with a massive black hole

We compute encounters of a realistic white dwarf model with a massive black hole in the regime where relativistic effects are important, using a three-dimensional, finite-difference, Eulerian, piecewise parabolic method (PPM) hydrodynamical code. Both disruptive and nondisruptive encounters are considered. We identify and discuss relativistic effects important for the problem: relativistic shift of the pericenter distance, time delay, relativistic precession, and the tensorial structure of the tidal forces. In the nondisruptive case, stripping of matter takes place. In the surface layers of the surviving core, complicated hydrodynamical phenomena are revealed. In both disruptive and nondispruptive encounters, material flows out in the form of two thin, S-shaped, supersonic jets. Our results provide realistic initial conditions for the subsequent investigation of the dynamics of the debris in the field of the black hole. We evaluate the critical conditions for complete disruption of the white dwarf, and compare our results with the corresponding results for nonrelativistic encounters.

Frolov, V. P.↗

The Use of Kirchhoff's Method in Jet Aeroacoustics

Supersonic jet aeroacoustics will be studied using computational techniques. In the study, a Kirchhoff method is used to predict flow generated noise in the mid- and far-fields. This type of method shows promise because it is based on surface integrals and not the volume integrals found in traditional acoustic prediction methods. The Kirchhoff method is dependent on accurate prediction of flow variables in the near-field. Here, computational fluid dynamics (CFD) programs are used for these predictions. Specifically, an existing large eddy simulation (LES) code will be modified for aeroacoustic applications. Issues involved in the implementation of the Kirchhoff method as well as the coupling with the CFD code will be discussed. Important physical noise parameters will be identified and investigated in the study.

Lyrintzis, Anastasios S.↗

Simulation of Conformal Spiral Slot Antennas on Composite Platforms

During the course of the grant, we wrote and distributed about 12 reports and an equal number of journal papers supported fully or in part by this grant. The list of reports (title & abstract) and papers are given in Appendices A and B. This grant has indeed been instrumental in developing a robust hybrid finite element method for the analysis of complex broadband antennas on doubly curved platforms. Previous to the grant, our capability was limited to simple printed patch antennas on mostly planar platforms. More specifically: (1) mixed element formulations were developed and new edge-based prisms were introduced; (2) these elements were important in permitting flexibility in geometry gridding for most antennas of interest; (3) new perfectly matched absorbers were introduced for mesh truncations associated with highly curved surfaces; (4) fast integral algorithms were introduced for boundary integral truncations reducing CPU time from O(N-2) down to O(N-1.5) or less; (5) frequency extrapolation schemes were developed for efficient broadband performance evaluations. This activity has been successfully continued by NASA researchers; (6) computer codes were developed and extensively tested for several broadband configurations. These include FEMA-CYL, FEMA-PRISM and FEMA-TETRA written by L. Kempel, T. Ozdemir and J. Gong, respectively; (7) a new infinite balun feed was designed nearly constant impedance over the 800-3000 MHz operational band; (8) a complete slot spiral antenna was developed, fabricated and tested at NASA Langley. This new design is a culmination of the projects goals and integrates the computational and experimental efforts. this antenna design resulted in a U.S. patent and was revised three times to achieve the desired bandwidth and gain requirements from 800-3000 MHz.

Volakis, J. L.↗

An Integrated Tool for the Coupled Thermal and Mechanical Analysis of Pyrolyzing Heatshield Materials

Materials that pyrolyze at elevated temperature have been commonly used as thermal protection materials in hypersonic flight, and advanced pyrolyzing materials for this purpose continue to be developed. Because of the large temperature gradients that can arise in thermal protection materials, significant thermal stresses can develop. Advanced applications of pyrolytic materials are calling for more complex heatshield configurations, making accurate thermal stress analysis more important, and more challenging. For non-pyrolyzing materials, many finite element codes are available and capable of performing coupled thermal-mechanical analyses. These codes do not, however, have a built-in capability to perform analyses that include pyrolysis effects. When a pyrolyzing material is heated, one or more components of the original virgin material pyrolyze and create a gas. This gas flows away from the pyrolysis zone to the surface, resulting in a reduction in surface heating. A porous residue, referred to as char, remains in place of the virgin material. While the processes involved can be complex, it has been found that a simple physical model in which virgin material reacts to form char and pyrolysis gas, will yield satisfactory analytical results. Specifically, the effects that must be modeled include: (1) Variation of thermal properties (density, specific heat, thermal conductivity) as the material composition changes; (2) Energy released or absorbed by the pyrolysis reactions; (3) Energy convected by the flow of pyrolysis gas from the interior to the surface; (4) The reduction in surface heating due to surface blowing; and (5) Chemical and mass diffusion effects at the surface between the pyrolysis gas and edge gas Computational tools for the one-dimensional thermal analysis these materials exist and have proven to be reliable design tools. The objective of the present work is to extend the analysis capabilities of pyrolyzing materials to axisymmetric configurations, and to couple thermal and mechanical analyses so that thermal stresses may be efficiently and accurately calculated.

Pronchick, Stephen W.↗

High-Area-Ratio Rocket Nozzle at High Combustion Chamber Pressure: Experimental and Analytical Validation

Experimental data were obtained on an optimally contoured nozzle with an area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 1800 to 2400 psia and mixture ratios of 3.89 to 6.15. This report compares the experimental performance, heat transfer, and boundary layer total pressure measurements with theoretical predictions of the current Joint Army, Navy, NASA, Air Force (JANNAF) developed methodology. This methodology makes use of the Two-Dimensional Kinetics (TDK) nozzle performance code. Comparisons of the TDK-predicted performance to experimentally attained thrust performance indicated that both the vacuum thrust coefficient and the vacuum specific impulse values were approximately 2.0-percent higher than the turbulent prediction for the 1025:1 configurations, and approximately 0.25-percent higher than the turbulent prediction for the 440:1 configuration. Nozzle wall temperatures were measured on the outside of a thin-walled heat sink nozzle during the test fittings. Nozzle heat fluxes were calculated front the time histories of these temperatures and compared with predictions made with the TDK code. The heat flux values were overpredicted for all cases. The results range from nearly 100 percent at an area ratio of 50 to only approximately 3 percent at an area ratio of 975. Values of the integral of the heat flux as a function of nozzle surface area were also calculated. Comparisons of the experiment with analyses of the heat flux and the heat rate per axial length also show that the experimental values were lower than the predicted value. Three boundary layer rakes mounted on the nozzle exit were used for boundary layer measurements. This arrangement allowed total pressure measurements to be obtained at 14 different distances from the nozzle wall. A comparison of boundary layer total pressure profiles and analytical predictions show good agreement for the first 0.5 in. from the nozzle wall; but the further into the core flow that measurements were taken, the more that TDK overpredicted the boundary layer thickness.

Jankovsky, Robert S.↗

Use of a Hybrid Edge Node-Centroid Node Approach to Thermal Modeling

A recent proposal submitted for an ESA mission required that models be delivered in ESARAD/ESATAN formats. ThermalDesktop was the preferable analysis code to be used for model development with a conversion done as the final step before delivery. However, due to some differences between the capabilities of the two codes, a unique approach was developed to take advantage of the edge node capability of ThermalDesktop while maintaining the centroid node approach used by ESARAD. In essence, two separate meshes were used: one for conduction and one for radiation. The conduction calculations were eliminated from the radiation surfaces and the capacitance and radiative calculations were eliminated from the conduction surfaces. The resulting conduction surface nodes were coincident with all nodes of the radiation surface and were subsequently merged, while the nodes along the edges remained free. Merging of nodes on the edges of adjacent surfaces provided the conductive links between surfaces. Lastly, all nodes along edges were placed into the subnetwork and the resulting supernetwork included only the nodes associated with radiation surfaces. This approach had both benefits and disadvantages. The use of centroid, surface based radiation reduces the overall size of the radiation network, which is often the most computationally intensive part of the modeling process. Furthermore, using the conduction surfaces and allowing ThermalDesktop to calculate the conduction network can save significant time by not having to manually generate the couplings. Lastly, the resulting GMM/TMM models can be exported to formats which do not support edge nodes. One drawback, however, is the necessity to maintain two sets of surfaces. This requires additional care on the part of the analyst to ensure communication between the conductive and radiative surfaces in the resulting overall network. However, with more frequent use of this technique, the benefits of this approach can far outweigh the additional effort.

Peabody, Hume L.↗

Use of a Hybrid Edge Node-Centroid Node Approach to Thermal Modeling

A recent proposal submitted for an ESA mission required that models be delivered in ESARAD/ESAT AN formats. ThermalDesktop was the preferable analysis code to be used for model development with a conversion done as the final step before delivery. However, due to some differences between the capabilities of the two codes, a unique approach was developed to take advantage of the edge node capability of ThermalDesktop while maintaining the centroid node approach used by ESARAD. In essence, two separate meshes were used: one for conduction and one for radiation. The conduction calculations were eliminated from the radiation surfaces and the capacitance and radiative calculations were eliminated from the conduction surfaces. The resulting conduction surface nodes were coincident with all nodes of the radiation surface and were subsequently merged, while the nodes along the edges remained free. Merging of nodes on the edges of adjacent surfaces provided the conductive links between surfaces. Lastly, all nodes along edges were placed into the subnetwork and the resulting supernetwork included only the nodes associated with radiation surfaces. This approach had both benefits and disadvantages. The use of centroid, surface based radiation reduces the overall size of the radiation network, which is often the most computationally intensive part of the modeling process. Furthermore, using the conduction surfaces and allowing ThermalDesktop to calculate the conduction network can save significant time by not having to manually generate the couplings. Lastly, the resulting GMM/TMM models can be exported to formats which do not support edge nodes. One drawback, however, is the necessity to maintain two sets of surfaces. This requires additional care on the part of the analyst to ensure communication between the conductive and radiative surfaces in the resulting overall network. However, with more frequent use of this technique, the benefits of this approach can far outweigh the additional effort.

Peabody, Hume L.↗

Computational Component Build-Up for the X-57 Distributed Electric Propulsion Aircraft

A computational study of the wing for the distributed electric propulsion X-57 Maxwell airplane configuration at cruise and takeoff/landing conditions was completed. Three unstructured-mesh, Navier-Stokes computational fluid dynamics methods, FUN3D, USM3D and Kestrel, were used to predict the performance buildup of components to the full X-57 configuration. The goal of the X-57 wing and distributed electric propulsion system design was to meet or exceed the required lift coefficient of 3.95 for a stall speed of 58 knots. The X-57 Maxwell airplane was designed with a small, high aspect ratio cruise wing that was designed for a high cruise lift coefficient of 0.75 at a cruise speed of 150 knots and altitude of 8,000 ft, with an angle of attack of approximately 0deg. The computational data indicates that the X-57 full aircraft drag would meet the cruise drag goal with a 25 count drag margin. The cruise configuration maximum lift coefficient is 2.07 and without including the stabilator is 1.86 at an angle of attack of 14 deg, predicted with the USM3D flow solver using the Spalart-Allmaras turbulence model. The maximum lift coefficient for the high-lift wing (with the 30deg flap deflection) without the stabilator contribution is 2.60 at an angle of attack of 13 deg. For high-lift blowing conditions with 13.7 hp/prop, the maximum lift coefficient excluding the stabilator is 4.426 at (alpha) = 13 deg. Therefore, the lift augmentation from the high-lift propellers is 1.7 and the total lift augmentation from the high-lift system (30 deg flap deflection and the high-lift blowing) is 2.38. The drag for the high-lift wing with 30 deg flap deflection is much higher than the cruise wing configuration, but the high-lift system is used only during a small portion of the flight envelope. The pitching moment is relatively constant for both blown and unblown conditions when the stabilator is excluded. Modeling the full geometry has indicated some adverse effects from the fuselage on the wing and stabilator. At high angles of attack, the solutions with the USM3D flow solver using the Spalart-Allmaras turbulence model indicates large flow separation on the wing upper surface between the two high-lift nacelles near the fuselage, and also a reduction in sectional lift on the stabilator in the first 50 percent of the stabilator semispan. However, the large flow separation near the fuselage is mostly eliminated in the solutions predicted with two codes, USM3D and Kestrel, using Hybrid Reynolds-averaged Navier Stokes/Large Eddy Simulation turbulence models.

Deere, Karen A.↗

Radiation Interaction with Particulates in Planetary Entry Flows

A hypersonic particle phase solver based on a continuum model is implemented in the NASA HyperSolve CFD code and used to predict particle phase conditions in a dusty Mars entry scenario. Black-body radiative emission from high temperature dust particles is predicted for a Mars relevant entry problem, and the emission integrated to an entry vehicle using a ray tracing approach to estimate surface heat flux. The surface heating due to this mechanism is found to be small relative to convective and radiation from the gas phase emission. Additionally, the effect of atmospheric haze particulates vaporizing in the Titan atmosphere is studied. The increase in carbon availability due to the composition of the haze particles is found to increase the gas radiative heating to the vehicle by up to 3% on the forebody of a representative Titan entry vehicle with approximately 0.15% haze in the freestream by mass.

Andrew Hinkle↗

Global Mapping of Underwater UV Irradiances and DNA-Weighted Exposures using TOMS and SeaWiFS Data Products

The global stratospheric ozone-layer depletion results In an increase in biologically harmful ultraviolet (UV) radiation reaching the surface and penetrating to ecologically significant depths in natural waters. Such an increase can be estimated on a global scale by combining satellite estimates of UV irradiance at the ocean surface from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument with the SeaWIFS satellite ocean-color measurements in the visible spectral region. In this paper we propose a model of seawater optical properties in the UV spectral region based on the Case I water model in the visible range. The inputs to the model are standard monthly SeaWiFS products: chlorophyll concentration and the diffuse attenuation coefficient at 490nm. Penetration of solar UV radiation to different depths in open ocean waters is calculated using the RT (radiative transfer) quasi-single scattering approximation (QSSA). The accuracy of the QSSA approximation in the water is tested using more accurate codes. The sensitivity study of the underwater UV irradiance to atmospheric and oceanic optical properties have shown that the main environmental parameters controlling the absolute levels of the UVB (280-320nm) and DNA-weighted irradiance underwater are: solar-zenith angle, cloud transmittance, water optical properties, and total ozone. Weekly maps of underwater UV irradiance and DNA-weighted exposure are calculated using monthly-mean SeaWiFS chlorophyll and diffuse attenuation coefficient products, daily SeaWiFS cloud fraction data, and the TOMS-derived surface UV irradiance daily maps. The final products include global maps of weekly-average UVB irradiance and DNA-weighted daily exposures at 3m and 10m, and depths where the UVB irradiance and DNA-weighted dose rate at local noon are equal to 10% of their surface values.

Vasilkov, Alexander↗