Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source term”

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 181 records · Page 10

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.

Maccormack, Robert W.↗

Modeling of Turbulence Generated Noise in Jets

A numerically calculated Green's function is used to predict jet noise spectrum and its far-field directivity. A linearized form of Lilley's equation governs the non-causal Green s function of interest, with the non-linear terms on the right hand side identified as the source. In this paper, contributions from the so-called self- and shear-noise source terms will be discussed. A Reynolds-averaged Navier-Stokes solution yields the required mean flow as well as time- and length scales of a noise-generating turbulent eddy. A non-compact source, with exponential temporal and spatial functions, is used to describe the turbulence velocity correlation tensors. It is shown that while an exact non-causal Green's function accurately predicts the observed shift in the location of the spectrum peak with angle as well as the angularity of sound at moderate Mach numbers, at high subsonic and supersonic acoustic Mach numbers the polar directivity of radiated sound is not entirely captured by this Green's function. Results presented for Mach 0.5 and 0.9 isothermal jets, as well as a Mach 0.8 hot jet conclude that near the peak radiation angle a different source/Green's function convolution integral may be required in order to capture the peak observed directivity of jet noise.

Khavaran, Abbas↗

Computation of three-dimensional viscous flows using a space-marching method

A space-marching method, developed to compute three-dimensional flows for internal geometries, has been utilized to predict viscous flows through a curved duct and over a swept wing. The Navier-Stokes equations have been posed as an initial value problem by neglecting the streamwise viscous diffusion terms and by treating the pressure gradient as a known source term. The resulting equations have been solved by a non-iterative (single pass) algorithm at each streamwise step. The results are compared with earlier computations (based on iterative methods) and the experimental data. The agreement between the present predictions, the experimental data, and the earlier predictions is good for the cases computed. The computation time is only a fraction of the iterative methods.

Murthy, K. N. S.↗

Determination of near and far field acoustics for advanced propeller configurations

A method has been studied for predicting the acoustic field of the SR-3 transonic propfan using flow data generated by two versions of the NASPROP-E computer code. Since the flow fields calculated by the solvers include the shock-wave system of the propeller, the nonlinear quadrupole noise source term is included along with the monopole and dipole noise sources in the calculation of the acoustic near field. Acoustic time histories in the near field are determined by transforming the azimuthal coordinate in the rotating, blade-fixed coordinate system to the time coordinate in a nonrotating coordinate system. Fourier analysis of the pressure time histories is used to obtain the frequency spectra of the near-field noise.

Korkan, K. D.↗

Plasma modification through boron particulate injection in the full tungsten environment of WEST

Recent experiments have confirmed the compatibility of extended boron particulate injections with high performance plasma discharges in the full tungsten (W) environment of WEST. Utilizing an impurity powder dropper (IPD) equipped with boron (B) powders a series of extended experimental programs providing controlled injections have quantified plasma response to varying levels of injection rate and total injection quantity. Calibration of injection quantities confirmed through post-situ testing of the IPD and cross-correlated with both high-speed camera illumination and spectroscopic measurement have allowed for the first time a fine scale determination of the effects of powder introduction on plasma performance. Plasma enhancement, consistent with turbulence reduction through profile modification, has been observed with sustained increases in the stored energy (WMHD), by 18%, electron temperature (T e ) by 35%, and neutron rate (N n ) by up to 200%, all of which scale positively with increasing powder injection rates. These injections have also resulted in both prompt and extended reductions in native impurity content, decreases in post injection radiated power, and strong decreases in divertor deuterium signatures signifying a reduction in recycling suggesting enhanced boron layer formation which provides a reduction of source terms and leads to enhanced gettering of main ion and impurity sources.

Lunsford, R. [Princeton Plasma Physics Laboratory ↗

Flux-difference split parabolized Navier-Stokes algorithm for non-equilibrium chemically reacting flows

A flux-difference split explicit finite-difference algorithm is presented for solving the parabolized form of the equations governing three-dimensional nonequilibrium chemically reacting flows. The algorithm is based on an explicit noniterative, upwind space-marching scheme developed by Korte, but differs in that the unsteady Riemann problem, rather than the steady Riemann problem, is solved. The algorithm allows either a second or an approximately third-order accurate upwind treatment of the convection terms by employing the unsteady approximate Riemann solver of Roe. The source terms of the species transport equations are treated in either an explicit or implicit manner, and the species diffusion terms are modeled with either a Fickian or a multicomponent model. A validation of the algorithm is performed by comparing computational results with the 2-D Mach 14, 15 degree compression-corner data of Holden. The three-dimensional capability of the algorithm is demonstrated by computing Mach 2.7 flow over a swept wedge scramjet fuel injector, and three-dimensional reacting flow capability is demonstrated by a computing a shock-jet interaction concept for mixing and combustion enhancement.

White, J. A.↗

On the Definition of Surface Potentials for Finite-Difference Operators

For a class of linear constant-coefficient finite-difference operators of the second order, we introduce the concepts similar to those of conventional single- and double-layer potentials for differential operators. The discrete potentials are defined completely independently of any notion related to the approximation of the continuous potentials on the grid. We rather use all approach based on differentiating, and then inverting the differentiation of a function with surface discontinuity of a particular kind, which is the most general way of introducing surface potentials in the theory of distributions. The resulting finite-difference "surface" potentials appear to be solutions of the corresponding continuous potentials. Primarily, this pertains to the possibility of representing a given solution to the homogeneous equation on the domain as a variety of surface potentials, with the density defined on the domain's boundary. At the same time the discrete surface potentials can be interpreted as one specific realization of the generalized potentials of Calderon's type, and consequently, their approximation properties can be studied independently in the framework of the difference potentials method by Ryaben'kii. The motivation for introducing and analyzing the discrete surface potentials was provided by the problems of active shielding and control of sound, in which the aforementioned source terms that drive the potentials are interpreted as the acoustic control sources that cancel out the unwanted noise on a predetermined region of interest.

Tsynkov, S. V.↗

A shock capturing technique for hypersonic, chemically relaxing flows

A fully coupled, shock capturing technique is presented for chemically reacting flows at high Mach numbers. The technique makes use of a total variation diminishing (TVD) dissipation operator which results in sharp, crisp shocks. The eigenvalues and eigenvectors of the fully coupled system, which includes species conversion equations in addition to the gas dynamics equations, are analytically derived for a general reacting gas. Species production terms for a model dissociating gas are introduced and are included in the algorithm. The convective terms are solved using a first-order TVD scheme while the source terms are solved using a fourth-order Runge-Kutta scheme to enhance stability. Results from one-dimensional numerical experiments are shown for a two species and a three species gas.

Eberhardt, S.↗

Measurements and Predictions of the Noise from Three-Stream Jets

An experimental and numerical investigation of the noise produced by high-subsonic and supersonic three-stream jets was conducted. The exhaust system consisted of externally-mixed-convergent nozzles and an external plug. Bypass- and tertiary-to-core area ratios between 1.0 and 2.5, and 0.4 and 1.0, respectively, were studied. Axisymmetric and offset tertiary nozzles were investigated for heated and unheated conditions. For axisymmetric configurations, the addition of the third stream was found to reduce peak- and high-frequency acoustic levels in the peak-jet-noise direction, with greater reductions at the lower bypass-to-core area ratios. For the offset configurations, an offset duct was found to decrease acoustic levels on the thick side of the tertiary nozzle relative to those produced by the simulated two-stream jet with up to 8 dB mid-frequency noise reduction at large angles to the jet inlet axis. Noise reduction in the peak-jet-noise direction was greater for supersonic core speeds than for subsonic core speeds. The addition of a tertiary nozzle insert used to divert the third-stream jet to one side of the nozzle system provided no noise reduction. Noise predictions are presented for selected cases using a method based on an acoustic analogy with mean flow interaction effects accounted for using a Green's function, computed in terms of its coupled azimuthal modes for the offset cases, and a source model previously used for round and rectangular jets. Comparisons of the prediction results with data show that the noise model predicts the observed increase in low-frequency noise with the introduction of a third, axisymmetric stream, but not the high-frequency reduction. For an offset third stream, the model predicts the observed trend of decreased sound levels on the thick side of the jet compared with the thin side, but the predicted azimuthal variations are much less than those seen in the data. Also, the shift of the spectral peak to lower frequencies with increasing polar angle is over-predicted. For an offset third stream with a heated core, it is shown that including the enthalpy-flux source terms in the acoustic analogy model improves predictions compared with those obtained using only the momentum flux.

Three-Stream Jets↗

Measurements and Predictions of the Noise from Three-Stream Jets

An experimental and numerical investigation of the noise produced by high-subsonic and supersonic three-stream jets was conducted. The exhaust system consisted of externally-mixed-convergent nozzles and an external plug. Bypass- and tertiary- to-core area ratios between 1.0 and 2.5, and 0.4 and 1.0, respectively, were studied. Axisymmetric and offset tertiary nozzles were investigated for heated and unheated conditions. For axisymmetric configurations, the addition of the third stream was found to reduce peak- and high-frequency acoustic levels in the peak-jet-noise direction, with greater reductions at the lower bypass-to-core area ratios. For the offset configurations, an offset duct was found to decrease acoustic levels on the thick side of the tertiary nozzle relative to those produced by the simulated two-stream jet with up to 8 dB mid-frequency noise reduction at large angles to the jet inlet axis. Noise reduction in the peak-jet-noise direction was greater for supersonic core speeds than for subsonic core speeds. The addition of a tertiary nozzle insert used to divert the third-stream jet to one side of the nozzle system provided no noise reduction. Noise predictions are presented for selected cases using a method based on an acoustic analogy with mean flow interaction effects accounted for using a Green's function, computed in terms of its coupled azimuthal modes for the offset cases, and a source model previously used for round and rectangular jets. Comparisons of the prediction results with data show that the noise model predicts the observed increase in low-frequency noise with the introduction of a third, axisymmetric stream, but not the high-frequency reduction. For an offset third stream, the model predicts the observed trend of decreased sound levels on the thick side of the jet compared with the thin side, but the predicted azimuthal variations are much less than those seen in the data. Also, the shift of the spectral peak to lower frequencies with increasing polar angle is over-predicted. For an offset third stream with a heated core, it is shown that including the enthalpy-flux source terms in the acoustic analogy model improves predictions compared with those obtained using only the momentum- flux.

Three-Stream Jets↗

Efficient yet accurate solution of the linear transport equation in the presence of internal sources - The exponential-linear-in-depth approximation

The present solutions to the linear transport equation pertain to monoenergetic particles' interaction with a multiple scattering/absorbing layered medium with a general anisotropic internal source term. Attention is given to a novel exponential-linear approximation to the internal source, as a function of scattering depth, which furnishes an at-once efficient and accurate solution to the linear transport equation through its reduction of the spatial mesh size. The great superiority of the proposed method is demonstrated by the numerical results obtained in the illustrative cases of (1) an embedded thermal source and (2) a rapidly varying beam pseudosource.

Kylling, Arve↗

Spectral Characteristics of Wake Vortex Sound During Roll-Up

This report presents an analysis of the sound spectra generated by a trailing aircraft vortex during its rolling-up process. The study demonstrates that a rolling-up vortex could produce low frequency (less than 100 Hz) sound with very high intensity (60 dB above threshold of human hearing) at a distance of 200 ft from the vortex core. The spectrum then drops o rapidly thereafter. A rigorous analytical approach has been adopted in this report to derive the spectrum of vortex sound. First, the sound pressure was solved from an alternative treatment of the Lighthill s acoustic analogy approach [1]. After the application of Green s function for free space, a tensor analysis was applied to permit the removal of the source term singularity of the wave equation in the far field. Consequently, the sound pressure is expressed in terms of the retarded time that indicates the time history and spacial distribution of the sound source. The Fourier transformation is then applied to the sound pressure to compute its spectrum. As a result, the Fourier transformation greatly simplifies the expression of the vortex sound pressure involving the retarded time, so that the numerical computation is applicable with ease for axisymmetric line vortices during the rolling-up process. The vortex model assumes that the vortex circulation is proportional to the time and the core radius is a constant. In addition, the velocity profile is assumed to be self-similar along the aircraft flight path, so that a benchmark vortex velocity profile can be devised to obtain a closed form solution, which is then used to validate the numerical calculations for other more realistic vortex profiles for which no closed form solutions are available. The study suggests that acoustic sensors operating at low frequency band could be profitably deployed for detecting the vortex sound during the rolling-up process.

Booth, Earl R., Jr.↗

Know Your Enemy - Implementation of Bioremediation within a Suspected DNAPL Source Zone Following High-Resolution Site Characterization at Contractors Road Heavy Equipment Area, Kennedy Space Center, Florida

The National Aeronautics and Space Administration (NASA), Kennedy Space Center (KSC), has adopted high-resolution site characterization (HRSC) sampling techniques during baseline sampling prior to implementation of remedies to confirm and refine the conceptual site model (CSM). HRSC sampling was performed at Contractors Road Heavy Equipment Area (CRHE) prior to bioremediation implementation to verify the extent of the trichloroethene (TCE) dense non-aqueous phase liquid (DNAPL) source area (defined as the area with TCE concentrations above 1% solubility) and its daughter product dissolved plume that had been identified during previous HRSC events. The results of HRSC pre-bioremediation implementation sampling suggested that the TCE source area was larger than originally identified during initial site characterization activities, leading to a design refinement to improve electron donor distribution and increase the likelihood of achieving remedial objectives. Approach/Activities: HRSC was conducted from 2009 through 2014 to delineate the vertical and horizontal extent of chlorinated volatile organic compounds (CVOCs) in the groundwater. Approximately 2,340 samples were collected from 363 locations using direct push technology (DPT) groundwater sampling techniques. Samples were collected from up to 14 depth intervals at each location using a 4-foot sampling screen. This HRSC approach identified a narrow (approx. 5 to 30 feet wide), approximately 3,000 square foot TCE DNAPL source area (maximum detected TCE concentration of 160,000 micrograms per liter [micro-g/L] at DPT sampling location DPT0225). Prior to implementation of a bioremediation interim measure, HRSC baseline sampling was conducted using DPT groundwater sampling techniques. Concentrations of TCE were an order of magnitude lower than previous reported (12,000 micro-g/L maximum at DPT sampling location DPT0225) at locations sampled adjacent to previous sampling locations. To further evaluate the variability in concentrations observed additional sampling was conducted in 2016. The results identified higher concentrations than originally detected within the previously defined source area and the presence of source zone concentrations upgradient of the previously defined source area (maximum concentration observed 570,000 micro-g/L). The HRSC baseline sampling data allowed for a revision of the bioremediation design prior to implementation. Bioremediation was implemented within the eastern portion of the source area in November and December 2016 and quarterly performance monitoring was completed in March and June 2017. Reductions in CVOC concentrations from baseline were observed at all performance monitoring wells in the treatment area, and by June 2017, an approximate 95% CVOC mass reduction was observed based on monitoring well sampling results. Results/Lessons Learned: The results of this project suggest that, due to the complexity of DNAPL source zones, HRSC during pre-implementation baseline sampling in the TCE source zone was an essential strategy for verifying the treatment area and depth prior to remedy implementation. If the upgradient source zone mass was not identified prior to bioremediation implementation, the mass would have served as a long-term source for the dissolved plume.

Chrest, Anne↗

Radiocarbon Fingerprinting Black Carbon Source History in the Himalayas

Abstract Black carbon (BC) is considered as an important contributor to the Himalayan glaciers melt in the past few decades. However, the long‐term source apportionment of BC remains unclear. Here we present the first radiocarbon ( 14 C)‐based annual variation of BC source apportionment in an ice core spanning the period of 1959–2012 drilled from the Southeastern Tibetan Plateau, a receptor site of South Asia outflow. We find fossil fuel combustion is a major contribution (73% ± 5%), yet the biomass burning fraction ( ƒ biomass ) has grown from 24% ± 4% to 30% ± 4% since 1990. Intriguingly, we further find the ƒ biomass demonstrating a robust correlation with South Asian wildfires linked to climate oscillations. Thus, for mitigating BC impacts on Himalayan glaciers, South Asia's transition from fossil fuels to clean energy is a more efficient and urgent strategy than reducing residential biomass burning.

Wang, Mo [State Key Laboratory of Tibetan Plateau ↗

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to develop and implement new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. We have just completed two (2) years of Phase I of this research. This annual report provides a brief and up-to-date summary of our activities during the period: September 1, 2000 through August 31, 2001. In the work within the past year, a methodology termed "velocity-scalar filtered density function" (VSFDF) is developed and implemented for large eddy simulation (LES) of turbulent flows. In this methodology the effects of the unresolved subgrid scales (SGS) are taken into account by considering the joint probability density function (PDF) of all of the components of the velocity and scalar vectors. An exact transport equation is derived for the VSFDF in which the effects of the unresolved SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source terms appear in closed form. The remaining unclosed terms in this equation are modeled. A system of stochastic differential equations (SDEs) which yields statistically equivalent results to the modeled VSFDF transport equation is constructed. These SDEs are solved numerically by a Lagrangian Monte Carlo procedure. The consistency of the proposed SDEs and the convergence of the Monte Carlo solution are assessed by comparison with results obtained by an Eulerian LES procedure in which the corresponding transport equations for the first two SGS moments are solved. The unclosed SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source in the Eulerian LES are replaced by corresponding terms from VSFDF equation. The consistency of the results is then analyzed for a case of two dimensional mixing layer.

Drozda, T. G.↗

An Improved Neutron Transport Algorithm for Space Radiation

A low-energy neutron transport algorithm for use in space radiation protection is developed. The algorithm is based upon a multigroup analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. This analysis is accomplished by solving a realistic but simplified neutron transport test problem. The test problem is analyzed by using numerical and analytical procedures to obtain an accurate solution within specified error bounds. Results from the test problem are then used for determining mean values associated with rescattering terms that are associated with a multigroup solution of the straight-ahead Boltzmann equation. The algorithm is then coupled to the Langley HZETRN code through the evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for a water and an aluminum-water shield-target configuration is then compared with LAHET and MCNPX Monte Carlo code calculations for the same shield-target configuration. The algorithm developed showed a great improvement in results over the unmodified HZETRN solution. In addition, a two-directional solution of the evaporation source showed even further improvement of the fluence near the front of the water target where diffusion from the front surface is important.

Heinbockel, John H.↗

Emission of Sound From Turbulence Convected by a Parallel Mean Flow in the Presence of a Confining Duct

An approximate method for calculating the noise generated by a turbulent flow within a semi-infinite duct of arbitrary cross section is developed. It is based on a previously derived high-frequency solution to Lilley's equation, which describes the sound propagation in transversely-sheared mean flow. The source term is simplified by assuming the turbulence to be axisymmetric about the mean flow direction. Numerical results are presented for the special case of a ring source in a circular duct with an axisymmetric mean flow. They show that the internally generated noise is suppressed at sufficiently large upstream angles in a hard walled duct, and that acoustic liners can significantly reduce the sound radiated in both the upstream and downstream regions, depending upon the source location and Mach number of the flow.

Goldstein, Marvin E.↗

Emission of Sound from Turbulence Convected by a Parallel Mean Flow in the Presence of a Confining Duct

An approximate method for calculating the noise generated by a turbulent flow within a semi-infinite duct of arbitrary cross section is developed. It is based on a previously derived high-frequency solution to Lilley's equation, which describes the sound propagation in a transversely-sheared mean flow. The source term is simplified by assuming the turbulence to be axisymmetric about the mean flow direction. Numerical results are presented for the special case of a ring source in a circular duct with an axisymmetric mean flow. They show that the internally generated noise is suppressed at sufficiently large upstream angles in a hard walled duct, and that acoustic liners can significantly reduce the sound radiated in both the upstream and downstream regions, depending upon the source location and Mach number of the flow.

Goldstein, Marvin E.↗