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 37 records · Page 2

Simulation of Supersonic Jet Noise with the Adaptation of Overflow CFD Code and Kirchhoff Surface Integral

An acoustic prediction capability for supersonic axisymmetric jets was developed on the basis of OVERFLOW Navier-Stokes CFD (Computational Fluid Dynamics) code of NASA Langley Research Center. Reynolds-averaged turbulent stresses in the flow field are modeled with the aid of Spalart-Allmaras one-equation turbulence model. Appropriate acoustic and outflow boundary conditions were implemented to compute time-dependent acoustic pressure in the nonlinear source-field. Based on the specification of acoustic pressure, its temporal and normal derivatives on the Kirchhoff surface, the near-field and the far-field sound pressure levels are computed via Kirchhoff surface integral, with the Kirchhoff surface chosen to enclose the nonlinear sound source region described by the CFD code. The methods are validated by a comparison of the predictions of sound pressure levels with the available data for an axisymmetric turbulent supersonic (Mach 2) perfectly expanded jet.

Kandula, Max↗

Simulation of Jet Noise with OVERFLOW CFD Code and Kirchhoff Surface Integral

An acoustic prediction capability for supersonic axisymmetric jets was developed on the basis of OVERFLOW Navier-Stokes CFD (Computational Fluid Dynamics) code of NASA Langley Research Center. Reynolds-averaged turbulent stresses in the flow field are modeled with the aid of Spalart-Allmaras one-equation turbulence model. Appropriate acoustic and outflow boundary conditions were implemented to compute time-dependent acoustic pressure in the nonlinear source-field. Based on the specification of acoustic pressure, its temporal and normal derivatives on the Kirchhoff surface, the near-field and the far-field sound pressure levels are computed via Kirchhoff surface integral, with the Kirchhoff surface chosen to enclose the nonlinear sound source region described by the CFD code. The methods are validated by a comparison of the predictions of sound pressure levels with the available data for an axisymmetric turbulent supersonic (Mach 2) perfectly expanded jet.

Kandula, M.↗

Surface acoustic wave coding for orthogonal frequency coded devices

Methods and systems for coding SAW OFC devices to mitigate code collisions in a wireless multi-tag system. Each device producing plural stepped frequencies as an OFC signal with a chip offset delay to increase code diversity. A method for assigning a different OCF to each device includes using a matrix based on the number of OFCs needed and the number chips per code, populating each matrix cell with OFC chip, and assigning the codes from the matrix to the devices. The asynchronous passive multi-tag system includes plural surface acoustic wave devices each producing a different OFC signal having the same number of chips and including a chip offset time delay, an algorithm for assigning OFCs to each device, and a transceiver to transmit an interrogation signal and receive OFC signals in response with minimal code collisions during transmission.

Malocha, Donald↗

ASHMET: A computer code for estimating insolation incident on tilted surfaces

A computer code, ASHMET, was developed by MSFC to estimate the amount of solar insolation incident on the surfaces of solar collectors. Both tracking and fixed-position collectors were included. Climatological data for 248 U. S. locations are built into the code. The basic methodology used by ASHMET is the ASHRAE clear-day insolation relationships modified by a clearness index derived from SOLMET-measured solar radiation data to a horizontal surface.

Elkin, R. F.↗

Computational procedure for three-dimensional boundary layers on aircraft and aerospace vehicles

An implicit numerical procedure is developed to solve the full three-dimensional laminar, boundary layer equations over realistic aircraft and spacecraft configurations. This procedure utilizes splined functions in the direction normal to the surface and finite differences in the cross flow direction. Calculation of off-pitch plane attachment and detachment lines can be calculated. The solution procedure permits the incorporation of a wide variety of turbulent models and boundary conditions including; entropy layer procedures, specified surface temperatures or heat flux distributions, surface catalysis, mass addition, and the coupling to both three-dimensional inviscid codes and surface ablation codes. In addition, reasonable computing times are obtained.

Kendall, R. M.↗

Capillary Pressure of a Liquid Between Uniform Spheres Arranged in a Square-Packed Layer

The capillary pressure in the pores defined by equidimensional close-packed spheres is analyzed numerically. In the absence of gravity the menisci shapes are constructed using Surface Evolver code. This permits calculation the free surface mean curvature and hence the capillary pressure. The dependences of capillary pressure on the liquid volume constructed here for a set of contact angles allow one to determine the evolution of basic capillary characteristics under quasi-static infiltration and drainage. The maximum pressure difference between liquid and gas required for a meniscus passing through a pore is calculated and compared with that for hexagonal packing and with approximate solution given by Mason and Morrow [l]. The lower and upper critical liquid volumes that determine the stability limits for the equilibrium capillary liquid in contact with square packed array of spheres are tabulated for a set of contact angles.

Alexader, J. Iwan D.↗

Optical Performance Modeling of FUSE Telescope Mirror

We describe the Metrology Data Processor (METDAT), the Optical Surface Analysis Code (OSAC), and their application to the image evaluation of the Far Ultraviolet Spectroscopic Explorer (FUSE) mirrors. The FUSE instrument - designed and developed by the Johns Hopkins University and launched in June 1999 is an astrophysics satellite which provides high resolution spectra (lambda/Delta(lambda) = 20,000 - 25,000) in the wavelength region from 90.5 to 118.7 nm The FUSE instrument is comprised of four co-aligned, normal incidence, off-axis parabolic mirrors, four Rowland circle spectrograph channels with holographic gratings, and delay line microchannel plate detectors. The OSAC code provides a comprehensive analysis of optical system performance, including the effects of optical surface misalignments, low spatial frequency deformations described by discrete polynomial terms, mid- and high-spatial frequency deformations (surface roughness), and diffraction due to the finite size of the aperture. Both normal incidence (traditionally infrared, visible, and near ultraviolet mirror systems) and grazing incidence (x-ray mirror systems) systems can be analyzed. The code also properly accounts for reflectance losses on the mirror surfaces. Low frequency surface errors are described in OSAC by using Zernike polynomials for normal incidence mirrors and Legendre-Fourier polynomials for grazing incidence mirrors. The scatter analysis of the mirror is based on scalar scatter theory. The program accepts simple autocovariance (ACV) function models or power spectral density (PSD) models derived from mirror surface metrology data as input to the scatter calculation. The end product of the program is a user-defined pixel array containing the system Point Spread Function (PSF). The METDAT routine is used in conjunction with the OSAC program. This code reads in laboratory metrology data in a normalized format. The code then fits the data using Zernike polynomials for normal incidence systems or Legendre-Fourier polynomials for grazing incidence systems. It removes low order terms from the metrology data, calculates statistical ACV or PSD functions, and fits these data to OSAC models for the scatter analysis. In this paper we briefly describe the laboratory image testing of FUSE spare mirror performed in the near and vacuum ultraviolet at John Hopkins University and OSAC modeling of the test setup performed at NASA/GSFC. The test setup is a double-pass configuration consisting of a Hg discharge source, the FUSE off-axis parabolic mirror under test, an autocollimating flat mirror, and a tomographic imaging detector. Two additional, small fold flats are used in the optical train to accommodate the light source and the detector. The modeling is based on Zernike fitting and PSD analysis of surface metrology data measured by both the mirror vendor (Tinsley) and JHU. The results of our models agree well with the laboratory imaging data, thus validating our theoretical model. Finally, we predict the imaging performance of FUSE mirrors in their flight configuration at far-ultraviolet wavelengths.

Saha, Timo T.↗

Extending a transonic small disturbance code to treat swept vertical surfaces

A flexible-swept vertical surface capability has been developed and implemented within the CAP-TSD transonic small disturbance (TSD) code. The new capability required a modification to the TSD equation and a grid transformation for swept vertical surfaces. Modifications to the vertical surface boundary conditions allow it to be treated as a flexible surface. The new capability extends the range of problems which the code can treat. In order to assess the accuracy of the modifications, calculations were performed for a rectangular T-tail configuration and an AGARD T-tail configuration. Unsteady forces and moments are presented for the rectangular T-tail oscillating in yaw for a range of reduced frequencies. Comparisons are presented with linear theory and experiment. Steady and unsteady surface pressures are presented for the AGARD T-tail along with generalized aerodynamic forces. Comparisons are made with linear theory. The comparisons demonstrate the accuracy of the vertical surface modifications.

Gibbons, Michael D.↗

Computer programs to predict induced effects of jets exhausting into a crossflow

A user's manual for two computer programs was developed to predict the induced effects of jets exhausting into a crossflow. Program JETPLT predicts pressures induced on an infinite flat plate by a jet exhausting at angles to the plate and Program JETBOD, in conjunction with a panel code, predicts pressures induced on a body of revolution by a jet exhausting normal to the surface. Both codes use a potential model of the jet and adjacent surface with empirical corrections for the viscous or nonpotential effects. This program manual contains a description of the use of both programs, instructions for preparation of input, descriptions of the output, limitations of the codes, and sample cases. In addition, procedures to extend both codes to include additional empirical correlations are described.

Perkins, S. C., Jr.↗

NASA Charging Analyzer Program - A computer tool that can evaluate electrostatic contamination

The data from the ATS-5 and -6 Auroral Particles Experiments have shown that the environment at geosynchronous orbit can charge spacecraft surfaces to appreciable negative values. Such surface charges could have undesirable effects on satellite equipment and the data obtained in satellite experiments. A computer program, the NASA Charging Analyzer Program (NASCAP), has, in this connection, been developed to determine the surface charging of spacecraft when encountering geomagnetic substorm conditions. Once a model of a spacecraft has been developed in the NASCAP code, the surface charging can be computed for any environmental input flux. On the basis of the obtained information regarding the surface charging, charged particle trajectories can be computed to determine if it is possible for a particle to return to a satellite surface. It is pointed out that the computer runs conducted with the NASCAP code represent the first attempts to use a self-consistent program to study the behavior of spacecraft surfaces experiencing a geomagnetic substorm condition.

Stevens, N. J.↗

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

Wing Weight Optimization Under Aeroelastic Loads Subject to Stress Constraints

A minimum weight optimization of the wing under aeroelastic loads subject to stress constraints is carried out. The loads for the optimization are based on aeroelastic trim. The design variables are the thickness of the wing skins and planform variables. The composite plate structural model incorporates first-order shear deformation theory, the wing deflections are expressed using Chebyshev polynomials and a Rayleigh-Ritz procedure is adopted for the structural formulation. The aerodynamic pressures provided by the aerodynamic code at a discrete number of grid points is represented as a bilinear distribution on the composite plate code to solve for the deflections and stresses in the wing. The lifting-surface aerodynamic code FAST is presently being used to generate the pressure distribution over the wing. The envisioned ENSAERO/Plate is an aeroelastic analysis code which combines ENSAERO version 3.0 (for analysis of wing-body configurations) with the composite plate code.

Kapania, Rakesh K.↗

Enhancement of surface definition and gridding in the EAGLE code

Algorithms for smoothing of curves and surfaces for the EAGLE grid generation program are presented. The method uses an existing automated technique which detects undesirable geometric characteristics by using a local fairness criterion. The geometry entity is then smoothed by repeated removal and insertion of spline knots in the vicinity of the geometric irregularity. The smoothing algorithm is formulated for use with curves in Beta spline form and tensor product B-spline surfaces.

Thompson, Joe F.↗

Overview of numerical codes developed for predicted electrothermal deicing of aircraft blades

An overview of the deicing computer codes that have been developed at the University of Toledo under sponsorship of the NASA-Lewis Research Center is presented. These codes simulate the transient heat conduction and phase change occurring in an electrothermal deicier pad that has an arbitrary accreted ice shape on its surface. The codes are one-dimensional rectangular, two-dimensional rectangular, and two-dimensional with a coordinate transformation to model the true blade geometry. All modifications relating to the thermal physics of the deicing problem that have been incorporated into the codes will be discussed. Recent results of reformulating the codes using different numerical methods to increase program efficiency are described. In particular, this reformulation has enabled a more comprehensive two-dimensional code to run in much less CPU time than the original version. The code predictions are compared with experimental data obtained in the NASA-Lewis Icing Research Tunnel with a UH1H blade fitted with a B. F. Goodrich electrothermal deicer pad. Both continuous and cyclic heater firing cases are considered. The major objective in this comparison is to illustrate which codes give acceptable results in different regions of the airfoil for different heater firing sequences.

Keith, Theo G.↗

The role of under-determined approximations in engineering and science application

There is currently a great deal of interest in using response surfaces in the optimization of aircraft performance. The objective function and/or constraint equations involved in these optimization problems may come from numerous disciplines such as structures, aerodynamics, environmental engineering, etc. In each of these disciplines, the mathematical complexity of the governing equations usually dictates that numerical results be obtained from large computer programs such as a finite element method program. Thus, when performing optimization studies, response surfaces are a convenient way of transferring information from the various disciplines to the optimization algorithm as opposed to bringing all the sundry computer programs together in a massive computer code. Response surfaces offer another advantage in the optimization of aircraft structures. A characteristic of these types of optimization problems is that evaluation of the objective function and response equations (referred to as a functional evaluation) can be very expensive in a computational sense. Because of the computational expense in obtaining functional evaluations, the present study was undertaken to investigate under-determinined approximations. An under-determined approximation is one in which there are fewer training pairs (pieces of information about a function) than there are undetermined parameters (coefficients or weights) associated with the approximation. Both polynomial approximations and neural net approximations were examined. Three main example problems were investigated: (1) a function of one design variable was considered; (2) a function of two design variables was considered; and (3) a 35 bar truss with 4 design variables was considered.

Carpenter, William C.↗

The design/analysis of flows through turbomachinery - A viscous/inviscid approach

A new design/analysis system for the flows through turbomachinery is currently being developed for studying turbomachinery problems with an axisymmetric viscous/inviscid 'average-passage' throughflow code. The advantage of this approach, compared to streamline curvature codes, is that the solutions obtained simulate some of the unsteadiness, compressibility and viscous effects of a multistage turbomachine. The design/analysis system consists of three elemental parts, the axisymmetric block grid generator, the blade surface element code, and the axisymmetric flow code. Each element of the system will be discussed and the flow solutions for three axisymmetric geometries will be shown compared to experimental data where available. The computations are shown to be in very good agreement with test data for SR7 spinner body and transonic boattail geometry obtained in the wind tunnels at NASA Lewis Research Center. The VIADAC Rotor 67 Fan results were compared to PARC2D calculated results and shown to be in very good agreement.

Miller, D. P.↗

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗