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 631 records · Page 35

Aperture amplitude and phase control of offset dual reflectors

An approximate synthesis method is presented for solving the problem of offset dual-shaped reflectors. It is shown that the equations obtained by this method are partial differential equations that can be integrated numerically, except that these equations are generally not total and therefore do not have a 'smooth' solution. However, it is demonstrated that these equations form a nearly total differential that enables one to generate a smooth subreflector by integrating the differential equations and then synthesizing a main reflector which gives an exact solution for the specified aperture phase distribution. Several important illustrative examples are treated, which validate the adequacy of the proposed method.

Galindo-Israel, V.↗

Super Smooth Optics for Extra-Solar Planet Detection

The goal of imaging planets around the nearby stars has important scientific significance but requires the use of advanced methods of controlling diffracted and scattered light. Over the last three years we have undertaken a study of coronagraphic methods of controlling diffracted light and of figuring hyper-contrast optics. Progress in these two general areas have led to a proposed space-based, 1.9 meter diameter coronagraphic telescope designed specifically for very high performance in the imaging of faint objects near bright sources. This instrument, called the Circumstellar Imaging Telescope (CIT), relies on a new high efficiency coronagraph design and the careful control of scattered light by extremely smooth optics. The high efficiency coronagraph uses focal plane apodization in order to concentrate diffracted light more efficiently in the pupil. This allows convenient removal of the diffracted light by masking off parts of the telescope pupil while not sacrificing the center of the field. Reductions of diffracted light by factors exceeding 1000 are not only possible but are required in order to detect extra-solar planets. Laboratory experiments with this new design have confirmed the theoretical diffraction reductions to the limits of the optics used (factors of about 300) . The extremely high efficiency of this coronagraph puts strong constraints on the narrow angle scattered light due to figure errors in the telescope mirror. Since planets orbiting nearby stars are expected at angular distances of about 1 arcsecond, it is in this small angular range in which scattering must be controlled. The figure errors responsible for scattering in this range come from mid-spatial frequencies corresponding to correlation lengths of about 10 cm on the primary mirror. A primary mirror about 15 times smoother than the Hubble Space Telescope mirror is required for the CIT. Laboratory experiments indicate that small test mirrors can be fabricated with existing technology which come within a factor of two of this requirement.

Terrile, Richard J.↗

High-order accurate finite difference discretisations on fully unstructured dual quadrilateral meshes

Here, we present a novel approach for high-order accurate numerical differentiation on unstructured meshes of quadrilateral elements. To differentiate a given function, an auxiliary function with greater smoothness properties is defined which when differentiated provides the derivatives of the original function. The method generalises traditional finite difference methods to meshes of arbitrary topology in any number of dimensions for any order of derivative and accuracy. We demonstrate the accuracy of the numerical scheme using dual quadrilateral meshes and a refinement method based on subdivision surfaces. The scheme is applied to the solution of a range of partial differential equations, including both linear and nonlinear, second and fourth order equations, and a time-dependent first order equation.

97 MATHEMATICS AND COMPUTING↗

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗

Quick-and-easy shear-load testing

Device for applying shear loads to test specimens can be attached and removed without damaging specimen surface. Because it is quick, clean, and inexpensive, method is expected to be useful in commercial testing laboratories for applying shear loads to smooth surfaces on which there is no provision for conventional attachments and which could be easily damaged.

Gustafson, J. A.↗

Stellar convection 3: Convection at large Rayleigh numbers

A three dimensional study of convection in a self gravitating sphere of Boussinesq fluid with a Rayleigh number of 10 to the 10th power and a Prandtl of 1 is presented. The velocity and temperature of the fluid are computed at the largest wavelengths using spectral methods. A confirmation that the fluid is anisotropic and that the energy spectra are not smooth functions of wavelength but have a large amount of fine structure is discussed. The parameterization of the transport properties of the unresolvable inertial subrange with eddy viscosities and diffusivities is described. The time dependent fluctuations in the energy spectra and how they cascade from large to small wavelengths is examined.

Marcus, P. S.↗

Remote sensing of soil moisture content over bare field at 1.4 GHz frequency

An algorithm for estimating moisture content of a bare soil from the observed brightness temperature at 1.4 GHz is discussed and applied to a limited data base. The method is based on a radiative transfer model calculation, which has been successfully used in the past to account for many observational results, with some modifications to take into account the effect of surface roughness. Besides the measured brightness temperatures, the three additional inputs required by the method are the effective soil thermodynamic temperature, the precise relation between moisture content and the smooth field brightness temperatures and a pair of parameters related to surface roughness. The procedures of estimating surface roughness parameters and of obtaining moisture content from observed brightness temperature are discussed. The algorithm is applied to observations from truck mounted and airborne radiometers. The estimated moisture contents compare favorably with the observations in the top 2 cm layer.

Wang, J. R.↗

Numerical generation of composite three dimensional grids by quasilinear elliptic systems

A three-dimensional grid generation technique designed to numerically construct a boundary-conforming grid within a three-dimensional region bounded by a closed surface is described. The boundary values are generated numerically by a modified elliptic system and are used to compute grid control parameters that are contained in the elliptic systems. The interior grid distribution is governed by the distribution of points on the boundary as well as by the boundary's geometric shape. The composite three-dimensional grid remains both continuous and smooth across the surface of juncture between any two adjoining subregions. The details of the method and its implementation are presented, displaying numerical results for three-dimensional composite grid for a wing-body combination and surface grids.

Thomas, P. D.↗

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without loss of derivatives.

Tadmor, E.↗

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without 'loss of derivatives'.

Tadmor, E.↗

Monte Carlo calculation of dynamical properties of the two-dimensional Hubbard model

A new method is introduced for analytically continuing imaginary-time data from quantum Monte Carlo calculations to the real-frequency axis. The method is based on a least-squares-fitting procedure with constraints of positivity and smoothness on the real-frequency quantities. Results are shown for the single-particle spectral-weight function and density of states for the half-filled, two-dimensional Hubbard model.

White, S. R.↗

High order filtering methods for approximating hyperbolic systems of conservation laws

The essentially nonoscillatory (ENO) schemes, while potentially useful in the computation of discontinuous solutions of hyperbolic conservation-law systems, are computationally costly relative to simple central-difference methods. A filtering technique is presented which employs central differencing of arbitrarily high-order accuracy except where a local test detects the presence of spurious oscillations and calls upon the full ENO apparatus to remove them. A factor-of-three speedup is thus obtained over the full-ENO method for a wide range of problems, with high-order accuracy in regions of smooth flow.

Lafon, F.↗

Data reduction using cubic rational B-splines

A geometric method is proposed for fitting rational cubic B-spline curves to data that represent smooth curves including intersection or silhouette lines. The algorithm is based on the convex hull and the variation diminishing properties of Bezier/B-spline curves. The algorithm has the following structure: it tries to fit one Bezier segment to the entire data set and if it is impossible it subdivides the data set and reconsiders the subset. After accepting the subset the algorithm tries to find the longest run of points within a tolerance and then approximates this set with a Bezier cubic segment. The algorithm uses this procedure repeatedly to the rest of the data points until all points are fitted. It is concluded that the algorithm delivers fitting curves which approximate the data with high accuracy even in cases with large tolerances.

Chou, Jin J.↗

A three dimensional multigrid multiblock multistage time stepping scheme for the Navier-Stokes equations

A general multiblock method for the solution of the three-dimensional, unsteady, compressible, thin-layer Navier-Stokes equations has been developed. The convective and pressure terms are spatially discretized using Roe's flux differencing technique while the viscous terms are centrally differenced. An explicit Runge-Kutta method is used to advance the solution in time. Local time stepping, adaptive implicit residual smoothing, and the Full Approximation Storage (FAS) multigrid scheme are added to the explicit time stepping scheme to accelerate convergence to steady state. Results for three-dimensional test cases are presented and discussed.

Elmiligui, Alaa↗

Quadrature-Free Implementation of the Discontinuous Galerkin Method for Hyperbolic Equations

A discontinuous Galerkin formulation that avoids the use of discrete quadrature formulas is described and applied to linear and nonlinear test problems in one and two space dimensions. This approach requires less computational time and storage than conventional implementations but preserves the compactness and robustness inherent to the discontinuous Galerkin method. Test problems include both linear and nonlinear one-dimensional scalar advection of botH smooth and discontinuous initial value problems, two-dimensional scalar advection of smooth initial value problems discretized by using unstructured grids with varying degrees of smoothness and regularity, and two-dimensional linear Euler solutions on unstructured grids.

Atkins, H. L.↗

Method for Ultrasonic Imaging and Device for Performing the Method

A method for ultrasonic imaging of interior structures and flaws in a test specimen with a smooth or irregular contact surfaces, in which an ultrasonic transducer is coupled acoustically to the contact surface via a plurality of ultrasonic wave guides with equal delay times. The wave guides are thin and bendable, so they adapt to variations in the distance between the transducer and different parts of the contact surface by bending more or less. All parts of the irregular contact surface accordingly receive sound waves that are in phase, even when the contact surface is irregular, so a coherent sound wave is infused in the test specimen. The wave guides can be arranged in the form of an ultrasonic brush, with a flat head for coupling to a flat transducer, and free bristles that can be pressed against the test specimen. By bevelling the bristle ends at a suitable angle, shear mode waves can be infused into the test specimen from a longitudinal mode transducer.

Madaras, Eric I.↗

Automated Approach to Very High-Order Aeroacoustic Computations

Computational aeroacoustics requires efficient, high-resolution simulation tools. For smooth problems, this is best accomplished with very high-order in space and time methods on small stencils. However, the complexity of highly accurate numerical methods can inhibit their practical application, especially in irregular geometries. This complexity is reduced by using a special form of Hermite divided-difference spatial interpolation on Cartesian grids, and a Cauchy-Kowalewski recursion procedure for time advancement. In addition, a stencil constraint tree reduces the complexity of interpolating grid points that am located near wall boundaries. These procedures are used to develop automatically and to implement very high-order methods (> 15) for solving the linearized Euler equations that can achieve less than one grid point per wavelength resolution away from boundaries by including spatial derivatives of the primitive variables at each grid point. The accuracy of stable surface treatments is currently limited to 11th order for grid aligned boundaries and to 2nd order for irregular boundaries.

Dyson, Rodger W.↗