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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

An analytic method to account for drag in the Vinti satellite theory

A quadrature algorithm is presented which employs analytical expressions for the variations of satellite orbital elements caused by air drag. The Hamiltonian is formally preserved and the Jacobi constants of the motion are advanced with time through the variational equations. The atmospheric density profile is written as a fitted exponential function of the eccentric anomaly, which adheres to tabulated data at all altitudes and simultaneously reduces the variational equations to definite integrals with closed form evaluations, whose limits are in terms of the eccentric anomaly. Results are given for two intense air drag satellites and indicate that the satellite ephemerides produced by this method in conjunction with the Vinti program are of very high accuracy.

Watson, J. S.↗
DOE OSTI2023

Operational control with set-points tuning - application to mobile robots

This paper proposes a novel method for the optimal tuning of set points for multiple-layered control system structure widely seen in robotics and other complex industrial processes composed of a number of subsystems. The terminal sliding mode control (SMC) is used as the low-level control strategy to ensure the stability of subsystems. When uncertainties exist, it can be shown that the deteriorated system performance will be improved by the outer loop with set points tuning. For this purpose, the learning of the new set point is designed to compensate for the effects caused by uncertainties during the system operation. At the same time, the system is proven to stay with the original set point when the compensation is introduced. A practical application to a holonomic mobile robot system is given to illustrate the presented method. Desired results have been obtained.

42 ENGINEERING↗
DOE OSTI2017

GEMPIC: geometric electromagnetic particle-in-cell methods

We present a novel framework for finite element particle-in-cell methods based on the discretization of the underlying Hamiltonian structure of the Vlasov–Maxwell system. We derive a semi-discrete Poisson bracket, which retains the defining properties of a bracket, anti-symmetry and the Jacobi identity, as well as conservation of its Casimir invariants, implying that the semi-discrete system is still a Hamiltonian system. In order to obtain a fully discrete Poisson integrator, the semi-discrete bracket is used in conjunction with Hamiltonian splitting methods for integration in time. Techniques from finite element exterior calculus ensure conservation of the divergence of the magnetic field and Gauss’ law as well as stability of the field solver. The resulting methods are gauge invariant, feature exact charge conservation and show excellent long-time energy and momentum behaviour. Due to the generality of our framework, these conservation properties are guaranteed independently of a particular choice of the finite element basis, as long as the corresponding finite element spaces satisfy certain compatibility conditions.

Physics↗
NASA NTRS1974

Hamilton/Jacobi perturbation methods applied to the rotational motion of a rigid body in a gravitational field

The formalism for studying perturbations of a triaxial rigid body within the Hamilton-Jacobi framework is developed. The motion of a triaxial artificial earth satellite about its center of mass is studied. Variables are found which permit separation, and the Euler angles and associated conjugate momenta are obtained as functions of canonical constants and time.

Fitzpatrick, P. M.↗
NASA NTRS1986

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗
NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗
NASA NTRS2005

Numerical Evaluation of P-Multigrid Method for the Solution of Discontinuous Galerkin Discretizations of Diffusive Equations

This paper describes numerical experiments with P-multigrid to corroborate analysis, validate the present implementation, and to examine issues that arise in the implementations of the various combinations of relaxation schemes, discretizations and P-multigrid methods. The two approaches to implement P-multigrid presented here are equivalent for most high-order discretization methods such as spectral element, SUPG, and discontinuous Galerkin applied to advection; however it is discovered that the approach that mimics the common geometric multigrid implementation is less robust, and frequently unstable when applied to discontinuous Galerkin discretizations of di usion. Gauss-Seidel relaxation converges 40% faster than block Jacobi, as predicted by analysis; however, the implementation of Gauss-Seidel is considerably more expensive that one would expect because gradients in most neighboring elements must be updated. A compromise quasi Gauss-Seidel relaxation method that evaluates the gradient in each element twice per iteration converges at rates similar to those predicted for true Gauss-Seidel.

Atkins, H. L.↗
DOE OSTI2021

The continuum and leading twist limits of parton distribution functions in lattice QCD

In this study, we present continuum limit results for the unpolarized parton distribution function of the nucleon computed in lattice QCD. This study is the first continuum limit using the pseudo-PDF approach with Short Distance Factorization for factorizing lattice QCD calculable matrix elements. Our findings are also compared with the pertinent phenomenological determinations. Inter alia, we are employing the summation Generalized Eigenvalue Problem (sGEVP) technique in order to optimize our control over the excited state contamination which can be one of the most serious systematic errors in this type of calculations. A crucial novel ingredient of our analysis is the parameterization of systematic errors using Jacobi polynomials to characterize and remove both lattice spacing and higher twist contaminations, as well as the leading twist distribution. This method can be expanded in further studies to remove all other systematic errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1996

Orbits Close to Asteroid 4769 Castalia

We use a radar-derived physical model of 4769 Castalia (1989 PB) to investigate close orbit dynamics around that kilometer- sized, uniformly rotating asteroid. Our methods of analysis provide a basis for systematic studies of particle dynamics close to any uniformly rotating asteroid. We establish that a Jacobi integral exists for particles orbiting this asteroid, examine the attendant zero-velocity surfaces, find families of periodic orbits, and determine their stability. All synchronous orbits and direct orbits within approx. 3 mean radii of Castalia are unstable and are subject to impact or escape from Castalia. Retrograde orbits are mostly stable and allow particles to orbit close to the asteroid surface. We derive a model which allows us to predict the escape conditions of a particle in orbit about Castalia and the (temporary) capture conditions for a hyperbolic interloper. Orbits within 1.5 km of Castalia are subject to immediate ejection from the system. Hyperbolic orbits with a V(sub infinity) less than 0.4 m/sec can potentially be captured by Castalia if their periapsis radius Is within approx. 2 km. For Castalia this capture region is small, but the results also apply to larger asteroids whose capture regions would also be larger. We determine bounds on ejecta speeds which either ensure ejecta escape or re-impact as functions of location on Castalia's surface. The speeds that ensure escape range from 0.28 to 0.84 m/sec and the speeds that ensure re-impact range from 0 to 0.18 m/sec. Speeds between these two bounds lead either to escape, re-impact, or potentially finite-time stable orbits. We develop a simple criterion which can establish whether a particle could have been ejected from the asteroid in the past and if it will impact the surface in the future.

Scheeres, D. J.↗
NASA NTRS1986

Fast algorithm for calculating chemical kinetics in turbulent reacting flow

This paper addresses the need for a fast batch chemistry solver to perform the kinetics part of a split operator formulation of turbulent reacting flows, with special attention focused on the solution of the ordinary differential equations governing a homogeneous gas-phase chemical reaction. For this purpose, a two-part predictor-corrector algorithm which incorporates an exponentially fitted trapezoidal method was developed. The algorithm performs filtering of ill-posed initial conditions, automatic step-size selection, and automatic selection of Jacobi-Newton or Newton-Raphson iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm, termed CREK1D (combustion reaction kinetics, one-dimensional), compared favorably with the code LSODE when tested on two representative problems drawn from combustion kinetics, and is faster than LSODE.

Radhakrishnan, K.↗
NASA NTRS1987

A full simulation of a vortex ring

A three-dimensional spectral method is developed for the solution to the incompressible Navier-Stokes equations in an unbounded domain. The spectral method relies on divergence-free basis functions as proposed by Leonard (1981). The basis functions are formed using vector spherical harmonics and Jacobi polynomials together with a mapping in the radial direction. An axisymmetric code was written and is verified using an exact solution of the Stokes equations. Preliminary results for the evolution of a vortex ring according to the Navier-Stokes equations are presented.

Stanaway, S. K.↗
DOE OSTI2014

PLANCK and WMAP constraints on generalised Hubble flow inflationary trajectories

We use the Hamilton-Jacobi formalism to constrain the space of possible single field, inflationary Hubble flow trajectories when compared to the WMAP and PLANK satellites Cosmic Microwave Background (CMB) results. This method yields posteriors on the space of Hubble Slow Roll (HSR) parameters that uniquely determine the history of the Hubble parameter during the inflating epoch. The trajectories are used to numerically determine the observable primordial power spectrum and bispectra that can then be compared to observations. Our analysis is used to infer the most likely shape of the inflaton potential V(φ) and also yields a prediction for, B, the dimensionless amplitude of the non-Gaussian bispectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems

Hamiltonian Neural Networks (HNNs) provide structure-preserving learning of Hamiltonian systems. In this paper, we extend HNNs to structure-preserving inversion of stochastic Hamiltonian systems (SHSs) from observational data. We propose the quadrature-based models according to the integral form of the SHSs’ solutions, where we denoise the loss-by-moment calculations of the solutions. The integral pattern of the models transforms the source of the essential learning error from the discrepancy between the modified Hamiltonian and the true Hamiltonian in the classical HNN models into that between the integrals and their quadrature approximations. This transforms the challenging task of deriving the relation between the modified and the true Hamiltonians from the (stochastic) Hamilton–Jacobi PDEs, into the one that only requires invoking results from the numerical quadrature theory. Meanwhile, denoising via moments calculations gives a simpler data fitting method than, e.g., via probability density fitting, which may imply better generalization ability in certain circumstances. Numerical experiments validate the proposed learning strategy on several concrete Hamiltonian systems. The experimental results show that both the learned Hamiltonian function and the predicted solution of our quadrature-based model are more accurate than that of the corrected symplectic HNN method on a harmonic oscillator, and the three-point Gaussian quadrature-based model produces higher accuracy in long-time prediction than the Kramers–Moyal method and the numerics-informed likelihood method on the stochastic Kubo oscillator as well as other two stochastic systems with non-polynomial Hamiltonian functions. Moreover, the Hamiltonian learning error εH arising from the Gaussian quadrature-based model is lower than that from Simpson’s quadrature-based model. These demonstrate the superiority of our approach in learning accuracy and long-time prediction ability compared to certain existing methods and exhibit its potential to improve learning accuracy via applying precise quadrature formulae.

Mathematics↗
NASA NTRS1988

A numerical study of viscous vortex rings using a spectral method

Viscous, axisymmetric vortex rings are investigated numerically by solving the incompressible Navier-Stokes equations using a spectral method designed for this type of flow. The results presented are axisymmetric, but the method is developed to be naturally extended to three dimensions. The spectral method relies on divergence-free basis functions. The basis functions are formed in spherical coordinates using Vector Spherical Harmonics in the angular directions, and Jacobi polynomials together with a mapping in the radial direction. Simulations are performed of a single ring over a wide range of Reynolds numbers (Re approximately equal gamma/nu), 0.001 less than or equal to 1000, and of two interacting rings. At large times, regardless of the early history of the vortex ring, it is observed that the flow approaches a Stokes solution that depends only on the total hydrodynamic impulse, which is conserved for all time. At small times, from an infinitely thin ring, the propagation speeds of vortex rings of varying Re are computed and comparisons are made with the asymptotic theory by Saffman. The results are in agreement with the theory; furthermore, the error is found to be smaller than Saffman's own estimate by a factor square root ((nu x t)/R squared) (at least for Re=0). The error also decreases with increasing Re at fixed core-to-ring radius ratio, and appears to be independent of Re as Re approaches infinity). Following a single ring, with Re=500, the vorticity contours indicate shedding of vorticity into the wake and a settling of an initially circular core to a more elliptical shape, similar to Norbury's steady inviscid vortices. Finally, we consider the case of leapfrogging vortex rings with Re=1000. The results show severe straining of the inner vortex core in the first pass and merging of the two cores during the second pass.

Stanaway, S. K.↗
DOE OSTI2017

Improvement of transport-corrected scattering stability and performance using a Jacobi inscatter algorithm for 2D-MOC

The MPACT code, being developed collaboratively by the University of Michigan and Oak Ridge National Laboratory, is the primary deterministic neutron transport solver being deployed within the Virtual Environment for Reactor Applications (VERA) as part of the Consortium for Advanced Simulation of Light Water Reactors (CASL). In many applications of the MPACT code, transport-corrected scattering has proven to be an obstacle in terms of stability, and considerable effort has been made to try to resolve the convergence issues that arise from it. Most of the convergence problems seem related to the transport-corrected cross sections, particularly when used in the 2D method of characteristics (MOC) solver, which is the focus of this work. Here in this paper, the stability and performance of the 2-D MOC solver in MPACT is evaluated for two iteration schemes: Gauss-Seidel and Jacobi. With the Gauss-Seidel approach, as the MOC solver loops over groups, it uses the flux solution from the previous group to construct the inscatter source for the next group. Alternatively, the Jacobi approach uses only the fluxes from the previous outer iteration to determine the inscatter source for each group. Consequently for the Jacobi iteration, the loop over groups can be moved from the outermost loop$-$as is the case with the Gauss-Seidel sweeper$-$to the innermost loop, allowing for a substantial increase in efficiency by minimizing the overhead of retrieving segment, region, and surface index information from the ray tracing data. Several test problems are assessed: (1) Babcock & Wilcox 1810 Core I, (2) Dimple S01A-Sq, (3) VERA Progression Problem 5a, and (4) VERA Problem 2a. The Jacobi iteration exhibits better stability than Gauss-Seidel, allowing for converged solutions to be obtained over a much wider range of iteration control parameters. Additionally, the MOC solve time with the Jacobi approach is roughly 2.0-2.5× faster per sweep. While the performance and stability of the Jacobi iteration are substantially improved compared to the Gauss-Seidel iteration, it does yield a roughly 8$-$10% increase in the overall memory requirement.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

The invariance of the restricted problem of three bodies

The restricted problem of three bodies in sidereal coordinates is investigated using Lie theoretical methods of extended groups. The application of the theory results in a single scalar generator with an associated invariant - the integral of the motion attributed to Jacobi.

Jezewski, D. J.↗
DOE OSTI2016

A Piecewise Deterministic Markov Toy Model for Traffic/Maintenance and Associated Hamilton–Jacobi Integrodifferential Systems on Networks

We study optimal control problems in infinite horizon whxen the dynamics belong to a specific class of piecewise deterministic Markov processes constrained to star-shaped networks (corresponding to a toy traffic model). We adapt the results in Soner (SIAM J Control Optim 24(6):1110–1122, 1986) to prove the regularity of the value function and the dynamic programming principle. Extending the networks and Krylov’s “shaking the coefficients” method, we prove that the value function can be seen as the solution to a linearized optimization problem set on a convenient set of probability measures. The approach relies entirely on viscosity arguments. As a by-product, the dual formulation guarantees that the value function is the pointwise supremum over regular subsolutions of the associated Hamilton–Jacobi integrodifferential system. This ensures that the value function satisfies Perron’s preconization for the (unique) candidate to viscosity solution.

97 MATHEMATICS AND COMPUTING↗