The radial heat polynomials and related functions.
Generalized heat equation solution expansion in terms of radial heat polynomials, Appel transforms and Laguerre polynomials
SEARCH · Engineering Papers
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.
Generalized heat equation solution expansion in terms of radial heat polynomials, Appel transforms and Laguerre polynomials
Radial-inflow turbine performance over range of specific speeds
Effect of Reynolds number on performance of 4.59-inch tip diameter radial inflow turbine
Radial heat equation and Laplace transforms, discussing representations alternative to two solutions
WKB approximation of quantum radial distribution function for neon intermolecular pair potential taking into account particle exchange effect
Approximate electron radial functions for uniformly charged nucleus /UCN/ used in analyzing boron 12 and nitrogen 12 beta decay spectra
Electromagnetically driven confined viscous vortex flow with high positive radial Reynolds numbers, solving energy equation for temperature distribution
NASA radial flow turbine research related to Brayton cycle power generator, investigating size, Reynolds number and speed effects on efficiency and flow
The use of radial basis function networks for fine pointing NASA's 70-meter deep space network antennas is described and evaluated.
The NEID spectrometer is an optical (380-930 nm), fiber-fed, precision Doppler spectrometer currently in development for the WIYN 3.5 m telescope at Kitt Peak National Observatory as part of the NN-EXPLORE partnership. Designed to achieve a radial velocity precision of < 30 cm/s, NEID will be sensitive enough to detect terrestrial-mass exoplanets around low-mass stars. Light from the target stars is focused by the telescope to a bent Cassegrain port at the edge of the primary mirror mechanical support. The specialized NEID "Port Adapter" system is mounted at this bent Cassegrain port and is responsible for delivering the incident light from the telescope to the NEID fibers. In order to provide stable, high-quality images to the science instrument, the Port Adapter houses several sub-components designed to acquire the target stars, correct for atmospheric dispersion, stabilize the light onto the science fibers, and calibrate the spectrometer by injecting known wavelength sources such as a laser frequency comb. Here we provide an overview of the overall opto-mechanical design and system requirements of the Port Adapter. We also describe the development of system error budgets and testplans to meet those requirements.
Two invariants associated with the parker model of the solar wind involved the radial, BR, and azimuthal, BT, components of the heliospheric magnetic field. These invariants have been investigated using Ulysses data obtained at high latitudes in both the north and south solar hemispheres.
Dipolar interactions are ever-present in supramolecular architectures, though their impact is typically revealed by making dipoles stronger. While it is also possible to assess the role of dipoles by altering their orientations by using synthetic design, doing so without altering the molecular shape is not straightforward. We have now done this by flipping one triazole unit in a rigid macrocycle, tricarb. The macrocycle is composed of three carbazoles (2 Debye) and three triazoles (5 Debye) defining an array of dipoles aligned radially but organized alternately in and out. These dipoles are believed to dictate edge-to-edge tiling and face-to-face stacking. We modified our synthesis to prepare isosteric macrocycles with the orientation of one triazole dipole rotated 40°. The new dipole orientation guides edge-to-edge contacts to reorder the stability of two surface-bound 2D polymorphs. The impact on dipole-enhanced π stacking, however, was unexpected. Our stacking model identified an unchanged set of short-range (3.4 Å) anti-parallel dipole contacts. Despite this situation, the reduction in self-association was attributed to long-range (~6.4 Å) dipolar repulsions between π-stacked macrocycles. This work highlights our ability to control the build-up and symmetry of macrocyclic skeletons by synthetic design, and the work needed to further our understanding of how dipoles control self-assembly.
We extract the relativistic classical radial action from scattering amplitudes, to all orders in perturbation theory, in the probe limit. Our sources include point charges and monopoles, as well as the Schwarzschild and pure-NUT gravitational backgrounds. A characteristic relativistic effect, that scattering trajectories may wind around these sources any number of times, can be recovered when all-order amplitudes are available. We show that the amplitude for scattering a probe off a pure NUT is given by the solution of a transcendental equation involving continued fractions, and explain how to solve this equation to any desired loop order.
We resolve subtleties in calculating the post-Minksowskian dynamics of binary systems, as a spin expansion, from massive scattering amplitudes of fixed finite spin. In particular, the apparently ambiguous spin Casimir terms can be fully determined from the gradient of the spin-diagonal part of the amplitudes with respect to S 2 = −s(s+1)ħ 2 , using an interpolation between massive amplitudes with different spin representations. From two-loop amplitudes of spin-0 and spin-1 particles minimally coupled to gravity, we extract the spin Casimir terms in the conservative scattering angle between a spinless and a spinning black hole at $ \mathcal{O} $(G 3 S 2 ), finding agreement with known results in the literature. This completes an earlier study [Phys. Rev. Lett. 130 (2023), 021601] that calculated the non-Casimir terms from amplitudes. We also illustrate our methods using a model of spinning bodies in electrodynamics, finding agreement between scattering amplitude predictions and classical predictions in a root-Kerr electromagnetic background up to $ \mathcal{O} $(α 3 S 2 ). For both gravity and electrodynamics, the finite part of the amplitude coincides with the two-body radial action in the aligned spin limit, generalizing the amplitude-action relation beyond the spinless case. Surprisingly, the two-loop amplitude displays a hidden spin-shift symmetry in the probe limit, which was previously observed at one loop. We conjecture that the symmetry holds to all orders in the coupling constant and is a consequence of integrability of Kerr orbits in the probe limit at the first few orders in spin.
Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. Further, the second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.
Functionally graded materials (FGMs) are of interest in multiple fields, yet many materials combinations are limited by coefficient of thermal expansion (CTE) mismatch. Here, a radially graded alumina/yttria-doped zirconia (Al 2 O 3 /8YZ) FGM is used to demonstrate processing strategies to mitigate CTE and sintering behavior differences between these oxides. FGM materials are especially sensitive to ink stability during printing, as all components (in this case, Al 2 O 3 and 8YZ) must be stabilized in the same dispersant or additive solution. Thus, this system is also ideal to demonstrate ink optimization best practices. Materials were characterized throughout processing to correlate the effects of common additives on both the ceramic particle suspensions and final sintered components. Aggregation observed in the initial additive-containing suspensions were present in the sintered component. The differences in sintering onset temperature and shrinkage rate resulted in internal stresses within the sintered component, which ultimately caused mechanical failure of the component under low stress. In conclusion, a processing strategy was recommended to mitigate the sintering behavior mismatch of alumina and 8 wt% yttria-stabilized zirconia.
A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.
Here, in this work, a multi-scale model of a radial flow fixed bed contactor packed with a carbon sorbent is developed and validated with laboratory-scale and pilot plant scale dynamic data. For the lab scale system, the model results were compared with low, medium and high gas and sweep flowrates, yielding root mean square error (RMSE) of 0.80, 0.63, 0.96 CO 2 mol%, respectively, for the outlet CO 2 concentration profile considering the entire A-D cycle. For the bed outer temperature profile, maximum RMSE was found to be 5.5 °C considering all flowrates and entire A-D cycles. An experimental campaign was developed and applied to a pilot plant at Technology Center Mongstad (TCM), Norway. Approaches were developed for pre-processing of data including consideration of the effect of gas mixing, measurement delay, and determination of cyclic steady-state conditions. Considering profiles during A-D cycles for all test runs, it was found that the maximum RMSE for pressure drop, temperature for the outer section of the bed, temperature for the middle section of the bed, and outlet CO 2 concentration profile remained less than 1.5 mbar, 3.5 °C, 2.8 °C, and 1.3 CO 2 mol%, respectively. The validated model was used to perform sensitivity studies on several key design operating variables for the adsorption-desorption cycle. It was found that the flow rate and concentration of flue gas have dominant nonlinear effects on the breakthrough time while the desorption time was strongly affected by the sweep gas flowrate for the specific sorbent being evaluated in this study.