Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete surface”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Variable-Domain Displacement Transfer Functions for Converting Surface Strains into Deflections for Structural Deformed Shape Predictions

Variable-Domain Displacement Transfer Functions were formulated for shape predictions of complex wing structures, for which surface strain-sensing stations must be properly distributed to avoid jointed junctures, and must be increased in the high strain gradient region. Each embedded beam (depth-wise cross section of structure along a surface strain-sensing line) was discretized into small variable domains. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. Through discretization, the embedded beam curvature equation can be piece-wisely integrated to obtain the Variable-Domain Displacement Transfer Functions (for each embedded beam), which are expressed in terms of geometrical parameters of the embedded beam and the surface strains along the strain-sensing line. By inputting the surface strain data into the Displacement Transfer Functions, slopes and deflections along each embedded beam can be calculated for mapping out overall structural deformed shapes. A long tapered cantilever tubular beam was chosen for shape prediction analysis. The input surface strains were analytically generated from finite-element analysis. The shape prediction accuracies of the Variable- Domain Displacement Transfer Functions were then determined in light of the finite-element generated slopes and deflections, and were fofound to be comparable to the accuracies of the constant-domain Displacement Transfer Functions

displacement transfer functions↗
NASA NTRS1977

Finite element analysis of aeroacoustic jet-flap flows

A computational analysis was performed on the steady, turbulent aerodynamic flowfields associated with a jet-blown flap. For regions devoid of flow separation, a parabolic approximation to the governing time-averaged Navier-Stokes equations was applied. Numerical results are presented for the symmetry plane flow of a slot-nozzle planar jet flap geometry, including prediction of flowfield evolution within the secondary mixing region immediately downstream of the trailing edge. Using a two equation turbulence kinetic energy closure model, rapid generation and decay of large spatial gradients in mean and correlated fluctuating velocity components within the immediate wake region were predicted. Modifications to the turbulent flow structure, as induced by porous surface treatment of the flap, were evaluated. The recirculating flow within a representative discrete slot in the surface was evaluated, using the two dimensional, time-averaged Navier-Stokes equations.

Baker, A. J.↗
NASA NTRS2008

Spline-Based Smoothing of Airfoil Curvatures

Constrained fitting for airfoil curvature smoothing (CFACS) is a splinebased method of interpolating airfoil surface coordinates (and, concomitantly, airfoil thicknesses) between specified discrete design points so as to obtain smoothing of surface-curvature profiles in addition to basic smoothing of surfaces. CFACS was developed in recognition of the fact that the performance of a transonic airfoil is directly related to both the curvature profile and the smoothness of the airfoil surface. Older methods of interpolation of airfoil surfaces involve various compromises between smoothing of surfaces and exact fitting of surfaces to specified discrete design points. While some of the older methods take curvature profiles into account, they nevertheless sometimes yield unfavorable results, including curvature oscillations near end points and substantial deviations from desired leading-edge shapes. In CFACS as in most of the older methods, one seeks a compromise between smoothing and exact fitting. Unlike in the older methods, the airfoil surface is modified as little as possible from its original specified form and, instead, is smoothed in such a way that the curvature profile becomes a smooth fit of the curvature profile of the original airfoil specification. CFACS involves a combination of rigorous mathematical modeling and knowledge-based heuristics. Rigorous mathematical formulation provides assurance of removal of undesirable curvature oscillations with minimum modification of the airfoil geometry. Knowledge-based heuristics bridge the gap between theory and designers best practices. In CFACS, one of the measures of the deviation of an airfoil surface from smoothness is the sum of squares of the jumps in the third derivatives of a cubicspline interpolation of the airfoil data. This measure is incorporated into a formulation for minimizing an overall deviation- from-smoothness measure of the airfoil data within a specified fitting error tolerance. CFACS has been extensively tested on a number of supercritical airfoil data sets generated by inverse design and optimization computer programs. All of the smoothing results show that CFACS is able to generate unbiased smooth fits of curvature profiles, trading small modifications of geometry for increasing curvature smoothness by eliminating curvature oscillations and bumps (see figure).

Li, W.↗
NASA NTRS2006

An Analytical Solution of Radiative Transfer in the Coupled Atmosphere-Ocean System with Rough Surface

Using the efficient discrete-ordinate method, we present an analytical solution for radiative transfer in the coupled atmosphere-ocean system with rough air-water interface. The theoretical formulations of the radiative transfer equation and solution are described. The effects of surface roughness on radiation field in the atmosphere and ocean are studied and compared with measurements. The results show that ocean surface roughness has significant effects on the upwelling radiation in the atmosphere and the downwelling radiation in the ocean. As wind speed increases, the angular domain of sunglint broadens, the surface albedo decreases, and the transmission to ocean increases. The downward radiance field in the upper ocean is highly anisotropic, but this anisotropy decreases rapidly as surface wind increases and as depth in ocean increases. The effects of surface roughness on radiation also depend greatly on both wavelength and angle of incidence (i.e., solar elevation); these effects are significantly smaller throughout the spectrum at high sun. The model-observation discrepancies may indicate that the Cox-Munk surface roughness model is not sufficient for high wind conditions.

Jin, Zhonghai↗
NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗
NASA NTRS2013

Large-Deformation Displacement Transfer Functions for Shape Predictions of Highly Flexible Slender Aerospace Structures

Large deformation displacement transfer functions were formulated for deformed shape predictions of highly flexible slender structures like aircraft wings. In the formulation, the embedded beam (depth wise cross section of structure along the surface strain sensing line) was first evenly discretized into multiple small domains, with surface strain sensing stations located at the domain junctures. Thus, the surface strain (bending strains) variation within each domain could be expressed with linear of nonlinear function. Such piecewise approach enabled piecewise integrations of the embedded beam curvature equations [classical (Eulerian), physical (Lagrangian), and shifted curvature equations] to yield closed form slope and deflection equations in recursive forms.

Deflection equations↗
NASA NTRS2016

Radiative Transfer Simulations of the Two-Dimensional Ocean Glint Reflectance and Determination of the Sea Surface Roughness

An optimized discrete-ordinate radiative transfer model (DISORT3) with a pseudo-two-dimensional bidirectional reflectance distribution function (BRDF) is used to simulate and validate ocean glint reflectances at an infrared wavelength (1036 nm) by matching model results with a complete set of BRDF measurements obtained from the NASA cloud absorption radiometer (CAR) deployed on an aircraft. The surface roughness is then obtained through a retrieval algorithm and is used to extend the simulation into the visible spectral range where diffuse reflectance becomes important. In general, the simulated reflectances and surface roughness information are in good agreement with the measurements, and the diffuse reflectance in the visible, ignored in current glint algorithms, is shown to be important. The successful implementation of this new treatment of ocean glint reflectance and surface roughness in DISORT3 will help improve glint correction algorithms in current and future ocean color remote sensing applications.

glint↗
DOE OSTI2022

Development of A Multi-Spectral Pyrometry Sensor for High-Speed Transient Surface-Temperature Measurements in Combustion-Relevant Harsh Environments

Accurate and high-speed transient surface-temperature measurements of combustion devices including internal combustion (IC) engines, gas turbines, etc., provide validation targets and boundary conditions for computational fluid dynamics models, and are broadly relevant to technology advancements such as performance improvement and emissions reduction. Development and demonstration of a multi-infrared-channel pyrometry-based optical instrument for high-speed surface-temperature measurement is described. The measurement principle is based on multi-spectral radiation thermometry (MRT) and uses surface thermal radiation at four discrete spectral regions and a corresponding emissivity model to obtain surface temperature via non-linear least squares (NLLS) optimization. Rules of thumb for specifying the spectral regions and considerations to avoid interference with common combustion products are developed; the impact of these along with linear and non-linear MRT analysis are assessed as a function of temperature and signal-to-noise ratio. A multi-start method to determine the MRT-solution global optimum is described and demonstrated. The resulting multi-channel transient pyrometry instrument is described along with practical considerations including optical-alignment drift, matching intra-channel transient response, and solution-confidence indicators. The instrument demonstrated excellent >97% accuracy and >99% 2-sigma precision over the 400–800 °C range, with ~20 µs (50 kHz, equivalent to 0.2 cad at 2000 RPM IC-engine operation) transient response in the bench validation.

47 OTHER INSTRUMENTATION↗
NASA NTRS2019

Turbulence Model Implementation and Verification in the SENSEI CFD Code

This paper outlines the implementation and verification of the negative Spalart-Allmaras turbulence model into the SENSEI CFD code. The SA-neg turbulence model is implemented in a flexible, object-oriented framework where additional turbulence models can be easily added. In addition to outlining the new turbulence modeling framework in SENSEI, an overview of the other general improvements to SENSEI is provided. The results for four 2D test cases are compared to results from CFL3D and FUN3D to verify that the turbulence models are implemented properly. Several differences in the results from SENSEI, CFL3D, and FUN3D are identified and are attributed to differences in the implementation and discretization order of the boundary conditions as well as the order of discretization of the turbulence model. When a solid surface is located near or intersects an inflow or outflow boundary, higher order boundary conditions should be used to limit their effect on the forces on the surface. When the turbulence equations are discretized using second order spatial accuracy, the edge of the eddy viscosity profile seems to be sharper than when a first order discretization is used. However, the discretization order of the turbulence equation does not have a significant impact on output quantities of interest, such as pressure and viscous drag, for the cases studied.

Jackson, Charles W.↗
NASA NTRS2003

Finite-Rate Ablation Boundary Conditions for Carbon-Phenolic Heat-Shield

A formulation of finite-rate ablation surface boundary conditions, including oxidation, nitridation, and sublimation of carbonaceous material with pyrolysis gas injection, has been developed based on surface species mass conservation. These surface boundary conditions are discretized and integrated with a Navier-Stokes solver. This numerical procedure can predict aerothermal heating, chemical species concentration, and carbonaceous material ablation rate over the heatshield surface of re-entry space vehicles. In this study, the gas-gas and gas-surface interactions are established for air flow over a carbon-phenolic heatshield. Two finite-rate gas-surface interaction models are considered in the present study. The first model is based on the work of Park, and the second model includes the kinetics suggested by Zhluktov and Abe. Nineteen gas phase chemical reactions and four gas-surface interactions are considered in the present model. There is a total of fourteen gas phase chemical species, including five species for air and nine species for ablation products. Three test cases are studied in this paper. The first case is a graphite test model in the arc-jet stream; the second is a light weight Phenolic Impregnated Carbon Ablator at the Stardust re-entry peak heating conditions, and the third is a fully dense carbon-phenolic heatshield at the peak heating point of a proposed Mars Sample Return Earth Entry Vehicle. Predictions based on both finite-rate gas- surface interaction models are compared with those obtained using B' tables, which were created based on the chemical equilibrium assumption. Stagnation point convective heat fluxes predicted using Park's finite-rate model are far below those obtained from chemical equilibrium B' tables and Zhluktov's model. Recession predictions from Zhluktov's model are generally lower than those obtained from Park's model and chemical equilibrium B' tables. The effect of species mass diffusion on predicted ablation rate is also examined.

Chen, Y.-K.↗
NASA NTRS1973

Position from gravity

Procedures for obtaining position from surface gravity observations are reviewed and their relevance assessed in the context of the application of modern geodetic techniques to programs of Earth and ocean physics. Solutions based on the use of surface layer techniques, the discrete value approach, and the development from Green's theorem are stated in summary, the latter being extended to order e cubed in the height anomaly. The representation of the surface gravity field which is required in order that this accuracy may be achieved is discussed. Interim techniques which could be used in the absence of such a representation are also outlined.

Mather, R. S.↗
DOE OSTI2021

The Mass-loss History of the Red Hypergiant VY CMa

Imaging and spectroscopy of the knots, clumps, and extended arcs in the complex ejecta of VY CMa confirm a record of high mass-loss events over the past few hundred years. Hubble Space Telescope/Space Telescope Imaging Spectrograph spectroscopy of numerous small knots close to the star allow us to measure their radial velocities from the strong K i emission and determine their separate motions, spatial orientations, and time since ejecta. Their ages concentrate around 70, 120, 200, and 250 yr ago. A K i emission knot only 50 mas from the star ejected as recently as 1985–1995 may coincide with an H{sub 2}O maser. Comparison with VY CMa’s historic light curve from 1800 to the present shows several knots with ejection times that correspond with extended periods of variability and deep minima. The similarity of this correspondence in VY CMa with the remarkable recent dimming of Betelgeuse and an outflow of gas is apparent. The evidence for similar outflows from the surface of a more typical red supergiant suggests that discrete ejections are more common and surface or convective activity is a major source of mass loss for red supergiants.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1976

A submerged singularity method for calculating potential flow velocities at arbitrary near-field points

A discrete singularity method has been developed for calculating the potential flow around two-dimensional airfoils. The objective was to calculate velocities at any arbitrary point in the flow field, including points that approach the airfoil surface. That objective was achieved and is demonstrated here on a Joukowski airfoil. The method used combined vortices and sources ''submerged'' a small distance below the airfoil surface and incorporated a near-field subvortex technique developed earlier. When a velocity calculation point approached the airfoil surface, the number of discrete singularities effectively increased (but only locally) to keep the point just outside the error region of the submerged singularity discretization. The method could be extended to three dimensions, and should improve nonlinear methods, which calculate interference effects between multiple wings, and which include the effects of force-free trailing vortex sheets. The capability demonstrated here would extend the scope of such calculations to allow the close approach of wings and vortex sheets (or vortices).

Maskew, B.↗
NASA NTRS1979

The computation of isentropic atmospheric trajectories using a 'discrete model' formulation

An explicit technique for calculating atmospheric trajectories is presented as an alternative method to the standard implicit scheme of Danielsen (1961). The technique uses the inviscid equations of motion and the discrete model formulation derived by Greenspan (1972, 1973) to compute trajectories on isentropic surfaces, assuming adiabatic flow. The discrete model formulation is designed specifically for a Lagrangian system and objectively accounts for the geostrophic departures, local psi-tendencies, and the subsequent accelerations along the entire length of the trajectory. Application of the discrete formulation to a diagnostic case study yielded favorable results.

Petersen, R. A.↗
NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗
NASA NTRS1976

Full-coverage film cooling heat transfer study: Summary of data for normal-hole injection and 30 deg slant-hole injection

Heat transfer to a full coverage film cooled turbulent boundary layer over a flat surface was studied. The surface consisted of a discrete hole test section containing 11 rows of holes spaced 5 diameters apart in a staggered array and an instrumented recovery region. Ten diameter spacing was also studied by plugging appropriate holes. Two test sections were used, one having holes normal to the surface and the other having holes angled 30 deg to the surface in the downstream direction. Stanton number data were obtained both in the full coverage region and in the downstream recovery region for a range of blowing ratios, or mass flux ratios, from 0 to 1.3. Initial conditions at the upstream edge of the blowing region were varied from 500 to 5000 for momentum thickness Reynolds number and from 100 to 1800 for enthalpy thickness Reynolds number. The range of Reynolds numbers based on hole diameter and mainstream velocity was 6000 to 22000. Initial boundary layer thicknesses range from 0.5 to 2.0 hole diameters. Air was used as the working fluid. The data were taken for the secondary injection temperature equal to the wall temperature and also equal to the mainstream temperature. Superposition was then used to obtain Stanton number as a continuous function of the injectant temperature. The heat transfer coefficient was defined on the basis of a mainstream-to-wall temperature difference. This definition permits direct comparison of performance between film cooling and transpiration cooling.

Crawford, M. E.↗
NASA NTRS1993

Estimating snow grain size using AVIRIS data

Estimates of snow grain size for the near-surface snow layer were calculated for the Tioga Pass region and Mammoth Mountain in the Sierra Nevada, California, using an inversion technique and data collected by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The Tioga Pass and Mammoth Mountain single-band AVIRIS radiance images were atmospherically corrected to obtain surface reflectance. A discrete-ordinate model was used to calculate directional reflectance as a function of snowpack grain size for a wide range of snow grain radii. The resulting radius vs. reflectance curves were each fit using a nonlinear least squares technique which provided a means of transforming surface reflectance in each AVIRIS image to optically equivalent grain size on a per-pixel basis. The model results and grain size estimates derived from the AVIRIS data show that, for solar incidence angles between 0 and 30, the technique provides good estimates of grain size. This work provides the first quantitative estimates for grain size using data acquired from an airborne remote sensing instrument and is an important step in improving our ability to retrieve snow physical properties independent of field measurements.

Nolin, Anne W.↗
NASA NTRS2006

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗