Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Exact two component”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

A Universal Tare Load Prediction Algorithm for Strain-Gage Balance Calibration Data Analysis

An algorithm is discussed that may be used to estimate tare loads of wind tunnel strain-gage balance calibration data. The algorithm was originally developed by R. Galway of IAR/NRC Canada and has been described in the literature for the iterative analysis technique. Basic ideas of Galway's algorithm, however, are universally applicable and work for both the iterative and the non-iterative analysis technique. A recent modification of Galway's algorithm is presented that improves the convergence behavior of the tare load prediction process if it is used in combination with the non-iterative analysis technique. The modified algorithm allows an analyst to use an alternate method for the calculation of intermediate non-linear tare load estimates whenever Galway's original approach does not lead to a convergence of the tare load iterations. It is also shown in detail how Galway's algorithm may be applied to the non-iterative analysis technique. Hand load data from the calibration of a six-component force balance is used to illustrate the application of the original and modified tare load prediction method. During the analysis of the data both the iterative and the non-iterative analysis technique were applied. Overall, predicted tare loads for combinations of the two tare load prediction methods and the two balance data analysis techniques showed excellent agreement as long as the tare load iterations converged. The modified algorithm, however, appears to have an advantage over the original algorithm when absolute voltage measurements of gage outputs are processed using the non-iterative analysis technique. In these situations only the modified algorithm converged because it uses an exact solution of the intermediate non-linear tare load estimate for the tare load iteration.

balance calibration analysis↗
NASA NTRS1981

The field and plasma configuration of a filament overlying a solar bipolar magnetic region

This paper presents an analytic model for a finite-size straight filament suspended horizontally in a steady state over a bipolar magnetic region. The equations of magnetostatic equilibrium are integrated exactly. The solution obtained illustrates the roles played by the electric current, magnetic field, pressure, and plasma weight in the balance of force everywhere in space. A specific example of a filament of diameter 50,000 km, with a density two orders of magnitude over the corona and supported by a magnetic field of about 4 gauss is included. The filament temperature can take values ranging from a small fraction to a few times the coronal temperature, depending on the internal electric current of the filament. To produce a cool filament, such as the quiescent prominence, the solution is required to have an internal field with a strong component along the filament, giving rise to helical structures. A hot filament such as the X-ray coronal loop can be produced as a twisted magnetic flux tube embedded in a strong background field aligned parallel to the filament and having lower density and temperature. The basic steps of construction can be used to develop models more realistic than the ones presented for their analytic simplicity.

Low, B. C.↗
NASA NTRS1992

Optical testing

Optical testing is one of the most vital elements in the process of preparing an optical instrument for launch. Without well understood, well controlled, and well documented test procedures, current and future mission goals will be jeopardized. We should keep in mind that the reason we test is to provide an opportunity to catch errors, oversights, and problems on the ground, where solutions are possible and difficulties can be rectified. Consequently, it is necessary to create tractable test procedures that truly provide a measure of the performance of all optical elements and systems under conditions which are close to those expected in space. Where testing is not feasible, accurate experiments are required in order to perfect models that can exactly predict the optical performance. As we stretch the boundaries of technology to perform more complex space and planetary investigations, we must expand the technology required to test the optical components and systems which we send into space. As we expand the observational wavelength ranges, so must we expand our range of optical sources and detectors. As we increase resolution and sensitivity, our understanding of optical surfaces to accommodate more stringent figure and scatter requirements must expand. Only with research and development in these areas can we hope to achieve success in the ever increasing demands made on optical testing by the highly sophisticated missions anticipated over the next two decades. Technology assessment and development plan for surface figure, surface roughness, alignment, image quality, radiometric quantities, and stray light measurement are presented.

Wyant, James↗
NASA NTRS2012

A Comparison of Filter-based Approaches for Model-based Prognostics

Model-based prognostics approaches use domain knowledge about a system and its failure modes through the use of physics-based models. Model-based prognosis is generally divided into two sequential problems: a joint state-parameter estimation problem, in which, using the model, the health of a system or component is determined based on the observations; and a prediction problem, in which, using the model, the stateparameter distribution is simulated forward in time to compute end of life and remaining useful life. The first problem is typically solved through the use of a state observer, or filter. The choice of filter depends on the assumptions that may be made about the system, and on the desired algorithm performance. In this paper, we review three separate filters for the solution to the first problem: the Daum filter, an exact nonlinear filter; the unscented Kalman filter, which approximates nonlinearities through the use of a deterministic sampling method known as the unscented transform; and the particle filter, which approximates the state distribution using a finite set of discrete, weighted samples, called particles. Using a centrifugal pump as a case study, we conduct a number of simulation-based experiments investigating the performance of the different algorithms as applied to prognostics.

prognostics↗
NASA NTRS1979

Fiber optic combiner and duplicator

The investigation of the possible development of two optical devices, one to take two images as inputs and to present their arithmetic sum as a single output, the other to take one image as input and present two identical images as outputs is described. Significant engineering time was invested in establishing precision fiber optics drawing capabilities, real time monitoring of the fiber size and exact measuring of fiber optics ribbons. Various assembly procedures and tooling designs were investigated and prototype models were built and evaluated that established technical assurance that the device was feasible and could be fabricated. Although the interleaver specification in its entirety was not achieved, the techniques developed in the course of the program improved the quality of images transmitted by fiber optic arrays by at least an order of magnitude. These techniques are already being applied to the manufacture of precise fiber optic components.

Source record↗
NASA NTRS1992

Distribution of Reynolds stress carried by mesoscale variability in the Antarctic Circumpolar Current

Satellite altimeter data collected by the Geosat Exact Repeat Mission were used to investigate turbulent stress resulting from the variability of surface geostrophic currents in the Antarctic Circumpolar Current. The altimeter measured sea level along the subsatellite track. The variability of the along-track slope of sea level is directly proportional to the variability of surface geostrophic currents in the cross-track direction. Because the grid of crossover points is dense at high latitudes, the satellite data could be used for mapping the temporal and spatial variability of the current. Two and a half years of data were used to compute the statistical structure of the variability. The statistics included the probability distribution functions for each component of the current, the time-lagged autocorrelation functions of the variability, and the Reynolds stress produced by the variability. The results demonstrate that stress is correlated with bathymetry. In some areas the distribution of negative stress indicate that eddies contribute to an acceleration of the mean flow, strengthening the hypothesis that baroclinic instability makes important contributions to strong oceanic currents.

Johnson, Thomas J.↗
NASA NTRS1974

The helium component of solar wind velocity streams

Systematic variations of the properties of the helium constituent of the solar wind in the velocity streams are described. It is found that the helium abundance varies by about a factor of 2 as the stream is crossed. The velocity of the helium differs from that of the hydrogen by a few kilometers per second throughout much of the stream structure. This velocity difference is greatest immediately after the proton density peak passes, the helium velocity being typically 20 km/sec faster than the protons at that position in the stream. A sharp dip in the helium to proton temperature ratio is centered on the proton density peak. Although it appears reasonable that at least the velocity and temperature effects are due to the dynamic interactions of the two streams, it is not yet clear exactly what physical processes are directly involved in producing the effects described here.

Hirshberg, J.↗
NASA NTRS2009

Evaluation of Carbon Dioxide Sensors for the Constellation Space Suit Life Support System for Surface Exploration

This paper presents the findings of the trade study to evaluate carbon dioxide (CO2) sensing technologies for the Constellation (Cx) space suit life support system for surface exploration. The trade study found that nondispersive infrared absorption (NDIR) is the most appropriate high Technology Readiness Level (TRL) technology for the CO2 sensor for the Cx space suit. The maturity of the technology is high, as it is the basis for the CO2 sensor in the Extravehicular Mobility Unit (EMU). The study further determined that while there is a range of commercial sensors available, the Cx CO2 sensor should be a new design. Specifically, there are light sources (e.g., infrared light emitting diodes) and detectors (e.g., cooled detectors) that are not in typical commercial sensors due to cost. These advanced technology components offer significant advantages in performance (weight, volume, power, accuracy) to be implemented in the new sensor. The exact sensor design (light source, transmitting optics, path length, receiving optics and detector) will be specific for the Cx space suit and will be determined by the performance requirements of the Cx space suit. The paper further identifies specifications for some of the critical performance parameters as well as discussing the engineering aspects of implementing the sensor into the Portable Life Support System (PLSS). The paper then presents testing results from three CO2 sensors with respect to issues important to Extravehicular Activity (EVA) applications; stability, humidity dependence and low pressure compatibility. The three sensors include two NDIR sensors, one commercial and one custom-developed by NASA (for a different purpose), and one commercial electrochemical sensor. The results show that both NDIR sensors have excellent stability, no dependence on ambient humidity (when the ambient temperature is above the dew point) and operate in low pressure conditions and after being exposed to a full vacuum. The commercial electrochemical sensor was not suitable for the Cx space suit for surface exploration. Finally, the paper identifies a number of techniques currently under development that offer significant advantages for EVA applications. These include miniaturized, room temperature, solid electrolyte systems and advanced optical detectors.

Dietrich, Daniel L.↗
DOE OSTI2022

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

Early Energetic Particle Irradiation of the HED Parent Body Regolith

Previous studies have shown that many individual grains within the dark phase of the Kapoeta howardite were irradiated with energetic particles while residing on the surface of the early HED regolith. Particle tracks in these grains vary in density by more than an order of magnitude and undoubtedly were formed by energetic heavy (Fe) ions associated with early solar flares. Early Irradiation of HED Regolith: Concentrations of excess Ne alone are not sufficient to decide between competing galactic and solar irradiation models. However, from recent studies of depth samples of oriented lunar rocks, we have shown that the cosmogenic 21-Ne/22-Ne ratio produced in feldspar differs substantially between Galactic Cosmic Radiation (GCR) and solar protons, and that this difference is exactly that predicted from cross-section data. Using Ne literature data and new isotopic data we obtained on acid-etched, separated feldspar from both the light and dark phases of Kapoeta, we derive 21-Ne/22-Ne = 0.80 for the recent GCR irradiation and 21-Ne/22-Ne = 0.68 for the early regolith irradiation. This derived ratio indicates that the early Ne production in the regolith occurred by both galactic and solar protons. If we adopt a likely one-component regolith model in which all grains were exposed to galactic protons but individual grains had variable exposure to solar protons, we estimate that this early GCR irradiation lasted for about 3-6 m.y. More complex two-component regolith models involving separate solar and galactic irradiation would permit this GCR age to be longer. Higher-energy solar protons would permit the GCR to be longer. Higher-energy solar protons would permit the GCR age to be shorter. Further, cosmogenic 126(Xe) in Kapoeta dark is no more than a factor of about 2 higher than that observed in Kapoeta light. Because 126(Xe) can only be formed by galactic protons and not solar protons, these data support a short GCR irradiation for the HED regolith. This would also be the maximum time peRiod for the solar irradiation. Various asteroidal regolith models, based on Monte Carlo modeling of impact rates as a function of size and on irradiation features of meteorites, predict surface exposure times of about 0.1 to 10 m.y., and depend on such factors as gravity, rock mechanical properties, and micrometeoroid flux. Because the depth at which solar Fe tracks are produced (is much less than 1 micrometer) is much less than the depth at which Solar Cosmic Rays (SCR) Ne is produced (about 1 cm), for a reasonably well-stirred HED regolith the "surface exposure time" for SCR 21-Ne production should be significantly longer than that for solar tracks and some other surface irradiation features. Enhanced Solar Proton Irradiation: For bulk samples of Kapoeta dark feldspar and a one-component regolith model, the derived ratio of 21-Ne/22-Ne = 0.68 implies that the early production ratio of SCR 21-Ne to GCR 21-Ne was about 0.5-1.5. This ratio is independent of any assumptions about the fraction of dark grains that are irradiated or of the variability in the degree of solar irradiation among grains. The 21-Ne SCR/GCR ratio indirectly derived from bulk Kapoeta pyroxene is somewhat larger, as is the ratio derived for simple two-component regolith models. Individual feldspar grains that were extensively solar irradiated would require even larger 21-Ne SCR/GCR production ratios. In contrast, the theoretical SCR/GCR production ratio for lunar feldspar with 0 g/CM2 shield ing is is less than or equal to 2, and the lowest ratio observed in near-surface samples of lunar anorthosites is less than or equal to 1. Considering the greater solar distance of Vesta (compared to the Moon), the likelihood that SCR 21-Ne was acquired under some shielding where production rates are lower, and the likelihood that the exposure time to galactic protons exceeded the exposure time to solar protons because of their very different penetration depths, the 21-Ne SCR/GCR production ratio on the HED parent body was probably < 0.1. The relatively large difference between the derived 21-Ne SCR/GCR ratio in Kapoeta dark feldspar and the estimated production ratio strongly indicates that the early solar irradiation involved a flux -20-50x the recent solar flux. This enhanced proton flux was probably associated with an overall greater solar activity in the first approximately 10(exp 7) to 10(exp 8) years of solar history.

Bogard, D. D.↗
NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Efforts in transferring computational work from the LRC computer to the IRIS Graphics Workstation at MSU are reported and the computation of a conservative solution of a simple hyperbolic equation on an overlapping grid is discussed. Several conclusions concerning computations on overlapping grids are apparent. Problems only occur when there is a major difference in grid spacing on the individual component grids. In the case of hyperbolic equations, it is necessary that both interpolation and extrapolation be applied at the grid boundaries. When interpolated values are used at outflow boundary points, excessive oscillations in the numerical solution may be the result. The same conclusions would be valid for more complicated systems of hyperbolic equations such as the Euler equations for inviscid flow. Some of the solution values would be extrapolated at the overlap boundary, the exact number depending on the number of characteristics pointing out of the overlap region. It is also possible that similar boundary conditions may be needed for some parabolic equations such as high Reynolds number viscous flow equations. Efforts were also expended on the development of three-dimensional conservative interpolation procedures. Finally, the investigation of grid smoothing procedures were initiated during this reporting period. It was decided that the first grid smoothing algorithms will be based on the concepts of variational grid generation.

Mastin, C. W.↗
DOE OSTI2023

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

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

Liquid lithium as divertor material to mitigate severe damage of nearby components during plasma transients

The successful operation of thermonuclear fusion reactors such as ITER, DEMO, and future commercial plants is mainly determined by the optimum choice of materials for various components. The objective of this work is to accurately and comprehensively simulate the entire device in 3D to predict pros and cons of various materials, e.g., liquid lithium in comparison to tungsten and carbon to predict future ITER-like and DEMO divertor performances. We used our comprehensive HEIGHTS simulation package to investigate ITER-like components response during transient events in exact 3D geometry. Starting from the lost hot core plasma particles through SOL, deposition on the divertor surface, and the generation of secondary plasma of divertor materials. Our simulations predicted significant reduction in the heat loading and damage to the divertor nearby and internal components in the case when lithium is used on the divertor plates. While if tungsten or carbon are used on the divertor plate, significant melting areas and vaporization spots can occur (less for carbon) on the reflector, dome, and stainless steel tubes, and even parts of the first walls can melt due to the high radiation power of the secondary divertor plasma. Lithium photon radiation deposition into the divertor and nearby surfaces was decreased by two orders of magnitude compared to tungsten and by one order of magnitude compared to carbon. This analysis showed that using liquid lithium for ITER-like surfaces and future DEMO can lead to significant enhancement in components lifetime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

The prediction of zenith range refraction from surface measurements of meteorological parameters

In the last two decades, increasingly sophisticated deep space missions have placed correspondingly stringent requirements on navigational accuracy. As part of the effort to increase navigational accuracy, and hence the quality of radiometric data, much effort has been expended in an attempt to understand and compute the tropospheric effect on range (and hence range rate) data. The general approach adopted has been that of computing a zenith range refraction, and then mapping this refraction to any arbitrary elevation angle via an empirically derived function of elevation. The prediction of zenith range refraction derived from surface measurements of meteorological parameters is presented. Refractivity is separated into wet (water vapor pressure) and dry (atmospheric pressure) components. The integration of dry refractivity is shown to be exact. Attempts to integrate wet refractivity directly prove ineffective; however, several empirical models developed by the author and other researchers at JPL are discussed. The best current wet refraction model is here considered to be a separate day/night model, which is proportional to surface water vapor pressure and inversely proportional to surface temperature. Methods are suggested that might improve the accuracy of the wet range refraction model.

Berman, A. L.↗
NASA NTRS1992

The effects of acid deposition on sulfate reduction and methane production in peatlands

Peatlands, as fens and bods, make up a large percentage of northern latitude terrestrial environments. They are organic rich and support an active community of anaerobic bacteria, such as methanogenic and sulfate-reducing bacteria. The end products of these microbial activities, methane and hydrogen sulfide, are important components in the global biogeochemical cycles of carbon and sulfur. Since these two bacterial groups compete for nutritional substrates, increases in sulfate deposition due to acid rain potentially can disrupt the balance between these processes leading to a decrease in methane production and emission. This is significant because methane is a potent greenhouse gas that effects the global heat balance. A section of Mire 239 in the Experimental Lakes Area, in Northwestern Ontario, was artificially acidified and rates of sulfate reduction and methane production were measured with depth. Preliminary results suggested that methane production was not affected immediately after acidification. However, concentrations of dissolved methane decreased and dissolved sulfide increased greatly after acidification and both took several days to recover. The exact mechanism for the decrease in methane was not determined. Analyses are under way which will be used to determine rates of sulfate reduction. These results will be available by Spring and will be discussed.

Murray, Georgia L.↗
NASA NTRSDate not supplied

SimUAM: A Comprehensive Microsimulation Toolchain to Evaluate the Impact of Urban Air Mobility in Metropolitan Areas

Over the past several years, Urban Air Mobility (UAM) has galvanized enthusiasm from investors and researchers, marrying expertise in aircraft design, transportation, logistics, artificial intelligence, battery chemistry, and broader policymaking. However, two significant questions remain unexplored: (1) What is the value of UAM in a region’s transportation network?, and (2) How can UAM be effectively deployed to realize and maximize this value to all stakeholders, including riders and local economies? To adequately understand the value proposition of UAM for metropolitan areas, we develop a holistic multi-modal toolchain, SimUAM, to model and simulate UAM and its impacts on travel behavior. This toolchain has several components: (1) MANTA: A fast, high-fidelity regional-scale traffic microsimulator, (2) VertiSim: A granular, discrete-event vertiport and pedestrian, (3) 3: A high-fidelity, trajectory-based aerial microsimulation. SimUAM, rooted in granular, GPU-based microsimulation, models millions of trips and their exact movements in the street network and in the air, producing interpretable and actionable performance metrics for UAM designs and deployments. The modularity, extensibility, and speed of the platform will allow for rapid scenario planning and sensitivity analysis, effectively acting as a detailed performance assessment tool. As a result, stakeholders in UAM can understand the impacts of critical infrastructure, and subsequently define policies, requirements, and investments needed to support UAM as a viable transportation mode.

Urban air mobility↗
DOE OSTI2020

Light cone dynamics in excitonic states of two-component Bose and Fermi gases

Herein we consider the non-equilibrium dynamics of two-component one dimensional quantum gases in the limit of extreme population imbalance where the minority species has but a single particle. We consider the situation where the gas is prepared in a state with a single spatially localized exciton: the single particle of the minority species is spatially localized while the density of the majority species in the vicinity of the minority particle sees a depression. Remarkably, we are able to consider cases where the gas contains on the order of $N$ = 100 particles, comparable to that studied in experiments on cold atomic gases. We are able to do so by exploiting the integrability of the gas together with the observation that the excitonic state can be constructed through a simple superposition of exact eigenstates of the gas. The number of states in this superposition, rather than being exponentially large in the number of particles, scales linearly with $N$. We study the evolution of such spatially localized states in both strongly interacting Bose and Fermi gases. The behavior of the light cones when the interaction strength and density of the gas is varied can be understood from exact results for the spin excitation spectrum in these systems. We argue that the light cone in both cases exhibits scaling collapse. However unique to the Bose gas, we show that the presence of gapped finite-momentum roton-like excitations provide the Bose gas dynamics with secondary light cones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

A computer program for anisotropic shallow-shell finite elements using symbolic integration

A FORTRAN computer program for anisotropic shallow-shell finite elements with variable curvature is described. A listing of the program is presented together with printed output for a sample case. Computation times and central memory requirements are given for several different elements. The program is based on a stiffness (displacement) finite-element model in which the fundamental unknowns consist of both the displacement and the rotation components of the reference surface of the shell. Two triangular and four quadrilateral elements are implemented in the program. The triangular elements have 6 or 10 nodes, and the quadrilateral elements have 4 or 8 nodes. Two of the quadrilateral elements have internal degrees of freedom associated with displacement modes which vanish along the edges of the elements (bubble modes). The triangular elements and the remaining two quadrilateral elements do not have bubble modes. The output from the program consists of arrays corresponding to the stiffness, the geometric stiffness, the consistent mass, and the consistent load matrices for individual elements. The integrals required for the generation of these arrays are evaluated by using symbolic (or analytic) integration in conjunction with certain group-theoretic techniques. The analytic expressions for the integrals are exact and were developed using the symbolic and algebraic manipulation language.

Andersen, C. M.↗