Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient 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 1,261 records · Page 70

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and complexity are studied for four nominally second-order accurate schemes: a node-centered scheme and three cell-centered schemes - a node-averaging scheme and two schemes with nearest-neighbor and adaptive compact stencils for least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Tests from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The tests of the second class are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes may degenerate on mixed grids, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to that of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping based on a distance function commonly available in practical schemes or modifying the scheme stencil to reflect the direction of strong coupling. The major conclusion is that accuracies of the node centered and the best cell-centered schemes are comparable at equivalent number of degrees of freedom.

Diskin, Boris↗

Implicit Large-Eddy Simulations of Zero-Pressure Gradient, Turbulent Boundary Layer

A set of direct simulations of zero-pressure gradient, turbulent boundary layer flows are conducted using various span widths (62-630 wall units), to document their influence on the generated turbulence. The FDL3DI code that solves compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Results are analyzed at two different Re values (500 and 1,400), and compared with spectral DNS data. They show that a minimum span width is required for the mere initiation of numerical turbulence. Narrower domains ((is) less than 100 w.u.) result in relaminarization. Wider spans ((is) greater than 600 w.u.) are required for the turbulent statistics to match reference DNS. The upper-wall boundary condition for this setup spawns marginal deviations in the mean velocity and Reynolds stress profiles, particularly in the buffer region.

Turblent boundary layer↗

Applications for Gradient Metal Alloys Fabricated Using Additive Manufacturing

Recently, additive manufacturing (AM) techniques have been developed that may shift the paradigm of traditional metal production by allowing complex net-shaped hardware to be built up layer-by-layer, rather than being machined from a billet. The AM process is ubiquitous with polymers due to their low melting temperatures, fast curing, and controllable viscosity, and 3D printers are widely available as commercial or consumer products. 3D printing with metals is inherently more complicated than with polymers due to their higher melting temperatures and reactivity with air, particularly when heated or molten. The process generally requires a high-power laser or other focused heat source, like an electron beam, for precise melting and deposition. Several promising metal AM techniques have been developed, including laser deposition (also called laser engineered net shaping or LENS® and laser deposition technology (LDT)), direct metal laser sintering (DMLS), and electron beam free-form (EBF). These machines typically use powders or wire feedstock that are melted and deposited using a laser or electron beam. Complex net-shape parts have been widely demonstrated using these (and other) AM techniques and the process appears to be a promising alternative to machining in some cases. Rather than simply competing with traditional machining for cost and time savings, the true advantage of AM involves the fabrication of hardware that cannot be produced using other techniques. This could include parts with "blind" features (like foams or trusses), parts that are difficult to machine conventionally, or parts made from materials that do not exist in bulk forms. In this work, the inventors identify that several AM techniques can be used to develop metal parts that change composition from one location in the part to another, allowing for complete control over the mechanical or physical properties. This changes the paradigm for conventional metal fabrication, which relies on an assortment of "post-processing" methods to locally alter properties (such as coating, heat treating, work hardening, shot peening, etching, anodizing, among others). Building the final part in an additive process allows for the development of an entirely new class of metals, so-called "functionally graded metals" or "gradient alloys." By carefully blending feedstock materials with different properties in an AM process, hardware can be developed with properties that cannot be obtained using other techniques but with the added benefit of the net-shaped fabrication that AM allows.

Hofmann, Douglas C.↗

The acoustics of aircraft engine-duct systems

Noise generated in aircraft engines is usually suppressed by acoustically treating the engine ducts. The optimization of this treatment requires an understanding of the transmission and attenuation of the acoustic waves. A critical review is presented of the state of the art regarding methods of determining the transmission and attenuation parameters and the effect on these parameters of (1) acoustic properties of liners, (2) the mean velocity, including uniform and shear profiles and nonparallel flow, (3) axial and transverse temperature gradients, (4) slowly and abruptly varying cross sections, and (5) finite-amplitude waves and nonlinear duct liners.

Nayfeh, A. H.↗

Computation of incompressible, three-dimensional turbulent boundary layers and comparison with experiment

Incompressible three dimensional, turbulent boundary layer (3DTBL) experiments were simulated numerically by integrating the boundary layer equations together with an algebraic eddy viscosity turbulence model. For the flow treated, the downstream portion, where the crossflow was large, was not predicted with the present computational method; the flow was significantly influenced by elliptic flow field effects. Departures from the boundary layer concept are indicated. Calculations agreed reasonably well with the mean flow development up to separation. In one experiment the normal pressure gradients were found to be neligible in regions with large skewing and allowed testing turbulence models using the boundary layer equations. The simulation of this flow compared favorably with the experimental data throughout the flow field and suggested the applicability of algebraic eddy viscosity models for 3DTBLs.

Mueller, U. R.↗

A Research Program in Active Control/aeroelasticity

Research activities concerning active control aeroelasticity are discussed. In aeroservoelastic system design of a flexible aircraft, it is often necessary to obtain specified steady state structural dynamic response and to maintain stability margins at both the plant (aircraft) input and output. The research during the present period is discussed. Steady state structural dynamic response constraints and gradients were formulated. Incorporation of the design software as an update to the PADLOCS synthesis software was examined. The stability margin improvement technique at both the plant input and output using singular value properties and constrained optimization method was validated.

Source record↗

CONDIF - A modified central-difference scheme with unconditional stability and very low numerical diffusion

This paper describes a method, called CONDIF, which retains the essential nature of the central difference scheme (CDS) but eliminates the over- and under-shoots. It modifies the CDS by introducing a controlled amount of numerical diffusion based on local gradients. The results obtained for a number of test problems show that the CONDIF is significantly more accurate than the hybrid scheme when the Peclet number is very high and the flow is at large angles to the grid. In the worst case, it yields results similar to those of the hybrid scheme. Though the CONDIF is unconditionally stable, under some conditions the rate of convergence deteriorates as the grid Peclet number is increased. Efforts are currently underway to develop a faster converging variation of CONDIF.

Runchal, Akshai K.↗

Progress in integrated analysis with adaptive unstructured meshing

Design of lightweight structures and thermal protection systems for hypersonic vehicles depend on accurate prediction of aerothermal loads, structural temperatures and their gradients, and structural deformations and stresses. Concentration is on an alternative meshing technique which generates an entirely new adaptive unstructured mesh based on the solution obtained from the earlier mesh. The technique combined with the finite element method has been shown to significantly improve the efficiency and accuracy of the fluid, thermal, and structural analyses. Current capability of the adaptive unstructured meshing technique for the integrated fluid-thermal-structural analysis is described first. The technique was extended to transient thermal analysis of structures with time-dependent adaptive meshing to capture the detailed temperature response with a minimum number of unknowns and computational cost. Both linear and higher-order finite elements are implemented to demonstrate the generality of the technique and to investigate their solution accuracy. Currently, the adaptive meshing technique is being developed for plane structures that can be modeled with membrane elements and built-up structures modeled with membrane and bending elements. The capability of the technique to these different disciplinary problems is demonstrated by several examples.

Dechaumphai, Pramote↗

Preliminary Exploration of Adaptive State Predictor Based Human Operator Modeling

Control-theoretic modeling of the human operator dynamic behavior in manual control tasks has a long and rich history. In the last two decades, there has been a renewed interest in modeling the human operator. There has also been significant work on techniques used to identify the pilot model of a given structure. The purpose of this research is to attempt to go beyond pilot identification based on collected experimental data and to develop a predictor of pilot behavior. An experiment was conducted to quantify the effects of changing aircraft dynamics on an operator s ability to track a signal in order to eventually model a pilot adapting to changing aircraft dynamics. A gradient descent estimator and a least squares estimator with exponential forgetting used these data to predict pilot stick input. The results indicate that individual pilot characteristics and vehicle dynamics did not affect the accuracy of either estimator method to estimate pilot stick input. These methods also were able to predict pilot stick input during changing aircraft dynamics and they may have the capability to detect a change in a subject due to workload, engagement, etc., or the effects of changes in vehicle dynamics on the pilot.

Trujillo, Anna C.↗

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗

Dynamic response of membranes with finite elements

A general set of finite-element equations is developed for treating the nonlinear dynamic response of membranes consisting of complex nonlinear materials. The approach proposed by Key (1974) for a most complete treatment of the finite-element method to the nonlinear dynamic analysis of solids is used in developing the nonlinear equations of motion for membranes with cubic quadrilateral finite elements. The computations are kept general with respect to material modeling by computing the deformation gradient and requiring a constitutive subroutine to supply the Cauchy stress. An example calculation of the vibration of a square membrane demonstrates the accuracy of the formulation for small deformation theory and shows the nonlinear effect of large deformations.

Benzley, S. E.↗

Nimbus 6 satellite vs rawinsonde capabilities

The paper establishes rms errors in Nimbus 6 satellite data from discrepancies between satellite and rawinsonde data, compares rms errors in Nimbus 6 with those in rawinsonde data, and examines the capability of Nimbus 6 satellite data to depict such quantities as horizontal gradients of basic variables and computed variables such as vorticity and lapse rate of temperature. The capability of Nimbus 6 data is illustrated by comparisons of rms errors in computed variables obtained through a propagation-of-error method with average and near-extreme values of the variables obtained from two of NASA's Atmospheric Variability Experiments (AVE II and AVE IV). It is shown that the rms errors in Nimbus 6 basic variables (temperature, mixing ratio, geopotential height, and wind speed) are between 2 and 12 times larger than the corresponding errors in rawinsonde data. For all variables examined, the magnitude of the rms errors are smaller than the near-extreme values of the variables determined from AVE II and AVE IV data.

Scoggins, J. R.↗

Ocean chlorophyll studies from a U-2 aircraft platform

Chlorophyll gradient maps of large ocean areas were generated from U-2/OCS data obtained over test sites in the Pacific and the Atlantic Oceans. The delineation of oceanic features using the upward radiant intensity relies on an analysis method which presupposes that radiation backscattered from the atmosphere and the ocean surface can be properly modeled by using a measurement made at 778 nm. The calculation of atmospheric radiance was performed by using a method developed by Dave. An estimation of the chlorophyll concentration is performed by properly ratioing radiances measured at 472 and 548 nm after removing the atmospheric effects. The correlation between the remotely sensed data and the in situ surface chlorophyll measurements has been validated in two sets of data. The results show that the correlation between the in situ measured chlorophyll and the derived quantity is a negative exponential function, and the correlation coefficient was calculated to be -0.965.

Kim, H. H.↗

The influence of actual and apparent geoid error on ocean analysis and prediction

The radar altimeter is the only satellite remote sensor with a proven capability for synoptically measuring an integral property of the dynamic ocean on a near global, all weather basis. Because acquisition of global, in situ ocean data with space/time resolution adequate to describe dynamically important ocean features is practically impossible, any attempt to develop a global ocean monitoring and forecasting system will rely heavily on altimetric data for initialization and updating. Maximizing useful information from the altimeter while minimizing error sources and developing methods for assimilating altimeter data into dynamical models are, therefore, vital areas for research and development. The limits imposed on ocean prediction by errors in the geoid or apparent errors associated with ground track variations near strong geopotential gradients are examined.

Thompson, J. D.↗

Development and testing of Parabolic Dish Concentrator No. 1

Parabolic Dish Concentrator No. 1 (PDC-1) is a 12-m-diameter prototype concentrator with low life-cycle costs for use with thermal-to-electric energy conversion devices. The concentrator assembly features panels made of a resin transfer molded balsa core/fiberglass sandwich with plastic reflective film as the reflective surface and a ribbed framework to hold the panels in place. The concentrator assembly tracks in azimuth and elevation on a base frame riding on a circular track. It is shown that the panels do not exhibit the proper parabolic contour. However, thermal gradients were discovered in the panels with daily temperature changes. The PDC-1 has sufficient optical quality to operate satisfactorily in a dish-electric system. The PDC-1 development provides the impetus for creating innovative optical testing methods and valuable information for use in designing and fabricating concentrators of future dish-electric systems.

Dennison, E. W.↗

Petrov-Galerkin methods for natural convection in directional solidification of binary alloys

A Petrov-Galerkin finite element method is presented for calculation of the steady, axisymmetric thermosolutal convection and interface morphology in a model for vertical Bridgman crystal growth of nondilute binary alloys. The Petrov-Galerkin method is based on the formulation for biquadratic elements developed by Heinrich and Zienkiewicz and is introduced into the calculation of the velocity, temperature and concentration fields. The algebraic system is solved simultaneously for the field variables and interface shape by Newton's method. The results of the Petrov-Galerkin method are compared critically with those of Galerkin's method using the same finite element grids. Significant improvements in accuracy are found with the Petrov-Galerkin method only when the mesh is refined and when the formulation of the residual equations is modified to account for the mixed boundary conditions that arise at the solidification interface. Calculations for alloys with stable and unstable solute gradients show the occurrence of classical flow transitions and morphological instabilities in the solidification system.

Adornato, Peter M.↗

Environmental noise characteristics of the MOD5-B (3.2 MW) wind turbine generator

Both narrow band and broad band acoustic data were obtained for the MOD5-B wind turbine for a range of wind speeds from 5.8 to 14.3 m/s; for a range of power outputs from 300 to 3100 kW; and for various azimuth angles and distances. Comparisons are made with those of other large machines and with predictions by available methods. The highest levels occur at the lower frequencies and generally decrease as the frequency increases. Low frequency rotational noise components were more intense than expected for an upwind machine and are believed to result from localized wind gradients across the rotor disk due to upwind terrain features. Predicted broad band spectra follow the general trends of the data but tend to underestimate the levels in the frequency range where the turbulent boundary layer-trailing edge interaction noise is expected to be significant.

Shepherd, Kevin P.↗

Conjugate Gradient Algorithms For Manipulator Simulation

Report discusses applicability of conjugate-gradient algorithms to computation of forward dynamics of robotic manipulators. Rapid computation of forward dynamics essential to teleoperation and other advanced robotic applications. Part of continuing effort to find algorithms meeting requirements for increased computational efficiency and speed. Method used for iterative solution of systems of linear equations.

Fijany, Amir↗