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 955 records · Page 53

Solution of the surface Euler equations for accurate three-dimensional boundary-layer analysis of aerodynamic configurations

The three-dimensional boundary-layer equations in the limit as the normal coordinate tends to infinity are called the surface Euler equations. The present paper describes an accurate method for generating edge conditions for three-dimensional boundary-layer codes using these equations. The inviscid pressure distribution is first interpolated to the boundary-layer grid. The surface Euler equations are then solved with this pressure field and a prescribed set of initial and boundary conditions to yield the velocities along the two surface coordinate directions. Results for typical wing and fuselage geometries are presented. The smoothness and accuracy of the edge conditions obtained are found to be superior to the conventional interpolation procedures.

Iyer, V.↗

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, Dimitri J.↗

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, D. J.↗

CFD-Predicted Tile Heating Bump Factors Due to Tile Overlay Repairs

A Computational Fluid Dynamics investigation of the Orbiter's Tile Overlay Repair (TOR) is performed to assess the aeroheating Damage Assessment Team's (DAT) existing heating correlation method for protuberance interference heating on the surrounding thermal protection system. Aerothermodynamic heating analyses are performed for TORs at the design reference damage locations body points 1800 and 1075 for a Mach 17.9 and a=39deg STS-107 flight trajectory point with laminar flow. Six different cases are considered. The computed peak heating bump factor on the surrounding tiles are below the DAT's heating bump factor values for smooth tile cases. However, for the uneven tiles cases the peak interference heating is shown to be considerably higher than the existing correlation prediction.

Lessard, Victor R.↗

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗

Eye movements during afterimage tracking under sinusoidal and random vestibular stimulation

The smooth portion of the horizontal vestibulo-ocular reflex was analyzed in terms of the frequency response, relating slow-phase eye velocity to angular velocity of a rotating chair under four different cases: sinusoidal rotation about a vertical axis in total darkness, sinusoidal rotation during afterimage tracking, pseudorandom head rotation in total darkness, and pseudorandom head rotation during afterimage tracking. Eye movements were recorded using a photoelectric limbus tracking method. The observation that the presence of an afterimage during vestibular stimulation increases the velocity of slow-phase eye movements is in support of the theory that such slow-phase movements are generated, at least in part, by the perceived velocity of the target. Since the target is immobilized on the retina, this perceived velocity is clearly not generated by retinal slip, but rather by a mechanism related to the eye movement such as corollary discharge.

Yasui, S.↗

Calculating climate attractor dimension from delta O-18 records by the Grassberger-Procaccia algorithm

The Grassberger-Procaccia method of calculating dimension from a time series is applied to 14 late Pleistocene delta O-18 records. A step-by-step sequence leading from data to the Grassberger-Procaccia dimension is outlined, and the problems encountered when dealing with observed (as opposed to theoretical) data are discussed; for the climatic proxy data these problems include situations where the time series is not very long, is noisy and/or smoothed, and is not sampled at a constant time interval. The delta O-18 records to be used are described, and the results are presented and compared with previously published dimension calculations. New dimension interpretations are assessed, and an example using a synthetic time series that illustrates the possible error due to inconsistencies in the time scale is analyzed.

Maasch, Kirk A.↗

A report on SHARP (Spacecraft Health Automated Reasoning Prototype) and the Voyager Neptune encounter

The development and application of the Spacecraft Health Automated Reasoning Prototype (SHARP) for the operations of the telecommunications systems and link analysis functions in Voyager mission operations are presented. An overview is provided of the design and functional description of the SHARP system as it was applied to Voyager. Some of the current problems and motivations for automation in real-time mission operations are discussed, as are the specific solutions that SHARP provides. The application of SHARP to Voyager telecommunications had the goal of being a proof-of-capability demonstration of artificial intelligence as applied to the problem of real-time monitoring functions in planetary mission operations. AS part of achieving this central goal, the SHARP application effort was also required to address the issue of the design of an appropriate software system architecture for a ground-based, highly automated spacecraft monitoring system for mission operations, including methods for: (1) embedding a knowledge-based expert system for fault detection, isolation, and recovery within this architecture; (2) acquiring, managing, and fusing the multiple sources of information used by operations personnel; and (3) providing information-rich displays to human operators who need to exercise the capabilities of the automated system. In this regard, SHARP has provided an excellent example of how advanced artificial intelligence techniques can be smoothly integrated with a variety of conventionally programmed software modules, as well as guidance and solutions for many questions about automation in mission operations.

Martin, R. G.↗

Effect of radiative transfer on convection in the deep photosphere of late-type dwarfs

A method is proposed to eliminate the compressional instability of a shallow layer in the upper part of stellar convective zones in standard mixing-length models. By equating the radiative cooling time of mixing eddies to their convective turnover time, the effective sizes of the eddies are assumed to be the smallest of those which are not eliminated by radiative transfer. Computations of the models with this assumption leads to smooth temperature profiles in the previously unstable layers and reductions of the convective velocity above its maximum value.

Fontenla, J. M.↗

Investigation of upwind, multigrid, multiblock numerical schemes for three dimensional flows. Volume 1: Runge-Kutta methods for a thin layer Navier-Stokes solver

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure.

Cannizzaro, Frank E.↗

Direct numerical simulation of incompressible axisymmetric flows

In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Loulou, Patrick↗

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

Propellant slosh is a potential source of disturbance critical to the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring-mass-damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. A Volume-Of-Fluid (VOF) based Computational Fluid Dynamics (CFD) program developed at MSFC was applied to extract slosh damping in the baffled tank from the first principle. First the experimental data using water with sub-scale smooth wall tank were used as the baseline validation. It is demonstrated that CFD can indeed accurately predict low damping values from the smooth wall at different fill levels. The damping due to a ring baffles at different depths from the free surface was then simulated, and fairly good agreement with experimental measurement was observed. Comparison with an empirical correlation of Miles equation is also made.

Yang, H. Q.↗

Technique for recovery of voice data from heat damaged magnetic tape

A method for conditioning, and thus enabling retrieval of intelligence from, magnetic tapes after damage from heat has caused the tape to wrinkle and curl severely thereby reducing tape width to less than one-half its original size. The damaged tape is superposed on a first piece of splicing tape with the oxide side of the magnetic tape in contact with the adhesive side of the splicing tape and then carefully smoothed by a special tool. A second piece of splicing tape is placed on the backing side of the magnetic tape then the resulting tape stack is trimmed to the original width of the magnetic tape. After the first piece of splicing tape is carefully removed from the oxide side of the damaged magnetic tape, the resulting magnetic tape is then ready to be placed into a recorder for playback.

Melugin, J. F.↗

Silicon ribbon growth by a capillary action shaping technique

The crystal growth method described is a capillary action shaping technique. Meniscus shaping for the desired ribbon geometry occurs at the vertex of a wettable die. As ribbon growth depletes the melt meniscus, capillary action supplies replacement material. A capillary die is so designed that the bounding edges of the die top are not parallel or concentric with the growing ribbon. The new dies allow a higher melt meniscus with concomitant improvements in surface smoothness and freedom from SiC surface particles, which can degrade perfection.

Schwuttke, G. H.↗

Improved vortex methods for three-dimensional flows

Robust numerical methods are developed for three-dimensional incompressible vortical flows, using Lagrangian vortex elements. A successful scheme must be able to handle regions of intense vortex stretching and vortex reconnection with reasonable accuracy (without diverging). Here, consideration is given to vortex particles, also commonly called vortons or vortex sticks. The following issues are discussed: (1) use of delta-function elements and weak solutions of the vorticity equation; (2) use of smoothed elements and the choice of the smoothing function; (3) representation of viscous effects and the redistribution of element strength; and (4) conservation laws (are they satisfied?). The various proposed schemes have been tested on flows involving a strong interaction between two vortex rings.

Winckelmans, G.↗

Pressure Drop in Radiator Air Tubes

This report describes a method for measuring the drop in static pressure of air flowing through a radiator and shows (1) a reason for the discrepancy noted by various observers between head resistance and drop in pressure; (2) a difference in degree of contraction of the jet in entering a circular cell and a square cell; (3) the ratio of internal frictional resistance to total head resistance for two representative types; (4) the effect of smoothness of surface on pressure gradient; and (5) the effects of supplying heat to the radiator on pressure gradient. The fact that the pressure gradients are found to be approximately proportional to the square of the rate of flow of air appears to indicate turbulent flow, even in the short tubes of the radiator. It was found that the drop in the static pressure in the air stream through a cellular radiator and the pressure gradient in the air tubes are practically proportional to the square of the air flow in a given air density; that the difference between the head resistance per unit area and the fall of static pressure through the air tubes in radiators is apparent rather than real; and that radiators of different types differ widely in the amount of contraction of the jet at entrance. The frictional resistance was found to vary considerably, and in one case to be two-thirds of the head resistance in the type using circular cells and one-half of the head resistance of the radiator type using square cells of approximately the same dimensions.

Parsons, S R↗

Effects of rotation on coolant passage heat transfer. Volume 2: Coolant passages with trips normal and skewed to the flow

An experimental program was conducted to investigate heat transfer and pressure loss characteristics of rotating multipass passages, for configurations and dimensions typical of modem turbine blades. This experimental program is one part of the NASA Hot Section Technology (HOST) Initiative, which has as its overall objective the development and verification of improved analysis methods that will form the basis for a design system that will produce turbine components with improved durability. The objective of this program was the generation of a data base of heat transfer and pressure loss data required to develop heat transfer correlations and to assess computational fluid dynamic techniques for rotating coolant passages. The experimental work was broken down into two phases. Phase 1 consists of experiments conducted in a smooth wall large scale heat transfer model. A detailed discussion of these results was presented in volume 1 of a NASA Report. In Phase 2 the large scale model was modified to investigate the effects of skewed and normal passage turbulators. The results of Phase 2 along with comparison to Phase 1 is the subject of this Volume 2 NASA Report.

Johnson, B. V.↗

Comparison of a two-dimensional shock impingement computation with experiment

Results of computations of two-dimensional viscous blunt-body flowfields with an impinging shock wave, with a time-dependent finite-difference method employed to solve the complete set of Navier-Stokes equations, are compared with experimental results. The experimental results were obtained in a 20-inch hypersonic tunnel with a planar shock impinging on the cylindrical leading edge of a fin, hence with the shock parallel to the centerline of the leading edge, so that type III and type IV interference patterns were generated. Close agreement is found. The overall effects of smoothing and grid size on the calculations are determined. A 31 x 51 mesh is adequate for wall pressure values (except in peaked regions).

Tannehill, J. C.↗