Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JACOBI INTEGRAL”

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 55 records · Page 3

Separation of variables in the special diagonal Hamilton-Jacobi equation: Application to the dynamical problem of a particle constrained on a moving surface

For a time-dependent, n-dimensional, special diagonal Hamilton-Jacobi equation a necessary and sufficient condition for the separation of variables to yield a complete integral of the form was established by specifying the admissible forms in terms of arbitrary functions. A complete integral was then expressed in terms of these arbitrary functions and also the n irreducible constants. As an application of the results obtained for the two-dimensional Hamilton-Jacobi equation, analysis was made for a comparatively wide class of dynamical problems involving a particle moving in Euclidean three-dimensional space under the action of external forces but constrained on a moving surface. All the possible cases in which this equation had a complete integral of the form were obtained and these are tubulated for reference.

Blanchard, D. L.↗

A guidance law for the aeroassisted plane change maneuver in the presence of atmospheric uncertainties

A stochastic feedback control law for a space vehicle performing an aeroassisted plane-change maneuver is developed. The stochastic control law is designed to minimize the energy loss while taking into consideration the uncertainty in the atmospheric density. The solution is based on expansion of the stochastic Hamilton-Jacobi-Bellman equation (or dynamic programming) about a zeroth-order known integrable solution. The resulting guidance law is expressed as a series expansion in the noise power spectral densities. A numerical example indicates the potential improvement of this method.

Mishne, D.↗

Optimal control of aeroassisted plane change maneuver using feedback expansions

A guidance law for an aeroassisted plane change maneuver is developed by an asymptotic expansion technique using a small parameter which essentially represents the ratio of the inertial forces to the atmospheric forces. This guidance law minimizes the energy loss while meeting terminal constraints on the altitude, flight path angle, and heading angle. By neglecting the inertial forces, the resulting optimization problem is integrable and can be determined in closed form. This zeroth-order solution is the first term in an asymptotic series solution of the Hamilton-Jacobi-Bellman equation. The remaining terms are determined from the solution of a first-order, linear partial differential equation whose solution requires only quadrature integration. Our initial results in using this guidance scheme are encouraging.

Mishne, D.↗

On the alleged collisional origin of the Kirkwood Gaps

This paper examines two proposed mechanisms whereby asteroidal collisions and close approaches may have given rise to the Kirkwood Gaps. The first hypothesis is that asteroids in near-resonant orbits have markedly increased collision probabilities and so are preferentially destroyed, or suffer decay in population density, within the resonance zones. A simple order-of-magnitude analysis shows that this hypothesis is untenable since it leads to conclusions which are either unrealistic or not in accord with present understanding of asteroidal physics. The second hypothesis is the Brouwer-Jefferys theory that collisions would smooth an asteroidal distribution function, as a function of Jacobi constant, thus forming resonance gaps. This hypothesis is examined by direct numerical integration of 50 asteroid orbits near the 2:1 resonance, with collisions simulated by random variables. No tendency to form a gap was observed.

Heppenheimer, T. A.↗

Asteroid collisional history - Effects on sizes and spins

The effects of asteroid collisional history on sizes and spins of present-day objects are discussed. Collisional evolution studies indicate that collisions have altered the spin-rates of small bodies, but that the largest asteroids may have retained their primordial rotation rates. Most asteroids larger than 100 km diam have probably been shattered, but have gravitationally recaptured their fragments to form a rubble-pile structure. Large angular momentum asteroids appear to have Maclaurian spheroidal or Jacobi-ellipsoid-like shapes; some of them may have fissioned into binaries. An integrated size and spin collisional evolution model is presented, with two critical parameters: one which determines the spin rates for small fragments resulting from a shattering collision, and the other determines the fraction of impact angular momentum that is retained by the target.

Davis, Donald R.↗

Redundancy of constraints in the classical and quantum theories of gravitation.

It is shown that in Dirac's version of the quantum theory of gravitation, the Hamiltonian constraints are greatly redundant. If the Hamiltonian constraint condition is satisfied at one point on the underlying, closed three-dimensional manifold, then it is automatically satisfied at every point, provided only that the momentum constraints are everywhere satisfied. This permits one to replace the usual infinity of Hamiltonian constraints by a single condition which may be taken in the form of an integral over the manifold. Analogous theorems are given for the classical Einstein Hamilton-Jacobi equations.

Moncrief, V.↗

An efficient technique for the computation of vector secondary patterns of offset paraboloid reflectors

A series approach for the rapid computation of the vector secondary pattern of offset paraboloid reflectors wherein the feed is displaced is presented. We show that the Jacobi polynomial series method, which has been demonstrated to provide an efficient means for evaluating the radiation integral of symmetric paraboloid reflectors, can be extended to the case of an offset paraboloid without compromising the ease or speed of computation. The analysis leading to the series formula is also useful for deriving an analytic expression for the optimum scan plane for the displacement of the feed. Representative numerical results illustrating the application of the method and the properties of the offset paraboloid are presented.

Mittra, R.↗

Capture of Planetesimals into a Circumterrestrial Swarm

The lunar origin model considered in this report involves processing of protolunar material through a circumterrestrial swarm of particles. Once such a swarm has formed, it can gain mass by capturing infalling planetesimals and ejecta from giant impacts on the Earth, although the angular momentum supply from these sources remains a problem. The first stage of formation of a geocentric swarm by capture of planetesimals from initially heliocentric orbits is examined. The only plausible capture mechanism that is not dependent on very low approach velocities is the mutual collision of planetesimals passing within Earth's sphere of influence. The dissipation of energy in inelastic collisions or accretion events changes the value of the Jacobi parameter, allowing capture into bound geocentric orbits. This capture scenario was tested directly by many body numerical integration of planetesimal orbits in near Earth space.

Weidenschilling, S. J.↗

Semiclassical theory and the Koopman-van Hove equation *

Abstract The phase space Koopman-van Hove (KvH) equation can be derived from the asymptotic semiclassical analysis of partial differential equations. Semiclassical theory yields the Hamilton–Jacobi equation for the complex phase factor and the transport equation for the amplitude. These two equations can be combined to form a nonlinear semiclassical version of the KvH equation in configuration space. There is a natural injection of configuration space solutions into phase space and a natural projection of phase space solutions onto configuration space. Hence, every solution of the configuration space KvH equation satisfies both the semiclassical phase space KvH equation and the Hamilton–Jacobi constraint. For configuration space solutions, this constraint resolves the paradox that there are two different conserved densities in phase space. For integrable systems, the KvH spectrum is the Cartesian product of a classical and a semiclassical spectrum. If the classical spectrum is eliminated, then, with the correct choice of Jeffreys–Wentzel–Kramers–Brillouin (JWKB) matching conditions, the semiclassical spectrum satisfies the Einstein–Brillouin–Keller quantization conditions which include the correction due to the Maslov index. However, semiclassical analysis uses different choices for boundary conditions, continuity requirements, and the domain of definition. For example, use of the complex JWKB method allows for the treatment of tunneling through the complexification of phase space. Finally, although KvH wavefunctions include the possibility of interference effects, interference is not observable when all observables are approximated as local operators on phase space. Observing interference effects requires consideration of nonlocal operations, e.g. through higher orders in the asymptotic theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Faster Computation Of Far Fields Of Dish Antennas

Modified Zernike-polynomial approach reduces number of terms in numerical integration. Yields greater or equal accuracy with fewer terms in series expansion and fewer sampling points across aperture and, less computation. Method based on Jacobi-Bessel-expansion concept.

Jamnejad, Vahraz↗

Fracture Mechanics Analysis of an Annular Crack in a Three-concentric-cylinder Composite Model

A boundary-value problem governing a three-phase concentric-cylinder model was analytically modeled to analyze annular interfacial crack problems with Love s strain functions in order to find the stress intensity factors (SIFs) and strain energy release rates (SERRs) at the tips of an interface crack in a nonhomogeneous medium. The complex form of a singular integral equation (SIE) of the second kind was formulated using Bessel s functions in the Fourier domain, and the SIF and total SERR were calculated using Jacoby polynomials. For the validity of the SIF equations to be established, the SIE of the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the results were compared with the previous results of Erdogan. A preliminary set of parametric studies was carried out to show the effect of interphase properties on the SERR. The method presented here provides insight about the effect of interphase properties on the crack driving force.

Kuguoglu, Latife H.↗

Optimal transfer between close near-circular orbits

A nonlinear analysis is carried out for the optimal transfer between close near-circular noncoplanar orbits with no limit on time. The degeneracy of the linear analyses of Marec and Edelbaum is, thereby, removed and with the aid of a generalized Jacobi test, it is established that the minimum-impulse transfer requires, at most, three impulses. A field of extremals is generated in an orbit element space starting with the initial state and anticipated adjoint variables and integrating these variables by using the optimal thrust control law that maximizes the variational Hamiltonian. The nonoptimality of any four-impulse strategy is determined numerically by observing the reflection of every extremal from either an envelope (conjugate point test) or a switching surface on or before the application of the fourth impulse. This test is carried out in Regions I, II, and III of Marec, as well as the special transition region of Breakwell near the boundary of Region III.

Kechichian, J.↗

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

An approximate atmospheric guidance law for aeroassisted plane change maneuvers

An approximate optimal guidance law for the aeroassisted plane change problem is presented which is based upon an expansion of the Hamilton-Jacobi-Bellman equation with respect to the small parameter of Breakwell et al. (1985). The present law maximizes the final velocity of the reentry vehicle while meeting terminal constraints on altitude, flight path angle, and heading angle. The integrable zeroth-order solution found when the small parameter is set to zero corresponds to a solution of the problem where the aerodynamic forces dominate the inertial forces. Higher order solutions in the expansion are obtained from the solution of linear partial differential equations requiring only quadrature integration.

Speyer, Jason L.↗

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗