Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kappa”

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 199 records · Page 11

Computation of Separated and Unsteady Flows with One- and Two-Equation Turbulence Models

The ability of one- and two-equation turbulence models to predict unsteady separated flows over airfoils is evaluated. An implicit, factorized, upwind-biased numerical scheme is used for the integration of the compressible, Reynolds averaged Navier-Stokes equations. The turbulent eddy viscosity is obtained from the computed mean flowfield by integration of the turbulent field equations. The two-equation turbulence models are discretized in space with an upwind-biased, second order accurate total variation diminishing scheme. One and two-equation turbulence models are first tested for a separated airfoil flow at fixed angle of incidence. The same models are then applied to compute the unsteady flowfields about airfoils undergoing oscillatory motion at low subsonic Mach numbers. Experimental cases where the flow has been tripped at the leading edge and where natural transition was allowed to occur naturally are considered. The more recently developed field-equation turbulence models capture the physics of unsteady separated flow significantly better than the standard kappa-epsilon and kappa-omega models. However, certain differences in the hysteresis effects are obtained. For an untripped high-Reynolds-number flow, it was found necessary to take into account the leading edge transitional flow region in order to capture the correct physical mechanism that leads to dynamic stall.

Ekaterinaris, John A.↗

Application of a Third Order Upwind Scheme to Viscous Flow over Clean and Iced Wings

A 3-D compressible Navier-Stokes solver has been developed and applied to 3-D viscous flow over clean and iced wings. This method uses a third order accurate finite volume scheme with flux difference splitting to model the inviscid fluxes, and second order accurate symmetric differences to model the viscous terms. The effects of turbulence are modeled using a Kappa-epsilon model. In the vicinity of the sold walls the kappa and epsilon values are modeled using Gorski's algebraic model. Sampling results are presented for surface pressure distributions, for untapered swept clean and iced wings made of NACA 0012 airfoil sections. The leading edge of these sections is modified using a simulated ice shape. Comparisons with experimental data are given.

Bangalore, A.↗

Nonguiding Center Motion and Substorm Effects in the Magnetotail

Thick and thin models of the middle magnetotail were developed using a consistent orbit tracing technique. It was found that currents carried near the equator by groups of ions with anisotropic distribution functions are not well approximated by the guiding center expressions. The guiding center equations fail primarily because the calculated pressure tensor is not magnetic field aligned. The pressure tensor becomes field aligned as one moves away from the equator, but here there is a small region in which the guiding center equations remain inadequate because the two perpendicular components of the pressure tensor are unequal. The significance of nonguiding center motion to substorm processes then was examined. One mechanism that may disrupt a thin cross-tail current sheet involves field changes that cause ions to begin following chaotic orbits. The lowest-altitude chaotic region, characterized by an adiabaticity parameter kappa approx. equal to 0.8, is especially important. The average cross-tail particle drift is slow, and we were unable to generate a thin current sheet using such ions. Therefore, any process that tends to create a thin current sheet in a region with kappa approaching 0.8 may cause the cross-tail current to get so low that it becomes insufficient to support the lobes. A different limit may be important in resonant orbit regions of a thin current sheet because particles reach a maximum cross-tail drift velocity. If the number of ions per unit length decreases as the tail is stretched, this part of the plasma sheet also may become unable to carry the cross-tail current needed to support the lobes. Thin sheets are needed for both resonant and chaotic orbit mechanisms because the distribution function must be highly structured. A description of current continuity is included to show how field aligned currents can evolve during the transition from a two-dimensional (2-D) to a 3-D configuration.

Kaufmann, Richard L.↗

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗

Detailed Comparison of Observed Dose-Time Profile of October 19-20, 1989 SPE on Mir with Model Calculations

The dose rate dynamics of the October 19-20,1989 solar energetic particle (SPE) event as observed by the Liulin instrument onboard the Mir orbital station was analyzed in light of new calculations of the geomagnetic cutoff and improved estimates of the less than 100 MeV energy spectra from the GOES satellite instrument. The new calculations were performed using the as-flown Mir orbital trajectory and includes time variations of the cutoff rigidity due to changes in the kappa (sub p) index. Although the agreement of total event integrated calculated dose to the measured dose is good, it results from some measured dose-time profile been higher and some lower than model calculations. They point to the need to include the diurnal variation of the geomagnetic cutoff and modifications of the cutoffs to variations in kappa (sub p) in model calculations. Understanding of such events in light of the upcoming construction of the International Space Station during the period of maximum solar activity needs to be vigorously pursued.

Badhwar, Gautam D.↗

First Year Wilkinson Microwave Anisotropy Probe(WMAP)Observations: The Angular Power Spectrum

We present the angular power spectrum derived from the first-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We study a variety of power spectrum estimation methods and data combinations and demonstrate that the results are robust. The data are modestly contaminated by diffuse Galactic foreground emission, but we show that a simple Galactic template model is sufficient to remove the signal. Point sources produce a modest contamination in the low frequency data. After masking approximately 700 known bright sources from the maps, we estimate residual sources contribute approximately 3500 mu sq Kappa at 41 GHz, and approximately 130 mu sq Kappa at 94 GHz, to the power spectrum [iota(iota + 1)C(sub iota)/2pi] at iota = 1000. Systematic errors are negligible compared to the (modest) level of foreground emission. Our best estimate of the power spectrum is derived from 28 cross-power spectra of statistically independent channels. The final spectrum is essentially independent of the noise properties of an individual radiometer. The resulting spectrum provides a definitive measurement of the CMB power spectrum, with uncertainties limited by cosmic variance, up to iota approximately 350. The spectrum clearly exhibits a first acoustic peak at iota = 220 and a second acoustic peak at iota approximately 540, and it provides strong support for adiabatic initial conditions. Researchers have analyzed the CT(sup Epsilon) power spectrum, and present evidence for a relatively high optical depth, and an early period of cosmic reionization. Among other things, this implies that the temperature power spectrum has been suppressed by approximately 30% on degree angular scales, due to secondary scattering.

Hinshaw, G.↗

Transition and Turbulence Modeling for Blunt-Body Wake Flows

This study attempts t o improve the modeling and computational prediction of high- speed transitional wake flows. The recently developed kappa - zeta (Enstrophy) turbulence model is coupled with a newly developed transition prediction method and implemented in an implicit flow solver well-suited to hypersonic flows. In this model, transition onset is determined as part of the solution. Results obtained using the new model for a 70- deg blunted cone/sting geometry demonstrate better agreement with experimental heat- transfer measurements when compared to laminar calculations as well as solutions using the kappa - omega model. Results are also presented for the situation where transition onset is preselected. It is shown that, in this case, results are quite sensitive to location of the transition point.

Nance, Robert P.↗

Electron Velocity Distribution Function in Magnetic Clouds in the Solar Wind

We present a study of the kinetic properties of the electron velocity distribution functions within magnetic clouds, since they are the dominant thermal component. The study is based on high time resolution data from the GSFC WIND/SWE electron spectrometer and the Berkeley 3DP electron plasma instruments. Recent studies on magnetic clouds have shown observational evidence of anti-correlation between the total electron density and electron temperature, which suggest a polytrope law P(sub e) = alpha(Nu(sub e) (sup gamma)) for electrons with the constant gamma approximates 0.5 < 1. This anti-correlation and small polytropic gamma-values is interpreted in the context of the presence of highly non-Maxwellian electron distributions (i.e. non-thermal) within magnetic clouds. These works suggested that the non-thermal electrons can contribute as much as 50% of the total electron pressure within magnetic clouds. We have revisited some of the magnetic cloud events previously studied and attempted to quantify the nature of the non-thermal electrons by modeling the electron velocity distribution function using a kappa distribution function to characterize the kinetic non-thermal effects. If non-thermal tail effects are the source for the anti-correlation between the moment electron temperature and density and if the kappa distribution is a reasonable representative model of non-thermal effects, then the electron velocity distribution within magnetic clouds should show indication for small K-values when gamma < 1.

Nieves-Chinchil, Teresa↗

Space Propulsion Design and Analysis

This software provides an improved methodology for predicting launcher base pressure and heat loads for RSRM (Reusable Solid Rocket Motor) launchers by accounting for complex anisotropic stress/strains and variable turbulent Prandtl and Schmidt numbers. A "building block" approach to turbulence model development, and validation has been applied for improved missile/launcher base region analysis. Modifications to existing kappa - epsilon turbulence models and application of scalar variance models are incorporated into a RANS-based method for aeropropulsive flow modeling, directly related to base flow methodology. (RANS stands for Reynolds-averaged Navier-Stokes.) The models are applied in a RANS solver framework and can improve analysis of other complex flow fields. The enhanced models provide a more accurate predictive capability for improving the design and analysis of RSRM launcher configuration. The kappa - epsilon model enhancements have been shown to improve the capability for predicting turbulence effects in base blow environments. The scalar variance models have been assessed over a wide range of flow configurations to improve prediction of turbulent scalar mixing.

Sinha, Neeraj↗

Crystal Lattice Controlled SiGe Thermoelectric Materials with High Figure of Merit

Direct energy conversion between thermal and electrical energy, based on thermoelectric (TE) effect, has the potential to recover waste heat and convert it to provide clean electric power. The energy conversion efficiency is related to the thermoelectric figure of merit ZT expressed as ZT=S(exp 2)(sigma)T/Kappa, T is temperature, S is the Seebeck coefficient, sigma is conductance and Kappa is thermal conductivity. For a lower thermal conductivity Kappa and high power factor (S(exp 2)(sigma)), our current strategy is the development of rhombohedrally strained single crystalline SiGe materials that are highly [111]-oriented twinned. The development of a SiGe "twin lattice structure (TLS)" plays a key role in phonon scattering. The TLS increases the electrical conductivity and decreases thermal conductivity due to phonon scattering at stacking faults generated from the 60 X rotated primary twin structure. To develop high performance materials, the substrate temperature, chamber working pressure, and DC sputtering power are controlled for the aligned growth production of SiGe layer and TLS on a c-plane sapphire. Additionally, a new elevated temperature thermoelectric characterization system, that measures the thermal diffusivity and Seebeck effect nondestructively, was developed. The material properties were characterized at various temperatures and optimized process conditions were experimentally determined. The present paper encompasses the technical discussions toward the development of thermoelectric materials and the measurement techniques.

Kim, Hyun-Jung↗

Using High Spatial Resolution Satellite Imagery to Map Forest Burn Severity Across Spatial Scales in a Pine Barrens Ecosystem

As a primary disturbance agent, fire significantly influences local processes and services of forest ecosystems. Although a variety of remote sensing based approaches have been developed and applied to Landsat mission imagery to infer burn severity at 30 m spatial resolution, forest burn severity have still been seldom assessed at fine spatial scales (less than or equal to 5 m) from very-high-resolution (VHR) data. We assessed a 432 ha forest fire that occurred in April 2012 on Long Island, New York, within the Pine Barrens region, a unique but imperiled fire-dependent ecosystem in the northeastern United States. The mapping of forest burn severity was explored here at fine spatial scales, for the first time using remotely sensed spectral indices and a set of Multiple Endmember Spectral Mixture Analysis (MESMA) fraction images from bi-temporal - pre- and post-fire event - WorldView-2 (WV-2) imagery at 2 m spatial resolution. We first evaluated our approach using 1 m by 1 m validation points at the sub-crown scale per severity class (i.e. unburned, low, moderate, and high severity) from the post-fire 0.10 m color aerial ortho-photos; then, we validated the burn severity mapping of geo-referenced dominant tree crowns (crown scale) and 15 m by 15 m fixed-area plots (inter-crown scale) with the post-fire 0.10 m aerial ortho-photos and measured crown information of twenty forest inventory plots. Our approach can accurately assess forest burn severity at the sub-crown (overall accuracy is 84% with a Kappa value of 0.77), crown (overall accuracy is 82% with a Kappa value of 0.76), and inter-crown scales (89% of the variation in estimated burn severity ratings (i.e. Geo-Composite Burn Index (CBI)). This work highlights that forest burn severity mapping from VHR data can capture heterogeneous fire patterns at fine spatial scales over the large spatial extents. This is important since most ecological processes associated with fire effects vary at the less than 30 m scale and VHR approaches could significantly advance our ability to characterize fire effects on forest ecosystems.

Meng, Ran↗

Virus Infection of Phytoplankton Increases Average Molar Mass and Reduces Hygroscopicity of Aerosolized Organic Matter

Viral infection of phytoplankton is a pervasive mechanism of cell death and bloom termination, which leads to the production of dissolved and colloidal organic matter that can be aerosolized into the atmosphere. Earth-observing satellites can track the growth and death of phytoplankton blooms on weekly time scales but the impact of viral infection on the cloud forming potential of associated aerosols is largely unknown. Here, we determine the influence of viral-derived organic matter, purified viruses, and marine hydrogels on the cloud condensation nuclei activity of their aerosolized solutions, compared to organic exudates from healthy phytoplankton. Dissolved organic material derived from exponentially growing and infected cells of well-characterized eukaryotic phytoplankton host-virus systems, including viruses from diatoms, coccolithophores and chlorophytes, was concentrated, desalted, and nebulized to form aerosol particles composed of primarily of organic matter. Aerosols from infected phytoplankton cultures resulted in an increase in critical activation diameter and average molar mass in three out of five combinations evaluated, along with a decrease in organic kappa (hygroscopicity) compared to healthy cultures and seawater controls. The infected samples also displayed evidence of increased surface tension depression at realistic cloud water vapor supersaturations. Amending the samples with xanthan gum to simulate marine hydrogels increased variability in organic kappa and surface tension in aerosols with high organic to salt ratios. Our findings suggest that the pulses of increased dissolved organic matter associated with viral infection in surface waters may increase the molar mass of dissolved organic compounds relative to surface waters occupied by healthy phytoplankton or low phytoplankton biomass.

Atmospheric science↗

Diffuse Auroral Intensities Produced By Whistler Mode and Electron Cyclotron Harmonic Waves

Whistler mode waves and ECH waves as observed by the THEMIS-D and THEMIS-E satellites have been analyzed. It is observed that ECH waves are very weak for pitch-angle diffusion whereas whistler mode waves are more efficient. Bounce-averaged pitch-angle diffusion rates at the edge of atmospheric loss-cone have been calculated for both waves. Further, these are used to obtain electron precipitation flux inside the loss-cone. The electron flux at the edge of the loss-cone is represented by the kappa distribution. Numerical calculation of precipitation flux is performed by varying the fitting parameters (which determine the shape of electron precipitation flux) appearing in the kappa distribution. Precipitation flux is used to obtain volume excitation rates and height-integrated volume excitation rates for the seven excitation states. Calculations have been performed using two standard atmosphere models corresponding to mean and high solar and geomagnetic activities. It has been observed from studies that the excitation rates produced by the ECH wave are several orders of magnitude smaller than the rates produced by whistler mode waves. We show that the scattering efficiency of whistler Mode is higher than that of ECH waves. The oblique whistler mode waves allow efficient scattering of electrons in diffuse aurora. Outcomes are discussed.

Arvind K. Tripathi↗

Estimating β eff from static, late-time Feynman neutron noise measurements

The effective delayed neutron fraction, β eff , is a kinematic parameter describing the contribution of delayed neutrons to the effective multiplication factor, $\kappa$ eff , and is of great interest to the nuclear criticality safety community and beyond. In fact, $\kappa$ eff presently must be inferred (as opposed to directly measured) since variables like β eff must be obtained from simulation, reference tables, or a dynamic or pulsed measurement where the assembly of fissile material is perturbed and not in steady state. Several recent works have proposed neutron noise techniques to estimate β eff from static measurements (such as benchmark experiments). This work experimentally investigates an approach based on looking at late-time (up to thousands of seconds) correlated neutrons in Feynman histograms. This work utilizes measurements of weapons grade plutonium and highly enriched uranium to construct late-time Feynman histograms and investigate a proposed methodology of estimating β eff from a fit of the distribution. Data is shown from both He-3 detectors and an organic scintillator array. This work is preliminary in nature, demonstrating some experimental feasibility and proposing improvements for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification of Upcoming MCNP Features For Estimating Nuclear Data Sensitivities in Fixed Source Simulations [Abstract]

Predictive simulation codes, like the Monte Carlo N-Particle (MCNP) transport code, are used throughout the nuclear community. These simulations are based on nuclear data. Maximizing the accuracy and precision of nuclear data maximizes the accuracy and precision of the overall simulation. This is imperative to applications that rely on simulations. For example, improving nuclear data for special nuclear material improves simulation accuracy in stockpile stewardship applications, which results in larger safety margins and decreased operational costs. The improvement and validation of nuclear data is completed through integral benchmark experiments. Past benchmarks have primarily been limited to focus on the effective multiplication factor ($\kappa$ eff ); broadening the purview of benchmarks beyond $\kappa$ eff -dependent nuclear data addresses nuclear data deficiencies. Different response types depend on different areas of nuclear data. This dependence is quantified as nuclear data sensitivity: the change in response due to perturbation of a contributing parameter. The larger the nuclear data sensitivity of a response, the more the experiment is influenced by the uncertainties of the nuclear data. The optimization of nuclear data sensitivities in future benchmarks would result in more detailed validation of lesser studied areas of nuclear data. Currently, direct sensitivity capabilities are not easily found for all experiment types and parameters. An MCNP tool to directly estimate the cross section sensitivities of tallied values is under development. Additionally, updates have been made to the perturbation feature of MCNP, which can be used in a less direct approach to estimating sensitivities. This work verifies these features to estimate nuclear data sensitivities in fixed source simulations of a 4.5-kg sphere of alpha- phase weapons-grade plutonium surrounded by differing amounts of copper and polyethylene. Integrated estimates made using MCNP’s tools were found to statistically agree with integrated estimates made from manual perturbation of nuclear data proving the validity of the MCNP tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

aosccnsmpskappa.c1

This product calculates the hygroscopicity parameter kappa from aerosol particle concentrations and size distributions

54 ENVIRONMENTAL SCIENCES↗

aosccnuhsaskappa.c1

Hygroscopicity parameter kappa derived from CCN and UHSAS data sources

54 ENVIRONMENTAL SCIENCES↗

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗