Search NASA⌕ Search

SEARCH · Search NASA

Results for “BESSEL FUNCTION”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Dynamics of Flexible Spinning Satellites with Radial Wire Antennas

A dynamic analysis is presented for a spin stabilized spacecraft employing four radial wire antennas with tip masses, a configuration first employed in the IMP-J spacecraft. The use of wires in place of the usual booms represents the ultimate in weight reduction at the expanse of flexibility. The satellite is modelled as a 14 degree of freedom system, and the linearized equations of motion are found. The lowest order vibrational modes and natural frequencies of the gyroscopically coupled system are then determined. Because the satellite spin rate is decreased by antenna deployment, a spin-up maneuver is needed. The response of the time varying mode equations during spin-up is found, for the planar modes, in terms of Bessel functions and a Struve function of order -1/4. Because tables of the latter are not readily available, the particular solution is expressed in various forms including an infinite series of Bessel functions and a particularly useful asymptotic expansion.

Longman, R. W.↗
NASA NTRS1992

Flare gamma ray continuum emission from neutral pion decay

We investigate, in detail, the production of solar flare gamma ray emission above 100 MeV via the interaction of high energy protons with the ambient solar atmosphere. We restrict our considerations to the broadband gamma ray spectrum resulting from the decay of neutral pions produced in p-H reactions. Thick-target calculations are performed to determine the photon fluences. However, proton transport is not considered. Inferences about the form of the proton spectrum at 10-100 MeV have already been drawn from de-excitation gamma ray lines. Our aim is to constrain the proton spectrum at higher energies. Thus, the injected proton spectrum is assumed to have the form of a Bessel Function, characteristics of stochastic energy at higher energies. The detailed shape of the gamma ray spectra around 100 MeV is found to have a strong dependence on the spectral index of the power law and on the turnover energy (from Bessel function to power law). As would be expected, the harder the photon spectrum the wider the 100 MeV feature. The photon spectra are to be compared with observations and used to place limits upon the number of particles accelerated and to constrain acceleration models.

Alexander, David↗
DOE OSTI2020

Beyond Limber: efficient computation of angular power spectra for galaxy clustering and weak lensing

Angular two-point statistics of large-scale structure observables are important cosmological probes. To reach the high accuracy required by the statistical precision of future surveys, some of these statistics may need to be computed without the commonly employed Limber approximation; the exact computation however requires integration over Bessel functions, and a brute-force evaluation is slow to converge. Here, we present a new method based on our generalized FFTLog algorithm for the efficient computation of angular power spectra beyond the Limber approximation. The new method significantly simplifies the calculation and improves the numerical speed and stability. It is easily extended to handle integrals involving derivatives of Bessel functions, making it equally applicable to numerically more challenging cases such as contributions from redshift-space distortions and Doppler effects. We implement our method for galaxy clustering and galaxy-galaxy lensing power spectra. We find that using the Limber approximation for galaxy clustering in future analyses like LSST Year 1 and DES Year 6 may cause significant biases in cosmological parameters, indicating that going beyond the Limber approximation is necessary for these analyses.

weak gravitational lensing↗
NASA NTRS1998

A High Frequency Model of Cascade Noise

Closed form asymptotic expressions for computing high frequency noise generated by an annular cascade in an infinite duct containing a uniform flow are presented. There are two new elements in this work. First, the annular duct mode representation does not rely on the often-used Bessel function expansion resulting in simpler expressions for both the radial eigenvalues and eigenfunctions of the duct. In particular, the new representation provides an explicit approximate formula for the radial eigenvalues obviating the need for solutions of the transcendental annular duct eigenvalue equation. Also, the radial eigenfunctions are represented in terms of exponentials eliminating the numerical problems associated with generating the Bessel functions on a computer. The second new element is the construction of an unsteady response model for an annular cascade. The new construction satisfies the boundary conditions on both the cascade and duct walls simultaneously adding a new level of realism to the noise calculations. Preliminary results which demonstrate the effectiveness of the new elements are presented. A discussion of the utility of the asymptotic formulas for calculating cascade discrete tone as well as broadband noise is also included.

Envia, Edmane↗
NASA NTRS1984

The energy spectra of solar energetic particles

A survey of recent results on the shapes and relative slopes of the spectra of various solar energetic particle populations is presented, with emphasis on the more extensive results currently available for protons, alphas and electrons. From previous work, it is found that proton spectra 0.8 to more than 400 MeV and alpha spectra 1.4 to 80 MeV/nucleon are best characterized, on average, by a functional form involving a Bessel function in momentum/nucleon. However, proton and alpha spectral slopes using this form are not equal, and there is significant variation from event to event. From other studies, electrons 0.02 to 20 MeV are also found to have curved spectra, but seem to be better fit with a double power law in energy. The spectral properties in both cases correlate with other measures of solar particle acceleration; e.g. gamma-ray line production, hard X-ray burst spectra and microwave fluxes.

Mcguire, R. E.↗
DOE OSTI2018

Systems of Functions Orthogonal Over the Domain and Their Application in Boundary-Value Problems of Mathematical Physics

We formulate a boundary-value problem for eigenvalues and eigenfunctions of the Helmholtz equation in a complex domain with the use of mutually conjugated complex variables. The obtained systems of functions are orthogonal in this domain and constructed by using the Bessel functions and powers of the conformal mappings of the analyzed domains onto a circle. The solutions of boundary-value problems for the principal equations of mathematical physics (hyperbolic, parabolic and elliptic types) are obtained in the form of the sums of series in the systems of functions orthogonal over the domain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fitting Matérn smoothness parameters using automatic differentiation

The Mat$\acute{e}$rn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $\nu$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Mat$\acute{e}$rn covariance function with respect to $\nu$ require derivatives of the modified second-kind Bessel function $K$ $\nu$ with respect to $\nu$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $\nu$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $\nu$ use finite difference estimates for $\partial$ $\nu$ $K$ $\nu$ . In this work, we introduce a new implementation of $K$$\nu$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. Here, we provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

A new linearized theory of laminar film condensation of two phase annular flow in a capillary pumped loop

Future large space based facilities, such as Space Station, will require energy management systems capable of transporting tens of kilowatts of heat over a hundred meters or more. This represents better than an order of magnitude improvement over current technology. Two-phase thermal systems are currently being developed to meet this challenge. Condensation heat transfer plays a very important role in this system. The present study attempts an analytic solution to the set of linearized partial differential equations. The axial velocity and temperature functions were found to be Bessel functions which have oscillatory behavior. This result agrees qualitatively with the experimental evidence from tests at both NASA Goddard Space Flight Center and elsewhere.

Hsu, Y. K.↗
NASA NTRS1991

Computer program for Bessel and Hankel functions

A set of FORTRAN subroutines for calculating Bessel and Hankel functions is presented. The routines calculate Bessel and Hankel functions of the first and second kinds, as well as their derivatives, for wide ranges of integer order and real or complex argument in single or double precision. Depending on the order and argument, one of three evaluation methods is used: the power series definition, an Airy function expansion, or an asymptotic expansion. Routines to calculate Airy functions and their derivatives are also included.

Kreider, Kevin L.↗
DOE OSTI2018

Sommerfeld effect in a constrained electromechanical system

In this paper, the dynamics of an electromechanical system is investigated. In this model, the phase of the system is a dynamical variable. The existence of four periodic orbits is proved. Two of them are asymptotically stable and the others are unstable. In the absence of viscous damping, when the parameters associated with the dimensionless voltages are adequately changed, each orbit is subject to a change of stability of the kind stable → unstable → stable with repetition of this pattern. This phenomenon is known as Sommerfeld effect. Moreover, the stability of these orbits depends on the zeros of the Bessel function J{sub 1}.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Properties of a potential energy matrix in oscillator basis

Matrix elements of potential energy are examined in detail. We consider a model problem — a particle in a central potential. The most popular forms of central potential are taken up, namely, square-well potential, Gaussian, Yukawa and exponential potentials. We study eigenvalues and eigenfunctions of the potential energy matrix constructed with oscillator functions. It is demonstrated that eigenvalues coincide with the potential energy in coordinate space at some specific discrete points. We establish approximate values for these points. It is also shown that the eigenfunctions of the potential energy matrix are the expansion coefficients of the spherical Bessel functions in a harmonic oscillator basis. We also demonstrate a close relation between the separable approximation and L{sup 2} basis (J-matrix) method for the quantum theory of scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Fokker-Planck equation in the presence of a uniform magnetic field

The Fokker-Planck equation in the presence of a uniform magnetic field is derived which has the same form as the case of no magnetic field but with different Fokker-Planck coefficients. The coefficients are calculated explicitly within the binary collision model, which are free from infinite sums of Bessel functions. They can be used to investigate relaxation and transport phenomena conveniently. The kinetic equation is also manipulated into the Landau form from which it is straightforward to compare with previous results and prove the conservation laws.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Dynamics of atoms in strong laser fields I: A quasi analytical model in momentum space based on a Sturmian expansion of the interacting nonlocal Coulomb potential

In this study we present a model that we have formulated in the momentum space to describe atoms interacting with intense laser fields. As a further step, it follows our recent theoretical approach in which the kernel of the reciprocal-space time-dependent Schrödinger equation (TDSE) is replaced by a finite sum of separable potentials, each of them supporting one bound state of atomic hydrogen (Tetchou Nganso et al. 2013). The key point of the model is that the nonlocal interacting Coulomb potential is expanded in a Coulomb Sturmian basis set derived itself from a Sturmian representation of Bessel functions of the first kind in the position space. As a result, this decomposition allows a simple spectral treatment of the TDSE in the momentum space. In order to illustrate the credibility of the model, we have considered the test case of atomic hydrogen driven by a linearly polarized laser pulse, and have evaluated analytically matrix elements of the atomic Hamiltonian and dipole coupling interaction. For various regimes of the laser parameters used in computations our results are in very good agreement with data obtained from other time-dependent calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Bifurcation of space-charge wave in a plasma waveguide including the wake potential effect

The wake potential effects on the propagation of the space-charge dust ion-acoustic wave are investigated in a cylindrically bounded dusty plasma with the ion flow. The results show that the wake potential would generate the double frequency modes in a cylindrically bounded dusty plasma. It is found that the upper mode of the wave frequency with the root of higher-order is smaller than that with the root of lower-order in intermediate wave number domains. However, the lower mode of the scaled wave frequency with the root of higher-order is found to be greater than that with the root of lower-order. It is found that the influence in the order of the root of the Bessel function on the wave frequency of the space-charge dust-ion-acoustic wave in a cylindrically confined dusty plasma decreases with an increase in the propagation wave number. It is also found that the double frequency modes increase with increasing Mach number due to the ion flow in a cylindrical dusty plasma. In addition, it is found that the upper mode of the group velocity decreases with an increase in the scaled radius of the plasma cylinder. However, it is shown that the lower mode of the scaled group velocity of the space-charge dust ion acoustic wave increases with an increase in the radius of the plasma cylinder. The variation of the space-charge dust-ion-acoustic wave due to the wake potential and geometric effects is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

On propagation of long waves in curved ducts.

A two-dimensional detailed study of the behavior of long waves in curved ducts and in junctions between straight and curved ducts will be given. The mathematical treatment of the problem utilizes the method of separation variables. Solutions and expressions for principal mode of the wave are obtained by using the linearized equation of motion solved for its characteristic values. The unavoidable approximations in the numerical solutions of the cylindrical functions are due to use of series expansion of Bessel functions and from restrictions necessary to solve infinite matrices.

Rostafinski, W.↗
NASA NTRS2022

Further Development of a Function Generator Technique

Further features and properties of functions generated by generalization of the process originally uncovered to describe the evolution of pressure in a rigid volume due to outgassing or desorption are explored. Properties presented include development of a general Maclaurin series associated with these functions, a limitation in using integration by parts to produce asymptotic expansions, and a general description of an implicit Adams-Moulton method for numerical quadrature of transformed functions. These developments are then applied to develop and explore functions that solve the Sievert integral and the modified Bessel function of the first kind, order zero. Features of the error function and the incomplete lower gamma function are also considered.

Michael S. Woronowicz↗