Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Comparisions of modeled N(+), O(+), H(+) and He(+) in the midlatitude ionospherewith mean densities and temperatures from Atmosphere Explorer

In this study, Atmosphere Explorer data and model results for the ion and electron temperature and the density of N(+), O(+), H(+), and He(+) between 120 and 1400 km altitude are compared for two midlatitude ranges (L=2 and L=4), noon and midnight local time, winter and summer, at solar minimum. The data for the heavy atomic ions (O(+) and N(+)) show that their densities are greater at noon than at midnight for a given season and greater in summer than winter for a given local time. There is only a weak latitudinal variation in the density of these ions. The data show that the light ion (H(+) and He(+)) densities are greater at midnight than at noon and are generally greater in winter than summer. There is a strong latitudinal variation of the light ion densities, with the densities decreasing with increasing latitude. The model densities are in good agreement with the AE densities for N(+), O(+), and H(+). Model He(+) densities are lower, by a factor of 2 or more, than the measured densities. Model ion and electron temperatures agree well with the measured temperatures with only a modest increase in plasmapheric heating.

Craven, P. D.↗
NASA NTRS1998

Dynamic Density: An Air Traffic Management Metric

The definition of a metric of air traffic controller workload based on air traffic characteristics is essential to the development of both air traffic management automation and air traffic procedures. Dynamic density is a proposed concept for a metric that includes both traffic density (a count of aircraft in a volume of airspace) and traffic complexity (a measure of the complexity of the air traffic in a volume of airspace). It was hypothesized that a metric that includes terms that capture air traffic complexity will be a better measure of air traffic controller workload than current measures based only on traffic density. A weighted linear dynamic density function was developed and validated operationally. The proposed dynamic density function includes a traffic density term and eight traffic complexity terms. A unit-weighted dynamic density function was able to account for an average of 22% of the variance in observed controller activity not accounted for by traffic density alone. A comparative analysis of unit weights, subjective weights, and regression weights for the terms in the dynamic density equation was conducted. The best predictor of controller activity was the dynamic density equation with regression-weighted complexity terms.

Laudeman, I. V.↗
NASA NTRS2001

Electron Densities Near Io from Galileo Plasma Wave Observations

This paper presents an overview of electron densities obtained near Io from the Galileo plasma wave instrument during the first four flybys of Io. These flybys were Io, which was a downstream wake pass that occurred on December 7, 1995; I24, which was an upstream pass that occurred on October 11, 1999; I25, which was a south polar pass that occurred on November 26, 1999; and I27, which was an upstream pass that occurred on February 22, 2000. Two methods were used to measure the electron density. The first was based on the frequency of upper hybrid resonance emissions, and the second was based on the low-frequency cutoff of electromagnetic radiation at the electron plasma frequency. For three of the flybys, Io, I25, and I27, large density enhancements were observed near the closest approach to Io. The peak electron densities ranged from 2.1 to 6.8 x 10(exp 4) per cubic centimeters. These densities are consistent with previous radio occultation measurements of Io's ionosphere. No density enhancement was observed during the I24 flyby, most likely because the spacecraft trajectory passed too far upstream to penetrate Io's ionosphere. During two of the flybys, I25 and I27, abrupt step-like changes were observed at the outer boundaries of the region of enhanced electron density. Comparisons with magnetic field models and energetic particle measurements show that the abrupt density steps occur as the spacecraft penetrated the boundary of the Io flux tube, with the region of high plasma density on the inside of the flux tube. Most likely the enhanced electron density within the Io flux tube is associated with magnetic field lines that are frozen to Io by the high conductivity of Io's atmosphere, thereby enhancing the escape of plasma along the magnetic field lines that pass through Io's ionosphere.

Gurnett, D. A.↗
NASA NTRS2004

Buoyancy Effects on Flow Structure and Instability of Low-Density Gas Jets

A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.

Pasumarthi, Kasyap Sriramachandra↗
NASA NTRS2009

Optimization of Layer Densities for Spacecraft Multilayered Insulation Systems

Numerous tests of various multilayer insulation systems have indicated that there are optimal densities for these systems. However, the only method of calculating this optimal density was by a complex physics based algorithm developed by McIntosh. In the 1970's much data were collected on the performance of these insulation systems with many different variables analyzed. All formulas generated included number of layers and layer density as geometric variables in solving for the heat flux, none of them was in a differentiable form for a single geometric variable. It was recently discovered that by converting the equations from heat flux to thermal conductivity using Fourier's Law, the equations became functions of layer density, temperatures, and material properties only. The thickness and number of layers of the blanket were merged into a layer density. These equations were then differentiated with respect to layer density. By setting the first derivative equal to zero, and solving for the layer density, the critical layer density was determined. Taking a second derivative showed that the critical layer density is a minimum in the function and thus the optimum density for minimal heat leak, this is confirmed by plotting the original function. This method was checked and validated using test data from the Multipurpose Hydrogen Testbed which was designed using McIntosh's algorithm.

Johnson, W. L.↗
NASA NTRS2017

Small-Scale Density Variations in the Lunar Crust Revealed by GRAIL

Data from the Gravity Recovery and Interior Laboratory (GRAIL) mission have revealed that approximately 98 percent of the power of the gravity signal of the Moon at high spherical harmonic degrees correlates with the topography. The remaining 2 percent of the signal, which cannot be explained by topography, contains information about density variations within the crust. These high-degree Bouguer gravity anomalies are likely caused by small-scale (10's of km) shallow density variations. Here we use gravity inversions to model the small-scale three-dimensional variations in the density of the lunar crust. Inversion results from three non-descript areas yield shallow density variations in the range of 100-200 kg/m3. Three end-member scenarios of variations in porosity, intrusions into the crust, and variations in bulk crustal composition were tested as possible sources of the density variations. We find that the density anomalies can be caused entirely by changes in porosity. Characteristics of density anomalies in the South Pole-Aitken basin also support porosity as a primary source of these variations. Mafic intrusions into the crust could explain many, but not all of the anomalies. Additionally, variations in crustal composition revealed by spectral data could only explain a small fraction of the density anomalies. Nevertheless, all three sources of density variations likely contribute. Collectively, results from this study of GRAIL gravity data, combined with other studies of remote sensing data and lunar samples, show that the lunar crust exhibits variations in density by plus or minus 10 percent over scales ranging from centimeters to 100’s of kilometers.

J C Jansen↗
DOE OSTI2016

Experimental study on magnetically insulated transmission line electrode surface evolution process under MA/cm current density

The design of high-current density magnetically insulated transmission line (MITL) is a difficult problem of current large-scale Z-pinch device. In particular, a thorough understanding of the MITL electrode surface evolution process under high current density is lacking. On the “QiangGuang-I” accelerator, the load area possesses a low inductance short-circuit structure with a diameter of 2.85 mm at the cathode, and three reflux columns with a diameter of 3 mm and uniformly distributed circumference at the anode. The length of the high density MITL area is 20 mm. A laser interferometer is used to assess and analyze the state of the MITL cathode and anode gap, and their evolution process under high current density. Experimental results indicate that evident current loss is not observed in the current density area at pulse leading edge, and peak when the surface current density reaches MA/cm. Analysis on electrode surface working conditions indicates that when the current leading edge is at 71.5% of the peak, the total evaporation of MITL cathode structure can be realized by energy deposition caused by ohmic heating. The electrode state changes, and diffusion conditions are reflected in the laser interferometer image. The MITL cathode area mainly exists in metal vapor form. The metal vapor density in the cathode central region is higher than the upper limit of laser penetration density (∼4 × 10{sup 21}/cm{sup 3}), with an expansion velocity of ∼0.96 km/s. The metal vapor density in the electrode outer area may lead to evident distortion of fringes, and its expansion velocity is faster than that in the center area (1.53 km/s).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Polyethylene density uncertainty in the simulation of neutron detectors

Simulations with the Monte Carlo N-Particle® (MCNP®) code are frequently used to study neutron coincidence counters for international safeguards. Reducing sources of uncertainty enables more accurate simulations, which reduces the need for expensive or impossible measurements. Accurately specifying the density of polyethylene is crucial for accurate simulations. The density of polyethylene ranges from 0.88 to 0.97 g/cm 3 and high-density polyethylene (HDPE) ranges from about 0.944 to 0.965 g/cm 3 . The density of a specific detector should be measured exactly for the best simulations. Here, the density of the blue/white Active Well Coincidence Counter (AWCC) bulk polyethylene was measured to be 0.9594 +- 0.0013 g/cm 3 at 1-sigma confidence. The density was calculated from repeated measurements of the various polyethylene dimensions for volume, and weighing the polyethylene for the mass. The AWCC neutron detection efficiency was simulated for each polyethylene density. The relative range of efficiencies (3σ) was 3.4% for unspecified polyethylene, 0.9% for HDPE, and 0.3% for the measured density. The impact of modeling small features such as screws was also studied and was found to be a negligible 0.2%. Measuring polyethylene density can reduce its role in uncertainty from one of the largest to one of the smallest contributors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

The non-thermal origin of the tokamak low-density stability limit

DIII-D plasmas at very low density exhibit onset of n=1 error field (EF) penetration (the `low-density locked mode') not at a critical density or EF, but instead at a critical level of runaway electron (RE) intensity. Raising the density during a discharge does not avoid EF penetration, so long as RE growth proceeds to the critical level. Penetration is preceded by non-thermalization of the electron cyclotron emission, anisotropization of the total pressure, synchrotron emission shape changes, as well as decreases in the loop voltage and bulk thermal electron temperature. The same phenomena occur despite various types of optimal EF correction, and in some cases modes are born rotating. Similar phenomena are also found at the low-density limit in JET. These results stand in contrast to the conventional interpretation of the low-density stability limit as being due to residual EFs and demonstrate a new pathway to EF penetration instability due to REs. Existing scaling laws for penetration project to increasing EF sensitivity as bulk temperatures decrease, though other possible mechanisms include classical tearing instability, thermo-resistive instability, and pressure-anisotropy driven instability. Regardless of first-principles mechanism, known scaling laws for Ohmic energy confinement combined with theoretical RE production rates allow rough extrapolation of the RE criticality condition, and thus, the low-density limit to other tokamaks. Furthermore, the extrapolated low-density limit by this pathway decreases with increasing machine size and is considerably below expected operating conditions for ITER. While likely unimportant for ITER, this effect can explain the low-density limit of existing tokamaks operating with small residual EFs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electrode biasing maintains the edge shear layer at high density in the J-TEXT tokamak

Collapse of the edge flow shear as the line-averaged density approaches the Greenwald density limit has been observed as a precursor to the enhanced edge particle flux characteristic of proximity to the density limit regime. Here, we report the use of a biased electrode to sustain the edge shear layer in high density discharges, in which the shear layer would otherwise collapse. A stable increase in line-averaged density is observed along with a strong increase in edge density. These experiments were carried out on the J-TEXT tokamak. The Reynolds stress at the edge is enhanced, and the zonal flow sustained, while density perturbation levels, the flux of turbulence internal energy (i.e., turbulence spreading), and particle and heat flux all decrease significantly. Electron adiabaticity increases, and bias voltage modulation experiments show that an increase in the edge shear leads the increase in adiabaticity. Furthermore, these results suggest that external edge E × B flow shear drive may be of interest for sustaining edge plasma states at high density, and support the hypothesis that collapse of the edge shear layer triggers the onset of the strong transport and turbulence characteristic of the density limit regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Disentangling core and edge mechanisms of the density limit in DIII-D negative triangularity plasmas

The density limit is investigated in the DIII-D negative triangularity plasmas which lack a standard H-mode edge. We find the limit may not be a singular disruptive boundary but a multifaceted density saturation phenomenon governed by distinct core and edge transport mechanisms. Sustained, non-disruptive operation is achieved at densities up to 1.8 times the Greenwald limit ($n_\mathrm{G}$) until the termination of auxiliary heating. Systematic power scans show distinct power scalings for the core ($n_e \propto P_\mathrm{SOL}^{0.27\pm0.03}$) and edge ($n_e \propto P_\mathrm{SOL}^{0.42\pm0.04}$) density limits. The edge density saturation is triggered by the onset of a non-disruptive, high-field side radiation front and the associated cooling, which clamps the edge density below . In contrast, the core density continues to rise until it saturates, a state characterized by enhanced core turbulence. Core transport evolves from a diffusive to an intermittent, avalanche-like state, as indicated by heavy-tailed probability density functions (kurtosis ≈ 6), increased Hurst exponents, and a $1/f$-type power spectrum. These findings suggest that the density limit in the low-confinement regime is determined by a combination of edge radiative cooling and core turbulent transport. This distinction provides separate targets for control strategies aimed at extending the operational space of future fusion devices.

density limit↗
DOE OSTI2019

Rethinking CO adsorption on transition-metal surfaces: Effect of density-driven self-interaction errors

Adsorption of the molecule CO on metallic surfaces is an important unsolved problem in Kohn-Sham density functional theory (KS-DFT). We present a detailed study of carbon monoxide adsorption on fcc (111) surfaces of 3d, 4d, and 5d metals using nonempirical semilocal density functionals for the exchange-correlation energy: the local-density approximation (LDA), two generalized gradient approximations or GGAs [Perdew-Burke-Ernzerhof (PBE) and PBE for solids (PBEsol)], and a meta-GGA [strongly constrained and appropriately normed (SCAN) functional]. The typical error pattern (as found earlier for free molecules and for free transition-metal surfaces), in which results improve from LDA to PBE or PBEsol to SCAN, due to the satisfaction of more exact constraints, is not found here. Instead, for CO adsorption on transition-metal surfaces, we find that, while SCAN overbinds much less than LDA, it overbinds slightly more than PBE. Moreover, the tested functionals often predict the wrong adsorption site, as first pointed out for LDA and GGA in the CO/Pt (111) puzzle. This abnormal pattern leads us to suspect that the errors of PBE and SCAN for this problem are density-driven self-interaction errors associated with incorrect charge transfer between molecule and metal surface. We point out that, by the variational principle, overbinding by an approximate functional would be reduced if that functional were applied not to its self-consistent density for the adsorbed system but to an exact or more correct density for that system. Finally, we show for CO on Pt(111) that the site preference is corrected and the adsorption energy is improved for the PBE functional by using not the self-consistent PBE density but a PBE+U density. The resulting correction to the PBE total energy is much larger for the adsorbed system than for its desorbed components, showing that the error is in the density of the adsorbed system. Furthermore, this seems to solve the Feibelman 2001 CO/Pt(111) puzzle, in principle if not fully in practice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Advanced High Torque Density Magnetically Geared Generator

The objective of this research project was to investigate the fundamental torque density limitations of multistage flux focusing magnetic gears and magnetically geared direct-drive generators for wind turbines. Magnetic gears create non-contact speed amplification by modulating the air-gap magnetic field harmonics created by magnetic rotors. By integrating a stator into the magnetic gear structure a very high torque density direct-drive generator can be created. Magnetic geared systems offer unique capabilities that are currently not available using existing mechanically geared systems, such as overload protection. Magnetic gears also offer the potential for higher reliability, higher efficiency and quieter operation than is currently achievable using a mechanically geared generator. By increasing reliability and design life magnetic gears could decrease the overall life-cycle costs when compared with equivalent mechanical or direct-drive drivetrains. This research project investigated and experimentally testing two sub-scale magnetic gearboxes (0.12 m diameter) and two large diameter magnetic gearboxes (0.56m and 0.63m diameter). The goal of the project was to demonstrate a torque density >300Nm/L for a multistage magnetic gearbox with integrated stator. A nested multistage magnetic gear structure was first proposed and studied. However, it was shown that the nested design greatly increased mechanical complexity and therefore a two-stage series typology was utilized. A multi-stage magnetic gear with a 59:1 gear ratio was successfully experimentally tested. A peak torque of 4.25 kN·m was measured and this resulted in an overall active region torque density of 141 N·m/L (when including support rings). The torque density was significantly reduced from what was computed using 2-D electromagnetic analysis. This reduction was due to the 3-D edge effects, increased power losses and the inclusion of additional mechanical support rings as well as other mechanical assembly design changes. At large diameters contiguous lamination parts could not be fabricated and so a segmented laminated rotor design had to be used. This reduced the mechanical strength of the design when compared to small diameter magnetic gears and resulted in the need for the addition of mechanical support rings to prevent deflection. The final experimentally tested design also had a high torque ripple and the power loss within the magnetic gear did not change significantly with load. This therefore resulted in the efficiency being high only at peak loads. Despite not experimentally achieving the 2-D calculated high torque density the analysis has shown that the problems encountered with the scaling can be overcome with careful design changes. For example, understanding the unique mechanical deflection requirements during the electromagnetic design of the large diameter magnetic gears is key to achieving both a high torque density based on the electromagnetic analysis and a high torque density based on the mechanical support needs. Furthermore, ensuring that end plates are made of non-conductive material and minimizing other conductive support elements within the magnetic structure is a key requirement for reducing eddy current losses. In summary, if the mechanical design can be made more robust and eddy current loss mitigation strategies are used in the mechanical design phase then the magnetic gear is capable of being scaled up whilst operating at a very high torque density. This will then enable it to be potentially performance and cost competitive with alternative direct-drive and mechanically geared wind turbine drivetrains.

17 WIND ENERGY↗
DOE OSTI2022

Development of a Direct Neutral Density Diagnostic for Fusion Edge Plasmas

As magnetically confined plasmas progress towards ignition and long pulse experiments, measurement and control of the neutral density in the plasma edge has become a critical issue for stability, formation of transport barriers, and fueling. Recent experiments by our research group have demonstrated that is possible to use two-photon absorption laser induced fluorescence (TA-LIF) to directly measure the density of neutral hydrogen in helicon sources and in the HIT-SI3 spheromak. While those experiments validated key elements of a diagnostic system that would enable similar measurements in tokamak plasmas, they did not fully address the issue of performance optimization in the presence of intense background light at the fluorescence wavelength and they also identified a number of other issues that must be resolved before successful TALIF neutral density measurements in a tokamak are likely to be achieved. Additional specific concerns raised during design reviews with the leadership of the DIII-D tokamak facility included: the minimum detection threshold for this TALIF diagnostic (currently ~ 5 x 10 15 m -3 ) and the rate at which neutral density measurements could be obtained. Improving the system performance to reduce the minimum detection threshold and demonstrating a faster rate of density determination are required to advance this diagnostic to the level where it could be considered for implementation on a major tokamak facility in the USA. The key issues that will be addressed through additional technological development are: validation of a new, chromatic aberration-free, xenon calibration scheme; increasing the output power of the laser; and suppression of background emission at the fluorescence wavelength through optimization of the collection optics and gating of the photomultiplier detector with a goal of obtaining a minimum detection threshold of 5 x 10 14 m -3 . Another key technological development, so far only tested in helicon source experiments, is the validation of Doppler-free TALIF as a means of obtaining higher speed, calibrated, neutral density measurements in tokamak-like conditions at the full cadence of the pulsed TALIF laser. In this work, we propose to complete these technological advancements through installation and testing of our prototype TALIF system on the proto-MPEX experiment at Oak Ridge National Laboratory (ORNL). The proto-MPEX facility will provide plasma conditions similar to the edge plasma of a major tokamak experiment but with pulse lengths and pulse repetition rates ideally suited for extensive development of a TALIF neutral density diagnostic. Measurement of the absolute neutral density in the edge of a magnetically confined plasma is necessary for plasma density control; calculation of charge-exchange power losses; control of plasma-wall interactions; determination of the braking of plasma flow; determination of the fuel mixture in deuterium-tritium plasma; and understanding the dynamics of the divertor region in the plasma edge. In terms of potential impact on the worldwide fusion program, we note that in burning plasma experiments, such as ITER, a 100 mJ/pulse Doppler-free TALIF diagnostic should be capable of directly measuring the deuterium/tritium (D/T) isotope ratio in the outer 0.5 to 1.0 m of the plasma radius. The D/T ratio is a critically important parameter for the control and optimization of burning plasmas. Therefore, development of the diagnostic system proposed here is also relevant to long-term US participation in diagnostic development for large tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Implementation and Investigation of a Weather- and Jam Density-tuned Network Perimeter Controller

This paper implements and evaluates the performance of a network fundamental diagram (NFD) proportional-integral perimeter controller (PC) using base tuned parameters (clear weather conditions and a base jam density of 160 veh/km/lane). The parameters were then re-tuned separately for different weather conditions and jam density values (reflecting different percentage of trucks), resulting in two new control methods: a weather-tuned perimeter controller (WTPC) and a jam density-tuned perimeter controller (JTPC). The WTPC was shown to outperform the no control strategy and the PC for different weather conditions. Specifically, the WTPC decreased the congestion inside the protected network (PN) and improved the overall performance of the full network (FN) by decreasing the average vehicle travel time, decreasing the vehicle total delay, increasing the average vehicle speed and decreasing the average vehicle fuel consumption. Alternatively, the JTPC was shown to perform similar to the PC for jam densities higher than the base-case jam density used in the PC (in our case 160 veh/km/lane). However, the JTPC outperformed the PC for smaller jam densities given that these smaller jam densities result in queues spilling back faster to upstream traffic signals. The results demonstrate the need to tune the controller to the actual jam density especially when the jam density decreases (e.g. trucks are introduced to the network). http://pubsindex.trb.org/

42 ENGINEERING↗
DOE OSTI2020

Variable-density effects in incompressible non-buoyant shear-driven turbulent mixing layers

The asymmetries that arise when a mixing layer involves two miscible fluids of differing densities are investigated using incompressible (low-speed) direct numerical simulations. The simulations are performed in the temporal configuration with very large domain sizes, to allow the mixing layers to reach prolonged states of fully turbulent self-similar growth. Imposing a mean density variation breaks the mean symmetry relative to the classical single-fluid temporal mixing layer problem. Unlike prior variable-density mixing layer simulations in which the streams are composed of the same fluids with dissimilar thermodynamic properties, the density variations are presently due to compositional differences between the fluid streams, leading to different mixing dynamics. Variable-density (non-Boussinesq) effects introduce strong asymmetries in the flow statistics that can be explained by the strongest turbulence increasingly migrating to the lighter fluid side as free-stream density difference increases. Interface thickness growth rates also reduce, with some thickness definitions particularly sensitive to the corresponding changes in alignment between density and streamwise velocity profiles. Additional asymmetries in the sense of statistical distributions of densities at a given position within the mixing layer reveal that fine scales of turbulence are preferentially sustained in lighter fluid, which also is where fastest mixing occurs. These effects influence statistics involving density fluctuations, which have important implications for mixing and more complicated phenomena that are sensitive to the mixing dynamics, such as combustion.

42 ENGINEERING↗
DOE OSTI2018

Effect of pore density on gas permeation through nanoporous graphene membranes

Pore density is an important factor dictating gas separations through one-atom-thin nanoporous membranes, but how it influences the gas permeation is not fully understood. We use molecular dynamics (MD) simulations to investigate gas permeation through nanoporous graphene membranes with the same pore size (3.0 Å × 3.8 Å in dimensions) but varying pore densities (from 0.01 to 1.28 nm -2 ). We find that higher pore density leads to higher permeation per unit area of membrane for both CO 2 and He, but the rate of the increase decreases greatly for CO 2 at high pore densities. As a result, the per-pore permeance decreases for CO 2 but remains relatively constant for He with the pore density, leading to a dramatic change in CO 2 /He selectivity. By separating the total flux into direct flux and surface flux, we find that He permeation is dominated by direct flux and hence the per-pore permeation rate is roughly constant with the pore density. In contrast, CO 2 permeation is dominated by surface flux and the overall decreasing trend of the per-pore permeation rate of CO 2 with the pore density can be explained by the decreasing per-pore coverage of CO 2 on the feed side with the pore density. Finally, our work now provides a complete picture of the pore-density dependence of gas permeation through one-atom-thin nanoporous membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Uncertainty quantification in density estimation from background-oriented Schlieren measurements

In this paper, we present an uncertainty quantification methodology for density estimation from background-oriented Schlieren (BOS) measurements, in order to provide local, instantaneous, a posteriori uncertainty bounds on each density measurement in the field of view. Displacement uncertainty quantification algorithms from cross-correlation-based particle image velocimetry are used to estimate the uncertainty in the dot pattern displacements obtained from cross-correlation for BOSs and assess their feasibility. In order to propagate the displacement uncertainty through the density integration procedure, we also develop a novel methodology via the Poisson solver using sparse linear operators. Testing the method using synthetic images of a Gaussian density field showed agreement between the propagated density uncertainties and the true uncertainty. Subsequently, the methodology is experimentally demonstrated for supersonic flow over a wedge, showing that regions with sharp changes in density lead to an increase in density uncertainty throughout the field of view, even in regions without these sharp changes. Lastly, the uncertainty propagation is influenced by the density integration scheme, and for the Poisson solver the density uncertainty on average increases on moving away from the regions where the Dirichlet boundary conditions are specified.

47 OTHER INSTRUMENTATION↗