Search NASASearch

SEARCH · Search NASA

Results for “Multipole expansion”

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

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.

Delta excitations and shell-model information in heavy-ion, charge-exchange reactions

We calculate total cross sections for coherent pion production using localized plane-wave approximations for the shell-structure of valence nucleons that are excited to delta particles in the intermediate state in the (12C, 12B) and (12C, 12N) charge-exchange, heavy-ion reactions. We find comparable agreement to projectile downshift data for 12C(12C, 12B)12N. Then we improve the formalism by replacing the localized plane wave bound states with harmonic oscillator states which are imbedded in a multipole expansion approach and calculate pion differential cross sections to test for the sensitivity of the spectra to the single-particle mass parameter.

NASA Program Space Medicine

Dusty Plasmas in Planetary Magnetospheres Award

This is my final report for the grant Dusty Plasmas in Planetary Magnetospheres. The funding from this grant supported our research on dusty plasmas to study: a) dust plasma interactions in general plasma environments, and b) dusty plasma processes in planetary magnetospheres (Earth, Jupiter and Saturn). We have developed a general purpose transport code in order to follow the spatial and temporal evolution of dust density distributions in magnetized plasma environments. The code allows the central body to be represented by a multipole expansion of its gravitational and magnetic fields. The density and the temperature of the possibly many-component plasma environment can be pre-defined as a function of coordinates and, if necessary, the time as well. The code simultaneously integrates the equations of motion with the equations describing the charging processes. The charging currents are dependent not only on the instantaneous plasma parameters but on the velocity, as well as on the previous charging history of the dust grains.

Horanyi, Mihaly

An Analysis of the Rapidly Rotating Bp star HD 133880

HD 133880 is a rapidly rotating chemically peculiar B-type (Bp) star (nu sin i approx = 103km/s) and is host to one of the strongest magnetic fields of any Ap/Bp star. A member of the Upper Centaurus Lupus association, it is a star with a well-determined age of 16 Myr. 12 new spectra, four of which are polarimetric, obtained from the FEROS, ESPaDOnS and HARPS instruments, provide sufficient material from which to re-evaluate the magnetic field and obtain a first approximation to the atmospheric abundance distributions of He, O, Mg, Si, Ti. Cr, Fe, Ni, Pr and Nd. An abundance analysis was carried out using ZEEMAN, a program which synthesizes spectral line profiles for stars with permeating magnetic fields. The magnetic field structure was characterized by a colinear multipole expansion from the observed variations of the longitudinal and surface fields with rotational phase. Both magnetic hemispheres are clearly visible during the stellar rotation, and thus a three-ring abundance distribution model encompassing both magnetic poles and magnetic equator with equal spans in colatitude was adopted. Using the new magnetic field measurements and optical photometry together with previously published data, we refine the period of HD 133880 to P = 0.877 476 +/- 0.000009 d. Our simple axisymmetric magnetic field model is based on a predominantly quadrupolar component that roughly describes the field variations. Using spectrum synthesis, we derived mean abundances for O, Mg, Si, Ti, Cr, Fe and Pr. All elements; except Mg, are overabundant compared to the Son. Mg appears to be approximately uniform over the stellar surface, while all other elements are more abundant in the negative magnetic hemisphere than in the positive magnetic hemisphere. In contrast to most Ap/Bp stars which show an underabundance in 0, in HD 133880 this element is clearly overabundant compared to the solar abundance ratio. In studying the Ha and Paschen lines in the optical spectra, we could not unambiguously detect information about the magnetosphere of HD 133880. However, radio emission data at both 3 and 6 cm suggest that the magnetospheric plasma is held in rigid rotation with the star by the magnetic field and further supported against collapse by the rapid rotation. Subtle differences in the shapes of the optically thick radio light curves at 3 and 6 cm suggest that the large-scale magnetic field is not fully axisymmetric at large distances from the star.

Bailey, J. D.

An Extension of the Truncated Ingard-Myers Impedance Boundary Condition for High Mach Number Grazing Flows

The Ingard-Myers impedance boundary condition is widely recognized to be prone to hydrodynamic Kelvin-Helmholtz-type instability, primarily due to its use of a vortex sheet to model the flow at the boundary. Recently, a stabilized approximation of the Ingard-Myers condition, called the Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC), has been introduced. This approximation offers a good representation of the Ingard-Myers condition for grazing flows characterized by low to mid subsonic Mach numbers. This paper explores an extension of the TIMIBC, referred to as TIMIBC-ext, designed for liners in grazing flows with high subsonic Mach numbers. The TIMIBC-ext introduces a tunable parameter and demonstrates that by judiciously selecting the parameter value, the accuracy of the TIMIBC can be improved for either upstream or downstream propagating waves in grazing flows with high subsonic Mach numbers. Consequently, the TIMIBC-ext represents an enhancement to the TIMIBC when the location of the source relative to the liner (upstream or downstream) is known. Time domain implementations of the TIMIBC-ext are also provided, employing a multipole expansion model for the impedance function. A numerical example is provided to assess the performance of the TIMIBC-ext, comparing computational results with measurements from a recent set of measurements in the NASA Langley Grazing Flow Impedance Tube.

time domain acoustic scattering

Elliptic multipoles and the modeling of narrow-gap bend magnets in accelerators

We highlight the virtues of 2D elliptic-multipole field expansions in modeling the magnetic fields of narrow-aperture, straight-axis bending magnets with parallel faces, addressing the limitations of the conventional circular multipole series when the beam-orbit sagitta exceeds the magnet's vertical half-gap. The elliptic multipoles provide a convenient way to represent the field in all aspects of the magnet development (design, particle-tracking simulations, measurements). We propose a numerically robust method of data analysis to determine the elliptic (or circular) multipoles from stretched-wire measurements with the wire moving on an arbitrary path.

Venturini, Marco

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Some improvements in DNA interaction calculations

Calculations are made on specific DNA-type complexes using refined expressions for electrostatic and polarization energies. Dispersion and repulsive terms are included in the evaluation of the total interaction energy. It is shown that the expansion of the electrostatic potential to include multipole moments up to octopole is necessary to achieve convergence of first-order energies. Polarization energies are not as sensitive to this expansion. The calculations also support the usefulness of the hard sphere model for DNA hydrogen bonds and indicate how stacking interactions are influenced by second-order energies.

Egan, J. T.

A Green’s function fast multipole method for computation of micromechanical fields in heterogeneous materials

Computation of micromechanical fields in heterogeneous materials is usually performed using either the finite element method or the Green’s function method based on FFTs. The finite element method allows for accurate discretization and for non-periodic boundary conditions but is computationally expensive. On the other hand, the FFT-based method is computationally efficient but requires discretization on a regular grid of hexahedral voxels. In this paper, a Green’s function method allowing for accurate discretization using tetrahedral elements and for non-periodic boundary conditions is proposed. The convolution is computed using the fast multipole method, which provides good accuracy even for low-order expansion due to the fast decay of interactions between elements. The proposed Green’s function fast multipole method is verified by comparison with analytical and FFT-based solutions. Furthermore, the computational time is analyzed and compared to the FFT-based method for non-periodic convolution. Finally, effective properties of an elastic polycrystalline microstructure containing thin intergranular cracks are computed and analyzed.

36 MATERIALS SCIENCE

Gravitational scattering of electromagnetic radiation

The scattering of electromagnetic radiation by linearized gravitational fields is studied to second order in a perturbation expansion. The incoming electromagnetic radiation can be of arbitrary multipole structure, and the gravitational fields are also taken to be advanced fields of arbitrary multipole structure. All electromagnetic multipole radiation is found to be scattered by gravitational monopole and time-varying dipole fields. No case has been found, however, in which any electromagnetic multipole radiation is scattered by gravitational fields of quadrupole or higher-order multipole structure. This lack of scattering is established for infinite classes of special cases, and is conjectured to hold in general. The results of the scattering analysis are applied to the case of electromagnetic radiation scattered by a moving mass. It is shown how the mass and velocity may be determined by a knowledge of the incident and scattered radiation.

Brooker, J. T.

On the Inversion for Mass (Re)Distribution from Global (Time-Variable) Gravity Field

The well-known non-uniqueness of the gravitational inverse problem states the following: The external gravity field, even if completely and exactly known, cannot Uniquely determine the density distribution of the body that produces the gravity field. This is an intrinsic property of a field that obeys the Laplace equation, as already treated in mathematical as well as geophysical literature. In this paper we provide conceptual insight by examining the problem in terms of spherical harmonic expansion of the global gravity field. By comparing the multipoles and the moments of the density function, we show that in 3-S the degree of knowledge deficiency in trying to inversely recover the density distribution from external gravity field is (n+l)(n+2)/2 - (2n+l) = n(n-1)/2 for each harmonic degree n. On the other hand, on a 2-D spherical shell we show via a simple relationship that the inverse solution of the surface density distribution is unique. The latter applies quite readily in the inversion of time-variable gravity signals (such as those observed by the GRACE space mission) where the sources over a wide range of the scales largely come from the Earth's Surface.

Chao, Benjamin F.

Exact formulas for multipole moments using Slater-type molecular orbitals

A triple infinite sum of formulas expressed as an expansion in Legendre polynomials is generated by use of computer algebra to represent the potential from the midpoint of two Slater-type orbitals; the charge density that determines the potential is given as the product of the two orbitals. An example using 1s orbitals shows that only a few terms are needed to obtain four-figure accuracy. Exact formulas are obtained for multipole moments by means of a careful study of expanded formulas, allowing an 'extrapolation to infinity'. This Loewdin alpha-function approach augmented by using a C matrix to characterize Slater-type orbitals can be readily generalized to all cases.

Jones, H. W.

Multipole moments of stellar oscillation modes

The oscillating mass 2(exp l)-pole moment, M(sub nl), of a star in a given (normalized) oscillation mode determines the energy that can be absorbed by the mode in a tidal interaction and the power radiated by the mode in gravitational waves, both of which are proportional to (absolute value of M(sub nl))(exp 2). The coefficients in the expansion of the vector fields del(r(exp l)Y(sub lm)(theta, phi)) in terms of the displacement fields of modes of given l and m are proportional to M(sub nl). This expansion leads to a sum rule sum over n(absolute value of M(sub nl))(exp 2) = constant. For stars of weak to moderate central condensation (such as neutron stars), the f-mode is well approximated by the vector field being expanded, and therefore it takes the lion's share of the sum. Thus the multipole moments of all other modes must be small. In there numerical evaluation, it is necessary to know the shape of the eigenfunctions quite precisly, since a small f-mode contamination can significantly increase the obtained values. This contamination occurs in some `hybrid' numerical computations of neutron star oscillations with relativistic equilibrium stars and Newtonian dynamics (e.g., McDermott et al. 1988). In this case, it is due to a slight inconsistency in the models and leads to a large overestimate of the power radiated in gravitational waves by modes other than the f-mode, although their oscillation periods are nearly unaffected.

Reisenegger, Andreas

Drift shell splitting by internal geomagnetic multipoles.

Computations on an 80-coefficient model of the earth's field illustrate the 'topography' of the magnetic equatorial surface and the geometry of the drift shells of geomagnetically trapped particles. Individual terms in the spherical harmonic expansion of the geomagnetic scalar potential V(r, theta, phi) are either even or odd in cos theta, where theta = 90 deg denotes the dipole equator. Terms that are even in cos theta tend to 'warp' the equatorial surface, but do not (in first order) distort particle drift shells radially nor split the drift shells of particles having different equatorial pitch angles. Azimuthally asymmetric terms that are odd in cos theta do cause shell splitting in first order. Shell splitting at large L values (neglecting deformation of the earth's field by the solar wind) is found to be dominated by the geomagnetic octupole. At L approximately equal to 1, shell splitting is strongly enhanced by the South American and South African anomalies. When combined with pitch angle diffusion caused by atmospheric scattering, these results may be able to account for anomalous radial diffusion of inner zone electrons.

Roederer, J. G.

Investigation of geomagnetic field forecasting and fluid dynamics of the core

Accomplishments to date were summarized in three parts submitted for publication. Goddard models and MAGSAT data were used heavily in the investigation which address: (1) the sensitivity of selected geomagnetic properties to truncation level of spherical harmonic expansions; (2) the pole strength of the Earth from MAGSAT, and magnetic determination of the core radius; and (3) frozen flux upper limits to the MAGSAT geomagnetic coefficients and relative multipole indices for Earth.

Benton, E. R.

The exterior source surface for force-free fields

Consideration is given to the exterior source surface for force-free fields. The spherical harmonic expansion is presented for boundary values on two concentric spheres. An upper limit on a constant which measures the strength of coronal currents is found to be a function of the lowest multipole moment of the prescribed boundary values. The solar atmosphere is in the class of magnetic fields for which the study is applicable.

Barbosa, D. D.

Equipotential transformation of multipole systems to dipole systems

The vector sum of fields produced by a set of magnetic dipoles with different magnetic moments that are both tilted from the planetary spin axis and offset from the planetary center by different amounts can be used to completely represent a planetary magnetic field in the conventional form of spherical harmonic expansion. The scalar equipotential transformation analytically results in 24 equations that may subsequently be solved for the 24 adjustable parameters in dipole systems with the predetermined main dipole. Attention is given to this method's application to the Jovian magnetic field.

Wang, W. X.

The Simons Observatory: forecasted constraints on primordial gravitational waves with the expanded array of Small Aperture Telescopes

We present updated forecasts for the scientific performance of the degree-scale (0.5 deg FWHM at 93 GHz), deep-field survey to be conducted by the Simons Observatory (SO). By 2027, the SO Small Aperture Telescope (SAT) complement will be doubled from three to six telescopes, including a doubling of the detector count in the 93 GHz and 145 GHz channels to 48,160 detectors. Combined with a planned extension of the survey duration to 2035, this expansion will significantly enhance SO's search for a B-mode signal in the polarisation of the cosmic microwave background, a potential signature of gravitational waves produced in the very early Universe. Assuming a 1/f noise model with knee multipole ℓ knee = 50 and a moderately complex model for Galactic foregrounds, we forecast a 1σ (or 68% confidence level) constraint on the tensor-to-scalar ratio r of σ r = 1.2 × 10 -3 , assuming no primordial B-modes are present. This forecast assumes that 70% of the B-mode lensing signal can ultimately be removed using high resolution observations from the SO Large Aperture Telescope (LAT) and overlapping large-scale structure surveys. For more optimistic assumptions regarding foregrounds and noise, and assuming the same level of delensing, this forecast constraint improves to σ r = 7 × 10 -4 . These forecasts represent a major improvement in SO's constraining power, being a factor of around 2.5 times better than what could be achieved with the originally planned campaign, which assumed the existing three SATs would conduct a five-year survey.

CMBR experiments