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 397 records · Page 22

Optimal Synthesis of Compliant Mechanisms using Subdivision and Commercial FEA (DETC2004-57497)

The field of distributed-compliance mechanisms has seen significant work in developing suitable topology optimization tools for their design. These optimal design tools have grown out of the techniques of structural optimization. This paper will build on the previous work in topology optimization and compliant mechanism design by proposing an alternative design space parameterization through control points and adding another step to the process, that of subdivision. The control points allow a specific design to be represented as a solid model during the optimization process. The process of subdivision creates an additional number of control points that help smooth the surface (for example a C(sup 2) continuous surface depending on the method of subdivision chosen) creating a manufacturable design free of some traditional numerical instabilities. Note that these additional control points do not add to the number of design parameters. This alternative parameterization and description as a solid model effectively and completely separates the design variables from the analysis variables during the optimization procedure. The motivation behind this work is to create an automated design tool from task definition to functional prototype created on a CNC or rapid-prototype machine. This paper will describe the proposed compliant mechanism design process and will demonstrate the procedure on several examples common in the literature.

Hull, Patrick V.↗

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗

NASA/ONERA Collaboration on Small Hovering Rotor Broadband Noise Prediction using Lattice-Boltzmann and Structured Navier-Stokes Solvers

This work compares two lattice-Boltzmann method solvers, PowerFLOW and ProLB, and two structured Navier-Stokes solvers, OVERFLOW2 and FAST, used by NASA and ONERA, respectively, for the broadband noise prediction of an ideally twisted rotor as part of Implementing Arrangement number FR-0685-0, ‘Comparing Computational Fluid Dynamics Solvers for Broadband Noise Prediction.’ Predicted results are evaluated against measured data from both smooth and rough blade sets acquired in the Small Hover Anechoic Chamber at the NASA Langley Research Center. Aerodynamic thrust predictions are seen to agree more favorably with the rough-blade measurements, whereas torque predictions agree better with the smooth-blade measurements. A tonal noise comparison shows better agreement to the measured data with the two structured Navier-Stokes solvers than with the two lattice-Boltzmann solvers, which is thought to be caused by the different geometric discretization associated with each solver paradigm. Broadband noise comparisons show that both lattice-Boltzmann method solvers trend well with the smooth-blade measurements, with the exception of an overprediction by ProLB between 4 kHz and 15 kHz. OVERFLOW2 is seen to capture the measured nondeterministic tonal content between 3 kHz and 8 kHz on a narrowband spectral basis and FAST agrees well with the rough blades on a one-third octave band basis.

Christopher S. Thurman↗

A time-accurate adaptive grid method and the numerical simulation of a shock-vortex interaction

A time accurate, general purpose, adaptive grid method is developed that is suitable for multidimensional steady and unsteady numerical simulations. The grid point movement is performed in a manner that generates smooth grids which resolve the severe solution gradients and the sharp transitions in the solution gradients. The temporal coupling of the adaptive grid and the PDE solver is performed with a grid prediction correction method that is simple to implement and ensures the time accuracy of the grid. Time accurate solutions of the 2-D Euler equations for an unsteady shock vortex interaction demonstrate the ability of the adaptive method to accurately adapt the grid to multiple solution features.

Bockelie, Michael J.↗

Performance Trades Study for Robust Airfoil Shape Optimization

From time to time, existing aircraft need to be redesigned for new missions with modified operating conditions such as required lift or cruise speed. This research is motivated by the needs of conceptual and preliminary design teams for smooth airfoil shapes that are similar to the baseline design but have improved drag performance over a range of flight conditions. The proposed modified profile optimization method (MPOM) modifies a large number of design variables to search for nonintuitive performance improvements, while avoiding off-design performance degradation. Given a good initial design, the MPOM generates fairly smooth airfoils that are better than the baseline without making drastic shape changes. Moreover, the MPOM allows users to gain valuable information by exploring performance trades over various design conditions. Four simulation cases of airfoil optimization in transonic viscous ow are included to demonstrate the usefulness of the MPOM as a performance trades study tool. Simulation results are obtained by solving fully turbulent Navier-Stokes equations and the corresponding discrete adjoint equations using an unstructured grid computational fluid dynamics code FUN2D.

Li, Wu↗

The multigrid method: Fast relaxation

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g. 4 x 2 mesh cells) to fine (e.g. 64 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with optimal iterations on the finest grid. To date the method was applied quite successfully to the solution of the transonic small-disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method has been studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g., 8 x 2 mesh cells) to fine (e.g., 128 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with, say, optimal SLOR iterations on the finest grid. The method was applied to the solution of the transonic small-disturbance equation for the velocity potential in the conservation form. Nonlifting transonic flow past a parabolic-arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗

Engine cyclic durability by analysis and material testing

The problem of calculating turbine engine component durability is addressed. Nonlinear, finite-element structural analyses, cyclic constitutive behavior models, and an advanced creep-fatigue life prediction method called strainrange partitioning were assessed for their applicability to the solution of durability problems in hot-section components of gas turbine engines. Three different component or subcomponent geometries are examined: a stress concentration in a turbine disk; a louver lip of a half-scale combustor liner; and a squealer tip of a first-stage high-pressure turbine blade. Cyclic structural analyses were performed for all three problems. The computed strain-temperature histories at the critical locations of the combustor linear and turbine blade components were imposed on smooth specimens in uniaxial, strain-controlled, thermomechanical fatigue tests of evaluate the structural and life analysis methods.

Kaufman, A.↗

Engine cyclic durability by analysis and material testing

The problem of calculating turbine engine component durability is addressed. Nonlinear, finite-element structural analyses, cyclic constitutive behavior models, and an advanced creep-fatigue life prediction method called strainrange partitioning were assessed for their applicability to the solution of durability problems in hot-section components of gas turbine engines. Three different component or subcomponent geometries are examined: a stress concentration in a turbine disk; a louver lip of a half-scale combustor linear; and a squealer tip of a first-stage high-pressure turbine blade. Cyclic structural analyses were performed for all three problems. The computed strain-temperature histories at the critical locations of the combustor linear and turbine blade components were imposed on smooth specimens in uniaxial, strain-controlled, thermomechanical fatigue tests of evaluate the structural and life analysis methods.

Kaufman, A.↗

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗

An analytical procedure for computing smooth transitions between two specified cross sections with applications to blended wing body configuration

An analytical procedure is described for designing smooth transition surfaces for blended wing-body configurations. Starting from two specified cross section shapes, the procedure generates a gradual transition from one cross section shape to the other as an analytic blend of the two shapes. The method utilizes a conformal mapping, with subsequent translation and scaling, to transform the specified and shapes to curves that can be combined more smoothly. A sample calculation is applied to a blended wing-body missile type configuration with a top mounted inlet.

Barger, R. L.↗

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.↗

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.↗