Search NASA⌕ Search

SEARCH · Search NASA

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 145 records · Page 8

Field-gradient partitioning for fracture and frictional contact in the material point method: Field-gradient partitioning for fracture and frictional contact in the material point method [Fracture and frictional contact in material point method using damage-field gradients for velocity-field partitioning]

Contact and fracture in the material point method require grid-scale enrichment or partitioning of material into distinct velocity fields to allow for displacement or velocity discontinuities at a material interface. We present a new method which a kernel-based damage field is constructed from the particle data. The gradient of this field is used to dynamically repartition the material into contact pairs at each node. Our approach avoids the need to construct and evolve explicit cracks or contact surfaces and is therefore well suited to problems involving complex 3-D fracture with crack branching and coalescence. A straightforward extension of this approach permits frictional ‘self-contact’ between surfaces that are initially part of a single velocity field, enabling more accurate simulation of granular flow, porous compaction, fragmentation, and comminution of brittle materials. Finally, numerical simulations of self contact and dynamic crack propagation are presented to demonstrate the accuracy of the approach.

97 MATHEMATICS AND COMPUTING↗

QUANTIFYING THE SIGNIFICANCE OF THE MAGNETIC FIELD FROM LARGE-SCALE CLOUD TO COLLAPSING CORE: SELF-SIMILARITY, MASS-TO-FLUX RATIO, AND STAR FORMATION EFFICIENCY

Dust polarization observational results are analyzed for the high-mass star formation region W51 from the largest parent cloud ({approx}2 pc, James Clerk Maxwell Telescope) to the large-scale envelope ({approx}0.5 pc, BIMA array) down to the collapsing core e2 ({approx}60 mpc, Submillimeter Array). Magnetic field and dust emission gradient orientations reveal a correlation which becomes increasingly more tight with higher resolution. The previously developed polarization-intensity-gradient method is applied in order to quantify the magnetic field significance. This technique provides a way to estimate the local magnetic field force compared to gravity without the need of any mass or field strength measurements, solely making use of measured angles which reflect the geometrical imprint of the various forces. All three data sets clearly show regions with distinct features in the field-to-gravity force ratio. Azimuthally averaged radial profiles of this force ratio reveal a transition from a field dominance at larger distances to a gravity dominance closer to the emission peaks. Normalizing these profiles to a characteristic core scale points toward self-similarity. Furthermore, the polarization-intensity-gradient method is linked to the mass-to-flux ratio, providing a new approach to estimate the latter one without mass and field strength inputs. A transition from a magnetically supercritical to a subcritical state as a function of distance from the emission peak is found for the e2 core. Finally, based on the measured radius-dependent field-to-gravity force ratio we derive a modified star formation efficiency with a diluted gravity force. Compared to a standard (free-fall) efficiency, the observed field is capable of reducing the efficiency down to 10% or less.

79 ASTRONOMY AND ASTROPHYSICS↗

Prediction method of basic domain structure in Fe3%Si(110) single crystal with grooved surface

This paper proposes the method to accurately predict the 180° basic domain width (D{sub w}) in demagnetized states of the grooved Fe3%Si(110) single crystal with the tilt angle of [001] out of the sheet surface (β). The evaluation of D{sub w} enables the estimation of the anomalous eddy current losses. In this paper, D{sub w} is optimized to minimize the magnetic Gibbs free energy represented by vector potentials using the finite element method and the conjugate gradient method. The μ*-method is adopted to approximate the magnetization relaxation. The stray field energy generated by the magnetic charges occurring on both grooves cross section and sheet surfaces is considered in our proposed method. The validity of the proposed method was confirmed by comparison with the observed D{sub w}. As a result, we could reveal the β dependence of D{sub w} against the groove depth. Moreover, the theoretical threshold of the domain refinement due to the grooves has been suggested.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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

Stochastic average model methods

We consider the solution of finite-sum minimization problems, such as those appearing in nonlinear least-squares or general empirical risk minimization problems. We are motivated by problems in which the summand functions are computationally expensive and evaluating all summands on every iteration of an optimization method may be undesirable. Here we present the idea of stochastic average model (SAM) methods, inspired by stochastic average gradient methods. SAM methods sample component functions on each iteration of a trust-region method according to a discrete probability distribution on component functions; the distribution is designed to minimize an upper bound on the variance of the resulting stochastic model. We present promising numerical results concerning an implemented variant extending the derivative-free model-based trust-region solver POUNDERS, which we name SAM-POUNDERS.

97 MATHEMATICS AND COMPUTING↗

Low-Field Nuclear Magnetic Resonance Characterization of Carbonate and Sandstone Reservoirs From Rock Spring Uplift of Wyoming

Laboratory measurements including gas (N 2 ) porosity and permeability, time-domain nuclear magnetic resonance, thin section, and scanning electron microscopy analysis were conducted to obtain petrographical and petrophysical descriptions of the Weber Sandstone and Madison Limestone at the Rock Spring Uplift, a potential carbon dioxide storage site in Southwestern Wyoming. The relationships between pore structures, such as pore geometry, pore-size distribution, pore network, and porosity/permeability are investigated. First, using thin sections combined with scanning electron microscopy for pore structures description, all samples are described in detail from the geological, petrographysical, and diagenetic viewpoint. Results show that within the Madison Limestone, pore types include intercrystalline, vuggy, moldic, or mixed (combination of all other pore types). Both moldic and vuggy pore types are associated with samples of high porosity and permeability. Nuclear magnetic resonance relaxation time distributions show either bimodal or multimodal distributions. Large relaxation time components are associated with samples with large pores, whereas small components are dominated by small pores. The T 2 geometric mean correlates well with gas permeability. Additionally, short-time diffusion coefficients (D) were measured by pulsed field gradient method using a series of gradient strengths. We found that diffusion coefficient distributions correlate with the corresponding T 2 distributions for macropores. By comparing the dominant peak position of T 2 distributions and their corresponding diffusion coefficient distributions, we predicted the surface relaxivity of different rock types. We found that surface relaxivities of Weber Sandstone samples can be well predicted, while for Madison Limestone samples, surface relaxivities are overestimated due to diffusive pore coupling effect.

58 GEOSCIENCES↗

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↗

An Efficient Numerical Algorithm for Solving Coupled Time-Dependent Ginzburg-Landau Equation for Superconductivity and Elasticity

A decoupled finite element algorithm is developed for simulating the vortex dynamics on an elastic superconductor which couples the time-dependent Ginzburg- Landau equation with the complex-valued superconducting order parameter and the vector-valued magnetic potential, and the elasticity equation. We present an iterative algorithm for the decoupled system arising from the time and spatial discretization using a combination of preconditioner, algebraic multigrid method (AMG) and preconditioned conjugate gradient method (PCG). The iterative algorithm allows us to perform large-scale three-dimensional simulations of mesoscale pattern formation during superconducting phase transitions with arbitrary elastic boundary conditions. Here, the performance and efficiency of the algorithm are numerically verified by several benchmark problems, exhibiting up to two orders of magnitude improvement depending on the scale of discrete system compared to the exact solver.

Efficiency↗

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