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,207 records · Page 67

Direct simulation of a turbulent boundary layer up to R sub theta = 1410

The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectral method with up to about 10 to the 7th grid points. Periodic spanwise and streamwise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Spalart, Philippe R.↗

Theoretical studies in support of the 3M-vapor transport (PVTOS-) experiments

Results are reported for a preliminary theoretical study of the coupled mass-, momentum-, and heat-transfer conditions expected within small ampoules used to grow oriented organic solid (OS-) films, by physical vapor transport (PVT) in microgravity environments. It is show that previous studies made restrictive assumptions (e.g., smallness of delta T/T, equality of molecular diffusivities) not valid under PVTOS conditions, whereas the important phenomena of sidewall gas creep, Soret transport of the organic vapor, and large vapor phase supersaturations associated with the large prevailing temperature gradients were not previously considered. Rational estimates are made of the molecular transport properties relevant to copper-phthalocyanine monomeric vapor in a gas mixture containing H2(g) and Xe(g). Efficient numerical methods have been developed and are outlined/illustrated here to making steady axisymmetric gas flow calculations within such ampoules, allowing for realistic realistic delta T/T(sub)w-values, and even corrections to Navier-Stokes-Fourier 'closure' for the governing continuum differential equations. High priority follow-on studies are outlined based on these new results.

Rosner, Daniel E.↗

Strain energy release rate analysis of delamination in a tapered laminate subjected to tension load

A tapered composite laminate subjected to tension load was analyzed using the finite-element method. The glass/epoxy laminate has a (+ or - 45)sub 3 group of plies dropped in three distinct steps, each 20 ply-thicknesses apart, thus forming a taper angle of 5.71 degrees. Steep gradients of interlaminar normal and shear stress on a potential delamination interface suggest the existence of stress singularities at the points of material and geometric discontinuities created by the internal plydrops. The delamination was assumed to initiate at the thin end of the taper on a -45/+45 interface and the delamination growth was simulated in both directions, i.e., along the taper and into the thin region. The strain-energy-release rate for a delamination growing into the thin laminate consisted predominantly of mode I (opening) component. For a delamination growing along the tapered region, the strain-energy-release rate was initially all mode I, but the proportion of mode I decreased with increase in delamination size until eventually total G was all mode II. The total G for both delamination tips increased with increase in delamination size, indicating that a delamination initiating at the end of the taper will grow unstably along the taper and into the thin laminate simultaneously.

Salpekar, S. A.↗

Model based estimation of image depth and displacement

Passive depth and displacement map determinations have become an important part of computer vision processing. Applications that make use of this type of information include autonomous navigation, robotic assembly, image sequence compression, structure identification, and 3-D motion estimation. With the reliance of such systems on visual image characteristics, a need to overcome image degradations, such as random image-capture noise, motion, and quantization effects, is clearly necessary. Many depth and displacement estimation algorithms also introduce additional distortions due to the gradient operations performed on the noisy intensity images. These degradations can limit the accuracy and reliability of the displacement or depth information extracted from such sequences. Recognizing the previously stated conditions, a new method to model and estimate a restored depth or displacement field is presented. Once a model has been established, the field can be filtered using currently established multidimensional algorithms. In particular, the reduced order model Kalman filter (ROMKF), which has been shown to be an effective tool in the reduction of image intensity distortions, was applied to the computed displacement fields. Results of the application of this model show significant improvements on the restored field. Previous attempts at restoring the depth or displacement fields assumed homogeneous characteristics which resulted in the smoothing of discontinuities. In these situations, edges were lost. An adaptive model parameter selection method is provided that maintains sharp edge boundaries in the restored field. This has been successfully applied to images representative of robotic scenarios. In order to accommodate image sequences, the standard 2-D ROMKF model is extended into 3-D by the incorporation of a deterministic component based on previously restored fields. The inclusion of past depth and displacement fields allows a means of incorporating the temporal information into the restoration process. A summary on the conditions that indicate which type of filtering should be applied to a field is provided.

Damour, Kevin T.↗

Hypersonic stability and transition

An evaluation is conducted of the development status of growing normal modes in the computation of flow fields over hypersonic vehicles, giving attention to the effects of surface cooling, pressure gradients, and suction on the growth rates of the normal modes. Also noted are the effects of bluntness on the stability properties of hypersonic vehicles, and the efficacy of the e exp N method of transition estimation at supersonic and hypersonic speeds. A survey is made of the precautions that must be taken to ensure the effectiveness of experimental and computational studies of hypersonic flow phenomena.

Reshotko, Eli↗

Study of natural convection and interface shape in directional solidification of succinonitrile

Flow visualization experiments in a Bridgman furnace at zero growth velocity have been performed. The experiments were intended to investigate if the photochromic dye activation method could be used in a crystal growth study. The results from this work have confirmed that with a carefully designed experimental setup, the photochromic-dye method permits qualitative and quantitative evaluation of the flow field near the interface in crystal growth experiments. For the horizontal orientation of a 6 mm x 6 mm ampoule and an axial temperature gradient of 26.5 C/cm, velocity profiles have been obtained accurately at various positions near the interface. The maximum velocity of 1.29 mm/sec was measured in the central vertical plane and the flow was symmetrical about that plane. A flow inversion point was also noted, above which the flow was towards the interface and below it, away from the interface. The results obtained are useful for validating 3-dimensional numerical models and establishing a link between the macroscopic processing conditions and the formation of crystalline defects.

Kawaji, M.↗

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

Interaction of Porosity with a Planar Solid/Liquid Interface

In this paper, an investigation of the interaction between gas porosity and a planar solid/liquid (SL) interface is reported. A two-dimensional numerical model able to accurately track sharp SL interfaces during solidification of pure metals and alloys is proposed. The finite difference method and a rectangular undeformed grid are used for computation. The SL interface is described through the points of intersection with the grid lines. Its motion is determined by the thermal and solute gradients at each particular point. Changes of the interface temperature because of capillarity or solute redistribution as well as any perturbation of the thermal and solute field produced by the presence of non-metallic inclusions can be computed. To validate the model, the dynamics of the interaction between a gas pore and a solidification front in metal alloys was observed using a state of the art X-ray Transmission Microscope. The experiments included observation of the distortion of the SL interface near a pore, real-time measurements of the growth rate and the change in shape of the porosity during interaction with an advancing SL interface in pure Al and Al-0.25 wt% Au alloy. In addition, porosity induced solute segregation patterns surrounding a pore were also quantified.

Catalina, Adrian V.↗

Interaction of Porosity with a Planar Solid/Liquid Interface

In this article, an investigation of the interaction between gas porosity and a planar solid/liquid (SL) interface is reported. A two-dimensional numerical model able to accurately track sharp SL interfaces during solidification of pure metals and alloys is proposed. The finite-difference method and a rectangular undeformed grid are used for computation. The SL interface is described through the points of intersection with the grid lines. Its motion is determined by the thermal and solute gradients at each particular point. Changes of the interface temperature because of capillarity or solute redistribution as well as any perturbation of the thermal and solute field produced by the presence of non-metallic inclusions can be computed. To validate the model, the dynamics of the interaction between a gas pore and a solidification front in metal alloys was observed using a state of the art X-ray transmission microscope (XTM). The experiments included observation of the distortion of the SL interface near a pore, real-time measurements of the growth rate, and the change in shape of the porosity during interaction with the SL interface in pure Al and Al-0.25 wt pct Au alloy. In addition, porosity-induced solute segregation patterns surrounding a pore were also quantified.

Catalina, Adrian V.↗

High-Fidelity Multidisciplinary Sensitivity Analysis Framework for Multipoint Rotorcraft Optimization

A multidisciplinary, gradient-based sensitivity-analysis methodology is evaluated for optimization of rotorcraft configurations. The tightly coupled discipline models include physics-based fluid dynamics and rotorcraft comprehensive analysis. A discretely-consistent adjoint method accounts for sensitivities of the unsteady flow and unstructured, dynamic, over set grids, while sensitivities of structural responses to aerodynamic loads are computed using a complex-variable method. The methodology is applied to optimize the shape of UH-60A Blackhawk helicopter blades for hover and forward flight conditions. The objective of the multipoint design is to simultaneously increase the rotorcraft figure of merit in a hover flight and reduce the rotor power in a forward flight. Trimmed loose-coupling solutions for the baseline configuration are used to initiate the tight-coupling multidisciplinary analysis. Target thrust and rolling and pitching moments are enforced as optimization constraints. The optimized configuration improves the optimization metrics at both design points. The improved performance and all constraints are maintained over many revolutions beyond the optimization interval, satisfying the required flight conditions. Computational cost of the optimization cycle is assessed in a high-performance computing environment and found affordable for design of rotorcraft in general level-flight conditions.

Wang, Li↗

Electronic heterodyne moire deflectometry: A method for transient and three dimensional density fields measurements

Effects of diffraction and nonlinear photographic emulsion characteristics on the performance of deferred electronic heterodyne moire deflectometry are investigated. The deferred deflectometry is used for measurements of nonsteady phase objects where it is difficult to complete the analysis of the field in real time. The sensitivity, accuracy and resolution of the system are calculated and it is shown that they are weakly affected by diffraction and by nonlinear recording. The feactures of the system are significantly improved compared with the conventional deferred intensity moire technique, and are comparable with the online heterodyne moire. The system was evaluated experimentally by deferred measurements of the refractive index gradients of a weak phase object consisting of a large KD*P crystal. This was done by photographing the phase object through a Ronchi grating and analyzing the tranparency with the electronic heterodyne readout system. The results are compared with the measurements performed on the same phase object with online heterodyne moire deflectometry and with heterodyne holographic interferometry methods. Some practical considerations for system improvement are discussed.

Stricker, Josef↗

Stimulus control in pigeons after extended discriminative training

The effects of amount of training on conditioned inhibition and on the degree of stimulus control were studied using pigeons. The ability of an S- associated with non-reinforcement of suppress positive reinforced behavior was acquired very rapidly during discriminative training. Increased S+, S- training appeared to weaken this conditioned inhibitory effect while at the same time more S+ training apparently increased the amount of external inhibition (non-conditioned inhibition) of positively reinforced behavior by a novel stimulus. Behavioral contrast and incremental generalization gradients along the S- dimension (inhibitory dimensional control) were absent at all stages of training. Behavioral contrast and inhibitory dimensional control are therefore not necessary concomitants of conditioned inhibition by an S-. A new method of assessing the suppressive effects of stimuli during generalization tests was described.

Yarczower, M.↗

New results and techniques in space radio astronomy.

The methods and results of early space radioastronomy experiments are reviewed, with emphasis on the RAE 1 spacecraft which was designed specifically and exclusively for radio astronomical studies. The RAE 1 carries two gravity-gradient-stabilized 229-m traveling-wave V-antennas, a 37-m dipole antenna, and a number of radiometer systems to provide measurements over the 0.2 to 9.2 MHz frequency range with a time resolution of 0.5 sec and an absolute accuracy of plus or minus 25%. Observations of solar bursts at frequencies down to 0.2 MHz provide new information on the density, plasma velocity, and dynamics of coronal streamers out to distances greater than 50 solar radii. New information on the distribution of the ionized component of the interstellar medium is being obtained from galactic continuum background maps at frequencies around 4 MHz. Cosmic noise background spectra measured down to 0.5 MHz produce new estimates on the interstellar flux of cosmic rays, on magnetic fields in the galactic halo, and on distant extragalactic radio sources.

Alexander, J. K.↗

A smoothing algorithm using cubic spline functions

Two algorithms are presented for smoothing arbitrary sets of data. They are the explicit variable algorithm and the parametric variable algorithm. The former would be used where large gradients are not encountered because of the smaller amount of calculation required. The latter would be used if the data being smoothed were double valued or experienced large gradients. Both algorithms use a least-squares technique to obtain a cubic spline fit to the data. The advantage of the spline fit is that the first and second derivatives are continuous. This method is best used in an interactive graphics environment so that the junction values for the spline curve can be manipulated to improve the fit.

Smith, R. E., Jr.↗

NASA-UVA Light Aerospace Alloy and Structures Technology Program (LA2ST)

The general objective of the Light Aerospace Alloy and Structures Technology (LA2ST) Program is to conduct interdisciplinary graduate student research on the performance of next generation, light weight aerospace alloys, composites, and associated thermal gradient structures in close collaboration with Langley researchers. Specific technical objectives are established for each research project. Relevant data and basic understanding of material behavior and microstructure, new monolithic and composite alloys, advanced processing methods, new solid and fluid mechanic analyses, measurement advances, and a pool of educated graduate students are sought.

Gangloff, Richard P.↗

Development of BEM for ceramic composites

It is evident that for proper micromechanical analysis of ceramic composites, one needs to use a numerical method that is capable of idealizing the individual fibers or individual bundles of fibers embedded within a three-dimensional ceramic matrix. The analysis must be able to account for high stress or temperature gradients from diffusion of stress or temperature from the fiber to the ceramic matrix and allow for interaction between the fibers through the ceramic matrix. The analysis must be sophisticated enough to deal with the failure of fibers described by a series of increasingly sophisticated constitutive models. Finally, the analysis must deal with micromechanical modeling of the composite under nonlinear thermal and dynamic loading. This report details progress made towards the development of a boundary element code designed for the micromechanical studies of an advanced ceramic composite. Additional effort has been made in generalizing the implementation to allow the program to be applicable to real problems in the aerospace industry.

Henry, D. P.↗

A Generalized Wall Function

The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended to more complex flows such as the flows with various pressure gradients, zero wall stress and rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to have the fine grids near the wall unless the near-wall flow structures are required to resolve. These solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both adverse and favorable pressure gradients on the surface flow is very significant. Such as unified wall function will be useful not only in deriving analytic expressions for surface flow properties but also bringing a great convenience for CFD methods to place accurate boundary conditions at any location away from the wall. The extended wall functions introduced in this paper can be used for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces.

Shih, Tsan-Hsing↗

Aerodynamic Design Optimization on Unstructured Grids with a Continuous Adjoint Formulation

A continuous adjoint approach for obtaining sensitivity derivatives on unstructured grids is developed and analyzed. The derivation of the costate equations is presented, and a second-order accurate discretization method is described. The relationship between the continuous formulation and a discrete formulation is explored for inviscid, as well as for viscous flow. Several limitations in a strict adherence to the continuous approach are uncovered, and an approach that circumvents these difficulties is presented. The issue of grid sensitivities, which do not arise naturally in the continuous formulation, is investigated and is observed to be of importance when dealing with geometric singularities. A method is described for modifying inviscid and viscous meshes during the design cycle to accommodate changes in the surface shape. The accuracy of the sensitivity derivatives is established by comparing with finite-difference gradients and several design examples are presented.

Anderson, W. Kyle↗