Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing methods”

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,063 records · Page 59

Toward Automatic Verification of Goal-Oriented Flow Simulations

We demonstrate the power of adaptive mesh refinement with adjoint-based error estimates in verification of simulations governed by the steady Euler equations. The flow equations are discretized using a finite volume scheme on a Cartesian mesh with cut cells at the wall boundaries. The discretization error in selected simulation outputs is estimated using the method of adjoint-weighted residuals. Practical aspects of the implementation are emphasized, particularly in the formulation of the refinement criterion and the mesh adaptation strategy. Following a thorough code verification example, we demonstrate simulation verification of two- and three-dimensional problems. These involve an airfoil performance database, a pressure signature of a body in supersonic flow and a launch abort with strong jet interactions. The results show reliable estimates and automatic control of discretization error in all simulations at an affordable computational cost. Moreover, the approach remains effective even when theoretical assumptions, e.g., steady-state and solution smoothness, are relaxed.

Simulations↗

Quasi-Lagrangian models of nascent thermals

Slab-symmetric and axisymmetric geometries are considered in the modeling of thermals by considering their boundaries as self-advecting vortex sheets, for the case of a neutral environment without turbulent processes. The initial accelerations for circular cross section thermals are obtained analytically, and a numerical method is developed to simulate the evolution of thermals by time-marching in which the vortex sheet is divided into finite segments whose positions are tracked in a quasi-Lagrangian fashion. Self-advection is considered, and a redistribution procedure is adopted to prevent the segments from bunching unmanageably. The induced field of motion is determined both inside and outside the thermal. It is shown that (1) the axisymmetric thermal rises more quickly than the slab-symmetric thermal, (2) vorticity maxima occur both within the concavities and on the trailing edges, and (3) the leading edges are very smooth.

Rambaldi, S.↗

Pair Identity and Smooth Variation Rules Applicable for the Spectroscopic Parameters of H2O Transitions Involving High-J States

Two basic rules (i.e. the pair identity and the smooth variation) applicable for H2O transitions involving high-J states have been discovered. The origins of these rules are the properties of the energy levels and wavefunctions of H2O states with the quantum number J above certain boundaries. As a result, for lines involving high-J states in individually defined groups, all their spectroscopic parameters (i.e. the transition wavenumber, intensity, pressure-broadened half-width, pressure-induced shift, and temperature exponent) must follow these rules. One can use these rules to screen spectroscopic data provided by databases and to identify possible errors. In addition, by using extrapolation methods within the individual groups, one is able to predict the spectroscopic parameters for lines in this group involving very high-J states. The latter are required in developing high-temperature molecular spectroscopic databases such as HITEMP.

Ma, Q.↗

A theory of supercritical wing sections, with computer programs and examples.

Mathematical methods for the design of supercritical wings, which depend on the numerical solution of the partial differential equations of two-dimensional gas dynamics, are developed. The main contribution is a computer program for the design of shockless transonic airfoils using the hodograph transformation and analytic continuation into the complex domain. The mathematical theory is described, and a manual for users of the programs is provided. Numerical examples are given and computational results are discussed, and the computer programs themselves are listed. The analysis routine can be used to ascertain whether the profiles behave well at off-design conditions, or to smooth coordinates and obtain a desirable shape more quickly when perfectly shockless flow is not essential.

Bauer, F.↗

Polarimeter measures sea state characteristics using emitted infrared radiation

An infrared polarimeter, capable of operating between 1 and 12 micrometers wavelength has been used to measure the polarization of emitted radiation from the sea. The observed polarization at 10.6 micrometers from a smooth sea was found to be positive, indicating the dominance of reflected infrared sky radiation over the emitted. With the appearance of waves, the percent polarization increased, as expected, for a zenith angle well above the Brewster angle for water. This is qualitatively in accordance with a model presented to explain the behavior. Initial analyses indicate that the polarized components of the sea's emitted and reflected radiation are affected by type and direction of waves, angle of viewing, and foam. The effects of variations in these parameters require further delineation. The infrared polarimetric technique appears to be a novel new passive method for remote monitoring of waves.

Egan, W. G.↗

Search for resonances in positron-atom systems

No one likes to see a scattering cross section curve that is too smooth; it is much more interesting to find bumps and wiggles and most interesting if it is possible to understand their cause. Several types of resonances were clearly established in positron-containing systems: those lying just below a degenerate threshold (like 2s-2p in hydrogenic atoms or ions) and those representing Coulomb bound states in a rearranged channel (like Ps + H reversing positron +H(-)). Recently, two new sorts of resonances have been reported for which the resonant mechanism is not clear. The first is a very low lying resonance in the e-Ps system (obtained by an adiabatic expansion method), and the second is a similarly low lying two channel resonance in the positron -H system (obtained by a close-coupling technique.) These developments encouraged the examination of such systems using the standard methods of stabilization and complex rotation. Most of the results are negative; the low lying resonances in either system are verified. Some indication of new resonances in the e(+)-He(+) system is found; this may be caused by the attraction between Ps in the n=2 state and the He(++) nucleus.

Bhatia, A. K.↗

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

Stiffener Jumping for Unitized Stitched Composite Preform Manufacturing

A single-sided stitching head is used to assemble unitized three-dimensional preforms using the NASA Langley Research Center (LaRC) Integrated Structural Assembly of Advanced Composites (ISAAC) robotic manufacturing system. Typically, assembly of a preform includes a “skin” that is stiffened by stiffener elements stitched to the skin. The single-sided stitching head has two needles that pass through the preform from the same side to form the stitch, typically from the smooth outer surface so the needles penetrate through the preform towards the stiffener side. The single-sided head cannot stitch across stiffening elements (e.g., stringer, frame, integral cap, etc.) because the needles would penetrate into the stiffener web, which would not allow the stitch form because the thread movement would be impeded. As a result, stitching seams cannot cross stiffening elements, and therefore are very short and must remain in the skin/flange connections within the space between the stiffening element webs. Each short segment requires a stop and cut process that takes a significant amount of time. A stiffener jumping stitching method was developed that eliminates the need to cut and restart seams when encountering a stiffener, thus providing a long, continuous seam with a significant reduction in stitching time by eliminating cutting operations and associated excess thread end trimming. This stiffener jumping process places what could be considered a very long stitch that bridges the stiffener to connect the end of one segment of the seam with the beginning of the next segment of the seam to provide a longer continuous seam. As a result, all intermediate thread cuts that would be done at the end of each individual segment of the seam are eliminated, with only a single cut at the end of the seam. The process and programming requirements for carrying out the stiffener jumping stitching approach are presented in this paper. Using the presented stiffener jumping approach to produce the long continuous seams, the stitching time for a representative commercial transport unitized stiffened wing cover panel can possibly be reduced by at least 25 percent.

stitched↗

The use of the AOA TTR-4P GPS receiver in operation at the BIPM for real-time restitution of GPS time

The Global Positioning System is an outstanding tool for the dissemination of time. Using mono-channel C/A-code GPS time receivers, the restitution of GPS time through the satellite constellation presents a peak-to-peak discrepancy of several tens of nanoseconds without SA but may be as high as several hundreds of nanoseconds with SA. As a consequence, civil users are more and more interested in implementing hardware and software methods for efficient restitution of GPS time, especially in the framework of the project of a real-time prediction of UTC (UTCp) which could be available in the form of time differences (UTCp - GPS time). Previous work, for improving the real-time restitution of GPS time with SA, to the level obtained without SA, focused on the implementation of a Kalman filter based on past data and updated at each new observation. An alternative solution relies upon the statistical features of the noise brought about by SA; it has already been shown that the SA noise is efficiently reduced by averaging data from numerous satellites observed simultaneously over a sufficiently long time. This method was successfully applied to data from a GPS time receiver, model AOA TTR-4P, connected to the cesium clock kept at the BIPM. This device, a multi-channel, dual frequency, P-code GPS time receiver, is one of the first TTR-4P units in operation in a civil laboratory. Preliminary comparative studies of this new equipment with conventional GPS time receivers are described in this paper. The results of an experimental restitution of GPS time, obtained in June 1993, are also detailed: 3 to 6 satellites were observed simultaneously with a sample interval of 15 s, an efficient smoothing of SA noise was realized by averaging data on all observed satellites over more than 1 hour. When the GPS system is complete in 1994, 8 satellites will be observable continuously from anywhere in the world and the same level of uncertainty will be obtained using a shorter averaging time.

Thomas, Claudine↗

NASA Tech Briefs, June 2003

Topics covered include: Nulling Infrared Radiometer for Measuring Temperature; The Ames Power Monitoring System; Hot Films on Ceramic Substrates for Measuring Skin Friction; Probe Without Moving Parts Measures Flow Angle; Detecting Conductive Liquid Leaking from Nonconductive Pipe; Adaptive Suppression of Noise in Voice Communications; High-Performance Solid-State W-Band Power Amplifiers; Microbatteries for Combinatorial Studies of Conventional Lithium-Ion Batteries; Correcting for Beam Aberrations in a Beam-Waveguide Antenna; Advanced Rainbow Solar Photovoltaic Arrays; Metal Side Reflectors for Trapping Light in QWIPs; Software for Collaborative Engineering of Launch Rockets; Software Assists in Extensive Environmental Auditing; Software Supports Distributed Operations via the Internet; Software Estimates Costs of Testing Rocket Engines; yourSky: Custom Sky-Image Mosaics via the Internet; Software for Managing Inventory of Flight Hardware; Lower-Conductivity Thermal-Barrier Coatings; Process for Smoothing an Si Substrate after Etching of SiO2; Flexible Composite-Material Pressure Vessel; Treatment to Destroy Chlorohydrocarbon Liquids in the Ground; Noncircular Cross Sections Could Enhance Mixing in Sprays; Small, Untethered, Mobile Roots for Inspecting Gas Pipes; Paint-Overspray Catcher; Preparation of Regular Specimens for Atom Probes; Inverse Tomo-Lithography for Making Microscopic 3D Parts; Predicting and Preventing Incipient Flameout in Combustors; MEMS-Based Piezoelectric/Electrostatic Inchworm Actuator; Metallized Capillaries as Probes for Raman Spectroscopy; Adaptation of Mesoscale Weather Models to Local Forecasting; Aerodynamic Design using Neural Networks; Combining Multiple Gyroscope Outputs for Increased Accuracy; and Improved Collision-Detection Method for Robotic Manipulator.

Source record↗

The composite load spectra project

Probabilistic methods and generic load models capable of simulating the load spectra that are induced in space propulsion system components are being developed. Four engine component types (the transfer ducts, the turbine blades, the liquid oxygen posts and the turbopump oxidizer discharge duct) were selected as representative hardware examples. The composite load spectra that simulate the probabilistic loads for these components are typically used as the input loads for a probabilistic structural analysis. The knowledge-based system approach used for the composite load spectra project provides an ideal environment for incremental development. The intelligent database paradigm employed in developing the expert system provides a smooth coupling between the numerical processing and the symbolic (information) processing. Large volumes of engine load information and engineering data are stored in database format and managed by a database management system. Numerical procedures for probabilistic load simulation and database management functions are controlled by rule modules. Rules were hard-wired as decision trees into rule modules to perform process control tasks. There are modules to retrieve load information and models. There are modules to select loads and models to carry out quick load calculations or make an input file for full duty-cycle time dependent load simulation. The composite load spectra load expert system implemented today is capable of performing intelligent rocket engine load spectra simulation. Further development of the expert system will provide tutorial capability for users to learn from it.

Newell, J. F.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars, and ongoing field programs such as DoE/ARM. These programs are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity.

Cahalan, Robert F.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.

Cahalan, Robert F.↗

Using Information From Prior Satellite Scans to Improve Cloud Detection Near the Day-Night Terminator

With geostationary satellite data it is possible to have a continuous record of diurnal cycles of cloud properties for a large portion of the globe. Daytime cloud property retrieval algorithms are typically superior to nighttime algorithms because daytime methods utilize measurements of reflected solar radiation. However, reflected solar radiation is difficult to accurately model for high solar zenith angles where the amount of incident radiation is small. Clear and cloudy scenes can exhibit very small differences in reflected radiation and threshold-based cloud detection methods have more difficulty setting the proper thresholds for accurate cloud detection. Because top-of-atmosphere radiances are typically more accurately modeled outside the terminator region, information from previous scans can help guide cloud detection near the terminator. This paper presents an algorithm that uses cloud fraction and clear and cloudy infrared brightness temperatures from previous satellite scan times to improve the performance of a threshold-based cloud mask near the terminator. Comparisons of daytime, nighttime, and terminator cloud fraction derived from Geostationary Operational Environmental Satellite (GOES) radiance measurements show that the algorithm greatly reduces the number of false cloud detections and smoothes the transition from the daytime to the nighttime clod detection algorithm. Comparisons with the Geoscience Laser Altimeter System (GLAS) data show that using this algorithm decreases the number of false detections by approximately 20 percentage points.

Yost, Christopher R.↗

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties↗

GIFTS SM EDU Data Processing and Algorithms

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) Sensor Module (SM) Engineering Demonstration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiances using a Fourier transform spectrometer (FTS). The GIFTS instrument employs three Focal Plane Arrays (FPAs), which gather measurements across the long-wave IR (LWIR), short/mid-wave IR (SMWIR), and visible spectral bands. The raw interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, which are further processed to obtain atmospheric profiles via retrieval algorithms. This paper describes the processing algorithms involved in the calibration stage. The calibration procedures can be subdivided into three stages. In the pre-calibration stage, a phase correction algorithm is applied to the decimated and filtered complex interferogram. The resulting imaginary part of the spectrum contains only the noise component of the uncorrected spectrum. Additional random noise reduction can be accomplished by applying a spectral smoothing routine to the phase-corrected blackbody reference spectra. In the radiometric calibration stage, we first compute the spectral responsivity based on the previous results, from which, the calibrated ambient blackbody (ABB), hot blackbody (HBB), and scene spectra can be obtained. During the post-processing stage, we estimate the noise equivalent spectral radiance (NESR) from the calibrated ABB and HBB spectra. We then implement a correction scheme that compensates for the effect of fore-optics offsets. Finally, for off-axis pixels, the FPA off-axis effects correction is performed. To estimate the performance of the entire FPA, we developed an efficient method of generating pixel performance assessments. In addition, a random pixel selection scheme is designed based on the pixel performance evaluation.

Tian, Jialin↗

Improved Discrete Approximation of Laplacian of Gaussian

An improved method of computing a discrete approximation of the Laplacian of a Gaussian convolution of an image has been devised. The primary advantage of the method is that without substantially degrading the accuracy of the end result, it reduces the amount of information that must be processed and thus reduces the amount of circuitry needed to perform the Laplacian-of- Gaussian (LOG) operation. Some background information is necessary to place the method in context. The method is intended for application to the LOG part of a process of real-time digital filtering of digitized video data that represent brightnesses in pixels in a square array. The particular filtering process of interest is one that converts pixel brightnesses to binary form, thereby reducing the amount of information that must be performed in subsequent correlation processing (e.g., correlations between images in a stereoscopic pair for determining distances or correlations between successive frames of the same image for detecting motions). The Laplacian is often included in the filtering process because it emphasizes edges and textures, while the Gaussian is often included because it smooths out noise that might not be consistent between left and right images or between successive frames of the same image.

Shuler, Robert L., Jr.↗

Comparison of central axis and jet ring coolant supply for turbine disk cooling on a SSME-HPOTP model

The test facility, test methods and results are presented for an experimental study modeling the cooling of turbine disks in the blade attachment regions with multiple impinging jets, in a configuration simulating the disk cooling method employed on the Space Shuttle Main Engine oxygen turbopump. The study's objective was to provide a comparison of detailed local convection heat transfer rates obtained for a single center-supply of disk coolant with those obtained with the present flight configuration where disk coolant is supplied through an array of 19 jets located near the disk outer radius. Specially constructed disk models were used in a program designed to evaluate possible benefits and identify any possible detrimental effects involved in employing an alternate disk cooling scheme. The study involved the design, construction and testing of two full scale rotating model disks, one plane and smooth for baseline testing and the second contoured to the present flight configuration, together with the corresponding plane and contoured stator disks. Local heat transfer rates are determined from the color display of encapsulated liquid crystals coated on the disk in conjunction with use of a computer vision system. The test program was composed of a wide variety of disk speeds, flowrates, and geometrical configurations, including testing for the effects of disk boltheads and gas ingestion from the gas path region radially outboard of the disk-cavity.

Kim, Y. W.↗