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 559 records · Page 31

A comparison of experiment and theory for sound propagation in variable area ducts

An experimental and analytical program has been carried out to evaluate sound suppression techniques in ducts that produce refraction effects due to axial velocity gradients. The analytical program employs a computer code based on the method of multiple scales to calculate the influence of axial variations due to slow changes in the cross-sectional area as well as transverse gradients due to the wall boundary layers. Detailed comparisons between the analytical predictions and the experimental measurements have been made. The circumferential variations of pressure amplitudes and phases at several axial positions have been examined in straight and variable area ducts, with hard walls and lined sections, and with and without a mean flow. Reasonable agreement between the theoretical and experimental results has been found.

Nayfeh, A. H.↗

Analytic gradients for coupled-cluster energies that include noniterative connected triple excitations - Application to cis- and trans-HONO

An efficient formulation of the analytic energy gradient for the single and double excitation coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations is presented. The formulation has a small computational cost, and the algebraic manipulations may be applied generally to the analytic gradient of Moller-Plesset perturbation theory energies. The new formulation has been implemented in an efficient set of programs that utilize highly vectorized algorithms and has been used to investigate the equilibrium structures, harmonic vibrational frequencies, IR intensities, and energy separation of cis- and trans-HONO.

Lee, Timothy J.↗

A mesh gradient technique for numerical optimization

A class of successive-improvement optimization methods in which directions of descent are defined in the state space along each trial trajectory are considered. The given problem is first decomposed into two discrete levels by imposing mesh points. Level 1 consists of running optimal subarcs between each successive pair of mesh points. For normal systems, these optimal two-point boundary value problems can be solved by following a routine prescription if the mesh spacing is sufficiently close. A spacing criterion is given. Under appropriate conditions, the criterion value depends only on the coordinates of the mesh points, and its gradient with respect to those coordinates may be defined by interpreting the adjoint variables as partial derivatives of the criterion value function. In level 2, the gradient data is used to generate improvement steps or search directions in the state space which satisfy the boundary values and constraints of the given problem.

Willis, E. A., Jr.↗

Controlling the Mobility of the Fluid Interface of Moving Gas Bubbles or Liquid Drops by Using Micellar Solutions of Surfactants

Microgravity processes must rely on mechanisms other than buoyancy to move bubbles or droplets from one region to another in a continuous liquid phase. One suggested method is thermocapillary migration in which a temperature gradient is applied to the continuous phase. A significant and as yet unresolved impediment to the use of thermocapillary migration to direct bubble or drop motion is that these migrations can be significantly retarded by the adsorption onto the fluid particle surface of surface active impurities unavoidably present in the continuous or (if the particle is a liquid) droplet phases. The focus of our research was to develop a theory for remobilizing fluid particle interfaces retarded by a surfactant impurity in an effort to make more viable the use of thermocapillary migrations for the management of bubbles and drops in microgravity. We postulated that a surfactant at high bulk concentration which kinetically exchanges rapidly with the surface can restore interface mobility. The scaling arguments along with a discussion of the previous literature is reviewed in the context of the scaling framework. The specific objectives of the research were twofold. The first was to prove the remobilization theory by studying a model problem. As the mechanism for remobilization is independent of the force which drives the particle, the fluid particle shape and the presence of fluid inertia, we chose the simplest model consisting of a spherical bubble rising steadily by buoyancy in creeping flow. We solved the hydrodynamic and surfactant transport equations for rapid kinetic exchange to demonstrate that as the concentration increases, the Marangoni retardation at first increases (the low k behavior) and then decreases (the high k behavior). The second objective was to develop a method to determine the kinetic rate constants of a surfactant molecule, since this information is necessary to select surfactants which will exchange rapidly enough relative to the convective rate in the thermocapillary process of interest. To measure the kinetic rate, we measure the dynamic tension change accompanying adsorption onto an initially clean interface, or the re-equilibration in tension when an equilibrium interface is compressed. The dynamic tension measurements are made by a pendant bubble method, in which surfactant adsorbs onto a pendant bubble, and the tension is measured by analyzing the shape change in the bubble. We conclude this report by detailing the publications, presentations and doctoral thesis completed under the auspices of this grant.

Maldarelli, Charles↗

A New Correction Technique for Strain-Gage Measurements Acquired in Transient-Temperature Environments

Significant strain-gage errors may exist in measurements acquired in transient-temperature environments if conventional correction methods are applied. As heating or cooling rates increase, temperature gradients between the strain-gage sensor and substrate surface increase proportionally. These temperature gradients introduce strain-measurement errors that are currently neglected in both conventional strain-correction theory and practice. Therefore, the conventional correction theory has been modified to account for these errors. A new experimental method has been developed to correct strain-gage measurements acquired in environments experiencing significant temperature transients. The new correction technique has been demonstrated through a series of tests in which strain measurements were acquired for temperature-rise rates ranging from 1 to greater than 100 degrees F/sec. Strain-gage data from these tests have been corrected with both the new and conventional methods and then compared with an analysis. Results show that, for temperature-rise rates greater than 10 degrees F/sec, the strain measurements corrected with the conventional technique produced strain errors that deviated from analysis by as much as 45 percent, whereas results corrected with the new technique were in good agreement with analytical results.

Richards, W. Lance↗

Turbulent atmospheric flow over a backward-facing step

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.

Kaul, U. K.↗

Simple methods for determining the virtual origin of turbulent boundary layers in hypersonic flow on sharp-edged flat plates

Two methods for determining the virtual origin of turbulent boundary layers in hypersonic flow are evaluated. The results of the analyses are restricted to wind-tunnel models having sharp-edged surfaces with zero or small pressure gradients. Virtual origin and skin friction estimates from these two methods are compared with values from a base method for which the virtual origin is calculated from the measured momentum thickness at a station downstream of boundary layer transition.

Hopkins, E. J.↗

Remobilizing the Interface of Thermocapillary Driven Bubbles Retarded By the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow that propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillary to direct the movement of bubbles in space is the fact that surfactant impurities, which are unavoidably present in the continuous phase, can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. (This can also be the case if kinetic desorption of surfactant at the back end of the bubble is much slower than convection.) For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like and bubbles translate as solid particles). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature-induced tension gradient, and can decrease to near zero the bubble's thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have be retarded by the adsorption onto the bubble surface of a surfactant impurity. Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble surface much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition, the rapid kinetic exchange and high bulk concentration maintain a saturated surface with uniform surface concentrations. This prevents retarding surface tension gradients and keeps the thermocapillary velocity high. In our reports over the first 2 years, we presented numerical simulations of the bubble motion and surfactant transport which verified theoretically the concept of remobilization, and the development of an apparatus to track and measure the velocity of rising bubbles in a glycerol/water surfactant solution. This year, we detail experimental observations of remobilization. Two polyethylene oxide surfactants were studied, C12E6 (CH3(CH2)11(OCH2)6OH) and C10E8 (CH3(CH2)4(OCH2CH2)8OH). Measurements of the kinetic exchange for these surfactants show that the one with the longer hydrophobe chain C12E6 has a lower rate of kinetic exchange. In addition, this surfactant is much less soluble in the glycerol/water mixture because of the shorter ethoxylate chain. As a result, we found that C12E6 had only a very limited ability to remobilize rising bubbles because of the limited kinetic exchange and reduced solubility. However, C10E8, with its higher solubility and more rapid exchange was found to dramatically remobilize rising bubbles. We also compared our theoretical calculations to the experimental measurements of velocity for both the non-remobilizing and remobilizing surfactants and found excellent agreement. We further observed that for C10E8 at high concentrations, which exceeded the critical micelle concentrations, additional remobilization was measured. In this case the rapid exchange of monomer between micelle and surfactant provides an additional mechanism for maintaining a uniform surface concentrations.

Palaparthi, Ravi↗

An adaptive grid, unsteady model for two-dimensional magnetohydrodynamic (MHD) flow

An adaptive grid finite difference method for solving multi-dimensional, time-dependent, magnetohydrodynamic (MHD) equations is developed. The method is capable of solving problems that include high gradients due to geometry, propagation of shock waves and unsteady boundary conditions. The grid generation technique is based on variational principles with direct control over grid concentration, smoothness and skewness. An example for a two-dimensional MHD simulation of the propagation of a solar-flare-generated shock wave in solar wind flow in the heliographic equatorial plane is selected for illustration of this method.

Panitchob, Supat↗

Data Images and Other Graphical Displays for Directional Data

Vectors, axes, and periodic phenomena have direction. Directional variation can be expressed as points on a unit circle and is the subject of circular statistics, a relatively new application of statistics. An overview of existing methods for the display of directional data is given. The data image for linear variables is reviewed, then extended to directional variables by displaying direction using a color scale composed of a sequence of four or more color gradients with continuity between sequences and ordered intuitively in a color wheel such that the color of the 0deg angle is the same as the color of the 360deg angle. Cross over, which arose in automating the summarization of historical wind data, and color discontinuity resulting from the use a single color gradient in computational fluid dynamics visualization are eliminated. The new method provides for simultaneous resolution of detail on a small scale and overall structure on a large scale. Example circular data images are given of a global view of average wind direction of El Nino periods, computed rocket motor internal combustion flow, a global view of direction of the horizontal component of earth's main magnetic field on 9/15/2004, and Space Shuttle solid rocket motor nozzle vectoring.

Morphet, Bill↗

Very Large Scale Optimization

The purpose of this research under the NASA Small Business Innovative Research program was to develop algorithms and associated software to solve very large nonlinear, constrained optimization tasks. Key issues included efficiency, reliability, memory, and gradient calculation requirements. This report describes the general optimization problem, ten candidate methods, and detailed evaluations of four candidates. The algorithm chosen for final development is a modern recreation of a 1960s external penalty function method that uses very limited computer memory and computational time. Although of lower efficiency, the new method can solve problems orders of magnitude larger than current methods. The resulting BIGDOT software has been demonstrated on problems with 50,000 variables and about 50,000 active constraints. For unconstrained optimization, it has solved a problem in excess of 135,000 variables. The method includes a technique for solving discrete variable problems that finds a "good" design, although a theoretical optimum cannot be guaranteed. It is very scalable in that the number of function and gradient evaluations does not change significantly with increased problem size. Test cases are provided to demonstrate the efficiency and reliability of the methods and software.

Vanderplaats, Garrett↗

A computer program for the calculation of laminar and turbulent boundary layer flows

The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.

Dwyer, H. A.↗

Geysering inhibitor for vertical cryogenic transfer piping

Geysering (i.e., the expulsion of boiling liquid and its vapor from a vertical tube) has been a problem for the missile industry in long vertical cryogenic propellant feed lines connecting the launch vehicle propellant tank with the rocket engines. A proposed novel method of inhibiting geysering and the associated pressure gradients provides a self-starting self-regulating action that is not dependent on other active systems or components. The inhibiting action is attained by incorporating a concentric tube within the main transfer tube to prevent constriction of natural convective flow.

Howard, F. S.↗

High-latitude truncation errors of box-type primitive equation models

The 'box-type' finite-difference method includes a weighted average of the pressure gradient with weights proportional to the surface of the grid walls. It is shown that this averaging introduces first-order truncation errors near the poles. An example is shown in which the relative error is of zero order and the scheme produces large distortions in the solution at high latitudes.

Kalnay-Rivas, E.↗

Bimetallic junctions

The formation of voids through interdiffusion in bimetallic welded structures exposed to high operating temperatures is inhibited by utilizing an alloy of the parent materials in the junction of the parent materials or by preannealing the junction at an ultrahigh temperature. These methods are also used to reduce the concentration gradient of a hardening agent.

Arcella, F. G.↗

Tropical cyclone outer surface winds derived from satellite microwave sounder data

Upper tropospheric temperature anomalies are detected in brightness temperature data from the Nimbus 6 Scanning Microwave Spectrometer (SCAMS). Brightness temperature anomalies are related to surface pressure anomalies through the radiative transfer and hydrostatic equation. Surface wind speeds at outer radii are then estimated using the gradient wind equation and a shearing parameter. The method is first tested using simulated satellite data constructed from temperature, pressure and height data recorded by aircraft reconnaissance of four hurricanes. Wind speeds in the 80-95 kPa region are estimated with 2-3 m/sec accuracy. Next, 55.45 GHz SCAMS data over eight typhoons during 1975 are used to estimate the radii of 15.4 m/sec (30 kt) and 27.5 m/sec (50 kt) winds. Accuracies of about + or - 80 and + or - 70 km, respectively, are found. It is suggested that the technique be further tested using data from the Microwave Sounding Unit on board the TIROS-N and NOAA 6 satellites.

Kidder, S. Q.↗

Optical observations of unidirectional solidification in microgravity

Optical interferometric, shadowgraph, and streak photographic methods are used to obtain temperature profiles, concentration gradients, and fluid velocities in NH4Cl-H2O metal-model solution during unidirectional solidification in microgravity. This study elucidates earlier low gravity suborbital rocket experiments and lays the groundwork for space shuttle experiments. The design and operation of the optical units are presented, and results are detailed to ground based experiments. The low gravity experiments were conducted using a KC-135 aircraft flying a parabolic trajectory.

Johnston, M. H.↗

Analysis of Doppler lidar wind measurements

New analysis methods were developed and refined for the synthesis of wind fields from multiple Doppler radar data. Considerable software was written and debugged for wind field synthesis. The analysis methods were applied to radar data collected during the CCOPE experiment in which lidar wind measurements were also made. The analysis method differs from existing methods in that interpolation (using filters with assigned weights) and advection of data, to account for time differential between data acquisition at different points, are not used. Instead, functional fits were used in space and time coordinates to transfer radial velocity data to a space-time grid. Subsequently the wind synthesis is performed. This method gives better space resolution in the synthesized winds, compared to other methods in common use, by better preserving the gradients of velocities. The wind fields show the cellular structure characteristic of boundary-layer convection. Also evident in the analysis are the sharp transitions and gradients in the wind field which were preserved by the synthesis.

Srivastava, R. C.↗