SEARCH · Engineering Papers

Results for “ACOUSTIC SIMULATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

NASA NTRS2021

Anti-Phase Vortex Reduction Control for Rotor Noise Suppression

An investigation is conducted in 2019 under the NASA Ames Center Innovation Funds (CIF) project entitled “Anti-Phase Vortex Reduction Control for Rotor Noise Suppression“ to experimentally validate an anti-phase rotor noise suppression concept. The first objective of the investigation is to conduct computational fluid dynamics (CFD) simulations to investigate the noise characteristics of several anti-phase rotor designs. The second objective is to conduct a series of acoustic tests of anti-phase rotors in an anechoic chamber at Pennsylvania State University (PSU) to evaluate the merit of the anti-phase rotor concept. The CFD investigation seeks to optimize the anti-phase alternating trailing edge patterns for rotor noise suppression. The design objective is to maximize the noise reduction while maintaining the aerodynamic thrust. The investigation is performed using a three-dimensional (3D) Unsteady Reynolds-Averaged Navier-Stokes (URANS) commercial solver STAR-CCM+ together with the Ffowcs-Williams and Hawkings (FW-H) formula to obtain the aerodynamic thrust and far-field noise level. An acoustic study is conducted for 13 anti-phase design candidates based on a proprietary rotor design. The CAD geometry of the rotor is furnished by PSU. These design candidates include different alternating trailing edge (TE) waveforms, TE segment lengths, TE deflection amplitudes, and transition characteristics. The best design candidate among those explored is an anti-phase rotor that has a four-period TE waveform which results in a reduction in far-field noise level of 2.1 dB in the hover condition and a reduction of 1.1 dB in the forward flight condition at 9.7 m/s. A further acoustic study is conducted for a different rotor manufactured by APC. Five APC rotor design candidates are simulated. The best anti-phase design candidate for the APC rotor results in a reduction in far-field noise level of 4.0 dB in the hover condition. An in-phase design candidate is also studied. This in-phase design provides a noise reduction of 2.5 dB. A series of acoustic experiments in the PSU anechoic chamber have been conducted in July 2019 and October 2019. Both anti-phase and in-phase rotors fabricated for the left-hand and right-hand rotations are tested. In the hover condition, all the rotors do not provide sufficient evidence of improved acoustic performance. However, the experimental data in the hover condition are deemed to be inconclusive due to the flow recirculation in the anechoic chamber caused by the rotor downwash. In the forward flight condition at 9.7 m/s, the anti-phase right-hand rotor produces a noise reduction by as much as 6.5 dB in the frequency range of 2000-4000 Hz, while the in-phase right-hand rotor produces a noise reduction by as much as 5 dB in the same frequency range. Both the anti-phase and in-phase left hand rotors offer no evidence of noise reduction. The difference in the acoustic performance of the left-hand and right-hand rotors could be explained by the location of the microphone array which is placed to the left side of the rotors. This microphone array location could create a bias in the sound pressure level in favor of the right-hand rotor. Using the average noise reduction values, the anti-phase 4H rotor could offer a noise reduction by as much as 3.25 dB while the in-phase 4I rotor could produce up to a 2.5 dB noise reduction. Both the computational and experimental results have provided sufficient evidence to support the noise suppression capability of the proposed anti-phase rotor concepts.

Rotor noise↗

NASA NTRS2019

Simulation-Based Assessment of a Full-Scale Installed Quiet Landing Gear

Full-scale simulations of a Gulfstream G-III aircraft, performed in support of the NASA Acoustic Research Measurements flights, are presented to complement results discussed in earlier studies. The flow solver employed in those studies, Dassault Systèmes’ lattice Boltzmann PowerFLOW®, was also used during this investigation to conduct time-dependent simulations of the entire aircraft in landing configuration with a fully dressed landing gear. The high-fidelity simulations, performed at a Mach number of 0.23 and a Reynolds number of 10.5 × 106 based on mean aerodynamic chord, captured all relevant airframe noise sources. The computations were used to assess the aeroacoustic performance of the main landing gear, with and without noise reduction fairings installed, of a G-III aircraft equipped with Adaptive Compliant Trailing Edge technology and conventional Fowler flaps. To facilitate comparison of predicted noise signatures with effective perceived noise levels obtained from flight test measurements, the “as-flown” nose landing gear geometry, missing in our earlier studies, was added to the simulated G-III aircraft configurations. The high fidelity, synthetic data were post-processed using a Ffowcs-Williams and Hawkings integral approach to estimate farfield acoustic behavior, with pressures on the model solid surface or pressure and velocity components on a permeable surface enveloping the acoustic near field used as input.

Duda, Benjamin↗

NASA NTRS1995

A test of a vortex method for the computation of flap side edge noise

Upon approach to landing, a major source location of airframe noise occurs at the side edges of the part span, trailing edge flaps. In the vicinity of these flaps, a complex arrangement of spanwise flow with primary and secondary tip vortices may form. Each of these vortices is observed to become fully three-dimensional. In the present study, a numerical model is developed to investigate the noise radiated from the side edge of a flap. The inherent three-dimensionality of this flow forces us to carefully consider a numerical scheme which will be both accurate in its prediction of the flow acoustics and also computationally efficient. Vortex methods have offered a fast and efficient means of simulating many two and three-dimensional, vortex dominated flows. In vortex methods, the time development of the flow is tracked by following exclusively the vorticity containing regions. Through the Biot-Savart law, knowledge of the vorticity field enables one to obtain flow quantities at any desired location during the flow evolution. In the present study, a numerical procedure has been developed which incorporates the Lagrangian approach of vortex methods into a calculation for the noise radiated by a flow-surface interaction. In particular, the noise generated by a vortex in the presence of a flat half plane is considered. This problem serves as a basic model of flap edge flow. It also permits the direct comparison between our computed results and previous acoustic analyses performed for this problem. In our numerical simulations, the mean flow is represented by the complex potential W(z) = Aiz(exp l/2), which is obtained through conformal mapping techniques. The magnitude of the mean flow is controlled by the parameter A. This mean flow has been used in the acoustic analysis by Hardin and is considered a reasonable model of the flow field in the vicinity of the edge and away from the leading and trailing edges of the flap. To represent the primary vortex which occurs near the flap, a point vortex is introduced just below the flat half plane. Using a technique from panel methods, boundary conditions on the flap surface are satisfied by the introduction of a row of stationary point vortices along the extent of the flap. At each time step in the calculation, the strength of these vortices is chosen to eliminate the normal velocity at intermediary collocation points. The time development of the overall flow field is then tracked using standard techniques from vortex methods. Vortex trajectories obtained through this computation are in good agreement with those predicted by the analytical solution given by Hardin, thus verifying the viability of this procedure for more complex flow arrangements. For the flow acoustics, the Ffowcs Williams-Hawkings equation is numerically integrated. This equation supplies the far field acoustic pressure based upon pressures occurring along the flap surface. With our vortex method solution, surface pressures may be obtained with exceptional resolution. The Ffowcs Williams-Hawkings equation is integrated using a spatially fourth order accurate Simpson's rule. Rational function interpolation is used to obtain the surface pressures at the appropriate retarded times. Comparisons between our numerical results for the acoustic pressure and those predicted by the Hardin analysis have been made. Preliminary results indicate the need for an improved integration technique. In the future, the numerical procedure developed in this study will be applied to the case of a rectangular flap of finite thickness and ultimately modified for application to the fully three-dimensional problem.

Martin, James E.↗

DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

FOAS4GEO - Fiber Optic Acoustic Sensors for Geophysical Applications (Final Report)

The motivation was to enhance interpretation of signals measured by distributed fiber optic seismic sensors, with primary emphasis on subsurface applications in energy production such as hydrocarbon production and geothermal energy although applications also exist in areas such as explosion monitoring. The work consisted of two major elements: 1) computational-based simulation of subsurface signals, with emphasis on hydraulic stimulation, as observed by distributed Fiber Optic Acoustic Sensors (FOAS), and 2) development of a capability to conduct lab-based measurements of distributed fiber optic acoustic sensors using a custom interrogator. Both goals were successfully achieved. For the first topic, simulations of hydraulic fracturing as measured by distributed fiber optic sensors were modeling using a computational rock mechanics code. Work on the second major element consisted of the design, construction, and testing of a distributed fiber optic interrogator and development of testbed to evaluate custom fibers. The results of this three-year effort have led to multiple conference presentations, one per-reviewed publication, and several new externally funded projects on topics directly related to the LDRD. Two post-doctoral researchers worked on the project.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗

DOE OSTI2019

Nalu-Wind and OpenFAST: A high-fidelity modeling and simulation environment for wind energy. Milestone ECP-Q2-FY19

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. This milestone represents the culmination of several parallel development activities towards the goal of establishing a full-physics simulation capability for modeling wind turbines operating in turbulent atmospheric inflow conditions. The demonstration simulation performed in this milestone is the first step towards the "ground truth" simulation and includes the following components: neutral atmospheric boundary layer inflow conditions generated using a precursor simulation, a hybrid RANS/LES simulation of the wall-resolved turbine geometry, hybridization of the turbulence equations using a blending function approach to transition from the atmospheric scales to the blade boundary layer scales near the turbine, fluid-structure interaction (FSI) that accounts for the complete set of blade deformations (bending, twisting and pitch motion, yaw and tower displacements) by coupling to a comprehensive turbine dynamics code (OpenFAST). The use of overset mesh methodology for the simulations in this milestone presents a significant deviation from the previous efforts where a sliding mesh approach was employed to model the rotation of the turbine blades. The choice of overset meshes was motivated by the need to handle arbitrarily large deformations of the blade and to allow for blade pitching in the presence of a controller and the ease of mesh generation compared to the sliding mesh approach. FSI and the new timestep algorithm used in the simulations were developed in partnership with the A2e High-Fidelity Modeling project. The individual physics components were verified and validated (V%V) through extensive code-to-code comparisons and with experiments where possible. The detailed V&V efforts provide confidence in the final simulation where these physics models were combined together even though no detailed experimental data is available to perform validation of the final configuration. Taken together, this milestone successfully demonstrated the most advanced simulation to date that has been performed with Nalu-Wind.

17 WIND ENERGY↗

DOE OSTI2019

High-order methods for wind turbine simulations and moving to next-generation platforms

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. This milestone represents an effort to increase the fidelity of Nalu-Wind at a fixed resolution through the implementation of a tensor-product based, matrix-free high order scheme. High order finite element methods have increased local work per datum communicated and have the potential to provide significantly more accurate solutions at a fixed number of degrees of freedom. Previous to this milestone, Nalu-Wind had an arbitrary order Control Volume Finite Element Method discretization as a solver option, but it required too much memory and was too slow to be of practical use. The work in this milestone addresses these issues by first implementing an implicit, high order solver that only partially assembles the global system. This reduces the memory footprint of the high-order scheme by orders of magnitude for higher polynomial orders. Second, a faster, tensor-product based method for evaluating the action of the left-hand side was implemented. This reduces the amount of computational work required by the scheme and dramatically enhanced the time-to-solution on example problems. Finally, this milestone is an evaluation of the value of high order methods in the wind application space. With the enhancements to memory and computational cost, accuracy vs. time-to-solution was evaluated for several resolutions on an under-resolved Taylor Green vortex test case. Results show that the high order scheme is cost-competitive with the production low-order schemes in Nalu-Wind, being moderately more expensive than the production edge-based vertex centered finite volume scheme. The evaluation of accuracy on the test case shows a potential benefit to high order at the highest resolution while not deteriorating accuracy on the lowest tested resolution. More work is needed to show value in the wind application, but positive strides have been made.

17 WIND ENERGY↗

NASA NTRS1990

Resonant instability of supersonic shear layers

A computer simulation of possible resonant instability of a supersonic shear layer is carried out. The resonance of two acoustic duct modes of the flow induced by periodic Mach waves generated by a wavy wall is sought. Results of the simulations are reported. Simulations are unable to document a resonant instability and the mixing characteristics remain unchanged. Possible weakness of the present simulations are discussed. A second set of simulations involving a mixing layer separating a supersonic and a subsonic stream were performed. A wavy wall placed adjacent to the supersonic stream to produce a set of periodic Mach waves terminating at the shear layer is modelled. The entire flow field is similar to that of an imperfectly expanded supersonic jet discharging into a subsonic coflowing stream for which enhanced mixing due to the onset of screech (feedback instability) is known to occur. The purpose of these simulations is to see if enhanced mixing and feedback instability would, indeed, take place. Some evidence of feedback oscillations is found in the simulated flow.

Tam, C. K. W.↗

DOE OSTI2019

Numerical Simulation and Experimental Validation of Nondendritic Structure Formation in Magnesium Alloy Under Oscillation and Ultrasonic Vibration

In this study, the formation of nondendritic structures in the primary phase of magnesium alloy solidified under oscillation and ultrasonic vibration was investigated by numerical simulation and experimentally. The growth and motion of a dendrite during solidification was simulated by a combination of the lattice Boltzmann method and the phase-field method. The simulation and experimental results indicated that higher oscillation amplitudes and acoustic streaming made the microstructures change from dendritic to nondendritic in the α-Mg primary phases. A sufficient shear stress and an appropriate flow time in the barrel should be satisfied when a given inclined angle is selected. The effects of the flow and the thermal field on the nucleation behavior, the constitutional undercooling, and the growth morphology are also discussed. It was found that a high shear rate and high turbulence can help homogenize the temperature and the concentration fields and collide and rotate the α-Mg primary phases.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Validation of the High-Resolution Salish Sea Tidal Hydrodynamic Model

In this study, a tidal hydrodynamic model was developed and validated to simulate tidal currents in Puget Sound, Washington, to support tidal energy resource characterization using the unstructured-grid, Finite Volume Community Ocean Model (FVCOM). The Salish Sea tidal hydrodynamic model was driven by tides along two open boundaries at the entrance of the Strait of Juan de Fuca and north end of Georgia Strait, and river flows from 19 major rivers in the Salish Sea. To simulate the tidal current in Puget Sound, a high-resolution model grid is required to accurately represent the complex coastlines and bathymetry. The spatial resolution of the model grid varies from ~10 m near river boundaries and ~30 m in small tidal channels and estuaries to near 1000 m inside Georgia Strait and at the open boundaries. Model validation was carried out by comparing simulated and observed water levels at 12 tidal stations and currents at 135 Acoustic Doppler Current Profiler stations in the model domain. A set of model performance metrics, including root mean square error, scatter index, bias, and linear correlation coefficient, were used to quantify the model skills in simulating the tidal hydrodynamics in Puget Sound. Error statistics showed an overall good agreement between simulated and observed tidal elevations and currents, which demonstrated that the Puget Sound tidal model can be used to accurately characterize the tidal stream energy resource in Puget Sound.

16 TIDAL AND WAVE POWER↗

NASA NTRS2019

Acoustic Phased Array Quantification of Quiet Technology Demonstrator 3 Advanced Inlet Liner Noise Component

Acoustic phased array flyover noise measurements were acquired as part of the Boeing 737 MAX-7 NASA Advanced Inlet Liner segment of the Quiet Technology Demonstrator 3 (QTD3) flight test program. This paper reports on the processes used for separating and quantifying the engine inlet, exhaust and airframe noise source components and provides sample phased array-based comparisons of the component noise source levels associated with the inlet liner treatment configurations. Full scale flyover noise testing of NASA advanced inlet liners was conducted as part of the Quiet Technology Demonstrator 3 flight test program in July and August of 2018. Details on the inlet designs and testing are provided in the companion paper of Reference 1. The present paper provides supplemental details relating to the acoustic phased array portion of the analyses provided in Ref. 1. In brief, the test article was a Boeing 737MAX-7 aircraft with a modified right hand (starboard side) engine inlet, which consisted of either a production inlet liner, a NASA designed inlet liner or a simulated hard wall configuration (accomplished by applying speed tape over the inlet acoustic treatment areas). In all three configurations, the engine forward fan case acoustic panel was replaced with a unperforated (hardwall) panel. No other modifications to any other acoustic treatment areas were made. The left hand (port side) engine was a production engine and was flown at idle thrust for all measurements in order to isolate the effects of the inlet liners to the right hand engine. As described in Ref. 1, the NASA inlet treatment consists of laterally cut slots (cut perpendicular to the flow direction) which are designed to reduce excrescence drag while maintaining or exceeding the liner acoustic noise reduction capabilities. The NASA inlet liner consists of a Multi-Degree of Freedom (MDOF) design with two breathable septum layers inserted into each honeycomb cell [1]. The aircraft noise measurements were acquired for both takeoff (flaps 1 setting, gear up) and approach (flaps 30 gear up and gear down) configurations. The inlet and flight test configurations are summarized in Table 1. Table 1: Inlet Treatment and Flight Configurations Inlet Forward Fan Case Aircraft Production Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down NASA Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down Hardwall Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down III.Test Description and Hardware The flight testing was conducted at the Grant County airport in Moses Lake, WA, between 27 July and 6 August 2018. The noise measurement instrumentation included 8 flush dish microphones arranged in a noise certification configuration as well as an 840 microphone phased array. The flush dish microphones were used to quantify the levels and differences in levels between the various inlet treatments. The phased array was used to separate and quantify the narrowband (tonal) and broadband noise component levels from the engine inlet/exhaust and from the airframe. Phased array extraction of the broadband component was critical to this study because it allowed for the separation of the inlet component from the total airplane level noise even when it was significantly below the total level. Figure 1 provides an overview of the phased array microphone layout as well as a detailed image of an individual phased array microphone mounted in a plate holder (the microphone sensor is the dot in the center of the plate). The ground plane ensemble array microphones (referred to as “ensemble array” in this paper) were mounted in plates with “flower petal” edges designed to minimize edge scattering effects. Fig. 1 Flyover test microphone layout. The phased array configuration was the result of a progressive development of concepts originally implemented in Ref. 2 and refined over the following years, consisting namely of multiple multi-arm logarithmic spiral subarrays designed to cover overlapping frequency ranges and optimized for various aircraft emission angles. For the present case, the signals from all 840 microphones were acquired on a single system. The 840 microphones were parsed into 11 primary subarray sets spanning from smallest to largest aperture size and labeled accordingly as “a, b, …, k”, where “a” corresponds to the smallest fielded subarray and “k” corresponds to the largest aperture subarray. The apertures ranged from approximately 10 ft to 427 ft in size (in the flight direction) with the subarrays consisting of between 215 and 312 microphones. Figure 2 shows three such subarrays, k, h and a. As done in Ref. 2, microphones were shared between subarrays in order to reduce total channel count. Fig. 2 Sample subarray sizes (20° from overhead – refer to Figure 3a discussion). In addition to the above, each of the 11 primary subarray sets consisted of four subarrays optimized to provide near equivalent array spatial resolution in both the flight and lateral directions within 30 degrees of overhead (i.e., airplane directly above the center of the array), namely, at angles  of 0, ±10, ±20 and ±30 degrees relative to overhead where angle  is defined as shown in Figure 3a. This allowed for optimized aircraft noise measurements from 60 to 120 degree emission angle.6 An example of this pletharray design is shown in Figure 3b for the k subarray. When the aircraft is at overhead, the microphones indicated by the blue markers are used for beamforming. When the aircraft is at angles ±10 degrees from overhead, both the blue and red colored microphones are used, and so on for the ±20 and ±30 degree aircraft locations. See Ref. 3 for extensive details on pletharray design for aeroacoustic phased array testing. 6 In the discussions that follow, emission angle values are used. These are the angles at the time sound is emitted relative to the engine axis and are calculated based on flight path angle, body aircraft body angle with respect to the relative wind direction, and engine axis angle relative to aircraft body angle.

Brusniak, Leon↗

DOE OSTI2018

A hybrid adaptive low-Mach number/compressible method: Euler equations

Flows in which the primary features of interest do not rely on high-frequency acoustic effects, but in which long-wavelength acoustics play a nontrivial role, present a computational challenge. Integrating the entire domain with low-Mach-number methods would remove all acoustic wave propagation, while integrating the entire domain with the fully compressible equations can in some cases be prohibitively expensive due to the CFL time step constraint. For example, simulation of thermoacoustic instabilities might require fine resolution of the fluid/chemistry interaction but not require fine resolution of acoustic effects, yet one does not want to neglect the long-wavelength wave propagation and its interaction with the larger domain. The present paper introduces a new multi-level hybrid algorithm to address these types of phenomena. In this new approach, the fully compressible Euler equations are solved on the entire domain, potentially with local refinement, while their low-Mach-number counterparts are solved on subregions of the domain with higher spatial resolution. The finest of the compressible levels communicates inhomogeneous divergence constraints to the coarsest of the low-Mach-number levels, allowing the low-Mach-number levels to retain the long-wavelength acoustics. Here, the performance of the hybrid method is shown for a series of test cases, including results from a simulation of the aeroacoustic propagation generated from a Kelvin-Helmholtz instability in low-Mach-number mixing layers. It is demonstrated that compared to a purely compressible approach, the hybrid method allows time-steps two orders of magnitude larger at the finest level, leading to an overall reduction of the computational time by a factor of 8.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1995

Simulation study of the role of ion kinetics in low-frequency wave train evolution

The evolution of uniform, parallel propagating, low-frequency (less than or approx. ion cyclotron) wave trains is followed with a one-dimensional hybrid numerical code with fluid electrons and particle ions. We show that moderate amplitude (delta B/B less than 1/2) wave trains give instabilities and saturated states which differ completely from pure fluid evolution. This is most clearly seen when beta greater than 1 and instability exists for wavenumbers both below and above the wavenumber of an inital, left-handed wave train or pump wave. For corresponding parameters a fluid theory gives only a narrow range of instability above the pump wavenumber where decay and beat instabilities can occur. in simulations wave energy inverse cascades to smaller wavenumbers and into a greater number of forward than backward going waves. In fluids energy by decay goes mostly to backward ones of smaller wavenumber, and energy by beat goes mostly to forward ones of larger wavenumber. Neither fluid instability explains simulation results. The instability is saturated by thermalizing ions and sometimes exciting small wavenumber electrostatic or acoustic modes. In contrast, saturation in fluids first occurs by generating the harmonics of the growing linear modes. Harmonic generation is mostly absent in simulations. Simulations are carried out to long times and mostly reach a limit beyond which no further significant evolution can occur. Application to Alfvenic fluctuations in the solar wind is discussed.

Vasquez, Bernard J.↗

NASA NTRS2016

First-Principles Modeling of IAT-Driven Anomalous Resistivity in Hollow Cathode Discharges II: Numerical Simulations and Comparisons with Experiments

We present a model that quantifies the magnitude of the ion-acoustic turbulence (IAT) in the plume of hollow cathodes and its effect on the resistivity and ion heating. The model takes the form of a partial differential equation (PDE) that can be solved concurrently with the equations of motion for a partially ionized plasma already included in our numerical code for the simulation of the plasma discharge in hollow cathodes, OrCa2D. We also determine that self-induced magnetic fields are not negligible in hollow cathodes operating at large discharge currents and implement in our code Ampere’s law and modifications to Ohm’s law that account for this effect. Numerical simulations that employed these models show large improvements in our agreement with experimental measurements with respect to a previous model, which assumed complete saturation of the IAT and did not account for the growth stage of the waves. In particular, the model is able to accurately predict the location and magnitude of the maximum resitivity to the electron current along the cathode centerline.

Lopez Ortega, Alejandro↗

NASA NTRS1976

Aero-acoustic performance characteristics of duct burning turbofan exhaust nozzles

A recent experimental investigation has identified the aero/acoustic characteristics of exhaust nozzles for duct heating turbofan engines over a range of simulated flow conditions. Jet noise and performance levels are summarized for a series of coannular nozzles representing both acoustically suppressed and unsuppressed designs operating in a static environment. The basic coannular nozzles were found to provide inherent noise suppression. Multi-element suppressor nozzles provided additional noise suppression, but with appreciable thrust loss. The impact of these results on the advanced supersonic transport studies is also presented, indicating potentially large reductions in take-off gross weight or community noise footprints.

Kozlowski, H.↗

DOE OSTI2014

Effects of small variations of speed of sound in optoacoustic tomographic imaging

Purpose: Speed of sound difference in the imaged object and surrounding coupling medium may reduce the resolution and overall quality of optoacoustic tomographic reconstructions obtained by assuming a uniform acoustic medium. In this work, the authors investigate the effects of acoustic heterogeneities and discuss potential benefits of accounting for those during the reconstruction procedure. Methods: The time shift of optoacoustic signals in an acoustically heterogeneous medium is studied theoretically by comparing different continuous and discrete wave propagation models. A modification of filtered back-projection reconstruction is subsequently implemented by considering a straight acoustic rays model for ultrasound propagation. The results obtained with this reconstruction procedure are compared numerically and experimentally to those obtained assuming a heuristically fitted uniform speed of sound in both full-view and limited-view optoacoustic tomography scenarios. Results: The theoretical analysis showcases that the errors in the time-of-flight of the signals predicted by considering the straight acoustic rays model tend to be generally small. When using this model for reconstructing simulated data, the resulting images accurately represent the theoretical ones. On the other hand, significant deviations in the location of the absorbing structures are found when using a uniform speed of sound assumption. The experimental results obtained with tissue-mimicking phantoms and a mouse postmortem are found to be consistent with the numerical simulations. Conclusions: Accurate analysis of effects of small speed of sound variations demonstrates that accounting for differences in the speed of sound allows improving optoacoustic reconstruction results in realistic imaging scenarios involving acoustic heterogeneities in tissues and surrounding media.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1994

Acoustic properties and durability of liner materials at non-standard atmospheric conditions

This report documents the results of an experimental study on how acoustic properties of certain absorbing liner materials are affected by nonstandard atmospheric conditions. This study was motivated by the need to assess risks associated with incorporating acoustic testing capability in wind tunnels with semicryogenic high Reynolds number aerodynamic and/or low pressure capabilities. The study consisted of three phases: 1) measurement of acoustic properties of selected liner materials at subatmospheric pressure conditions, 2) periodic cold soak and high pressure exposure of liner materials for 250 cycles, and 3) determination of the effect of periodic cold soak on the acoustic properties of the liner materials at subatmospheric conditions and the effect on mechanical resiliency. The selected liner materials were Pyrell foam, Fiberglass, and Kevlar. A vacuum facility was used to create the subatmospheric environment in which an impedance tube was placed to measure acoustic properties of the test materials. An automated cryogenic cooling system was used to simulate periodic cold soak and high pressure exposure. It was found that lower ambient pressure reduced the absorption effectiveness of the liner materials to varying degrees. Also no significant change in the acoustic properties occurred after the periodic cold soak. Furthermore, mechanical resiliency tests indicated no noticeable change.

Ahuja, K. K.↗