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

Viscous-fluid reaction to a torsionally oscillating spheroid Linearized steady-state weak-separation solution

For both the linearized prolate and oblate problems, perturbation series for all orders of eccentricity are obtained. The 'weak separation' method of solving partial-differential-equation boundary-value problems assumes that each linearly independent solution is a sum of products of single-variable functions. The method is applied to a hierarchy of partial differential equations. The weak-separation solutions are constructed by generalizing only the Bessel functions of the homogeneous solution.

Tompkins, D. R., Jr.↗

A note of expansions of functions of velocity in the two-body problem

Fourier expansions based on both the true anomaly and the mean anomaly are obtained for the functions of velocity in the two-body problem; the series of coefficients is written from classical formulae involving associated Legendre polynomials, Gegenbauer polynomials, or Bessel functions. The Fourier expansions are compared with the expansions in powers of eccentricity developed by Broucke (1974) through use of computerized Poisson series manipulation.

Kinoshita, H.↗

Radiation Losses Due to Tapering of a Double-Core Optical Waveguide

The theoretical model we designed parameterizes the power losses as a function of .the profile shape for a tapered, single mode, optical dielectric coupler. The focus of this project is to produce a working model that determines the power losses experienced by the fibers when light crosses a taper region. This phenomenon can be examined using coupled mode theory. The optical directional coupler consists of a parallel, dual-channel, waveguide with minimal spacing between the channels to permit energy exchange. Thus, power transfer is essentially a function of the taper profile. To find the fields in the fibers, the approach used was that of solving the Helmholtz equation in cylindrical coordinates involving Bessel and modified Bessel functions depending on the location.

Lyons, Donald R.↗

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Generalized Airy functions for use in one-dimensional quantum mechanical problems

The solution of the one dimensional, time independent, Schroedinger equation in which the energy minus the potential varies as the nth power of the distance is obtained from proper linear combinations of Bessel functions. The linear combinations called generalized Airy functions, reduce to the usual Airy functions Ai(x) and Bi(x) when n equals 1 and have the same type of simple asymptotic behavior. Expressions for the generalized Airy functions which can be evaluated by the method of generalized Gaussian quadrature are obtained.

Eaves, J. O.↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electrostatic approximation

The plasma weak turbulence theory is a perturbative nonlinear theory, which has been proven to be quite valid in a number of applications. However, the standard weak turbulence theory found in the literature is fully developed for highly idealized unmagnetized plasmas. As many plasmas found in nature and laboratory are immersed in a background static magnetic field, it is necessary to extend the existing discussions to include the effects of ambient magnetic field. Such a task is quite formidable, however, which has prevented fundamental and significant progresses in the subject matter. The central difficulty lies in the formulation of the complete nonlinear response functions for magnetized plasmas. The present paper derives the nonlinear susceptibilities for weakly turbulent magnetized plasmas up to the third order nonlinearity, but in doing so, a substantial reduction in mathematical complexity is achieved by the use of Bessel function addition theorem (or sum rule). The present paper also constructs the weak turbulence wave kinetic equation in a formal sense. For the sake of simplicity, however, the present paper assumes the electrostatic interaction among plasma particles. Fully electromagnetic generalization is a subject of a subsequent paper.

Physics↗

The classical solution of a definite integral occurring in the satellite theory in the extended phase space

A method of integrating the functions defined by Scheifele and Graf (1974), arising in the elimination of time from the arguments of satellite theory using Delaunay elements in the extended phase space, is proposed. By repeated applications of an identity of Stiefel and Scheifele (1971), neglecting certain terms, and assuming a solution in the form of a product of Bessel functions, a solution expanded in the eccentricity is obtained.

Bond, V. R.↗

Some MACSYMA program for solving difference equations

A set of MACSYMA programs are described for finding closed form solutions to linear recurrence relations in equations having either constant or variable coefficients. In the homogenous case, a polymonial equation is obtained and the solution to the recurrence relation can be written as a linear combination of the roots of the polynomial. Exponential generating functions are used to solve variable coefficient relations. Taking successive derivatives and using the recurrence relation, an ordinary differential equation is obtained. Expanding the solution to the differential equation in a Taylor series, shows that the nth term of the series is the solution to the recurrence relation. For second order recurrences, a check is made for those that can be solved in terms of Bessel functions.

Ivie, J.↗

Solution to the backward-Kolmogorov equation for a nonstationary oscillation problem

The transition probability density function of a Markovian approximation of the response amplitude of an oscillator under nonstationary excitation is determined in an analytical form. A solution is presented for the associated backward-Kolmogorov equation by transforming the equation into a form amenable to solution by the method of the separation of variables. This procedure results in a boundary value problem which is then solved by using an infinite series of Laguerre polynomials. It is found that the infinite series solution is equivalent to a closed-form solution involving a Bessel function.

Solomos, G. P.↗

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗

Non-symmetric two-stream instability

A theoretical investigation is performed concerning the instability spectrum of quasi-electrostatic waves at shifted half-odd-integer values of the cyclotron frequency due to nonsymmetric counterstreaming electron beams. The beam velocities parallel and perpendicular to the static magnetic field are represented by double Dirac delta functions with no imposition of parameter limitations. Systematic consideration is given to the coupling between plasma modes and cyclotron modes as well as the coupling between cyclotron modes of the two beams that result in shifted half-odd-integer multiples of the cyclotron frequency. The general dispersion equation for quasi-electrostatic waves is analyzed, and plasma-cyclotron coupling in the nonsymmetric case is treated by deriving approximate analytical expressions for maximum growth rates and marginal stability. Exact numerical solutions in both frequency and wavenumber space are obtained and compared with the analytical expressions. Cyclotron-cyclotron coupling modes are treated in the same way, and the results for both types of coupling are compared. It is found that certain modes may be weakened or completely suppressed when the Bessel functions for specific plasma parameters vanish.

Cuperman, S.↗

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determining all thermodynamic transport coefficients for an interacting large N quantum field theory

Thermodynamic transport coefficients can be calculated directly from quantum field theory without requiring analytic continuation to real time. We determine all second-order thermodynamic transport coefficients for the uncharged N-component massless (critical) scalar field theory with quartic interaction in the large N limit, for any value of the coupling. We find that in the large N limit, all thermodynamic transport coefficients for the interacting theory can be expressed analytically in terms of the in-medium mass and sums over modified Bessel functions. We expect our technique to allow a similar determination of all thermodynamic transport coefficients for all theories that are solvable in the large N limit, including certain gauge theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

Microscopic-macroscopic level densities for low excitation energies

Level density ρ(E,Q) is derived within the micro-macroscopic approximation (MMA) for a system of strongly interacting Fermi particles with the energy E and additional integrals of motion Q , in line with several topics of the universal and fruitful activity of A. S. Davydov. Within the extended Thomas Fermi and semiclassical periodic orbit theory beyond the Fermi-gas saddle-point method, we obtain ρ ∝ I ν (S)/S ν , where I ν (S) is the modified Bessel function of the entropy S . For small shell-structure contribution, one finds ν = κ/2 + 1, where κ is the number of additional integrals of motion. This integer number is a dimension of Q, Q = { N, Z , …} for the case of two-component atomic nuclei, where N and Z are the numbers of neutrons and protons, respectively. For much larger shell structure contributions, one obtains ν = κ /2 + 2. The MMA level density ρ reaches the well-known Fermi gas asymptote for large excitation energies and the finite micro-canonical combinatoric limit for low excitation energies. Further, the additional integrals of motion can also be the projection of the angular momentum of a nuclear system for nuclear rotations of deformed nuclei, number of excitons for collective dynamics, and so on. Fitting the MMA total level density ρ( E , Q) for a set of the integrals of motion Q = { N, Z }, to experimental data on a long nuclear isotope chain for low excitation energies, one obtains the results for the inverse level-density parameter K , which differs significantly from those of neutron resonances due to shell, isotopic asymmetry, and pairing effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electromagnetic formalism

This is a companion paper to the previous work [P. H. Yoon, Phys. Plasmas 31, 032309 (2024)] in which the nonlinear susceptibilities of weakly turbulent magnetized plasma are derived under a simplifying assumption of electrostatic interaction. The present paper extends the analysis to a general situation of electromagnetic interaction. The main novelty of the previous and present papers is that by employing the Bessel function addition theorem, the mathematical definitions for the susceptibilities are substantially simplified, a procedure that has not been discussed in the existing literature. In the present paper, a full set of Maxwell’s equations are considered in conjunction with the nonlinear Vlasov equation, which is solved by a perturbative method. The result is a fully general nonlinear susceptibility, given in tensorial form, which is applicable for weakly turbulent magnetized plasmas.

Materials Science↗

Semiclassical shell-structure micro-macroscopic approach for the level density

Level density ρ(E,A) is derived for a one-component nucleon system with a given energy E and particle number A within the mean-field semiclassical periodic-orbit theory beyond the saddle-point method of the Fermi gas model. We obtain ρ∝I ν (S)/S ν , with I ν (S) being the modified Bessel function of the entropy S. Within the micro-macro-canonical approximation (MMA), for a small thermal excitation energy U, with respect to rotational excitations E rot , one obtains ν = 3/2 for ρ(E,A). In the case of excitation energy U larger than E rot but smaller than the neutron separation energy, one finds a larger value of ν = 5/2. A role of the fixed spin variables for rotating nuclei is discussed. The MMA level density ρ reaches the well-known grand-canonical ensemble limit (Fermi gas asymptote) for large S related to large excitation energies, and also reaches the finite micro-canonical limit for small combinatorial entropy S at low excitation energies (the constant “temperature” model). Fitting the ρ(E,A) of the MMA to the experimental data for low excitation energies, taking into account shell and, qualitatively, pairing effects, one obtains for the inverse level density parameter K a value which differs essentially from that parameter derived from data on neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shell-structure and asymmetry effects in level densities

Level density [Formula: see text] is derived for a nuclear system with a given energy [Formula: see text], neutron [Formula: see text], and proton [Formula: see text] particle numbers, within the semiclassical extended Thomas–Fermi and periodic-orbit theory beyond the Fermi-gas saddle-point method. We obtain [Formula: see text], where [Formula: see text] is the modified Bessel function of the entropy [Formula: see text], and [Formula: see text] is related to the number of integrals of motion, except for the energy [Formula: see text]. For small shell structure contribution one obtains within the micro–macroscopic approximation (MMA) the value of [Formula: see text] for [Formula: see text]. In the opposite case of much larger shell structure contributions one finds a larger value of [Formula: see text]. The MMA level density [Formula: see text] reaches the well-known Fermi gas asymptote for large excitation energies, and the finite micro-canonical limit for low excitation energies. Fitting the MMA [Formula: see text] to experimental data on a long isotope chain for low excitation energies, due mainly to the shell effects, one obtains results for the inverse level density parameter [Formula: see text], which differs significantly from that of neutron resonances.

Physics↗