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 91 records · Page 5

Induced electric currents in the Alaska oil pipeline measured by gradient, fluxgate, and SQUID magnetometers

The field gradient method for observing the electric currents in the Alaska pipeline provided consistent values for both the fluxgate and SQUID method of observation. These currents were linearly related to the regularly measured electric and magnetic field changes. Determinations of pipeline current were consistent with values obtained by a direct connection, current shunt technique at a pipeline site about 9.6 km away. The gradient method has the distinct advantage of portability and buried- pipe capability. Field gradients due to the pipe magnetization, geological features, or ionospheric source currents do not seem to contribute a measurable error to such pipe current determination. The SQUID gradiometer is inherently sensitive enough to detect very small currents in a linear conductor at 10 meters, or conversely, to detect small currents of one amphere or more at relatively great distances. It is fairly straightforward to achieve imbalance less than one part in ten thousand, and with extreme care, one part in one million or better.

Campbell, W. H.↗

Comparison of solution approaches for minimum-fuel, low-thrust, power-limited orbital transfers

An initial assessment of the feasibility of a function space gradient method for computing solutions to minimum-fuel power-limited transfers encompassing a wide range of thrust to weight ratios is conducted. Three transfers between coplanar ellipses are used as test cases. The gradient method performs best at the high end of the thrust to weight ratio range. At the lower end, there is reduced sensitivity of the fuel consumption to the control profiles. The minimum fuel consumption and the trajectory are computed quite accurately but the control profiles are in error. An approximate analytical solution, obtained by Edelbaum using the method of averaging, is discussed.

Mease, Kenneth D.↗

Microstratigraphy of the lunar regolith and compaction ages of lunar breccias

Fossil track analyses of impregnated core sections from the Apollo 15 and 16 deep drill cores were used to study two types of distinct layers of about 1 mm and about several mm thick in the lunar regolith. The data supports a mixing model hypothesis for the origin of the sampled soil column dominated by large, fresh feldspar fragments. The boundary track method and the track gradient method were used to study the compaction ages of lunar breccias.

Goswami, J. N.↗

A finite element conjugate gradient FFT method for scattering

An extension of a two dimensional formulation is presented for a three dimensional body of revolution. With the introduction of a Fourier expansion of the vector electric and magnetic fields, a coupled two dimensional system is generated and solved via the finite element method. An exact boundary condition is employed to terminate the mesh and the fast fourier transformation (FFT) is used to evaluate the boundary integrals for low O(n) memory demand when an iterative solution algorithm is used. By virtue of the finite element method, the algorithm is applicable to structures of arbitrary material composition. Several improvements to the two dimensional algorithm are also described. These include: (1) modifications for terminating the mesh at circular boundaries without distorting the convolutionality of the boundary integrals; (2) the development of nonproprietary mesh generation routines for two dimensional applications; (3) the development of preprocessors for interfacing SDRC IDEAS with the main algorithm; and (4) the development of post-processing algorithms based on the public domain package GRAFIC to generate two and three dimensional gray level and color field maps.

Collins, Jeffery D.↗

A finite element conjugate gradient FFT method for scattering

Validated results are presented for the new 3D body of revolution finite element boundary integral code. A Fourier series expansion of the vector electric and mangnetic fields is employed to reduce the dimensionality of the system, and the exact boundary condition is employed to terminate the finite element mesh. The mesh termination boundary is chosen such that is leads to convolutional boundary operatores of low O(n) memory demand. Improvements of this code are discussed along with the proposed formulation for a full 3D implementation of the finite element boundary integral method in conjunction with a conjugate gradiant fast Fourier transformation (CGFFT) solution.

Collins, Jeffery D.↗

Bulk Crystal Growth of Nonlinear Optical Organic Materials Using Inverted Vertical Gradient Freeze Method

A new process for producing large bulk single crystals of benzil (C6H5COCOC6H5) is reported in this paper. Good quality crystals have been successfully grown using this approach to crystal growth. This method seems to be very promising for other thermally stable NLO organic materials also. The entire contents vycor crucible 1.5 inch in diameter and 2 inch deep was converted to single crystal. Purity of the starting growth material is also an important factor in the final quality of the grown crystals. The entire crystal can be very easily taken out of the crucible by simple maneuvering. Initial characterization of the grown crystals indicated that the crystals are as good as other crystals grown by conventional Bridgman Stockbarger technique.

Choi, J.↗

Coherent gradient sensing method and system for measuring surface curvature

A system and method for determining a curvature of a specularly reflective surface based on optical interference. Two optical gratings are used to produce a spatial displacement in an interference field of two different diffraction components produced by one grating from different diffraction components produced by another grating. Thus, the curvature of the surface can be determined.

Rosakis, Ares J.↗

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗

Algorithms for parallel and vector computations

This is a final report on work performed under NASA grant NAG-1-1112-FOP during the period March, 1990 through February 1995. Four major topics are covered: (1) solution of nonlinear poisson-type equations; (2) parallel reduced system conjugate gradient method; (3) orderings for conjugate gradient preconditioners, and (4) SOR as a preconditioner.

Ortega, James M.↗

Schwarz Methods: To Symmetrize or not to Symmetrize

A preconditioning theory for Schwarz methods is presented. The theory establishes sufficient conditions for multiplicative and additive Schwarz algorithms to yield self-adjoint positive definite preconditioners. It allows for the analysis and use of non-variational and non-convergent linear methods as preconditioners for conjugate gradient methods, and it is applied to domain decomposition and multigrid. This paper illustrates why symmetrizing may be a bad idea for linear methods. Numerical examples are presented for a test problem.

Holst, Michael↗

A displacement gradient BEM for accurate stress computation near boundaries in 2-D anisotropic problems

A displacement gradient method for 2D anisotropic elasticity problems is presented, which effectively minimizes the boundary layer effect through a two-step procedure. First, the boundary integral equations are solved for the unknown boundary displacements and tractions. Second, a direct integral equation for displacement gradients is developed in terms of boundary tractions. Three methods based on different evaluation procedures and locations for determining the displacement gradients are proposed. In the first method the displacement gradients are averaged at nodes common to adjacent elements. The second method stores the gradients element-wise. In the third method, the gradients are evaluated at the nodes of discontinuous elements. The three methods are applied to near-isotropic plates with circular and elliptic cutouts. It is concluded that all three methods can yield accurate stress distributions.

Sistla, R.↗

Reanalysis procedure for large structural systems

Global-basis-vector approximate reanalysis techniques for use in automated structural optimization schemes are developed and demonstrated. The vibrational response of a modified structure (MS) is estimated by lumping the design variables into a single tracing parameter, applying an operator splitting procedure to express the FEM equations of the MS in terms of the original-structure equations plus correction terms, and reducing the MS equations via a classical Bubnov-Galerkin scheme. The sensitivity of the vibrational response to structural modifications is evaluated, and the relationship between this method and the preconditioned conjugate-gradient method (Noor and Peters, 1988) is explored. Numerical results for linear static and free vibration problems involving beamlike lattices, double-layered hexahedral grids, and structural gridworks are presented in tables and graphs and briefly characterized.

Noor, Ahmed K.↗

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.↗

Detection of Sub-Surface Water on Mars by Controlled and Natural Source Electromagnetic Induction

Detection of subsurface liquid water on Mars is a leading scientific objective for Mars exploration in this decade. We describe electromagnetic induction (EM) methods that are both uniquely well suited for detection of subsurface liquid water on Mars and practical within the context of a Mars exploration program. EM induction methods are ideal for detection of more highly conducting (liquid water bearing) soils and rock beneath a more resistive overburden. A combined natural source and controlled source method offers an efficient and unambiguous characterization of the depth to liquid water and the extent of the aqueous region. The controlled source method employs an ac vertical dipole source (horizontal loop) to probe the depth to the conductor and a natural source method (gradient sounding) to characterize its conductivity-thickness product. These methods are proven in geophysical exploration and can be tailored to cope with any reasonable Mars crustal electrical conductivity. We describe a practical experiment and discuss experiment optimization to address the range of material properties likely encountered in the Mars crust.

Connerney, J. E. P.↗

Spectral ordering techniques for incomplete LU preconditoners for CG methods

The effectiveness of an incomplete LU (ILU) factorization as a preconditioner for the conjugate gradient method can be highly dependent on the ordering of the matrix rows during its creation. Detailed justification for two heuristics commonly used in matrix ordering for anisotropic problems is given. The bandwidth reduction and weak connection following heuristics are implemented through an ordering method based on eigenvector computations. This spectral ordering is shown to be a good representation of the heuristics. Analysis and test cases in two and three dimensional diffusion problems demonstrate when ordering is important, and when an ILU decomposition will be ordering insensitive. The applicability of the heuristics is thus evaluated and placed on a more rigorous footing.

Clift, Simon S.↗

Minimum impulse three-body trajectories.

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

D'Amario, L.↗