Search NASA⌕ Search

SEARCH · Search NASA

Results for “Jacobi Method”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗
NASA NTRS1991

Unstructured grid generation using the distance function

A new class of methods for obtaining level sets to generate unstructured grids is presented. The consecutive grid levels are computed using the distance functions, which corresponds to solving the Hamilton-Jacobi equations representing the equations of motion of fronts propagating with curvature-dependent speed. The relationship between the distance function and the governing equations will be discussed as well as its application to generating grids. Multi-ply connected domains and complex geometries are handled naturally, with a straightforward generalization to several space dimensions. The grid points for the unstructured grid are obtained simultaneously with the grid levels. The search involved in checking for overlapping triangles is minimized by triangulating the entire domain one level at a time.

Bihari, Barna L.↗
NASA NTRS1993

Uniform convergence of multigrid V-cycle iterations for indefinite and nonsymmetric problems

In this paper, we present an analysis of a multigrid method for nonsymmetric and/or indefinite elliptic problems. In this multigrid method various types of smoothers may be used. One type of smoother which we consider is defined in terms of an associated symmetric problem and includes point and line, Jacobi, and Gauss-Seidel iterations. We also study smoothers based entirely on the original operator. One is based on the normal form, that is, the product of the operator and its transpose. Other smoothers studied include point and line, Jacobi, and Gauss-Seidel. We show that the uniform estimates for symmetric positive definite problems carry over to these algorithms. More precisely, the multigrid iteration for the nonsymmetric and/or indefinite problem is shown to converge at a uniform rate provided that the coarsest grid in the multilevel iteration is sufficiently fine (but not depending on the number of multigrid levels).

Bramble, James H.↗
NASA NTRS2005

Actively Cooled SLMS(TM) Technology for HEL Applications

Mr. Jacoby is the Chief Scientist for Schafer's Lightweight Optical Systems business area with twenty four years experience in laser and optical systems for space and military applications. He and colleague Dr. Goodman conceived and developed Silicon Lightweight Mirrors (SLMS(TM)) technologies for space applications from the extreme UV to FAR IR wavelengths. Schafer has demonstrated two different methods for actively cooling our Silicon Lightweight Mirrors (SLMS(TM)) technology. Direct internal cooling was accomplished by flowing liquid nitrogen through the continuous open cell core of the SLMS(TM) mirror. Indirect external cooling was accomplished by flowing liquid nitrogen through a CTE matched Cesic square-tube manifold that was bonded to the back of the mirror in the center. Testing was done in the small 4-foot thermal/vacuum chamber located at the NASA/MSFC X-Ray Calibration Facility. Seven thermal diodes were located over the front side of the 5 inch diameter mirror and one was placed on the outlet side of the Cesic manifold. Results indicate that the mirror reaches steady state at 82K in less than four minutes for both cooling methods. The maximum temperature difference of the eight diodes was less than 200 mK when the mirror was internally cooled and covered with MLI to insulate it from the large 300 K aluminum plate that was used to mount it.

Jacoby, Marc T.↗
DOE OSTI2015

Trajectory approach to the Schrödinger–Langevin equation with linear dissipation for ground states

The Schrödinger–Langevin equation with linear dissipation is integrated by propagating an ensemble of Bohmian trajectories for the ground state of quantum systems. Substituting the wave function expressed in terms of the complex action into the Schrödinger–Langevin equation yields the complex quantum Hamilton–Jacobi equation with linear dissipation. We transform this equation into the arbitrary Lagrangian–Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation is simultaneously integrated with the trajectory guidance equation. Then, the computational method is applied to the harmonic oscillator, the double well potential, and the ground vibrational state of methyl iodide. The excellent agreement between the computational and the exact results for the ground state energies and wave functions shows that this study provides a synthetic trajectory approach to the ground state of quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Algorithm For Optimal Control Of Large Structures

Cost of computation appears competitive with other methods. Problem to compute optimal control of forced response of structure with n degrees of freedom identified in terms of smaller number, r, of vibrational modes. Article begins with Hamilton-Jacobi formulation of mechanics and use of quadratic cost functional. Complexity reduced by alternative approach in which quadratic cost functional expressed in terms of control variables only. Leads to iterative solution of second-order time-integral matrix Volterra equation of second kind containing optimal control vector. Cost of algorithm, measured in terms of number of computations required, is of order of, or less than, cost of prior algoritms applied to similar problems.

Salama, Moktar A.↗
NASA NTRS1990

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
NASA NTRS1991

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
DOE OSTI2014

Calculation of state-to-state differential and integral cross sections for atom-diatom reactions with transition-state wave packets

A recently proposed transition-state wave packet method [R. Welsch, F. Huarte-Larrañaga, and U. Manthe, J. Chem. Phys. 136, 064117 (2012)] provides an efficient and intuitive framework to study reactive quantum scattering at the state-to-state level. It propagates a few transition-state wave packets, defined by the eigenfunctions of the low-rank thermal flux operator located near the transition state, into the asymptotic regions of the reactant and product arrangement channels separately using the corresponding Jacobi coordinates. The entire S-matrix can then be assembled from the corresponding flux-flux cross-correlation functions for all arrangement channels. Since the transition-state wave packets can be defined in a relatively small region, its transformation into either the reactant or product Jacobi coordinates is accurate and efficient. Furthermore, the grid/basis for the propagation, including the maximum helicity quantum number K, is much smaller than that required in conventional wave packet treatments of state-to-state reactive scattering. This approach is implemented for atom-diatom reactions using a time-dependent wave packet method and applied to the H + D{sub 2} reaction with all partial waves. Excellent agreement with benchmark integral and differential cross sections is achieved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Solution of D dimensional Dirac equation for coulombic potential using NU method and its thermodynamics properties

The analytical solution of Ddimensional Dirac equation for Coulombic potential is investigated using Nikiforov-Uvarov method. The D dimensional relativistic energy spectra are obtained from relativistic energy eigenvalue equation by using Mat Lab software.The corresponding D dimensional radial wave functions are formulated in the form of generalized Jacobi and Laguerre Polynomials. In the non-relativistic limit, the relativistic energy equation reduces to the non-relativistic energy which will be applied to determine some thermodynamical properties of the system. The thermodynamical properties of the system are expressed in terms of error function and imaginary error function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

RIACS

Topics considered include: high-performance computing; cognitive and perceptual prostheses (computational aids designed to leverage human abilities); autonomous systems. Also included: development of a 3D unstructured grid code based on a finite volume formulation and applied to the Navier-stokes equations; Cartesian grid methods for complex geometry; multigrid methods for solving elliptic problems on unstructured grids; algebraic non-overlapping domain decomposition methods for compressible fluid flow problems on unstructured meshes; numerical methods for the compressible navier-stokes equations with application to aerodynamic flows; research in aerodynamic shape optimization; S-HARP: a parallel dynamic spectral partitioner; numerical schemes for the Hamilton-Jacobi and level set equations on triangulated domains; application of high-order shock capturing schemes to direct simulation of turbulence; multicast technology; network testbeds; supercomputer consolidation project.

Oliger, Joseph↗
DOE OSTI2013

A quasi-Bohmian approach for a homogeneous spherical solid body based on its geometric structure

In this paper we express the space of rotation as a Riemannian space and try to generalize the classical equations of motion of a homogeneous spherical solid body in the domain of quantum mechanics. This is done within Bohm's view of quantum mechanics, but we do not use the Schrödinger equation. Instead, we assume that in addition to the classical potential there is an extra potential and try to obtain it. In doing this, we start from a classical picture based on Hamilton-Jacobi formalism and statistical mechanics but we use an interpretation which is different from the classical one. Then, we introduce a proper action and extremize it. This procedure gives us a mathematical identity for the extra potential that limits its form. The classical mechanics is a trivial solution of this method. In the simplest cases where the extra potential is not a constant, a mathematical identity determines it uniquely. In fact the first nontrivial potential, apart from some constant coefficients which are determined by experiment, is the usual Bohmian quantum potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

An unstructured multigrid method for elliptic problems

A multigrid algorithm for implementation on unstructured meshes is proposed. The algorithm uses a sequence of unnested grids and requires the development of efficient inter-grid interpolation procedures. It is demonstrated how elliptic problems can be solved in this fashion by using Jacobi smoothers.

Lohner, R.↗
DOE OSTI2014

Analytical solutions of a generalized non-central potential in N-dimensions

We present that N-dimensional non-relativistic wave equation for the generalized non-central potential with arbitrary angular momentum is analytically solvable in the hyperspherical coordinates. Asymptotic iteration method as a different approach is applied to obtain N-dimensional energy eigenvalues and the corresponding eigenfunctions. In hyperspherical coordinates, the wave function solutions are obtained in terms of hypergeometric functions and Jacobi polynomials. The bound states of quantum systems under consideration for some special cases, such as Hartmann and Makarov potentials, have been discussed in N-dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Stability regions for quasiperiodic motion in the restricted problem of three bodies

Surfaces of section, plotted in configuration space, have been computed for the motion of the massless particle in the restricted problem of three bodies. Nine mass ratios and a wide variety of Jacobi constants were investigated; over four thousand orbits were computed, about half of which were finally used. The plots of surface of section have been reduced to plots of stability regions, following a method due to Henon (1965a, 1965b, 1966b, 1969). Sample surfaces of section are also given. The complete set of 276 surfaces of section has been published as a report (Jefferys, 1971) and is available upon request from the author.

Jefferys, W. H.↗
NASA NTRS2003

Compressed Semi-Discrete Central-Upwind Schemes for Hamilton-Jacobi Equations

We introduce a new family of Godunov-type semi-discrete central schemes for multidimensional Hamilton-Jacobi equations. These schemes are a less dissipative generalization of the central-upwind schemes that have been recently proposed in series of works. We provide the details of the new family of methods in one, two, and three space dimensions, and then verify their expected low-dissipative property in a variety of examples.

Bryson, Steve↗
DOE OSTI2022

Weighted relaxation for multigrid reduction in time

Current trends in computer architectures now mean that faster computation speed must come primarily from increased concurrency, not faster clock speeds, which are stagnating. Thus, this situation creates bottlenecks for serial algorithms, including the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in time is with multigrid. Here, in this work, we consider multigrid-reduction-in-time (MGRIT), a multilevel method applied to the time dimension that computes multiple time-steps in parallel. Like all multigrid methods, MGRIT relies on the complementary relationship between relaxation on a fine-grid and a correction from the coarse grid to solve the problem. All current MGRIT implementations are based on unweighted-Jacobi relaxation; here we introduce the concept of weighted relaxation to MGRIT. We derive new convergence bounds for weighted relaxation, and use this analysis to guide the selection of relaxation weights. Numerical results then demonstrate that by choosing appropriate non-unitary relaxation weights, one can achieve faster convergence rates and lower iteration counts for MGRIT when compared with unweighted relaxation. In most cases, weighted relaxation yields a 10%–20% saving in iterations, which is significant when using large high-performance computers. For A-stable integration schemes, results also illustrate that under-relaxation can restore convergence in some cases where unweighted relaxation is not convergent.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Next-to-$\mathrm{MHV}$ Yang-Mills kinematic algebra

Kinematic numerators of Yang-Mills scattering amplitudes possess a rich Lie algebraic structure that suggest the existence of a hidden infinite-dimensional kinematic algebra. Explicitly realizing such a kinematic algebra is a longstanding open problem that only has had partial success for simple helicity sectors. In past work, we introduced a framework using tensor currents and fusion rules to generate BCJ numerators of a special subsector of NMHV amplitudes in Yang-Mills theory. Here we enlarge the scope and explicitly realize a kinematic algebra for all NMHV amplitudes. Master numerators are obtained directly from the algebraic rules and through commutators and kinematic Jacobi identities other numerators can be generated. Inspecting the output of the algebra, we conjecture a closed-form expression for the master BCJ numerator up to any multiplicity. We also introduce a new method, based on group algebra of the permutation group, to solve for the generalized gauge freedom of BCJ numerators. It uses the recently introduced binary BCJ relations to provide a complete set of NMHV kinematic numerators that consist of pure gauge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗