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 919 records · Page 51

The reduction correction in North America

An inverse Poisson integral technique was used to determine a gravity field on the geoid which, when continued by analytic free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages, each stage refining the previous solution using data at progressively increasing resolution (1 x 1 deg, 5 x 5', 5/8 x 5/8') from a decreasing area of integration. Reduction corrections are computed at 5/8 x 5/8' granularity by differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5' intervals. The 1 x 1 deg averages of the reduction corrections thus obtained for 172 1 x 1 deg squares in western North America are discussed. The 1 x 1 deg mean reduction corrections are predominantly positive, varying from -3 to +15 mgal, with values in excess of 5 mgal for 26 squares. Their mean and rms values are +2.4 and 3.6 mgal respectively and they correlate well with the mean terrain corrections. The mean and rms contributions from the three stages of computation are: 1 x 1 deg stage +0.15 and 0.7 mgal; 5 x 5' stage + 1.0 and 1.6 mgal; and 5/8 x 5/8' stage +1.3 and 1.8 mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification, due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface being systematically more positive than those on the topography.

Martzen, P. D.↗

Image feature extraction using Gabor-like transform

Noisy and highly textured images were operated on with a Gabor-like transform. The results were evaluated to see if useful features could be extracted using spatio-temporal operators. The use of spatio-temporal operators allows for extraction of features containing simultaneous frequency and orientation information. This method allows important features, both specific and generic, to be extracted from images. The transformation was applied to industrial inspection imagery, in particular, a NASA space shuttle main engine (SSME) system for offline health monitoring. Preliminary results are given and discussed. Edge features were extracted from one of the test images. Because of the highly textured surface (even after scan line smoothing and median filtering), the Laplacian edge operator yields many spurious edges.

Finegan, Michael K., Jr.↗

New approaches to optimization in aerospace conceptual design

Aerospace design can be viewed as an optimization process, but conceptual studies are rarely performed using formal search algorithms. Three issues that restrict the success of automatic search are identified in this work. New approaches are introduced to address the integration of analyses and optimizers, to avoid the need for accurate gradient information and a smooth search space (required for calculus-based optimization), and to remove the restrictions imposed by fixed complexity problem formulations. (1) Optimization should be performed in a flexible environment. A quasi-procedural architecture is used to conveniently link analysis modules and automatically coordinate their execution. It efficiently controls a large-scale design tasks. (2) Genetic algorithms provide a search method for discontinuous or noisy domains. The utility of genetic optimization is demonstrated here, but parameter encodings and constraint-handling schemes must be carefully chosen to avoid premature convergence to suboptimal designs. The relationship between genetic and calculus-based methods is explored. (3) A variable-complexity genetic algorithm is created to permit flexible parameterization, so that the level of description can change during optimization. This new optimizer automatically discovers novel designs in structural and aerodynamic tasks.

Gage, Peter J.↗

Evolutionary optimization of a Genetically Refined Truss

Structural optimization is a field of research that has experienced noteworthy growth for many years. Researchers in this area have developed optimization tools to successfully design and model structures, typically minimizing mass while maintaining certain deflection and stress constraints. Numerous optimization studies have been performed to minimize mass, deflection and stress on a benchmark cantilever truss problem. Predominantly traditional optimization theory is applied to this problem. The cross-sectional area of each member is optimized to minimize the aforementioned objectives. This paper will present a structural optimization technique that has been previously applied to compliant mechanism design. This technique demonstrates a method that combines topology optimization, geometric refinement, finite element analysis, and two forms of evolutionary computation: Genetic Algorithms and Differential Evolution to successfully optimize a benchmark structural optimization problem. An non-traditional solution to the benchmark problem is presented in this paper, specifically a geometrically refined topological solution. The design process begins with an alternate control mesh formulation, multilevel geometric smoothing operation, and an elastostatic structural analysis. The design process is wrapped in an evolutionary computing optimization toolset.

Hull, Patrick V.↗

Aeroheating Test of CEV Entry Vehicle at Turbulent Conditions

An investigation of the aeroheating environment of the Project Orion Crew Entry Vehicle has been performed in the Arnold Engineering Development Center Tunnel 9. Data were measured on a approx. 3.5% scale model (0.1778m/7-inch diam.) of the vehicle using coaxial thermocouples in the Mach 8 and Mach 10 nozzles of Tunnel 9. Runs were performed at free stream Reynolds numbers of 1 106/ft to 20 10(exp 6)/ft in the Mach 10 nozzle and 8 10(exp 6)/ft to 48 10(exp 6)/ft in the Mach 8 nozzle. The test gas in Tunnel 9 is pure N2, which at these operating conditions remains un-dissociated and may be treated as a perfect gas. At these conditions, laminar, transitional, and turbulent flow was produced on the model at Mach 10, and transitional and turbulent conditions were produced on the model at Mach 8. The majority of runs were made on a clean, smooth-surface model configuration and a limited number of runs were made in which inserts with varying boundary-layer trips configurations were used to force the occurrence of transition. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins for the computational method. Data from both the wind tunnel test and the computations are presented herein. Figure 1 shows a schematic of the thermocouple locations on the model and figures 2 and 3 show a photo and schematic of the AEDC Hypervelocity Tunnel 9. Figure 4 shows a typical grid used in the computations. From the comparisons shown in figures 5 through 8 it was concluded that for perfect-gas conditions, the computations could predict either fully-laminar or full-turbulent flow to within +/-10% of the experimental data. The experimental data showed that transition began on the leeside of the heatshield at a free stream Reynolds number of 9 10(exp 6)/ft in the Mach 10 nozzle and fully-developed turbulent flow was produced at 20 10(exp 6)/ft. In the Mach 8 nozzle, transition on the leeside of the heat-shield was observed for all test conditions, and full-developed turbulent flow occurred at a free stream Reynolds number of 18 10(exp 6)/ft. On the aftbody of the vehicle no evidence of turbulence was detected at Mach 10 conditions, and at Mach 8 conditions, transition appeared to begin on the windside of the aftbody at free stream Reynolds number of 18x10(exp 6)/ft with fully-developed turbulent flow occurring only at the highest test condition of 48x10(exp 6)/ft.

Hollis, Brian R.↗

An unconditionally stable Runge-Kutta method for unsteady flows

A quasi-three-dimensional analysis was developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream surface thickness. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing was employed to accelerate the solution of the time accurate computations. The scheme is described and accuracy analyses are given. Results are shown for a supersonic through-flow fan designed for NASA Lewis. The rotor:stator blade ratio was taken as 1:1. Results are also shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Here the blade ratio is 2:3. Implicit residual smoothing was used to increase the time step limit of the unsmoothed scheme by a factor of six with negligible differences in the unsteady results. It is felt that the implicitly smoothed Runge-Kutta scheme is easily competitive with implicit schemes for unsteady flows while retaining the simplicity of an explicit scheme.

Jorgenson, Philip C. E.↗

Sensitivity of Spaceborne and Ground Radar Comparison Results to Data Analysis Methods and Constraints

With the availability of active weather radar observations from space from the Precipitation Radar (PR) on board the Tropical Rainfall Measuring Mission (TR.MM) satellite, numerous studies have been performed comparing PR reflectivity and derived rain rates to similar observations from ground-based weather radars (GR). These studies have used a variety of algorithms to compute matching PR and GR volumes for comparison. Most studies have used a fixed 3-dimensional Cartesian grid centered on the ground radar, onto which the PR and GR data are interpolated using a proprietary approach and/or commonly available GR analysis software (e.g., SPRINT, REORDER). Other studies have focused on the intersection of the PR and GR viewing geometries either explicitly or using a hybrid of the fixed grid and PR/GR common fields of view. For the Dual-Frequency Precipitation Radar (DPR) of the upcoming Global Precipitation Measurement (GPM) mission, a prototype DPR/GR comparison algorithm based on similar TRMM PR data has been developed that defines the common volumes in terms of the geometric intersection of PR and GR rays, where smoothing of the PR and GR data are minimized and no interpolation is performed. The PR and GR volume-averaged reflectivity values of each sample volume are accompanied by descriptive metadata, for attributes including the variability and maximum of the reflectivity within the sample volume, and the fraction of range gates in the sample average having reflectivity values above an adjustable detection threshold (typically taken to be 18 dBZ for the PR). Sample volumes are further characterized by rain type (Stratiform or Convective), proximity to the melting layer, underlying surface (land/water/mixed), and the time difference between the PR and GR observations. The mean reflectivity differences between the PR and GR can differ between data sets produced by the different analysis methods; and for the GPM prototype, by the type of constraints and categorization applied to the data. In this paper, we will show results comparing the 3-D gridded analysis "black box" approach to the GPM prototype volume-matching approach, using matching TRMM PR and WSR-88D ground radar data. The affects of applying data constraints and data categorizations on the volume-matched data to the results will be shown, and explanations of the differences in terms of data and analysis algorithm characteristics will be presented. Implications of the differences to the determination of PR/DPR calibration differences and use of ground radar data to evaluate the PR and DPR attenuation correction algorithms will be discussed.

Morris, Kenneth R.↗

Flux-vector splitting and Runge-Kutta methods for the Euler equations

Runge-Kutta schemes have been used as a method of solving the Euler equations exterior to an airfoil. In the past this has been coupled with central differences and an artificial vesocity in space. In this study the Runge-Kutta time-stepping scheme is coupled with an upwinded space approximation based on flux-vector splitting. Several acceleration techniques are also considered including a local time step, residual smoothing and multigrid.

Turkel, E.↗

Interpolation methods for shaped reflector analysis

The diffraction analysis of reflector surfaces which are described only at a discrete set of locations usually leads to the requirement of an interpolation to determine the surface characteristics over a continuum of locations. Two methods of interpolation, the global and the local methods, are presented. The global interpolation representation is a closed-form or series expression valid over the entire surface. The coefficients of a series expression are found by an integration of all of the raw data. Since the number of coefficients used to describe the surface is much smaller than the number of raw data points, the integration effectively provides a smoothing of the raw data. The local interpolation provides a closed-form expression for only a small area of the reflector surface. The subreflector is divided into sectors each of which has constant discretized data. Each area segment is then locally described by a two-dimensional quadratic surface. The second derivative data give the desired smoothed values.

Galindo-Israel, Victor↗

Enhancing Adhesion: Relative Merits of Different Approaches

Adhesive performance is improved mainly by manipulation of the bimaterials interface zone, which is only a few molecules thick. There are three approaches to enhancement of interfacial adhesion at the molecular level. They are 1) changing the nonchemically bonded interactions across the interface from weak ones to strong ones, 2) making the true interfacial area much larger than the simple geometric area, and 3) inducing chemical bonding between the two materials forming the interface. Our goal this summer was to question some of the built-in assumptions contained within these approaches and to determine the most promising approach, both theoretically and practically, for enhancing adhesion in NASA structures. Our computations revealed that all three of these approaches have, in theory, the potential to enhance molecular adhesion approximately ten-fold. Experiments, however, revealed that this excellent level of enhancement is not likely to be reached in practice. Each approach was found to be severely limited by practical problems. In addition, some of the built-in assumptions associated with these approaches were found to be insufficient or inadequate. The first approach, changing the nonchemically bonded interactions from weak to strong, Is an example of one containing inadequate assumptions. The extensive literature on intermolecular interactions, based on solution studies, shows that certain functional group pairs interact much more strongly than others. It has always been assumed that these data can be reliably extended to systems where only one member of the pair is in solution and the other Is contained in a solid surface. Our experiments this summer demonstrated that solution data do not adequately predict the strength of functional group interaction at the solid-liquid interface. Furthermore, the strong solvents needed to dissolve the monomers or polymers to which the functional groups of interest are attached compete successfully with the solid surface for the functional group. As a result, functional groups in solution cannot pair with the complementary groups in the solid surface, and the expected enhancement of nonchemically bonded interactions is not realized. The second approach, increasing the true interfacial area, is an example of one containing inadequate assumptions and suffering from numerous practical problems. First, practitioners have assumed that material removal, such as bead blasting or etching, increases true surface area (and therefore interfacial area) in a meaningful way. Our geometric analysis demonstrated that removal methods increase area by a factor of two at most. To increase interfacial area by an order of magnitude or more. a thin layer of high porosity must be added to the substrate surface prior to application of the adhesive phase. Consistent with this finding, we attempted to create a thin layer of rigid, highly porous glass on the surface of our smooth glass substrate by means of sol-gel technology. We were unable to surmount a wide variety of practical problems and obtained only collapsed, nonporous layers. Thus this approach, appealing in principle, would require long term development and is not promising in the near term. The third approach, inducing chemical bonding at the interface, is an example of one having neither inadequate assumptions nor insurmountable practical problems. When silicate glass is the substrate, there are only a few chemical reactions that can be successfully conducted to create these chemical bonds, and these reactions usually involve silicon-containing reagents. We compared the silazane reagents to the silane reagents and found through experiment that the silazanes react with the glass surface much more readily, and under milder conditions, than the silanes. The functional groups attached to the glass surface by silazane reactions were not able to be removed by solvent extraction, elevated temperature exposure, or mechanical action. This clearly indicates that the formation of chemical bonds at the interface is the most effective approach for enhancing molecular adhesion.

Penn, L. S.↗

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗

A multigrid method for the transonic full potential equation discretized with finite elements on an arbitrary body fitted mesh

A multigrid method for the acceleration of transonic potential flow calculations based on a Galerkin finite element approach is described. In order to allow the use of arbitrary body fitted meshes it is necessary to introduce nonuniform interpolation and residual weighting. Emphasis is put on the construction of these operators consistent with the finite element approximation, while standard successive line overrelaxation is used as a smoothing step. Substantial convergence acceleration is obtained and results are presented for different transonic flow configurations including shocks.

Deconinck, H.↗

Convergence Acceleration of the Navier-Stokes Equations Through Time-Derivative Preconditioning

Chorin's method of artificial compressibility is extended to both compressible and incompressible fluids by using physical arguments to define artificial fluid properties that make up a local preconditioning matrix. In particular, perturbation expansions are used to provide appropriate temporal derivatives for the equations of motion at both low speeds and low Reynolds numbers. These limiting forms are then combined into a single function that smoothly merges into the physical time derivatives at high speeds so that the equations are left unchanged at transonic, high Reynolds number conditions. The effectiveness of the resulting preconditioning procedures for the Navier-Stokes equations is demonstrated for a wide speed and Reynolds number ranges by means of stability results and computational solutions. Nevertheless, the preconditioned equations sometimes fail to provide a solution for applications for which the non-preconditioned equations converge. Often this is because the reduced dissipation in the preconditioned equations results in an unsteady solution while the more dissipative non-preconditioned equations result in a steady state. Problems of this type represent a computational challenge; it is important to distinguish between non-convergence of algorithms, and the non-existence of steady state solutions.

Merkle, Charles L.↗

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Tanks

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involve the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations and knowledge were mainly carried out by extensive experimental studies. A Volume-Of-Fluid (VOF) based CFD program developed at NASA MSFC was applied to extract slosh damping in a baffled tank from the first principle. First, experimental data using water with subscale smooth wall tank were used as the baseline validation. CFD simulation was demonstrated to be capable of accurately predicting natural frequency and very low damping value from the smooth wall tank at different fill levels. The damping due to a ring baffle at different liquid fill levels from barrel section and into the upper dome was then investigated to understand the slosh damping physics due to the presence of a ring baffle. Based on this study, the Root-Mean-Square error of our CFD simulation in estimating slosh damping was less than 4.8%, and the maximum error was less than 8.5%. Scalability of subscale baffled tank test using water was investigated using the validated CFD tool, and it was found that unlike the smooth wall case, slosh damping with baffle is almost independent of the working fluid and it is reasonable to apply water test data to the full scale LOX tank when the damping from baffle is dominant. On the other hand, for the smooth wall, the damping value must be scaled according to the Reynolds number. Comparison of experimental data, CFD, with the classical and modified Miles equations for upper dome was made, and the limitations of these semi-empirical equations were identified.

Yang, H. Q.↗

Precision heat forming of tetrafluoroethylene tubing

An invention that provides a method of altering the size of tetrafluoroethylene tubing which is only available in limited combination of wall thicknesses and diameter are discussed. The method includes the steps of sliding the tetrafluoroethylene tubing onto an aluminum mandrel and clamping the ends of the tubing to the mandrel by means of clamps. The tetrafluorethylene tubing and mandrel are then placed in a supporting coil which with the mandrel and tetrafluorethylene tubing are then positioned in a insulated steel pipe which is normally covered with a fiber glass insulator to smooth out temperature distribution therein. The entire structure is then placed in an event which heats the tetrafluorethylene tubing which is then shrunk by the heat to the outer dimension of the aluminum mandrel. After cooling the aluminum mandrel is removed from the newly sized tetrafluorethylene tubing by a conventional chemical milling process.

Ruiz, W. V.↗

Method for machining holes in composite materials

A method for boring well defined holes in a composite material such as graphite/epoxy is discussed. A slurry of silicon carbide powder and water is projected onto a work area of the composite material in which a hole is to be bored with a conventional drill bit. The silicon carbide powder and water slurry allow the drill bit, while experiencing only normal wear, to bore smooth, cylindrical holes in the composite material.

Daniels, Julia G.↗

The Small Angular Oscillations of Airplanes in Steady Flight

This investigation was carried out by the National Advisory Committee for Aeronautics at the request of the Army Air Service to provide data concerning the small angular oscillations of several types of airplanes in steady flight under various atmospheric conditions. The data are of use in the design of bomb sights and other aircraft instruments. The method used consisted in flying the airplane steadily in one direction for at least one minute, while recording the angle of the airplane with the sun by means of a kymograph. The results show that the oscillations differ but little for airplanes of various types, but that the condition of the atmosphere is an important factor. The average angular excursion from the mean in smooth air is 0.8 degrees in pitch, 1.4 degrees in roll, and 0.9 degrees in yaw, without special instruments to aid the pilot in holding steady conditions. In bumpy air the values given above are increased about 50 per cent. (author)

Norton, F H↗

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗