Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HARMONIC POTENTIAL”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Third-order solution of an artificial-satellite theory

A third-order solution is developed for the motions of artificial satellites moving in the gravitational field of the earth, whose potential includes the second-, third-, and fourth-order zonal harmonics. Third-order periodic perturbations with fourth-order secular perturbations are derived by Hori's perturbations method. All quantities are expanded into power series of the eccentricity, but the solution is obtained so as to be closed with respect to the inclination. A comparison with the results of numerical integration of the equations of motion indicates that the solution can predict the position of a close-earth satellite with a small eccentricity with an accuracy of better than 1 cm over 1 month.

Kinoshita, H.↗
NASA NTRS1981

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.↗
NASA NTRS1983

Geomagnetism of earth's core

Instrumentation, analytical methods, and research goals for understanding the behavior and source of geophysical magnetism are reviewed. Magsat, launched in 1979, collected global magnetometer data and identified the main terrestrial magnetic fields. The data has been treated by representing the curl-free field in terms of a scalar potential which is decomposed into a truncated series of spherical harmonics. Solutions to the Laplace equation then extend the field upward or downward from the measurement level through intervening spaces with no source. Further research is necessary on the interaction between harmonics of various spatial scales. Attempts are also being made to analytically model the main field and its secular variation at the core-mantle boundary. Work is also being done on characterizing the core structure, composition, thermodynamics, energetics, and formation, as well as designing a new Magsat or a tethered satellite to be flown on the Shuttle.

Benton, E. R.↗
NASA NTRS1989

Optimal trajectories for the Aeroassisted Flight Experiment. Part 1: Equations of motion in an Earth-fixed system

The determination of optimal trajectories for the aeroassisted flight experiment (AFE) is discussed. The AFE refers to the study of the free flight of an autonomous spacecraft, shuttle-launched and shuttle-recovered. Its purpose is to gather atmospheric entry environmental data for use in designing aeroassisted orbital transfer vehicles (AOTV). It is assumed that: (1) the spacecraft is a particle of constant mass; (2) the Earth is rotating with constant angular velocity; (3) the Earth is an oblate planet, and the gravitational potential depends on both the radial distance and the latitude (harmonics of order higher than four are ignored); and (4) the atmosphere is at rest with respect to the Earth. Under these assumptions, the equations of motion for hypervelocity atmospheric flight (which can be used not only for AFE problems, but also for AOT problems and space shuttle problems) are derived in an Earth-fixed system. Transformation relations are supplied which allow one to pass from quantities computed in an Earth-fixed system to quantities computed in an inertial system, and vice versa.

Miele, A.↗
NASA NTRS1989

Optimal trajectories for the Aeroassisted Flight Experiment. Part 2: Equations of motion in an inertial system

The determination of optimal trajectories for the aeroassisted flight experiment (AFE) is discussed. The AFE refers to the study of the free flight of an autonomous spacecraft, shuttle-launched and shuttle-recovered. Its purpose is to gather atmospheric entry environmental data for use in designing aeroassisted orbital transfer vehicles (AOTV). It is assumed that: (1) the spacecraft is a particle of constant mass; (2) the Earth is rotating with constant angular velocity; (3) the Earth is an oblate planet, and the gravitational potential depends on both the radial distance and the latitude (harmonics of order higher than four are ignored); and (4) the atmosphere is at rest with respect to the Earth. Under these assumptions, the equations of motion for hypervelocity atmospheric flight (which can be used not only for AFE problems, but also for AOT problems and space shuttle problems) are derived in an inertial system. Transformation relations are supplied which allow one to pass from quantities computed in an inertial system to quantities computed in an Earth-fixed system and vice versa.

Miele, A.↗
DOE OSTI2020

Will LISA Detect Harmonic Gravitational Waves from Galactic Cosmic String Loops?

Rapid advancement in the observation of cosmic strings has been made in recent years placing increasingly stringent constraints on their properties, with $G\mu\lesssim 10^{-11}$ from Pulsar Timing Array (PTA). Cosmic string loops with low string tension clump in the Galaxy due to slow loop decay and low gravitational recoil, resulting in great enhancement to loop abundance in the Galaxy. With an average separation of down to just a fraction of a kpc, and the total power of gravitational wave (GW) emission dominated by harmonic modes spanning a wide angular scale, resolved loops located in proximity are powerful, persistent, and highly monochromatic sources of GW with a harmonic signature not replicated by any other sources, making them prime targets for direct detection by the upcoming Laser Interferometer Space Antenna (LISA), whose frequency range is well-matched. Unlike detection of bursts where the detection rate scales with loop abundance, the detection rate for harmonic signal is the result of a complex interplay between the strength of GW emission, loop abundance, and other sources of noise, and is most suitably studied through numerical simulations. We develop a robust and flexible framework for simulating loops in the Galaxy for predicting direct detection of harmonic signal from resolved loops by LISA. Our simulation reveals that the most accessible region in the parameter space consists of large loops $\alpha=0.1$ with low tension $10^{-21}\lesssim G\mu\lesssim 10^{-19}$. Direct detection of field theory cosmic strings is unlikely, with the detection probability $p_{\mathrm{det}}\lesssim 2\%$ for a 1-year mission. An extension suggests that direct detection of cosmic superstrings with a low intercommutation probability is very promising. Searching for harmonic GW signal from resolved loops through LISA observations will potentially place physical constraints on string theory.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

An efficient finite element technique for sound propagation in axisymmetric hard wall ducts carrying high subsonic Mach number flows

The general non-linear three-dimensional equation for acoustic potential is derived by using a perturbation technique. The linearized axisymmetric equation is then solved by using a finite element algorithm based on the Galerkin formulation for a harmonic time dependence. The solution is carried out in complex number notation for the acoustic velocity potential. Linear, isoparametric, quadrilateral elements with non-uniform distribution across the duct section are implemented. The resultant global matrix is stored in banded form and solved by using a modified Gauss elimination technique. Sound pressure levels and acoustic velocities are calculated from post element solutions. Different duct geometries are analyzed and compared with experimental results.

Tag, I. A.↗
NASA NTRS1996

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Kreider, Kevin L.↗
NASA NTRS2023

Spacecraft Floating Potential Measurements for the Wind Spacecraft

Analysis of 8,804,545 electron velocity distribution functions, observed by the Wind spacecraft near 1 au between 2005 January 1 and 2022 January 1, was performed to determine the spacecraft floating potential, ∅ sc . Wind was designed to be electrostatically clean, which helps keep the magnitude of ∅ sc small (i.e., ∼5–9 eV for nearly all intervals) and the potential distribution more uniform. We observed spectral enhancements of ∅ sc at frequencies corresponding to the inverse synodic Carrington rotation period with at least three harmonics. The two-dimensional histogram of ∅ sc versus time also shows at least two strong peaks, with a potential third, much weaker peak. These peaks vary in time, with the intensity correlated with solar maximum. Thus, the spectral peaks and histogram peaks are likely due to macroscopic phenomena like coronal mass ejections (solar cycle dependence) and stream interaction regions (Carrington rotation dependence). The values of ∅ sc are summarized herein and the resulting data set is discussed.

Lynn B Wilson III↗
NASA NTRS2004

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗
DOE OSTI2025

Exploring the onset of collectivity approaching 𝑁=40 with masses of manganese isotopes

Isotopes in the region of the nuclear chart below 68 Ni have been the subject of intense experimental and theoretical effort due to the potential onset of a new “island of inversion” when crossing the harmonic oscillator subshell closure at 𝑁=40. Here, we have measured the masses of 64−68 Mn using TITAN's multiple-reflection time-of-flight mass spectrometer, resulting in the first precision mass measurements of 67 Mn and 68 Mn. These results are compared to ab initio calculations and modern shell model calculations and show an increase in collectivity approaching 𝑁=40.

binding energy & masses↗
NASA NTRS1992

A submillimeter tripler using a quasi-waveguide structure

A new type of frequency multiplier structure is being developed which is suitable for application at frequencies above 1 THz. This structure preserves some of the properties of waveguide for mode control, yet is not truly single mode. The device resembles a sectoral horn, with a varactor diode mounted near the throat. Input and output coupling are through the same aperture, requiring a quasi-optical diplexer. Initial tests are directed at building a tripler at 500 GHz, for comparison with waveguide structures. The diplexer is a blazed diffraction grating with appropriate focusing optics. Model studies show that the impedance match to a varactor should be good, and initial tests of the beam patterns of the prototype indicate that optical coupling efficiency should be very high. The structure also has the potential for use as a fundamental mixer, or as a third harmonic mixer.

Erickson, Neal R.↗
NASA NTRS1992

High resolution gravity models combining terrestrial and satellite data

Spherical harmonic expansions to degree 360 have been developed that combine satellite potential coefficient information, terrestrial gravity data, satellite altimeter information as a direct tracking data type and topographic information. These models define improved representations of the Earth's gravitational potential beyond that available from just satellite or terrestrial data. The development of the degree 360 models, however, does not imply a uniform accuracy in the determination of the gravity field as numerous geographic areas are devoid of terrestrial data or the resolution of such data is limited to, for example, 100 km. This paper will consider theoretical and numerical questions related to the combination of the various data types. Various models of the combination process are discussed with a discussion of various correction terms for the different models. Various sources of gravity data will be described. The new OSU91 360 model will be discussed with comparisons made to previous 360 models and to other potential coefficient models that are complete to degree 50. Future directions in high degree potential coefficient models will be discussed.

Rapp, Richard H.↗
NASA NTRS1994

Method for confining the magnetic field of the cross-tail current inside the magnetopause

A method is presented for analytically representing the magnetic field due to the cross-tail current and its closure on the magnetopause. It is an extension of a method used by Tsyganenko (1989b) to confine the dipole field inside an ellipsoidal magnetopause using a scalar potential. Given a model of the cross-tail current, the implied net magnetic field is obtained by adding to the cross-tail current field a potential field B = - del gamma, which makes all field lines divide into two disjoint groups, separated by the magnetopause (i.e., the combined field is made to have zero normal component with the magnetopause). The magnetopause is assumed to be an ellipsoid of revolution (a prolate spheroid) as an approximation to observations (Sibeck et al., 1991). This assumption permits the potential gamma to be expressed in spheroidal coordinates, expanded in spheroidal harmonics and its terms evaluated by performing inversion integrals. Finally, the field outside the magnetopause is replaced by zero, resulting in a consistent current closure along the magnetopause. This procedure can also be used to confine the modeled field of any other interior magnetic source, though the model current must always flow in closed circuits. The method is demonstrated on the T87 cross-tail current, examples illustrate the effect of changing the size and shape of the prescribed magnetopause and a comparison is made to an independent numerical scheme based on the Biot-Savart equation.

Sotirelis, T.↗
NASA NTRS1976

On the computation of the transonic perturbation flow field around two- and three-dimensional oscillating wings

A finite difference method for solving the unsteady flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. It is determined that there is a limit on reduced frequency, which is a function of Mach number and size of mesh region, above which the relaxation procedures will not converge. It is found that row line relaxation is more efficient than column relaxation and results are presented for a rectangular wing in harmonic pitch.

Weatherill, W. H.↗
DOE OSTI2022

Spatiotemporally resolved measurements of electric field around a piezoelectric transformer using electric-field induced second harmonic (E-FISH) generation

When a piezoelectric transformer (PT) is actuated at its second harmonic frequency by a low input voltage, the generated electric field at the distal end can be sufficient to breakdown the surrounding gas, making them attractive power sources for non-equilibrium plasma generation. Understanding the potential and electric fields produced in the surrounding medium by the PT is important for effectively designing and using PT plasma devices. In this work, the spatiotemporally resolved characteristics of the electric field generated by a PT operating in open air have been investigated using the femtosecond electric field-induced second harmonic generation (E-FISH) method. Here, electric field components were determined by simultaneously conducting E-FISH measurements with the incident laser polarized in two orthogonal directions relative to the PT crystal. Results of this work demonstrate the spatial distribution of electric field around the PT’s output distal end and how it evolves as a function of time. Notably, the strongest electric field appears on the face of the PT’s distal surface, near the top and bottom edges and decreases by approximately 70% over 3 mm. The time delay between the PT’s input voltage and measured electric field indicates that there is an about 0.45 phase difference between the PT’s input voltage and output signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗