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 235 records · Page 13

Numerical approximations to nonlinear conservation laws with locally varying time and space grids

Numerical approximations to the initial value problem for nonlinear systems of conservation laws are considered. The considered system is said to be hyperbolic when all eigenvalues of every real linear combination of the Jacobian matrices are real. Solutions may develop discontinuities in finite time, even when the initial data are smooth. In the investigation, explicit finite difference methods which use locally varying time grids are considered. The global CFL restriction is replaced by a local restriction. The numerical flux function is studied from a finite volume viewpoint, and a differencing technique is developed at interface points between regions of distinct time increments.

Osher, S.↗

Impact cratering calculations

In the course of carrying out finite difference calculations, it was discovered that for large craters, a previously unrecognized type of crater (diameter) growth occurred which was called lip wave propagation. This type of growth is illustrated for an impact of a 1000 km (2a) silicate bolide at 12 km/sec (U) onto a silicate half-space at earth gravity (1 g). The von Misses crustal strength is 2.4 kbar. The motion at the crater lip associated with this wave type phenomena is up, outward, and then down, similar to the particle motion of a surface wave. It is shown that the crater diameter has grown d/a of approximately 25 to d/a of approximately 4 via lip propagation from Ut/a = 5.56 to 17.0 during the time when rebound occurs. A new code is being used to study partitioning of energy and momentum and cratering efficiency with self gravity for finite-sized objects rather than the previously discussed planetary half-space problems. These are important and fundamental subjects which can be addressed with smoothed particle hydrodynamic (SPH) codes. The SPH method was used to model various problems in astrophysics and planetary physics. The initial work demonstrates that the energy budget for normal and oblique impacts are distinctly different than earlier calculations for silicate projectile impact on a silicate half space. Motivated by the first striking radar images of Venus obtained by Magellan, the effect of the atmosphere on impact cratering was studied. In order the further quantify the processes of meteor break-up and trajectory scattering upon break-up, the reentry physics of meteors striking Venus' atmosphere versus that of the Earth were studied.

Ahrens, Thomas J.↗

Calculation of compressible boundary layer flow about airfoils by a finite element/finite difference method

Preliminary results are presented of a finite element/finite difference method (semidiscrete Galerkin method) used to calculate compressible boundary layer flow about airfoils, in which the group finite element scheme is applied to the Dorodnitsyn formulation of the boundary layer equations. The semidiscrete Galerkin (SDG) method promises to be fast, accurate and computationally efficient. The SDG method can also be applied to any smoothly connected airfoil shape without modification and possesses the potential capability of calculating boundary layer solutions beyond flow separation. Results are presented for low speed laminar flow past a circular cylinder and past a NACA 0012 airfoil at zero angle of attack at a Mach number of 0.5. Also shown are results for compressible flow past a flat plate for a Mach number range of 0 to 10 and results for incompressible turbulent flow past a flat plate. All numerical solutions assume an attached boundary layer.

Strong, Stuart L.↗

Effects of magnetometer calibration and maneuvers on accuracies of magnetometer-only attitude-and-rate determination

Two different algorithms - a deterministic magnetic-field-only algorithm and a Kalman filter for gyroless spacecraft - are used to estimate the attitude and rates of the Rossi X-Ray Timing Explorer (RXTE) using only measurements from a three-axis magnetometer. The performance of these algorithms is examined using in-flight data from various scenarios. In particular, significant enhancements in accuracies are observed when' the telemetered magnetometer data are accurately calibrated using a recently developed calibration algorithm. Interesting features observed in these studies of the inertial-pointing RXTE include a remarkable sensitivity of the filter to the numerical values of the noise parameters and relatively long convergence time spans. By analogy, the accuracy of the deterministic scheme is noticeably lower as a result of reduced rates of change of the body-fixed geomagnetic field. Preliminary results show the filter-per-axis attitude accuracies ranging between 0.1 and 0.5 deg and rate accuracies between 0.001 deg/sec and 0.005 deg./sec, whereas the deterministic method needs a more sophisticated techniques for smoothing time derivatives of the measured geomagnetic field to clearly distinguish both attitude and rate solutions from the numerical noise. Also included is a new theoretical development in the deterministic algorithm: the transformation of a transcendental equation in the original theory into an 8th-order polynomial equation. It is shown that this 8th-order polynomial reduces to quadratic equations in the two limiting cases-infinitely high wheel momentum, and constant rates-discussed in previous publications.

Challa, M.↗

The European Russia Drought Atlas (1400-2016 CE)

We present the European Russia Drought Atlas (ERDA) that covers the East European Plain to the Ural Mountains from 1400–2016 CE. Like the Old World Drought Atlas (OWDA) for the Euro-Mediterranean region, the ERDA is a one-half degree gridded reconstruction of summer Palmer Drought Severity Indices estimated from a network of annual tree-ring chronologies. Ensemble point-by-point regression is used to generate the ERDA with the identical protocols used for developing the OWDA. Split calibration/validation tests of the ERDA indicate that it has significant skill over most of its domain and is much more skillful than the OWDA where they overlap in the western part of ERDA domain. Comparisons to historical droughts over European Russia additionally support the ERDA’s overall validity. The ERDA has been spatially smoothed and infilled using a local regression method to yield a spatially complete drought atlas back to 1400 CE. EOF analysis indicates that there are three principal modes of hydroclimatic variability in the ERDA. After Varimax rotation, these modes correlate significantly with independent climate data sets extending back to the late nineteenth century in a physically interpretable way and relate to atmospheric circulation dynamics of droughts and heatwaves over European Russia based on more recent instrumental data.

Cook, Edward R.↗

Uniform hydrogen fuel layers for inertial fusion targets by microgravity

A critical concern in the fabrication of targets for inertial confinement fusion (ICF) is ensuring that the hydrogenic (D(sub 2) or DT) fuel layer maintains spherical symmetry. Solid layered targets have structural integrity, but lack the needed surface smoothness. Liquid targets are inherently smooth, but suffer from gravitationally induced sagging. One method to reduce the effective gravitational field environment is freefall insertion into the target chamber. Another method to counterbalance field gravitational force is to use an applied magnetic field combined with a gradient field to induce a magnetic dipole force on the liquid fuel layer. Based on time dependent calculations of the dynamics of the liquid fuel layer in microgravity environments, we show that it may be possible to produce a liquid layered ICF target that satisfies both smoothness and symmetry requirements.

Parks, P. B.↗

Aerodynamic design of turbomachinery blading in three-dimensional flow - An application to radial inflow turbines

A computational method, based on a theory for turbomachinery blading design in three-dimensional inviscid flow, is applied to a parametric design study of a radial inflow turbine wheel. As the method requires the specification of swirl distribution, a technique for its smooth generation within the blade region is proposed. Excellent agreements have been obtained between the computed results from this design method and those from direct Euler computations, demonstrating the correspondence and consistency between the two. The computed results indicate the sensitivity of the pressure distribution to a lean in the stacking axis and a minor alteration in the hub/shroud profiles. Analysis based on Navier-Stokes solver shows no breakdown of flow within the designed blade passage and agreement with that from design calculation; thus the flow in the designed turbine rotor closely approximates that of an inviscid one. These calculations illustrate the use of a design method coupled to an analysis tool for establishing guidelines and criteria for designing turbomachinery blading.

Yang, Y. L.↗

On the practical use of high-order methods for hyperbolic systems

The paper tests a number of high order methods on a variety of dynamic problems in one, two, and three space dimensions. The problems covered include wave propagation phenomena as well as an asymptotic approach to a steady state. Consideration is given to both smooth and shocked flows. It is shown that the methods compared require only minor modifications of many existing second-order schemes. Further, the results show that significant gains can be expected from the use of fourth-order methods. Finally, spectral methods are also considered for some of the problems presented.

Turkel, E.↗

Global Geochemical Variation on the Lunar Surface: A Three-Element Approach

We present a method for displaying the relative abundances of three important elements (Th, Fe, and Ti) on the same map projection of the lunar surface. Using Th-, Fe-, and Ti-elemental abundances from orbital geochemical data and assigning each element a primary color, a false-color map of the lunar surface was created. This approach is similar to the ternary diagram approach presented by Davis and Spudis with some important differences, discussed later. For the present maps, Th abundances were measured by the Lunar Prospector (LP) Gamma-Ray Spectrometer(GRS).The new LPGRS low-altitude dataset was used in this analysis. Iron and Ti weight percentages were based on Clementine spectral reflectance data smoothed to the LP low altitude footprint. This method of presentation was designed to aid in the location and recognition of three principal lunar compositions: ferroan anorthosite (FAN), mare basalts (MB), and the Mg suite/ KREEP-rich rocks on the lunar surface, with special emphasis on the highlands and specific impact basins. In addition to the recognition of these endmember rock compositions, this method is an attempt to examine the relationship between elemental compositions that do not conform readily to previously accepted or observed endmember rocks in various specific regions of interest, including eastern highlands regions centered on 150 deg longitude, and a northern highlands Th-rich region observed. The LP low-altitude data has full width at half-maximum spatial resolution of about 40 km. The Clementine spectral reflectance datasets were adapted using an equal-area, gaussian smoothing routine to this footprint. In addition, these datasets, reported in weight percent of FeO and of Ti02, were adjusted to Fe and Ti weight percentages. Each dataset was then assigned one of the three primary colors: blue for Th, red for Fe, and green for Ti. For each element, the data range was normalized to represent the ratio of each point to the maximum in the dataset. (To view the color table, go to http://cass.jsc.nasa.gov/meetings/moon99/pdf/8033.pdf.) The full range of lunar longitudes is represented, but due to the lack of coverage of the Clementine data for latitudes > 70 deg and <-70 deg, the data for these regions is excluded. The differences between this approach and the ternary diagram approach of Davis and Spudis eliminate some of the uncertainty and ambiguity of the ternary diagram approach. Rather than using a ratio of Th to Ti normalized to CI chondritic ratios, and a ternary diagram with ternary apexes located at specific endmember compositional values, elemental compositions were used independently, eliminating the errors resulting from dividing numbers that can have high uncertainties, especially at low concentration. The three elements used in this method of presentation were chosen for several reasons. One reason for the inclusion of Th in this study is that it is an accurate indicator of KREEP. Iron and Ti concentrations are both low in highland regolith, causing any small fluctuations in Th to stand out very well. In addition, Fe and Ti are good compositional indicators of different mare basalts. Mixed with red for Fe, the green for Ti produces a yellow signal in high-Ti basalts. While universally high in Fe relative to the surrounding highlands, mare basalts have a diverse range of Ti values, making Ti concentration a valuable asset to the classification and identification of different basalt types. Finally, an important constraint in element selection is the availability of the global data, both from LP and Clementine results. Data for Th, Fe, and Ti are among the highest quality of existing lunar remote-sensing data. In addition, LP data for Fe and Ti will become available, enabling these data to be incorporated into the analysis. Using upper-limit values for end member rock compositions calculated from Korotev et al., attempts were made to locate the different endmember compositions of terranes on this diagram. Most strikingly, ferroan anorthosite (Th < and = 0.37 micro g/g; Fe (wt%)< and =2.29; Ti (wt%) < and = 0.22), which should appear as an almost black, reddish color, does not appear on the diagram at any noticeable frequency. Based on this analysis, the suggestion of extensive FAN regions on the lunar surface is not strong, especially at the presently accepted values for Fe and Th. However, to make sure this effect is not due to systematic errors, a thorough investigation of the precision, accuracy, and uncertainties of the Fe, Ti, and Th abundances needs to be carried out, especially at low concentrations. A particular region of interest is an area of high Th concentrations relative to Fe and Ti content north and east of Humboldtianum Crater. First observed by Lawrence et al., this region does not coincide with any visible impact structure and comprises one of the closest approximations to pure blue (high Th, very low Ti and Fe) on the lunar surface. Such an elemental composition does not lend itself readily to classification, and presents something of an anomaly. More detailed analysis of this region is needed to understand its structure and origin. There seems to be a longitudinal asymmetry in the Th concentrations of the highlands regolith. High-Th, low-Ti, and Fe regions are located between 135 deg and 180 deg longitude and between -30 deg and +30 deg latitude. While the Th levels are not high enough to attract attention in a single elemental display, the variation in the abundance of Th relative to Fe and Ti abundances can be clearly seen. The composition that these data suggest is not well represented in the sample return suite. In addition, these regions were largely missed by the Apollo orbital ground tracks, which only covered the outer edge of the areas of interest. The LP orbital Th data represent the first information about the Th concentrations in these regions of the highlands. Additional information contained in original.

Thomsen, D. R.↗

High order filtering methods for approximating hyberbolic systems of conservation laws

In the computation of discontinuous solutions of hyperbolic systems of conservation laws, the recently developed essentially non-oscillatory (ENO) schemes appear to be very useful. However, they are computationally costly compared to simple central difference methods. A filtering method which is developed uses simple central differencing of arbitrarily high order accuracy, except when a novel local test indicates the development of spurious oscillations. At these points, the full ENO apparatus is used, maintaining the high order of accuracy, but removing spurious oscillations. Numerical results indicate the success of the method. High order of accuracy was obtained in regions of smooth flow without spurious oscillations for a wide range of problems and a significant speed up of generally a factor of almost three over the full ENO method.

Lafon, F.↗

Design of optimally smoothing multi-stage schemes for the Euler equations

In this paper, a method is developed for designing multi-stage schemes that give optimal damping of high-frequencies for a given spatial-differencing operator. The objective of the method is to design schemes that combine well with multi-grid acceleration. The schemes are tested on a nonlinear scalar equation, and compared to Runge-Kutta schemes with the maximum stable time-step. The optimally smoothing schemes perform better than the Runge-Kutta schemes, even on a single grid. The analysis is extended to the Euler equations in one space-dimension by use of 'characteristic time-stepping', which preconditions the equations, removing stiffness due to variations among characteristic speeds. Convergence rates independent of the number of cells in the finest grid are achieved for transonic flow with and without a shock. Characteristic time-stepping is shown to be preferable to local time-stepping, although use of the optimally damping schemes appears to enhance the performance of local time-stepping. The extension of the analysis to the two-dimensional Euler equations is hampered by the lack of a model for characteristic time-stepping in two dimensions. Some results for local time-stepping are presented.

Van Leer, Bram↗

Multigrid calculation of three-dimensional turbomachinery flows

Research was performed in the general area of computational aerodynamics, with particular emphasis on the development of efficient techniques for the solution of the Euler and Navier-Stokes equations for transonic flows through the complex blade passages associated with turbomachines. In particular, multigrid methods were developed, using both explicit and implicit time-stepping schemes as smoothing algorithms. The specific accomplishments of the research have included: (1) the development of an explicit multigrid method to solve the Euler equations for three-dimensional turbomachinery flows based upon the multigrid implementation of Jameson's explicit Runge-Kutta scheme (Jameson 1983); (2) the development of an implicit multigrid scheme for the three-dimensional Euler equations based upon lower-upper factorization; (3) the development of a multigrid scheme using a diagonalized alternating direction implicit (ADI) algorithm; (4) the extension of the diagonalized ADI multigrid method to solve the Euler equations of inviscid flow for three-dimensional turbomachinery flows; and also (5) the extension of the diagonalized ADI multigrid scheme to solve the Reynolds-averaged Navier-Stokes equations for two-dimensional turbomachinery flows.

Caughey, David A.↗

Experience in grid optimization

Two optimization methods for solving a variational problem in grid generation are described and evaluated. The smoothness, cell volumes, and orthogonality of the variational integrals are examined. The Jacobi-Newton iterative method is compared to the Fletcher-Reeves conjugate gradient method. It is observed that a combination of the Jacobi-Newton iteration and the direct solution of the variational problem produces an algorithm which is easy to program and requires less storage and computer time/iteration than the conjugate gradient method.

Mastin, C. W.↗

Computing Viscoplastic Behavior Of A Material

Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.

Arya, V. K.↗

Use of a genetic algorithm for the analysis of eye movements from the linear vestibulo-ocular reflex

It is common in vestibular and oculomotor testing to use a single-frequency (sine) or combination of frequencies [sum-of-sines (SOS)] stimulus for head or target motion. The resulting eye movements typically contain a smooth tracking component, which follows the stimulus, in which are interspersed rapid eye movements (saccades or fast phases). The parameters of the smooth tracking--the amplitude and phase of each component frequency--are of interest; many methods have been devised that attempt to identify and remove the fast eye movements from the smooth. We describe a new approach to this problem, tailored to both single-frequency and sum-of-sines stimulation of the human linear vestibulo-ocular reflex. An approximate derivative is used to identify fast movements, which are then omitted from further analysis. The remaining points form a series of smooth tracking segments. A genetic algorithm is used to fit these segments together to form a smooth (but disconnected) wave form, by iteratively removing biases due to the missing fast phases. A genetic algorithm is an iterative optimization procedure; it provides a basis for extending this approach to more complex stimulus-response situations. In the SOS case, the genetic algorithm estimates the amplitude and phase values of the component frequencies as well as removing biases.

Non-NASA Center↗

Two-Point Turbulence Closure Applied to Variable Resolution Modeling

Variable resolution methods have become frontline CFD tools, but in order to take full advantage of this promising new technology, more formal theoretical development is desirable. Two general classes of variable resolution methods can be identified: hybrid or zonal methods in which RANS and LES models are solved in different flow regions, and bridging or seamless models which interpolate smoothly between RANS and LES. This paper considers the formulation of bridging methods using methods of two-point closure theory. The fundamental problem is to derive a subgrid two-equation model. We compare and reconcile two different approaches to this goal: the Partially Integrated Transport Model, and the Partially Averaged Navier-Stokes method.

Girimaji, Sharath S.↗

Airplane Mesh Development with Grid Density Studies

Automatic Grid Generation Wish List Geometry handling, including CAD clean up and mesh generation, remains a major bottleneck in the application of CFD methods. There is a pressing need for greater automation in several aspects of the geometry preparation in order to reduce set up time and eliminate user intervention as much as possible. Starting from the CAD representation of a configuration, there may be holes or overlapping surfaces which require an intensive effort to establish cleanly abutting surface patches, and collections of many patches may need to be combined for more efficient use of the geometrical representation. Obtaining an accurate and suitable body conforming grid with an adequate distribution of points throughout the flow-field, for the flow conditions of interest, is often the most time consuming task for complex CFD applications. There is a need for a clean unambiguous definition of the CAD geometry. Ideally this would be carried out automatically by smart CAD clean up software. One could also define a standard piece-wise smooth surface representation suitable for use by computational methods and then create software to translate between the various CAD descriptions and the standard representation. Surface meshing remains a time consuming, user intensive procedure. There is a need for automated surface meshing, requiring only minimal user intervention to define the overall density of mesh points. The surface mesher should produce well shaped elements (triangles or quadrilaterals) whose size is determined initially according to the surface curvature with a minimum size for flat pieces, and later refined by the user in other regions if necessary. Present techniques for volume meshing all require some degree of user intervention. There is a need for fully automated and reliable volume mesh generation. In addition, it should be possible to create both surface and volume meshes that meet guaranteed measures of mesh quality (e.g. minimum and maximum angle, stretching ratios, etc.).

Cliff, Susan E.↗

Automatic cloud tracking applied to GOES and Meteosat observations

An improved automatic processing method for the tracking of cloud motions as revealed by satellite imagery is presented and applications of the method to GOES observations of Hurricane Eloise and Meteosat water vapor and infrared data are presented. The method is shown to involve steps of picture smoothing, target selection and the calculation of cloud motion vectors by the matching of a group at a given time with its best likeness at a later time, or by a cross-correlation computation. Cloud motion computations can be made in as many as four separate layers simultaneously. For data of 4 and 8 km resolution in the eye of Hurricane Eloise, the automatic system is found to provide results comparable in accuracy and coverage to those obtained by NASA analysts using the Atmospheric and Oceanographic Information Processing System, with results obtained by the pattern recognition and cross correlation computations differing by only fractions of a pixel. For Meteosat water vapor data from the tropics and midlatitudes, the automatic motion computations are found to be reliable only in areas where the water vapor fields contained small-scale structure, although excellent results are obtained using Meteosat IR data in the same regions. The automatic method thus appears to be competitive in accuracy and coverage with motion determination by human analysts.

Endlich, R. M.↗