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 127 records · Page 7

An investigation of rotor harmonic noise by the use of small scale wind tunnel models

Noise measurements of small scale helicopter rotor models were compared with noise measurements of full scale helicopters to determine what information about the full scale helicopters could be derived from noise measurements of small scale helicopter models. Comparisons were made of the discrete frequency (rotational) noise for 4 pairs of tests. Areas covered were tip speed effects, isolated rotor, tandem rotor, and main rotor/tail rotor interaction. Results show good comparison of noise trends with configuration and test condition changes, and good comparison of absolute noise measurements with the corrections used except for the isolated rotor case. Noise measurements of the isolated rotor show a great deal of scatter reflecting the fact that the rotor in hover is basically unstable.

Sternfeld, H., Jr.↗

Discrete rotor noise.

Description of a single unified theory to account generally for the discrete radiation from the whole family of heavily loaded axial flow rotors, including helicopter rotors, propellers, and compressors. The approach used is to analyze, in a Fourier way, the periodic blade loading variation (azimuth profile) into blade loading harmonics (BLH) and then sum the total radiation from each BLH. The basic differences between high and low solidity rotor spectrums are considered, emphasizing the difference between subsonic and supersonic rotor noise. The special effect of helicopter blade slap and compressor rotor-stator interaction is discussed, and measured rotor spectrums have been included to help illustrate the theory.

Wright, S. E.↗

Rotation of the photospheric magnetic fields: A north-south asymmetry

During most of solar cycle 21 the large-scale photospheric field rotated more rapidly in the Northern Hemisphere than in the southern. The large-scale northern field rotated with a 26.9 day period (synodic), was centered at 15 degress N, and covered a latitude zone about 24 degrees wide. The large-scale southern field rotated with a periodicity of 28.1 days, was centered at 26 degrees S, and covered a latitude zone about 32 degrees wide. Our analysis showed rotational power at only a few discrete latitudes and frequencies in each hemisphere. The center of each peak lies near the sunspot differential rotation curve. The largest scale field contributes to the configuration of the coronal and interplanetary magnetic field (IMF). The strength of the first harmonic of the northern field suggests that this structure may be related to the 4-sector pattern observed in the IMF polarity. The southern field had much lower power at the first harmonic of the solar rotation rate and so would contribute only to a 2-sector structure in the IMF. These results were discovered in Fourier analysis of photospheric synoptic charts obtained at the Wilcox Solar Observatory from 1976 to 1986 and confirmed in higher resolution maps from the National Solar Observatory. Mt. Wilson magnetic field measurements from solar cycle 20 show a similar north-south asymmetry.

Antonucci, E.↗

A Higher Harmonic Optimal Controller to Optimise Rotorcraft Aeromechanical Behaviour

Three methods to optimize rotorcraft aeromechanical behavior for those cases where the rotorcraft plant can be adequately represented by a linear model system matrix were identified and implemented in a stand-alone code. These methods determine the optimal control vector which minimizes the vibration metric subject to constraints at discrete time points, and differ from the commonly used non-optimal constraint penalty methods such as those employed by conventional controllers in that the constraints are handled as actual constraints to an optimization problem rather than as just additional terms in the performance index. The first method is to use a Non-linear Programming algorithm to solve the problem directly. The second method is to solve the full set of non-linear equations which define the necessary conditions for optimality. The third method is to solve each of the possible reduced sets of equations defining the necessary conditions for optimality when the constraints are pre-selected to be either active or inactive, and then to simply select the best solution. The effects of maneuvers and aeroelasticity on the systems matrix are modelled by using a pseudo-random pseudo-row-dependency scheme to define the systems matrix. Cases run to date indicate that the first method of solution is reliable, robust, and easiest to use, and that it was superior to the conventional controllers which were considered.

Leyland, Jane Anne↗

An unsteady lifting surface method for single rotation propellers

The mathematical formulation of a lifting surface method for evaluating the steady and unsteady loads induced on single rotation propellers by blade vibration and inflow distortion is described. The scheme is based on 3-D linearized compressible aerodynamics and presumes that all disturbances are simple harmonic in time. This approximation leads to a direct linear integral relation between the normal velocity on the blade (which is determined from the blade geometry and motion) and the distribution of pressure difference across the blade. This linear relation is discretized by breaking the blade up into subareas (panels) on which the pressure difference is treated as approximately constant, and constraining the normal velocity at one (control) point on each panel. The piece-wise constant loads can then be determined by Gaussian elimination. The resulting blade loads can be used in performance, stability and forced response predictions for the rotor. Mathematical and numerical aspects of the method are examined. A selection of results obtained from the method is presented. The appendices include various details of the derivation that were felt to be secondary to the main development in Section 1.

Williams, Marc H.↗

Rotation of the photospheric magnetic fields - A north-south asymmetry

During most of solar cycle 21 the large-scale photospheric field rotated more rapidly in the Northern Hemisphere than in the southern. The large-scale northern field rotated with a 26.9 day period (synodic), was centered at 15 degrees N, and covered a latitude zone about 24 degrees wide. The large-scale southern field rotated with a periodicity of 28.1 days, was centered at 26 degrees S, and covered a latitude zone about 32 degrees wide. Our analysis showed rotational power at only a few discrete latitudes and frequencies in each hemisphere. The center of each peak lies near to the sunspot differential rotation curve. The largest scale field contributes to the configuration of the coronal and interplanetary magnetic field (IMF). The strength of the first harmonic of the northern field suggests that this structure may be related to the 4-sector pattern observed in the IMF polarity. The southern field had much lower power at the first harmonic of the solar rotation rate and so would contribute only to a 2-sector structure in the IMF. These results were discovered in Fourier analysis of photospheric synoptic charts obtained at the Wilcox Solar Observatory from 1976 to 1986 and confirmed in higher resolution maps from the National Solar Observatory. Mt. Wilson magnetic field measurements from solar cycle 20 show a similar north-south asymmetry.

Antonucci, E.↗

VLF line radiation in the earth's magnetosphere and its association with power system radiation

In a recent experiment, discrete VLF emissions from the magnetosphere were triggered by a transmitter at Siple Station in Antarctica. Spectrograms of these signals as received at the conjugate point, Roberval, Quebec, showed changes in slope, entrainments, and cutoffs at frequencies (several kilohertz) close to the harmonic induction lines from the local 60-Hz power system. This observation led to the suggestion that harmonic radiation from the power system enters the magnetosphere and interacts with the triggered emissions. New evidence supporting this suggestion has been found in spectrograms of simultaneous recordings made at Roberval and at Siple Station in Antarctica. It is shown that line radiation, near harmonics of 60 Hz, travels along the earth's magnetic field in the whistler mode and is received in the conjugate hemisphere at Siple Station. Echoing of the line radiation between Siple and Roberval is often observed. The magnetospheric lines are usually shifted in frequency by 20-30 Hz with respect to the adjacent induction line, but their spacings are near 120 Hz. They may trigger and cut off emissions as do signals from VLF transmitters.

Helliwell, R. A.↗

Stability analysis of large space structure control systems with delayed input

Large space structural systems, due to their inherent flexibility and low mass to area ratio, are represented by large dimensional mathematical models. For implementation of the control laws for such systems a finite amount of time is required to evaluate the control signals; and this time delay may cause instability in the closed loop control system that was previously designed without taking the input delay into consideration. The stability analysis of a simple harmonic oscillator representing the equation of a single mode as a function of delay time is treated analytically and verified numerically. The effect of inherent damping on the delay is also analyzed. The control problem with delayed input is also formulated in the discrete time domain.

Reddy, A. S. S. R.↗

Stability analysis of large space structure control systems with delayed input

Large space structural systems, due to their inherent flexibility and low mass to area ratio, are represented by large dimensional mathematical models. For implementation of the control laws for such systems a finite amount of time is required to evaluate the control signals; and this time delay may cause instability in the closed loop control system that was previously designed without taking the input delay into consideration. The stability analysis of a simple harmonic oscillator representing the equation of a single mode as a function of delay time is analyzed analytically and verified numerically. The effect of inherent damping on the delay is also analyzed. The control problem with delayed input is also formulated in the discrete time domain.

Reddy, A. S. S. R.↗

Conical Euler solution for a highly-swept delta wing undergoing wing-rock motion

Modifications to an unsteady conical Euler code for the free-to-roll analysis of highly-swept delta wings are described. The modifications involve the addition of the rolling rigid-body equation of motion for its simultaneous time-integration with the governing flow equations. The flow solver utilized in the Euler code includes a multistage Runge-Kutta time-stepping scheme which uses a finite-volume spatial discretization on an unstructured mesh made up of triangles. Steady and unsteady results are presented for a 75 deg swept delta wing at a freestream Mach number of 1.2 and an angle of attack of 30 deg. The unsteady results consist of forced harmonic and free-to-roll calculations. The free-to-roll case exhibits a wing rock response produced by unsteady aerodynamics consistent with the aerodynamics of the forced harmonic results. Similarities are shown with a wing-rock time history from a low-speed wind tunnel test.

Lee, Elizabeth M.↗

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗

A numerical method based on the Fourier-Fourier transform approach for modeling 1-D electron plasma evolution

A numerical method is presented for studying one-dimensional electron plasma evolution under typical interplanetary conditions. The method applies the Fourier-Fourier transform approach to a plasma model that is a generalization of the electrostatic Vlasov-Poisson system of equations. Conservation laws that are modified to include the plasma model generalization and also the boundary effects of nonperiodic solutions are given. A new conservation law for entropy in the transformed space is then introduced. These conservation laws are used to verify the numerical solutions. A discretization error analysis is presented. Two numerical instabilities and the methods used for their suppression are treated. It is shown that in interplanetary plasma conditions, the bump-on-tail instability produces significant excitation of plasma oscillations at the Bohm-Gross frequency and its second harmonic. An explanation of the second harmonic excitation is given in terms of wave-wave coupling during the growth phase of the instability.

Klimas, A. J.↗

Quadratic Optimization in the Problems of Active Control of Sound

We analyze the problem of suppressing the unwanted component of a time-harmonic acoustic field (noise) on a predetermined region of interest. The suppression is rendered by active means, i.e., by introducing the additional acoustic sources called controls that generate the appropriate anti-sound. Previously, we have obtained general solutions for active controls in both continuous and discrete formulations of the problem. We have also obtained optimal solutions that minimize the overall absolute acoustic source strength of active control sources. These optimal solutions happen to be particular layers of monopoles on the perimeter of the protected region. Mathematically, minimization of acoustic source strength is equivalent to minimization in the sense of L(sub 1). By contrast. in the current paper we formulate and study optimization problems that involve quadratic functions of merit. Specifically, we minimize the L(sub 2) norm of the control sources, and we consider both the unconstrained and constrained minimization. The unconstrained L(sub 2) minimization is certainly the easiest problem to address numerically. On the other hand, the constrained approach allows one to analyze sophisticated geometries. In a special case, we call compare our finite-difference optimal solutions to the continuous optimal solutions obtained previously using a semi-analytic technique. We also show that the optima obtained in the sense of L(sub 2) differ drastically from those obtained in the sense of L(sub 1).

Loncaric, J.↗

Electromagnetic Scattering by Discrete Random Media. I: The Dispersion Equation and the Configuration-Averaged Exciting Field

We consider the scattering of a plane electromagnetic wave obliquely incident on a plane-parallel layer of discrete random medium with non-scattering boundaries. We solve the Lax integral equation for the conditional configuration-averaged exciting field coefficients by assuming a special-form solution, that is, by representing the conditional configuration-averaged exciting field coefficients as a linear combination of the coefficients corresponding to an up-going and a down-going wave. This solution representation is supposed to be valid within the whole domain occupied by the particles, even in the close proximity of the boundaries. By balancing the waves with different propagation directions and wavenumbers we derive two homogeneous systems of equations corresponding to the generalized Lorenz–Lorentz law and two inhomogeneous systems of equations corresponding to the generalized Ewald–Oseen extinction theorem. It is shown that (i) the two homogeneous systems of equations of the generalized Lorenz–Lorentz law reduce to a single homogeneous system of equations corresponding to a semi-infinite discrete random medium at normal incidence; (ii) the dispersion equation is direction and polarization independent; and (iii) the two inhomogeneous systems of equations of the generalized Ewald–Oseen extinction theorem can be reduced to two scalar equations by means of the addition theorem for vector spherical harmonics. It is also shown that the same dispersion equation can be obtained without assuming a special-form solution representation in the proximity of the boundaries.

Adrian Doicu↗

Applications of the modified Hulthén-Kohn method for bound and scattering states

We adapt the Hulthén–Kohn method suggested by Efros [Phys. Rev. C 99, 034620 (2019)] for calculating various observables in the continuum and discrete spectrum using two-body interactions in single- and coupled channel systems. We explore the convergence of phase shifts and wave functions as well as the location of S-matrix poles which enables obtaining both resonance and bound state parameters. We find that employing a harmonic oscillator basis, together with an interaction smoothing scheme introduced by Gyarmati et al. [Nucl. Phys. A 326, 119 (1979)], and adopting approximate bound-state solutions for the short-range components of basis wave functions lead to good convergence even with restricted oscillator quanta accessible for modern no-core shell model codes. The adapted Efros method will facilitate ab initio many-body nuclear structure applications.

Nuclear reactions↗

Multi-Port Autonomous Reconfigurable Solar Power Plant

As the penetration level of power electronics increases and remote photovoltaic (PV) generation is integrated into the alternating current (ac) grid, the short-circuit ratio (SCR) at the point of interconnection of a hybrid PV–energy storage system (ESS) plant may be low. Additionally, the inertia of the alternating current (ac) grid may be low. The low SCRs and inertias can lead to reliability challenges in the power grid. These operating conditions require additional reinforcements, such as synchronous condensers, static var compensators, static synchronous compensators, and high-voltage direct current (HVdc) links/grids. HVdc links or grids may also provide the additional capability of connecting the plant to asynchronous grids and/or connecting asynchronous grids, among others. This scenario leads to discrete development of solar inverters, energy storage inverters, HVdc converters, and several transformers. Some of the problems associated with this discrete development and inverter-based generation include increased cost, lower reliability, and reduced efficiency associated with duplication of power electronics (PEs); competing controls of several individual discrete inverter-based generators due to the presence of multiple PEs, which leads to derating of the system; and transient stability problems arising from inverter-based generation, such as voltage/frequency events leading inverter shutdowns, voltage instability and control interactions in the formed weak grid, and harmonics caused by resonances of multiple inverters.

14 SOLAR ENERGY↗

Conical Euler methodology for unsteady vortical flows about rolling delta wings

A conical Euler methodology was developed to study unsteady vortex-dominated flows about rolling highly swept delta wings undergoing either forced or free-to-roll motions. The flow solver of the code involves a multistage Runge-Kutta time-stepping scheme which uses a finite-volume spatial discretization of the Euler equations on an unstructured grid of triangles. Results are presented for a 75-deg swept sharp-leading-edge delta wing at a freestream Mach number of 1.2 and at alpha = 10, 20, and 30 deg. At the 10 and 20 deg, forced harmonic analyses indicate that the rolling moment coefficients provide a positive damping which is verified by free-to-roll calculations. In contrast, at 30 deg, a forced harmonic analysis indicates that the rolling moment coefficient provides a negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation.

Lee, Elizabeth M.↗

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Here, we propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for nonlinear hyperbolic PDEs based on Burgers’ and Euler’s equations.

97 MATHEMATICS AND COMPUTING↗