Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The planetary nebula system and dynamics of NGC 5128. I - Planetary nebulae as standard candles

We present the result of a planetary nebula (PN) survey of the nearby giant elliptical galaxy NGC 5128 performed with CCD cameras at the prime focus of the CTIO 4 m telescope. By comparing CCD images centered on the characteristic emission line forbidden O III 5007 A and on the adjacent continuum, we identify a total of 785 PNs in areas extending 20 kpc along the photometric major axis and covering the whole galaxy to 10 kpc. From these data, we form a complete sample of 224 PNs extending to a dereddened limiting magnitude of m5007 = 24.8, which extends 1.5 mag down the PN luminosity function (PNLF). Adopting a foreground extinction of E(B-V) = 0.1, we derive a distance to the galaxy of 3.5 +/- 0.2 Mpc, in excellent agreement with the surface brightness fluctuation method. No population effect on the bright cutoff of PNLF is observed in the isophotal radius range of 2-16 kpc, but the luminosity specific PN density seems to increase with radius inside of 7 kpc, in agreement with the alpha(2.5)-color relation observed for other galaxies.

Hui, Xiaohui↗
NASA NTRS1985

Secondary pattern computation of an arbitrarily shaped main reflector

The secondary pattern of a perfectly conducting offset main reflector being illuminated by a point feed at an arbitrary location is studied. The method of analysis is based upon the application of the Fast Fourier Transform (FFT) to the aperture fields obtained using geometrical optics (GO) and geometrical theory of diffraction (GTD). Key features of the present work are (1) the reflector surface is completely arbitrary, (2) the incident field from the feed is most general with arbitrary polarization and location, and (3) the edge diffraction is calculated by either UAT or by UTD. Comparison of this technique for an offset parabolic reflector with the Jacobi-Bessel and Fourier-Bessel techniques shows good agreement. Near field, far field, and scan data of a large refelctor are presented.

Lee, S. W.↗
NASA NTRS2007

Experimental Design for the Evaluation of Detection Techniques of Hidden Corrosion Beneath the Thermal Protective System of the Space Shuttle Orbiter

The detection of corrosion beneath Space Shuttle Orbiter thermal protective system is traditionally accomplished by removing the Reusable Surface Insulation tiles and performing a visual inspection of the aluminum substrate and corrosion protection system. This process is time consuming and has the potential to damage high cost tiles. To evaluate non-intrusive NDE methods, a Proof of Concept (PoC) experiment was designed and test panels were manufactured. The objective of the test plan was three-fold: establish the ability to detect corrosion hidden from view by tiles; determine the key factor affecting detectability; roughly quantify the detection threshold. The plan consisted of artificially inducing dimensionally controlled corrosion spots in two panels and rebonding tile over the spots to model the thermal protective system of the orbiter. The corrosion spot diameter ranged from 0.100" to 0.600" inches and the depth ranged from 0.003" to 0.020". One panel consisted of a complete factorial array of corrosion spots with and without tile coverage. The second panel consisted of randomized factorial points replicated and hidden by tile. Conventional methods such as ultrasonics, infrared, eddy current and microwave methods have shortcomings. Ultrasonics and IR cannot sufficiently penetrate the tiles, while eddy current and microwaves have inadequate resolution. As such, the panels were interrogated using Backscatter Radiography and Terahertz Imaging. The terahertz system successfully detected artificially induced corrosion spots under orbiter tile and functional testing is in-work in preparation for implementation.

Kemmerer, Catherine C.↗
DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.

abundances↗
DOE OSTI2022

Smoothed aggregation for difficult stretched mesh and coefficient variation problems

Abstract Four adaptations of the smoothed aggregation algebraic multigrid (SA‐AMG) method are proposed with an eye toward improving the convergence and robustness of the solver in situations when the discretization matrix contains many weak connections. These weak connections can cause higher than expected levels of fill‐in within the coarse discretization matrices and can also give rise to suboptimal smoothing within the prolongator smoothing phase. These smoothing drawbacks are due to the relatively small size of some diagonal entries within the filtered matrix that one obtains after dropping the weak connections. The new algorithms consider modifications to the Jacobi‐like step that defines the prolongator smoother, modifications to the filtered matrix, and also direct modifications to the resulting grid transfer operators. Numerical results are given illustrating the potential benefits of the proposed adaptations.

Hu, Jonathan J.↗
NASA NTRS1981

A new look at Fresnel field computation using the Jacobi-Bessel series

Computational procedures that would be useful in finding the Fresnel field from a knowledge of the Jacobi-Bessel expansion of the far field are considered. The range of validity of the Fresnel approximation is carefully examined by comparing it with the exact closed form solution for the uniform circular aperture. Also investigated numerically, and in great detail, is the range of validity (over theta) of the Fresnel small angle (FSA) approximation. For moderate sized apertures as small as 10 wavelengths, it is found that the FSA approximation is very accurate to angles as wide as four or more sidelobes (as seen in the far zone). A very efficient computational method is shown to exist for the radiation integral in the form of a single series expansion that is analytically continuous and convergent for a wide range of observation points in three-dimensional space.

Galindo-Israel, V.↗
DOE OSTI2013

YNOGK: A NEW PUBLIC CODE FOR CALCULATING NULL GEODESICS IN THE KERR SPACETIME

Following the work of Dexter and Agol, we present a new public code for the fast calculation of null geodesics in the Kerr spacetime. Using Weierstrass's and Jacobi's elliptic functions, we express all coordinates and affine parameters as analytical and numerical functions of a parameter p, which is an integral value along the geodesic. This is the main difference between our code and previous similar ones. The advantage of this treatment is that the information about the turning points does not need to be specified in advance by the user, and many applications such as imaging, the calculation of line profiles, and the observer-emitter problem, become root-finding problems. All elliptic integrations are computed by Carlson's elliptic integral method as in Dexter and Agol, which guarantees the fast computational speed of our code. The formulae to compute the constants of motion given by Cunningham and Bardeen have been extended, which allow one to readily handle situations in which the emitter or the observer has an arbitrary distance from, and motion state with respect to, the central compact object. The validation of the code has been extensively tested through applications to toy problems from the literature. The source FORTRAN code is freely available for download on our Web site http://www1.ynao.ac.cn/{approx}yangxl/yxl.html.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

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↗
DOE OSTI2016

MOC Efficiency Improvements Using a Jacobi Inscatter Approximation

In recent weeks, attention has been given to resolving the convergence issues encountered with TCP 0 by trying a Jacobi (J) inscatter approach when group sweeping, where the inscatter source is constructed using the previous iteration flux. This is in contrast to a Gauss-Seidel (GS) approach, which has been the default to-date, where the scattering source uses the most up-to-date flux values. The former is consistent with CASMO, which has no issues with TCP 0 convergence. Testing this out on a variety of problems has demonstrated that the Jacobi approach does indeed provide substantially more stability, though can take more outer iterations to converge. While this is not surprising, there are improvements that can be made to the MOC sweeper to capitalize on the Jacobi approximation and provide substantial speedup. For example, the loop over groups, which has traditionally been the outermost loop in MPACT, can be moved to the interior, avoiding duplicate modular ray trace and coarse ray trace setup (mapping coarse mesh surface indexes), which needs to be performed repeatedly when group is outermost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

On an asymptotic viscosity solution property of solutions of discrete Hamilton–Jacobi–Bellman equations

In this paper, we show that a proper limit of solutions of discrete Hamilton–Jacobi–Bellman (dHJB) equations in a random walk model becomes a viscosity solution of a Hamilton–Jacobi–Bellman (HJB) variational inequality in a continuous-time geometric Brownian model. HJB variational inequalities are used to analyze singular stochastic control problems in mathematical finance. By our result, with the help of dHJB equations, we can obtain viscosity solutions of HJB variational inequalities which are usually identified with the value functions of the singular stochastic control problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Generating functionals and Lagrangian partial differential equations

The main goal of this paper is to derive an alternative characterization of the multisymplectic form formula for classical field theories using the geometry of the space of boundary values. We review the concept of Type-I/II generating functionals defined on the space of boundary data of a Lagrangian field theory. On the Lagrangian side, we define an analogue of Jacobi's solution to the Hamilton–Jacobi equation for field theories, and we show that by taking variational derivatives of this functional, we obtain an isotropic submanifold of the space of Cauchy data, described by the so-called multisymplectic form formula. As an example of the latter, we show that Lorentz's reciprocity principle in electromagnetism is a particular instance of the multisymplectic form formula. We also define a Hamiltonian analogue of Jacobi's solution, and we show that this functional is a Type-II generating functional. We finish the paper by defining a similar framework of generating functions for discrete field theories, and we show that for the linear wave equation, we recover the multisymplectic conservation law of Bridges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Using virtual menus in a virtual environment

Virtual environment interfaces to computer programs in several diverse application areas are currently being developed. The users of virtual environments will require many different methods to interact with the environments and the objects in them. This paper reports on our use of virtual menus as a method of interacting with virtual environments. Several aspects of virtual environments make menu interactions different from interactions with conventional menus. We review the relevant aspects of conventional menus and virtual environments, in order to provide a frame of reference for the design of virtual menus. We discuss the features and interaction methodologies of two different versions of virtual menus which have been developed and used in our lab. We also examine the problems associated with our original version, and the enhancements incorporated into our current version.

Jacoby, Richard H.↗
DOE OSTI2018

Fundamental Solutions for Two-Point Boundary Value Problems in Orbital Mechanics

We consider a two-point boundary value problem (TPBVP) in orbital mechanics involving a small body (e.g., a spacecraft or asteroid) and N larger bodies. The least action principle TPBVP formulation is converted into an initial value problem via the addition of an appropriate terminal cost to the action functional. The latter formulation is used to obtain a fundamental solution, which may be used to solve the TPBVP for a variety of boundary conditions within a certain class. In particular, the method of convex duality allows one to interpret the least action principle as a differential game, where an opposing player maximizes over an indexed set of quadratics to yield the gravitational potential. In the case where the time duration is less than a specific bound, there exists a unique critical point for the resulting differential game, which yields the fundamental solution given in terms of the solutions of associated Riccati equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1998

Sensor Spatial Distortion, Visual Latency, and Update Rate Effects on 3D Tracking in Virtual Environments

Several common defects that we have sought to minimize in immersing virtual environments are: static sensor spatial distortion, visual latency, and low update rates. Human performance within our environments during large amplitude 3D tracking was assessed by objective and subjective methods in the presence and absence of these defects. Results show that 1) removal of our relatively small spatial sensor distortion had minor effects on the tracking activity, 2) an Adapted Cooper-Harper controllability scale proved the most sensitive subjective indicator of the degradation of dynamic fidelity caused by increasing latency and decreasing frame rates, and 3) performance, as measured by normalized RMS tracking error or subjective impressions, was more markedly influenced by changing visual latency than by update rate.

Ellis, S. R.↗
NASA NTRS1988

Hypercube matrix computation task

A major objective of the Hypercube Matrix Computation effort at the Jet Propulsion Laboratory (JPL) is to investigate the applicability of a parallel computing architecture to the solution of large-scale electromagnetic scattering problems. Three scattering analysis codes are being implemented and assessed on a JPL/California Institute of Technology (Caltech) Mark 3 Hypercube. The codes, which utilize different underlying algorithms, give a means of evaluating the general applicability of this parallel architecture. The three analysis codes being implemented are a frequency domain method of moments code, a time domain finite difference code, and a frequency domain finite elements code. These analysis capabilities are being integrated into an electromagnetics interactive analysis workstation which can serve as a design tool for the construction of antennas and other radiating or scattering structures. The first two years of work on the Hypercube Matrix Computation effort is summarized. It includes both new developments and results as well as work previously reported in the Hypercube Matrix Computation Task: Final Report for 1986 to 1987 (JPL Publication 87-18).

Calalo, Ruel H.↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Approximate optimal guidance for the advanced launch system

A real-time guidance scheme for the problem of maximizing the payload into orbit subject to the equations of motion for a rocket over a spherical, non-rotating earth is presented. An approximate optimal launch guidance law is developed based upon an asymptotic expansion of the Hamilton - Jacobi - Bellman or dynamic programming equation. The expansion is performed in terms of a small parameter, which is used to separate the dynamics of the problem into primary and perturbation dynamics. For the zeroth-order problem the small parameter is set to zero and a closed-form solution to the zeroth-order expansion term of Hamilton - Jacobi - Bellman equation is obtained. Higher-order terms of the expansion include the effects of the neglected perturbation dynamics. These higher-order terms are determined from the solution of first-order linear partial differential equations requiring only the evaluation of quadratures. This technique is preferred as a real-time, on-line guidance scheme to alternative numerical iterative optimization schemes because of the unreliable convergence properties of these iterative guidance schemes and because the quadratures needed for the approximate optimal guidance law can be performed rapidly and by parallel processing. Even if the approximate solution is not nearly optimal, when using this technique the zeroth-order solution always provides a path which satisfies the terminal constraints. Results for two-degree-of-freedom simulations are presented for the simplified problem of flight in the equatorial plane and compared to the guidance scheme generated by the shooting method which is an iterative second-order technique.

Feeley, T. S.↗
DOE OSTI2019

Analog of Hamilton-Jacobi theory for the time-evolution operator

In this paper we develop an analog of Hamilton-Jacobi theory for the time-evolution operator of a quantum many-particle system. The theory offers a useful approach to develop approximations to the time-evolution operator, and also provides a unified framework and starting point for many well-known approximations to the time-evolution operator. In the important special case of periodically driven systems at stroboscopic times, we find relatively simple equations for the coupling constants of the Floquet Hamiltonian, where a straightforward truncation of the couplings leads to a powerful class of approximations. Using our theory, we construct a flow chart that illustrates the connection between various common approximations, which also highlights some missing connections and associated approximation schemes. These missing connections turn out to imply an analytically accessible approximation that is the “inverse” of a rotating frame approximation and thus has a range of validity complementary to it. Here, we numerically test the various methods on the one-dimensional Ising model to confirm the ranges of validity that one would expect from the approximations used. The theory provides a map of the relations between the growing number of approximations for the time-evolution operator. We describe these relations in a table showing the limitations and advantages of many common approximations, as well as the approximations introduced in this paper.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗