Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete harmonics”

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 91 records · Page 5

Analysis of Satellite-to-Satellite Tracking (SST) and altimetry data from GEOS-C

Radar altimetry and satellite-to-satellite (SST) range and range rate tracking measurements were used to infer the exterior gravitational field of the earth and the structure of the geoid from GEOS-C metric data. Under the SST analysis, a direct point-by-point estimate of gravity disturbance by means of a recursive filter with backward smoothing was attempted but had to be forsaken because of poor convergence. The adopted representation consists of a more or less uniform grid of discrete masses at a depth of approximately 400 km from the earth's surface. The layer is superimposed on a spherical harmonics model. The procedure for smoothing the altimetry and inferring the fine-structured gravity field over the Atlantic test area is described. The local disturbances are represented by means of a density layer. The altimeter height biases were first estimated by a least squares adjustment at orbital crossover points. After taking out the bias, long wavelength contributions from GEM-6 as well as a calibration correction were subtracted. The residual heights were then represented by a mass distribution beneath the earth's surface.

Wong, L.↗

Improvements to the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.

97 MATHEMATICS AND COMPUTING↗

Characterizing Aeroelastic Systems Using Eigenanalysis, Explicitly Retaining The Aerodynamic Degrees of Freedom

Discrete time aeroelastic models with explicitly retained aerodynamic modes have been generated employing a time marching vortex lattice aerodynamic model. This paper presents analytical results from eigenanalysis of these models. The potential of these models to calculate the behavior of modes that represent damped system motion (noncritical modes) in addition to the simple harmonic modes is explored. A typical section with only structural freedom in pitch is examined. The eigenvalues are examined and compared to experimental data. Issues regarding the convergence of the solution with regard to refining the aerodynamic discretization are investigated. Eigenvector behavior is examined; the eigenvector associated with a particular eigenvalue can be viewed as the set of modal participation factors for that particular mode. For the present formulation of the equations of motion, the vorticity for each aerodynamic element appears explicitly as an element of each eigenvector in addition to the structural dynamic generalized coordinates. Thus, modal participation of the aerodynamic degrees of freedom can be assessed in M addition to participation of structural degrees of freedom.

Heeg, Jennifer↗

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Annoyance response to simulated advanced turboprop aircraft interior noise containing tonal beats

A study is done to investigate the effects on subjective annoyance of simulated advanced turboprop (ATP) interior noise environments containing tonal beats. The simulated environments consisted of low-frequency tones superimposed on a turbulent-boundary-layer noise spectrum. The variables used in the study included propeller tone frequency (100 to 250 Hz), propeller tone levels (84 to 105 dB), and tonal beat frequency (0 to 1.0 Hz). Results indicated that propeller tones within the simulated ATP environment resulted in increased annoyance response that was fully predictable in terms of the increase in overall sound pressure level due to the tones. Implications for ATP aircraft include the following: (1) the interior noise environment with propeller tones is more annoying than an environment without tones if the tone is present at a level sufficient to increase the overall sound pressure level; (2) the increased annoyance due to the fundamental propeller tone frequency without harmonics is predictable from the overall sound pressure level; and (3) no additional noise penalty due to the perception of single discrete-frequency tones and/or beats was observed.

Leatherwood, Jack D.↗

Structural response to discrete and continuous gusts of an airplane having wing bending flexibility and a correlation of calculated and flight results

An analysis is made of the structural response to gusts of an airplane having the degrees of freedom of vertical motion and wing bending flexibility and basic parameters are established. A convenient and accurate numerical solution of the response equations is developed for the case of discrete-gust encounter, an exact solution is made for the simpler case of continuous-sinusoidal-gust encounter, and the procedure is outlined for treating the more realistic condition of continuous random atmospheric turbulence, based on the methods of generalized harmonic analysis. Correlation studies between flight and calculated results are then given to evaluate the influence of wing bending flexibility on the structural response to gusts of two twin-engine transports and one four-engine bomber. It is shown that calculated results obtained by means of a discrete-gust approach reveal the general nature of the flexibility effects and lead to qualitative correlation with flight results. In contrast, calculations by means of the continuous-turbulence approach show good quantitative correlation with flight results and indicate a much greater degree of resolution of the flexibility effects.

Houbolt, John C↗

Sub Coulomb barrier d+ 208 Pb scattering in the time-dependent basis function approach

We employ the non-perturbative time-dependent basis function (tBF) approach to study the scattering of the deuteron on 208 Pb below the Coulomb barrier. We obtain the bound and discretized scattering states of the projectile, which form the basis representation of the tBF approach, by diagonalizing a realistic Hamiltonian in a large harmonic oscillator basis. We find that the higher-order inelastic scattering effects are noticeable for sub barrier scatterings with the tBF method. We have successfully reproduced experimental sub Coulomb barrier elastic cross section ratios with the tBF approach by considering only the electric dipole (E1) component of the Coulomb interaction between the projectile and the target during scatterings. Here, we find that the correction of the polarization potential to the Rutherford trajectory is dominant in reproducing the data at very low bombarding energies, whereas the role of internal transitions of the deuteron projectile induced by the E1 interaction during the scattering becomes increasingly significant at higher bombarding energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interference Control of Fano Resonances and Dynamical Imaging of an Electron Wave Packet

The Fano line shape, arising from the interference of pathways for the excitation of discrete and continuum states, plays a fundamental role in many branches of physics, chemistry, and materials science. Exciting the resonance with a high harmonic provides naturally a phase delay between the pathways leading to a complex asymmetry parameter. We demonstrate that its amplitude and phase can be controlled on the femtosecond and attosecond time scales, respectively. With our high-energy-resolution (10-meV) experiment, we dynamically image a resonance-enhanced electron wave packet during its temporal evolution, extracting both the amplitude and the phase. Calculations reproduce our experimental results. Our approach constitutes a method for measuring the photoionization delays of a resonance and enables the reconstruction of the electron wave packet in the time domain. This concept of an interference-controlled Fano line shape is a step toward attosecond quantum optics with potential ramifications into nanoscience and next-generation optical materials.

Han, Meng [ETH Zürich (Switzerland); Kansas State ↗

Vectorization and parallelization of the finite strip method for dynamic Mindlin plate problems

The finite strip method is a semi-analytical finite element process which allows for a discrete analysis of certain types of physical problems by discretizing the domain of the problem into finite strips. This method decomposes a single large problem into m smaller independent subproblems when m harmonic functions are employed, thus yielding natural parallelism at a very high level. In this paper we address vectorization and parallelization strategies for the dynamic analysis of simply-supported Mindlin plate bending problems and show how to prevent potential conflicts in memory access during the assemblage process. The vector and parallel implementations of this method and the performance results of a test problem under scalar, vector, and vector-concurrent execution modes on the Alliant FX/80 are also presented.

Chen, Hsin-Chu↗

Spectral Analysis of Vibrational Harmonic Motion by use of a Continuous-Wave CO2 Doppler Lidar

Vibrational motion of a harmonic oscillator was investigated using a focused continuous wave CO2 Doppler lidar at 9.1 microns wavelength. A continuum of frequencies along with many discrete, equally spaced, resonant frequency modes was observed. The frequency modes are similar in structure to the oscillatory longitudinal modes of a laser cavity and arise because of interference of the natural resonant frequency of the oscillator with specific frequencies within the continuum. The spectra revealed departures from linear motion for vigorous vibrations of the oscillator. Each consecutive resonant frequency mode occurred for a movement of the oscillator much less than the wavelength of incident lidar radiation.

Jarzembski, Maurice A.↗

Formation of cyclotron lines in gamma-ray burst spectra

A transmission model of gamma-ray burst sources is studied using the relativistic QED magnetic-resonant opacities including multiple photon scattering, incorporated into a discrete-ordinate radiative-transport scheme. The physics of the cyclotron line-producing region is discussed in general, and the expected line profiles, relative harmonic strengths, and polarizations are indicated under various conditions. The calculated spectra for these models show good agreement with the spectra reported from Ginga for GB 880205 and GB 870303.

Alexander, S. G.↗

A numerical modeling study of the interaction between the tides and the circulation forced by high-latitude plasma convection

A spectral, time-varying thermospheric general circulation model has been used to study the nonlinear interaction at high latitudes between the tides propagating into the thermosphere from below and the circulation induced by magnetospheric forcing and in situ solar heating. The model is discrete in the vertical with 27 layers spaced by half a scale height. In the horizontal, the fields are expanded in a series of spherical harmonics using a triangular truncation at wave number 31, equivalent to a homogeneous global resolution with a minimum wavelength of 1270 km. A hypothetical uniform grid point model would require a horizontal spacing of 417 km to describe the same minimum wavelength. In the high-latitude F region the tides affect the dusk vortex of the neutral flow very little, but the dawn vortex is either suppressed or amplified dependent upon the universal time and tidal phase. In the E region neutral flow, both the dusk and dawn vortices are shifted in local time by the tides, again as a function of universal time and tidal phase. At dusk a nonlinear amplification of the sunward winds occurs for certain combinations of parameters, and at dawn the winds may be completely suppressed. Below 120 km altitude the magnetospheric forcing creates a single cyclonic vortex which is also sensitive to the high-latitude tidal structure.

Mikkelsen, I. S.↗

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.↗

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

On the nonlinear development of the most unstable Goertler vortex mode TLSP - Final report

The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

A New Model of the Crustal Magnetic Field of Mars Using MGS and MAVEN

While devoid of an active magnetic dynamo field today, Mars possesses a remanent magnetic field that may reach several thousand nanoteslas locally. The exact origin and the events that have shaped the crustal magnetization remain largely enigmatic. Three magnetic field data sets from two spacecraft collected over 13 cumulative years have sampled the Martian magnetic field over a range of altitudes from 90 up to 6,000 km: (a) Mars Global Surveyor (MGS) magnetometer (1997–2006), (b) MGS Electron Reflectometer (1999–2006), and (c) Mars Atmosphere and Volatile EvolutioN (MAVEN) magnetometer (2014 to today). In this paper we combine these complementary data sets for the first time to build a new model of the Martian internal magnetic field. This new model improves upon previous ones in several aspects: comprehensive data coverage, refined data selection scheme, modified modeling scheme, discrete-to-continuous transformation of the model, and increased model resolution. The new model has a spatial resolution of ∼160 km at the surface, corresponding to spherical harmonic degree 134. It shows small scales and well-defined features, which can now be associated with geological signatures.

Magnetic↗

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Tbeodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modem three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Theodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modern three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗