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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Complexity of Kronecker Operations on Sparse Matrices with Applications to the Solution of Markov Models

We present a systematic discussion of algorithms to multiply a vector by a matrix expressed as the Kronecker product of sparse matrices, extending previous work in a unified notational framework. Then, we use our results to define new algorithms for the solution of large structured Markov models. In addition to a comprehensive overview of existing approaches, we give new results with respect to: (1) managing certain types of state-dependent behavior without incurring extra cost; (2) supporting both Jacobi-style and Gauss-Seidel-style methods by appropriate multiplication algorithms; (3) speeding up algorithms that consider probability vectors of size equal to the "actual" state space instead of the "potential" state space.

Buchholz, Peter↗
NASA NTRS1987

Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-Jacobi formulations

New numerical algorithms are devised (PSC algorithms) for following fronts propagating with curvature-dependent speed. The speed may be an arbitrary function of curvature, and the front can also be passively advected by an underlying flow. These algorithms approximate the equations of motion, which resemble Hamilton-Jacobi equations with parabolic right-hand-sides, by using techniques from the hyperbolic conservation laws. Non-oscillatory schemes of various orders of accuracy are used to solve the equations, providing methods that accurately capture the formation of sharp gradients and cusps in the moving fronts. The algorithms handle topological merging and breaking naturally, work in any number of space dimensions, and do not require that the moving surface be written as a function. The methods can be used also for more general Hamilton-Jacobi-type problems. The algorithms are demonstrated by computing the solution to a variety of surface motion problems.

Osher, Stanley↗
DOE OSTI2016

State-resolved differential and integral cross sections for the Ne + H{sub 2}{sup +} (v = 0–2, j = 0) → NeH{sup +} + H reaction

State-to-state quantum dynamic calculations for the proton transfer reaction Ne + H{sub 2}{sup +} (v = 0–2, j = 0) are performed on the most accurate LZHH potential energy surface, with the product Jacobi coordinate based time-dependent wave packet method including the Coriolis coupling. The J = 0 reaction probabilities for the title reaction agree well with previous results in a wide range of collision energy of 0.2-1.2 eV. Total integral cross sections are in reasonable agreement with the available experiment data. Vibrational excitation of the reactant is much more efficient in enhancing the reaction cross sections than translational and rotational excitation. Total differential cross sections are found to be forward-backward peaked with strong oscillations, which is the indication of the complex-forming mechanism. As the collision energy increases, state-resolved differential cross section changes from forward-backward symmetric peaked to forward scattering biased. This forward bias can be attributed to the larger J partial waves, which makes the reaction like an abstraction process. Differential cross sections summed over two different sets of J partial waves for the v = 0 reaction at the collision energy of 1.2 eV are plotted to illustrate the importance of large J partial waves in the forward bias of the differential cross sections.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Non-gaussian signatures of general inflationary trajectories

We carry out a numerical calculation of the bispectrum in generalised trajectories of canonical, single-field inflation. The trajectories are generated in the Hamilton-Jacobi (HJ) formalism based on Hubble Slow Roll (HSR) parameters. The calculation allows generally shape and scale dependent bispectra, or dimensionless f{sub NL}, in the out-of-slow-roll regime. The distributions of f{sub NL} for various shapes and HSR proposals are shown as an example of how this procedure can be used within the context of Monte Carlo exploration of inflationary trajectories. We also show how allowing out-of-slow-roll behaviour can lead to a bispectrum that is relatively large for equilateral shapes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Numerical simulation of boundary-layer transition

The transition to turbulence in boundary layers was investigated by direct numerical solution of the nonlinear, three-dimensional, incompressible Navier-Stokes equations in the half-infinite domain over a flat plate. Periodicity was imposed in the streamwise and spanwise directions. A body force was applied to approximate the effect of a nonparallel mean flow. The numerical method was spectra, based on Fourier series and Jacobi polynomials, and used divergence-free basis functions. Extremely rapid convergence was obtained when solving the linear Orr-Sommerfeld equation. The early nonlinear and three-dimensional stages of transition, in a boundary layer disturbed by a vibrating ribbon, were successfully simulated. Excellent qualitative agreement was observed with either experiments or weakly nonlinear theories. In particular, the breakdown pattern was staggered or nonstaggered depending on the disturbance amplitude.

Spalart, P. R.↗
NASA NTRS1989

Nonlinear longitudinal control of a supermaneuverable aircraft

A technique is described which can be used for design of feedback controllers for high-performance aircraft operating in flight conditions in which nonlinearities significantly affect performance. Designs are performed on a mathematical model of the longitudinal dynamics of a hypothetical aircraft similar to proposed supermaneuverable flight test vehicles. Nonlinear controller designs are performed using truncated solutions of the Hamilton-Jacobi-Bellman equation. Preliminary results show that the method yields promising results.

Garrard, William L.↗
DOE OSTI2025

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗
DOE OSTI2014

Physics Integration KErnels (PIKE)

Pike is a software library for coupling and solving multiphysics applications. It provides basic interfaces and utilities for performing code-to-code coupling. It provides simple "black-box" Picard iteration methods for solving the coupled system of equations including Jacobi and Gauss-Seidel solvers. Pike was developed originally to couple neutronics and thermal fluids codes to simulate a light water nuclear reactor for the Consortium for Simulation of Light-water Reactors (CASL) DOE Energy Innovation Hub. The Pike library contains no physics and just provides interfaces and utilities for coupling codes. It will be released open source under a BSD license as part of the Trilinos solver framework (trilinos.org) which is also BSD. This code provides capabilities similar to other open source multiphysics coupling libraries such as LIME, AMP, and MOOSE.

Pawlowski, Roger↗
NASA NTRS1988

Problems in fluid dynamics

A scheme was developed for the parametric differentiation and integration of gas dynamics equations. A numerical integration of the gas dynamics equations is necessarily performed for a specific set of parameter values. The linear variational equations are obtained by differentiating the exact equations with respect to each of the relevant parameters. The resulting matrix of flow quantities is referred to as the Jacobi matrix. The subsequent procedure is then straightforward. The method was tested for two dimensional supersonic flow past an airfoil, with airfoil thickness, camber, and angle of attack varied. This approach has great potential value for rapidly assessing the effect of design changes. The other focus of the work was on problems in fluid stability, bifurcations, and turbulence.

Sirovich, Lawrence↗
DOE OSTI2014

Physics Integration KErnels (PIKE)

Pike is a software library for coupling and solving multiphysics applications. It provides basic interfaces and utilities for performing code-to-code coupling. It provides simple "black-box" Picard iteration methods for solving the coupled system of equations including Jacobi and Gauss-Seidel solvers. Pike was developed originally to couple neutronics and thermal fluids codes to simulate a light water nuclear reactor for the Consortium for Simulation of Light-water Reactors (CASL) DOE Energy Innovation Hub. The Pike library contains no physics and just provides interfaces and utilities for coupling codes. It will be released open source under a BSD license as part of the Trilinos solver framework (trilinos.org) which is also BSD. This code provides capabilities similar to other open source multiphysics coupling libraries such as LIME, AMP, and MOOSE.

Pawlowski, Roger↗
DOE OSTI2018

An updated Lagrangian discontinuous Galerkin hydrodynamic method for gas dynamics

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for gas dynamics. The new method evolves conserved unknowns in the current configuration, which obviates the Jacobi matrix that maps the element in a reference coordinate system or the initial coordinate system to the current configuration. The density, momentum, and total energy (ρ, ρu, E ) are approximated with conservative higher-order Taylor expansions over the element and are limited toward a piecewise constant field near discontinuities using a limiter. Two new limiting methods are presented for enforcing the bounds on the primitive variables of density, velocity, and specific internal energy (ρ, u, e). The nodal velocity, and the corresponding forces, are calculated by solving an approximate Riemann problem at the element nodes. An explicit second-order method is used to temporally advance the solution. This new Lagrangian DG hydrodynamic method conserves mass, momentum, and total energy. 1D Cartesian coordinates test problem results are presented to demonstrate the accuracy and convergence order of the new DG method with the new limiters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Hamilton–Jacobi analysis for three dimensional gravity without dynamics

Highlights: •We report the HJ. analysis for three dimensional gravity without dynamics. •We report the complete set of Hamiltonians. •The characteristic equations are reported. •The canonical covariant method is performed. •A closed and gauge invariant symplectic structure is constructed. -- Abstract: The Hamilton–Jacobi analysis for gravity without dynamics is performed. We report a detailed analysis where the complete set of Hamilton–Jacobi constraints, the characteristic equations and the gauge transformations of the theory are found. We compare our results with those reported in the literature where alternative approaches are used. In addition, we complete our work by performing the canonical covariant analysis by constructing a gauge invariant symplectic structure, and we find a full consistency between the results obtained from both approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Full-sky formulae for weak lensing power spectra from total angular momentum method

We systematically derive full-sky formulae for the weak lensing power spectra generated by scalar, vector and tensor perturbations from the total angular momentum (TAM) method. Based on both the geodesic and geodesic deviation equations, we first give the gauge-invariant expressions for the deflection angle and Jacobi map as observables of the CMB lensing and cosmic shear experiments. We then apply the TAM method, originally developed in the theoretical studies of CMB, to a systematic derivation of the angular power spectra. The TAM representation, which characterizes the total angular dependence of the spatial modes projected along a line-of-sight, can carry all the information of the lensing modes generated by scalar, vector, and tensor metric perturbations. This greatly simplifies the calculation, and we present a complete set of the full-sky formulae for angular power spectra in both the E-/B-mode cosmic shear and gradient-/curl-mode lensing potential of deflection angle. Based on the formulae, we give illustrative examples of non-vanishing B-mode cosmic shear and curl-mode of deflection angle in the presence of the vector and tensor perturbations, and explicitly compute the power spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Large Deviations and Importance Sampling for Systems of Slow-Fast Motion

In this paper we develop the large deviations principle and a rigorous mathematical framework for asymptotically efficient importance sampling schemes for general, fully dependent systems of stochastic differential equations of slow and fast motion with small noise in the slow component. We assume periodicity with respect to the fast component. Depending on the interaction of the fast scale with the smallness of the noise, we get different behavior. We examine how one range of interaction differs from the other one both for the large deviations and for the importance sampling. We use the large deviations results to identify asymptotically optimal importance sampling schemes in each case. Standard Monte Carlo schemes perform poorly in the small noise limit. In the presence of multiscale aspects one faces additional difficulties and straightforward adaptation of importance sampling schemes for standard small noise diffusions will not produce efficient schemes. It turns out that one has to consider the so called cell problem from the homogenization theory for Hamilton-Jacobi-Bellman equations in order to guarantee asymptotic optimality. We use stochastic control arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Variable-Position Acoustic Levitation

Method of acoustic levitation supports objects at positions other than acoustic nodes. Acoustic force is varied so it balances gravitational (or other) force, thereby maintaining object at any position within equilibrium range. Levitation method applicable to containerless processing. Such objects as table-tennis balls, hollow plastic spheres, and balsa-wood spheres levitated in laboratory by new method.

Barmatz, M. B.↗
NASA NTRS2013

Robust Maneuvering Envelope Estimation Based on Reachability Analysis in an Optimal Control Formulation

This paper discusses an algorithm for estimating the safe maneuvering envelope of damaged aircraft. The algorithm performs a robust reachability analysis through an optimal control formulation while making use of time scale separation and taking into account uncertainties in the aerodynamic derivatives. Starting with an optimal control formulation, the optimization problem can be rewritten as a Hamilton- Jacobi-Bellman equation. This equation can be solved by level set methods. This approach has been applied on an aircraft example involving structural airframe damage. Monte Carlo validation tests have confirmed that this approach is successful in estimating the safe maneuvering envelope for damaged aircraft.

optimal control↗