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 37 records · Page 2

Tilt rotor hover aeroacoustics

The methodology, results, and conclusions of a study of tilt rotor hover aeroacoustics and aerodynamics are presented. Flow visualization and hot wire velocity measurement were performed on a 1/12-scale model of the XV-15 Tilt Rotor Aircraft in hover. The wing and fuselage below the rotor cause a complex recirculating flow. Results indicate the physical dimensions and details of the flow including the relative unsteadiness and turbulence characteristics of the flow. Discrete frequency harmonic thickness and the loading noise mechanism were predicted using WOPWOP for the standard metal blades and the Advanced Technology Blades. The recirculating flow created by the wing below the rotor is a primary sound mechanism for a hovering tilt rotor. The effects of dynamic blade response should be included for fountain flow conditions which produce impulsive blade loading. Broadband noise mechanisms were studied using Amiet's method with azimuthally varying turbulence characteristics derived from the measurements. The recirculating fountain flow with high turbulence levels in the recirculating zone is the dominant source of broadband noise for a hovering rotor. It is shown that tilt rotor hover aeroacoustic noise mechanisms are now understood. Noise predictions can be made based on reasonably accurate aerodynamic models developed here.

Coffen, Charles David↗

Advances in tilt rotor noise prediction

The two most serious tilt rotor external noise problems, hover noise and blade-vortex interaction noise, are studied. The results of flow visualization and inflow velocity measurements document a complex, recirculating highly unsteady and turbulent flow due to the rotor-wing-body interactions characteristic of tilt rotors. The wing under the rotor is found to obstruct the inflow, causing a deficit in the inflow velocities over the inboard region of the rotor. Discrete frequency harmonic thickness and loading noise mechanisms in hover are examined by first modeling tilt rotor hover aerodynamics and then applying various noise prediction methods using the WOPWOP code. The analysis indicates that the partial ground plane created by the wing below the rotor results in a primary sound source for hover.

George, A. R.↗

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗

Comparison between solar wind latitude distribution derived from Lyman-alpha observations and Ulysses measurements

The all-sky interplanetary Lyman-alpha pattern is sensitive to the latitude distribution of the solar wind because of destruction of neutral H by charge-exchange with solar wind protons. Lyman-alpha intensities recorded by Prognoz 5 and 6 in 1976 in a few parts of the sky were demonstrating a decrease of solar wind mass flux by about 30 % from equator to pole, when assuming a sinusoidal variation of this mass flux (harmonic distribution). A new analysis with a discrete variation with latitude has shown a decrease from 0 to 30 deg and then a plateau of constant mass flux up to the pole. This distribution bears a striking resemblance with Ulysses in-situ measurements, showing a clear similarity at 19 years interval. The Ulysses measurements were then used as a model input to calculate an all-sky Lyman-alpha pattern, either with a discrete model or with a harmonic solar wind variation with the same Ulysses equator-to-pole variation. There are conspicuous differences between the two Lyman-alpha patterns, in particular in the downwind region which are discussed in the context of future all-sky measurements with SWAN experiment on SOHO.

Quemarais, E.↗

Spectral decontamination of a real-time helicopter simulation

Nonlinear mathematical models of a rotor system, referred to as rotating blade-element models, produce steady-state, high-frequency harmonics of significant magnitude. In a discrete simulation model, certain of these harmonics may be incompatible with realistic real-time computational constraints because of their aliasing into the operational low-pass region. However, the energy is an aliased harmonic may be suppressed by increasing the computation rate of an isolated, causal nonlinearity and using an appropriate filter. This decontamination technique is applied to Sikorsky's real-time model of the Black Hawk helicopter, as supplied to NASA for handling-qualities investigations.

Mcfarland, R. E.↗

Comparison of measured and calculated sound pressure levels around a large horizontal axis wind turbine generator

Results are reported from a large number of simultaneous acoustic measurements around a large horizontal axis downwind configuration wind turbine generator. In addition, comparisons are made between measurements and calculations of both the discrete frequency rotational harmonics and the broad band noise components. Sound pressure time histories and noise radiation patterns as well as narrow band and broadband noise spectra are presented for a range of operating conditions. The data are useful for purposes of environmental impact assessment.

Shepherd, Kevin P.↗

Discrete modes and continuous spectra in supersonic boundary layers

The disturbance field induced due to a harmonic point source consists of discrete eigenmodes and a continuous spectrum; these are studied by using generalized Fourier transform techniques. For a supersonic boundary layer, there exist seven branches of the continuous spectrum in the complex wavenumber space, four of which (two acoustic waves, one vorticity wave and one entropy wave) contribute to the flowfield downstream of the source. The discrete eigenmodes spring off from these branches at some critical Reynolds numbers. The results for Mach 2 and 4.5 boundary layers show that the receptivity coefficients for the stable discrete modes are much larger than that for the unstable mode. Therefore, the flow very near the source is dominated by the continuous spectrum and the stable discrete modes. However, the unstable mode takes over sufficiently far away from the source. It is shown that it is only necessary to consider the first few discrete modes to construct the solution. Calculations also show that, in a supersonic boundary layer, upstream influence from a localized disturbance is minimal.

Balakumar, P.↗

Determination of rotor harmonic blade loads from acoustic measurements

The magnitude of discrete frequency sound radiated by a rotating blade is strongly influenced by the presence of a nonuniform distribution of aerodynamic forces over the rotor disk. An analytical development and experimental results are provided for a technique by which harmonic blade loads are derived from acoustic measurements. The technique relates, on a one-to-one basis, the discrete frequency sound harmonic amplitudes measured at a point on the axis of rotation to the blade-load harmonic amplitudes. This technique was applied to acoustic data from two helicopter types and from a series of test results using the NASA-Langley Research Center rotor test facility. The inferred blade-load harmonics for the cases considered tended to follow an inverse power law relationship with harmonic blade-load number. Empirical curve fits to the data showed the harmonic fall-off rate to be in the range of 6 to 9 db per octave of harmonic order. These empirical relationships were subsequently used as input data in a compatible far field rotational noise prediction model. A comparison between predicted and measured off-axis sound harmonic levels is provided for the experimental cases considered.

Kasper, P. K.↗

Comparison of Radiation Models for a Turbulent Piloted Methane/Air Jet Flame: A Frozen-Field Study

Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.

David, Chloe↗

SCEPTRE 2.1 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general purpose C++ code for solving the Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates and spherical harmonics. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL's are provided. The TPL's needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autoconf build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

97 MATHEMATICS AND COMPUTING↗

SCEPTRE 2.2 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general- purpose C++ code for solving the li near Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates and spherical harmonics. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Part y Libraries (TPL) to be available, and example scripts for building these TPLs are provided. The TPLs needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autotools build system , and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Extending radiative transfer models by use of Bayes rule

This paper presents a procedure that extends some existing radiative transfer modeling techniques to problems in atmospheric science where curvature and layering of the medium and dynamic range and angular resolution of the signal are important. Example problems include twilight and limb scan simulations. Techniques that are extended include successive orders of scattering, matrix operator, doubling, Gauss-Seidel iteration, discrete ordinates and spherical harmonics. The procedure for extending them is based on Bayes' rule from probability theory.

Whitney, C.↗

N-stream approximations to radiative transfer

Schuster's two-stream approximation (1905) is first derived from Chandrasekhar's radiative transfer equation (1950), and then extended to an arbitrary number of streams. The resulting technique for solving the transfer function similar to the discrete ordinate and spherical harmonic methods, is useful for modeling atmospheres with complicated phase functions and moderate optical depths. The resulting n coupled linear differential equations are simple and consume less computer time than other approximations, yet have the same required accuracy. The approximation is also flexible with respect to the choice of patch functions, and no approximations are made on the form of the phase function, other than its expansion into Legendre polynomials. A four-stream approximation is evaluated for a Henyey-Greenstein phase function with an asymmetry factor equal to 0.5.

Acquista, C.↗

Kravchuk functions for the finite oscillator approximation

Kravchuk orthogonal functions - Kravchuk polynomials multiplied by the square root of the weight function - simplify the inversion algorithm for the analysis of discrete, finite signals in harmonic oscillator components. They can be regarded as the best approximation set. As the number of sampling points increases, the Kravchuk expansion becomes the standard oscillator expansion.

Atakishiyev, Natig M.↗

SCEPTRE 2.3 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general- purpose C++ code for solving the linear Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates (Sn) and spherical harmonics (Pn). Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL are provided. The TPL needed by SCEPTRE are Trilinos, Boost, and Netcdf. SCEPTRE uses an autotools build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Helicopter crew/passenger vibration sensitivity -

Helicopter crew and passenger vibration sensitivity are presented. Pilot subjective ratings are established for discrete frequencies and the impact of combinations of harmonic frequencies is examined. A passenger long term comfort level and a short term limit are defined for discrete frequencies and compared with pilot ratings. The results show reasonable agreement between pilot and passenger. Subjective comfort levels obtained for mixed frequency environments clearly demonstrate the need for a multi-frequency criterion.

Gabel, R.↗

Adaptive Angular Quadrature Scheme for a backwards-in-time Method of Characteristics Solution to the Radiative Transfer Equation [Slides]

Radiative transfer/radiation transport are important problems to solve in astrophysics and high energy density physics. Various methods exist to solve radiation transport, such as Monte Carlo (MC), Discrete Ordinates (S N ), Method of Characteristics (MOC), and the spherical harmonics (P N ) method. Method of Characteristics requires “launching” of rays in discrete directions. Unresolved details of angular mesh create ray effects and can miss sources in the domain. Ray effects can lead to unphysical “stepping” in solution and incorrect energy deposition. Adaptive quadrature schemes can be used to detect and mitigate these effects. The Method of Characteristics (MOC) is a common method for solving hyperbolic PDEs in radiation transport and supersonic flow problems. Generally in MOC for radiation transport, virtual particles are tracked from birth to the end of a timestep. This requires interpolation to go from final location to cell averaged or corner values of angular intensity. Backwards-in-Time (BIT) particle tracking avoids this by prescribing the final position of the virtual particle at the cell nodes/corners. Angular intensities are computed at time k + 1 by launching ray back to previous timestep(s), or t = 0. Scheme allows solution to be computed as the characteristic ray is traced backwards in time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗