Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic separation”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Sr-Nd-Pb isotope systematics of Australasian tektites: Implications for the nature and composition of target materials and possible volatile loss of Pb

The Australasian tektite (AAT) strewn field is the largest strewn field on the Earth with about ~10–30% coverage, both land and ocean, but a clearly identified source impact crater is absent despite the young age of AAT of ca. 790 ka. A genetic link between the Australasian tektites and their unequivocal parental materials is therefore largely impossible to establish. Nevertheless, the nature of the parental materials and the extent of volatilization can be constrained using the splash form tektites, carrying the chemical signatures of high-temperature processes, and the layered (so-called Muong Nong-type) tektites, which are less chemically homogenized and exceptionally abundant in the AAT field. New high-precision Sr, Nd and Pb isotopic measurements were obtained for a chemically and petrographically well-characterized suite of AAT, which included the Muong Nong-type (MN-AAT) with precisely known field locations in Laos and splash forms (SF-AAT) from different parts of the strewn field. Additionally, optically dark and light zones of the MN-AAT were also separately analyzed. Homogeneous ε Nd values from -11.8 to -11.2, combined with a narrow range of two-stage Nd model ages from 1.67 to 1.72 Ga for the entire AAT suite, point to a well-mixed source, in terms of REE, of the crustal segment from which the sedimentary material for tektites was ultimately derived. The Sr isotopic data largely overlap for SF-AAT and MN-AAT ( 87 Sr/ 86 Sr = 0.71636–0.72021) and indicate Paleozoic to Mesozoic sedimentary parentage. However, late Neogene to early Quaternary re-deposition and formation of a thick silt-sized sedimentary section with vertical stratification is required to comply with 10Be data. Lead isotope systematics documents at least three different components which can perhaps be represented by different mineral phases, such as feldspar, zircon, organic matter adsorbed on young sediments etc., sorted during fluvial transport and final deposition. In addition, the SF-AAT have systematically lower Pb contents than the MN-AAT, and generally show isotopically heavier Pb isotopic ratios. This is theoretically consistent with a preferential volatilization of lighter Pb isotopes during evaporation and considerably larger Pb loss from SF-AAT when compared to MN-AAT. Nevertheless, further experimental work would be necessary to unambiguously distinguish kinetic fractionation from source mixing.

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

In-situ thermodynamics measurements at metal oxides-solution interfaces using Flow Adsorption Microcalorimetry.

Mineral-fluid interfaces are the principal sites of geochemical processes near Earth’s surface, hosting chemical reactions that play a fundamental role in (bio)geochemical cycles and in the fate and transport of anthropogenic contaminants, and that effectively control the compositions of soil and water environments. These complex interfaces are critical for our energy and environmental future. The mineral-fluid interface has been studied in an unprecedented level of detail with both experimental and computational approaches, separately and in combination. However, conspicuously missing from studies of the mineral-fluid interface are direct measurements of the energies of ion sorption and exchange. The literature on energetics and enthalpies of exchange, adsorption, dissolution, precipitation, and surface protonation reactions, especially those directly supported by experimental data, remains scarce despite their fundamental nature. The overarching goal of this project is to complete a systematic study of the thermodynamics properties of interfacial reactions at four MO surfaces (Rutile (α-TiO2), Quartz (SiO2), boehmite (γ-AlOOH) and goethite (α-FeOOH)) through the application and construction of novel flow adsorption microcalorimetry techniques and instrumentations. These unique and specialized microcalorimeters will operate at various temperatures and solution chemical compositions allowing for in-situ measurements across metal oxides and ligands of various characteristics. The overall research goal will be accomplished by completing the following three specific objectives (O): O1) Determine the energetics of surface protonation and deprotonation, ion exchange and ligand sorption reactions; O2) Investigate the surface charge thermodynamic properties under a range of temperature and solution chemical compositions; and O3) Develop predictive trends of the interplay between MO structure, surface coverage and surface reactivity. In addition to key thermodynamics parameters, calorimetric measurements provide a wealth of mechanistic information about reactions energetics and kinetics, surface charge characteristics, and structure-reactivity or selectivity relationships, all obtained in-situ and in real-time. This report includes science highlights from various projects completed over the performance period of the project. Also listed are dissemination opportunities, people supported on the grant and the impact on available physical resources and the discipline as a whole.

58 GEOSCIENCES↗
NASA NTRS2012

Can We Infer Ocean Dynamics from Altimeter Wavenumber Spectra?

The wavenumber spectra of sea surface height (SSH) and kinetic energy (KE) have been used to infer the dynamics of the ocean. When quasi-geostrophic dynamics (QG) or surface quasi-geostrophic (SQG) turbulence dominate and an inertial subrange exists, a steep SSH wavenumber spectrum is expected with k-5 for QG turbulence and a flatter k-11/3 for SQG turbulence. However, inspection of the spectral slopes in the mesoscale band of 70 to 250 km shows that the altimeter wavenumber slopes typically are much flatter than the QG or SQG predictions over most of the ocean. Comparison of the altimeter wavenumber spectra with the spectra estimated from the output of an eddy resolving global ocean circulation model (the Hybrid Coordinate Ocean Model, HYCOM, at 1/25 resolution), which is forced by high frequency winds and includes the astronomical forcing of the sun and the moon, suggests that the flatter slopes of the altimeter may arise from three possible sources, the presence of internal waves, the lack of an inertial subrange in the 70 to 250 km band and noise or submesoscales at small scales. When the wavenumber spectra of SSH and KE are estimated near the internal tide generating regions, the resulting spectra are much flatter than the expectations of QG or SQG theory. If the height and velocity variability are separated into low frequency (periods greater than 2 days) and high frequency (periods less than a day), then a different pattern emerges with a relatively flat wavenumber spectrum at high frequency and a steeper wavenumber spectrum at low frequency. The stationary internal tides can be removed from the altimeter spectrum, which steepens the spectral slopes in the energetic internal wave regions. Away from generating regions where the internal waves

Richman, James↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Flow Characteristics of UHV Nitrogen Arcjet

Flow characteristics in the nozzle of a nitrogen arcjet have been simulated by a model. Electrons, ions, and neutral atoms and molecules are represented as separate species. Dissociation, ionization, and recombination are treated as separate reactions. Thermal non-equilibrium is represented by a two-temperature model. Energy input to plasma from the arc is modeled as a source determined by local current and electrical conductivity. Momentum and energy losses and recombinations at the nozzle wall are included as source/sink determined by fluxes to the wall. The results show essentially frozen gas-phase reactions and thermal nonequilibrium due to Joule heating of electrons and in the expansion part of the nozzle. Wall interactions have strong effects on the results, indicating that they play important roles in the flow due to very high surface area relative to the volume of the plasma. For example, plasma velocity decreases from 6000 m per s without wall interactions to 3000 m per s with modest amount of wall interactions. These results also show reasonable agreement with optical emission measurements, which confirms that the arcjet plasma is far from LTE. The spectra suggest nitrogen dissociation levels of 0.3% - 9%, depending on nitrogen flow rate and arc plasma current, which also determine the relative amounts of excited atomic and molecular nitrogen. Langmuir probe studies of the source show that electron and ion fluxes increase with increasing power, and that the ion energy distribution shifts to lower energies. Typical ion fluxes were on the order of 4E-9 A/cm**2 with a maximum ion kinetic energy of 3.5eV. The median electron energy was 1eV, with a maximum of less than 4eV.

Chang, C. H.↗
NASA NTRS2003

High-Time-Resolution Study of Magnetic Holes in the Solar Wind

The objectives of this investigation are to determine the internal plasma structure of kinetic-scale and larger scale magnetic holes, and to determine their stability, their source mechanism(s), and their spatial extent. It is also of importance to determine the relationship between kinetic-scale holes and long-duration holes. As smaller and smaller magnetic depressions are investigated in order to make this a complete study, a robust criterion is necessary for distinguishing magnetic holes from random or unresolvable fluctuations in the interplanetary magnetic field. In order to resolve this ambiguity, we obtained from the MFI experiments magnetic field measurements from the WIND spacecraft at a time resolution of 46 to 184 ms over certain periods. We have also devised a measure of certainty for magnetic hole detections. The certainty factor, q, is defined as the difference between the mean magnetic field in the hole and the local magnetic field, in units of the local standard deviation of the field strength. For fullest generality, it is necessary to calculate this q over the range of available scales of interest, from 60 ms up to 300 s. This technique establishes a two dimensional matrix of relative probabilities that a hole of some duration (d) might exist in the data set at a given time (t). In identifying q-peaks in time and duration, we also come upon a natural method for distinguishing holes with internal structure from multiple holes in close proximity or holes nested inside of others. If two q-peaks are more than a half-width apart, they are simply said to be separate events.

Lazarus, Alan↗
DOE OSTI2020

Radiation-Induced Effects on the Extraction Properties of Hexa-n-octylnitrilo-triacetamide (HONTA) Complexes of Americium and Europium

The candidate An(III)/Ln(III) separation ligand hexa-n-octylnitrilo-triacetamide (HONTA) was irradiated under envisioned SELECT (Solvent Extraction from Liquid waste using Extractants of CHON-type for Transmutation) process conditions (n-dodecane/0.1 M HNO 3 ) using a solvent test loop in conjunction with cobalt-60 gamma irradiation. The extent of HONTA radiolysis and complimentary degradation product formation was quantified by HPLC-ESI-MS/MS. Further, the impact of HONTA radiolysis on process performance was evaluated by measuring the change in 243 Am and 154 Eu distribution ratios as a function of absorbed gamma dose. HONTA was found to decay exponentially with increasing dose, affording a dose coefficient of d = (4.48 ± 0.19) × 10 –3 kGy –1 . Multiple degradation products were detected by HPLC-ESI-MS/MS with dioctylamine being the dominant quantifiable species. Both 243 Am and 154 Eu distribution ratios exhibited an induction period of ~70 kGy for extraction (0.1 M HNO 3 ) and back-extraction (4.0 M HNO 3 ) conditions, after which both values decreased with absorbed dose. The decrease in distribution ratios was attributed to a combination of the destruction of HONTA and ingrowth of dioctylamine, which is capable of interfering in metal ion complexation. Furthermore, the loss of HONTA with absorbed gamma dose was predominantly attributed to its reaction with the n-dodecane radical cation (R˙ + ). These R˙ + reaction kinetics were measured for HONTA and its 241 Am and 154 Eu complexes using picosecond pulsed electron radiolysis techniques. All three second-order rate coefficients (k) were essentially diffusion limited in n-dodecane indicating a significant reaction pathway: k(HONTA + R˙ + ) = (7.6 ± 0.8) × 10 9 M –1 s –1 , k(Am(HONTA) 2 + R˙ + ) = (7.1 ± 0.7) × 10 10 M –1 s –1 , and k(Eu(HONTA) 2 + R˙ + ) = (9.5 ± 0.5) × 10 10 M –1 s –1 . HONTA-metal ion complexation afforded an order-of-magnitude increase in rate coefficient. Additional nanosecond time-resolved measurements showed that both direct and indirect HONTA radiolysis yielded the short-lived (< 100 ns) HONTA radical cation and a second long-lived (µs) species identified as the HONTA triplet excited state. The latter was confirmed by a series of oxygen quenching picosecond pulsed electron measurements, affording a quenching rate coefficient of k( 3 [HONTA]* + O 2 ) = 2.2 × 10 8 M –1 s –1 . Overall, both the HONTA radical cation and triplet excited state are important precursors to the suite of measured HONTA degradation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1991

Computer-based video digitizer analysis of surface extension in maize roots: kinetics of growth rate changes during gravitropism

We used a video digitizer system to measure surface extension and curvature in gravistimulated primary roots of maize (Zea mays L.). Downward curvature began about 25 +/- 7 min after gravistimulation and resulted from a combination of enhanced growth along the upper surface and reduced growth along the lower surface relative to growth in vertically oriented controls. The roots curved at a rate of 1.4 +/- 0.5 degrees min-1 but the pattern of curvature varied somewhat. In about 35% of the samples the roots curved steadily downward and the rate of curvature slowed as the root neared 90 degrees. A final angle of about 90 degrees was reached 110 +/- 35 min after the start of gravistimulation. In about 65% of the samples there was a period of backward curvature (partial reversal of curvature) during the response. In some cases (about 15% of those showing a period of reverse bending) this period of backward curvature occurred before the root reached 90 degrees. Following transient backward curvature, downward curvature resumed and the root approached a final angle of about 90 degrees. In about 65% of the roots showing a period of reverse curvature, the roots curved steadily past the vertical, reaching maximum curvature about 205 +/- 65 min after gravistimulation. The direction of curvature then reversed back toward the vertical. After one or two oscillations about the vertical the roots obtained a vertical orientation and the distribution of growth within the root tip became the same as that prior to gravistimulation. The period of transient backward curvature coincided with and was evidently caused by enhancement of growth along the concave and inhibition of growth along the convex side of the curve, a pattern opposite to that prevailing in the earlier stages of downward curvature. There were periods during the gravitropic response when the normally unimodal growth-rate distribution within the elongation zone became bimodal with two peaks of rapid elongation separated by a region of reduced elongation rate. This occurred at different times on the convex and concave sides of the graviresponding root. During the period of steady downward curvature the elongation zone along the convex side extended farther toward the tip than in the vertical control. During the period of reduced rate of curvature, the zone of elongation extended farther toward the tip along the concave side of the root. The data show that the gravitropic response pattern varies with time and involves changes in localized elongation rates as well as changes in the length and position of the elongation zone. Models of root gravitropic curvature based on simple unimodal inhibition of growth along the lower side cannot account for these complex growth patterns.

NASA Discipline Plant Biology↗
NASA NTRS1989

Numerical investigation of an internal layer in turbulent flow over a curved hill

The development of an internal layer in a turbulent boundary layer flow over a curved hill is investigated numerically. The turbulence field of the boundary layer flow over the curved hill is compared with that of a turbulent flow over a symmetric airfoil (which has the same geometry as the curved hill except that the leading and trailing edge plates were removed) to study the influence of the strongly curved surface on the turbulence field. The turbulent flow equations are solved by a control-volume based finite difference method. The turbulence is described by a multiple-time-scale turbulence model supplemented with a near-wall turbulence model. Computational results for the mean flow field (pressure distributions on the walls, wall shearing stresses and mean velocity profiles), the turbulence structure (Reynolds stress and turbulent kinetic energy profiles), and the integral parameters (displacement and momentum thicknesses) compared favorably with the measured data. Computational results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. Development of the internal layer was more obviously observed in the Reynolds stress profiles and in the turbulent kinetic energy profiles than in the mean velocity profiles. In this regard, the internal layers is significantly different from wall-bounded simple shear layers in which the mean velocity profile characterizes the boundary layer most distinguishably. Development of such an internal layer, characterized by an intense turbulence field, is attributed to the enormous mean flow strain rate caused by the streamline curvature and the strong pressure gradient. In the turbulent flow over the curved hill, the internal layer begin to form near the forward corner of the hill, merges with the external boundary layer, and develops into a new fully turbulent boundary layer as the fluid flows in the downstream direction. For the flow over the symmetric airfoil, the boundary layer began to form from almost the same location as that of the curved hill, grew in its strength, and formed a fully turbulent boundary layer from mid-part of the airfoil and in the downstream region. Computational results also show that the detailed turbulence structure in the region very close to the wall of the curved hill is almost the same as that of the airfoil in most of the curved regions except near the leading edge. Thus the internal layer of the curved hill and the boundary layer of the airfoil were also almost the same. Development of the wall shearing stress and separation of the boundary layer at the rear end of the curved hill mostly depends on the internal layer and is only slightly influenced by the external boundary layer flow.

Kim, S-W.↗
NASA NTRS2001

The catecholamine response to spaceflight: role of diet and gender

Compared with men, women appear to have a decreased sympathetic nervous system (SNS) response to stress. The two manifestations where the sexual dimorphism has been the most pronounced involve the response of the SNS to fluid shifts and fuel metabolism during exercise. The objectives of this study were to investigate whether a similar sexual dimorphism was found in the response to spaceflight. To do so, we compared catecholamine excretion by male and female astronauts from two similar shuttle missions, Spacelab Life Sciences 1 (SLS1, 1991) and 2 (SLS2, 1993) for evidence of sexual dimorphism. To evaluate the variability of the catecholamine response in men, we compared catecholamine excretion from the two SLS missions against the 1996 Life and Microgravity Sciences Mission (LMS) and the 1973 Skylab missions. RESULTS: No gender- or mission-dependent changes were found with epinephrine. Separating out the SLS1/2 data by gender shows that norepinephrine excretion was essentially unchanged with spaceflight in women (98 +/- 10%; n = 3) and substantially decreased with the men (41 +/- 9%; n = 4, P < 0.05). Data are a percentage of mean preflight value +/- SE. Comparisons among males demonstrated significant mission effects on norepinephrine excretion. After flight, there was a transient increase in norepinephrine but no evidence of any gender-specific effects. We conclude that norepinephrine excretion during spaceflight is both mission and gender dependent. Men show the greater response, with at least three factors being involved, a response to microgravity, energy balance, and the ratio of carbohydrate to fat in the diet.

NASA Discipline Regulatory Physiology↗
DOE OSTI2021

Recent Advances in Nonclassical Crystallization: Fundamentals, Applications, and Challenges

We report crystals are ubiquitous and play an important role in our daily life. Recent studies on properties and formation of crystals have proven the attachment of particle ranging from ion pairs to well-crystallized nanocrystal is one of the most important ways for crystallization. Compared with the monomer-by-monomer (i.e., atoms, ions, or molecules) addition described in longstanding classical theories, so-called nonclassical crystallization pathways are more complex as the non-specific interaction of variation in system free energy and reaction kinetics. In view of the detailed insights into nonclassical crystallization pathways has broadened the scope of research across various disciplines such as geological events, biomineralization mechanisms, environmental remediation, and the development of functional materials, this paper reviews the current state of the art in nonclassical nucleation (the prenucleation clusters pathway, aggregation pathway, etc.) and growth pathways (oriented attachment and random attachment) based on several spotlights. In addition, the observed specific phenomena on crystallization in the last 10 years as well as recent advancements and achievements made on the synthesis of mesocrystals are also elucidated. With implications and perspectives of nonclassical crystallization researches proposed in this paper, this review aims to summarize research advancement regarding crystallizing systems. Emphasis will be placed on novel insights into nonclassical crystallization of materials, which are expected to provide completely new perspectives on the underlying crystallization mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗
NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗
DOE OSTI2020

Overcoming Fundamental Limitations in Adsorbent Design: Alkene Adsorption by Non-porous Copper(I) Complexes

Purifying alkenes from alkanes requires cryogenic distillation. This consumes energy equivalent to countries of ca. 5 million people. Replacing distillation with adsorption processes would significantly increase energy efficiency. Trade-offs between kinetics, selectivity, capacity, and heat of adsorption have prevented production of an optimal adsorbent. Here, we report adsorbents that overcome these trade-offs. [Cu-Br] 3 and [Cu-H] 3 are air-stable trinuclear complexes that undergo reversible solid-state inter-molecular rearrangements to produce dinuclear [Cu-Br∙(alkene)] 2 and [Cu-H∙(alkene)] 2 . The reversible solid-state rearrangement, confirmed in situ using powder X-ray diffraction, allows adsorbent design trade-offs to be overcome, coupling low heat of adsorption (-10 to -17 kJ mol -1 alkene ), high alkene:alkane selectivity (47; 29), and uptake capacity (>2.5 molalkene mol -1 Cu3 ). Most remarkably, [Cu-H] 3 displays fast uptake and regenerates capacity within 10 minutes.

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

Colloidal Clusters and Networks Formed by Oppositely Charged Nanoparticles with Varying Stiffnesses

Colloidal clusters and gels are ubiquitous in science and technology. Particle softness has a strong effect on interparticle interactions; however, our understanding of the role of this factor in the formation of colloidal clusters and gels is only beginning to evolve. Here, we report the results of experimental and simulation studies of the impact of particle softness on the assembly of clusters and networks from mixtures of oppositely charged polymer nanoparticles (NPs). Experiments were performed below or above the polymer glass transition temperature, at which the interaction potential and adhesive forces between the NPs were significantly varied. Hard NPs assembled in fractal clusters that subsequently organized in a kinetically arrested colloidal gel, while soft NPs formed dense precipitating aggregates, due to the NP deformation and the decreased interparticle distance. Importantly, interactions of hard and soft NPs led to the formation of discrete precipitating NP aggregates at a relatively low volume fraction of soft NPs. A phenomenological model was developed for interactions of oppositely charged NPs with varying softnesses. Furthermore, the experimental results were in agreement with molecular dynamics simulations based on the model. This work provides insight on interparticle interactions before, during, and after the formation of hard–hard, hard–soft, and soft–soft contacts and has impact for numerous applications of reversible colloidal gels, including their use as inks for additive manufacturing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2018

Detonation to Deflagration-Mode Transitions in Pulsed Plasma Accelerators

PULSED plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas. The current in such discharges will ionize the gas and produce a strong magnetic field, which interacts with the flowing current to accelerate the plasma through the Lorentz body force. For the present work, two plasma accelerator types employing this general scheme are of interest: the gas-fed pulsed plasma thruster (PPT) 1 and the quasi-steady magnetoplasmadynamic (MPD) accelerator. 2 The gas-fed pulsed plasma accelerator is generally understood as a completely transient device discharging in ∼1-10 μs. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j × B body force, entraining and accelerating propellant it encounters. This process is sometimes referred to in literature as 'snowplowing' the propellant or accelerating the gas in a detonation-mode because the current sheet representation approximates that of a strong detonation shockwave propagating through the gas. For these devices, acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new 'crowbar' discharge at the breech and depriving the initial sheet of additional acceleration. In general, PPTs typically claim thrust efficiencies (ratio of jet kinetic energy to input electrical energy) registering in the teens or lower. 3 In the quasi-steady MPD accelerator, the pulse is lengthened to ∼1 ms or longer and maintained at an approximately constant level during discharge through the use of a pulse-forming network (PFN) of capacitors. After an initial transient discharge, which is typically short relative to the overall discharge period, the plasma assumes a relatively steady-state configuration, known as 'quasi-steady' MPD operation. 2 In this state, ionized gas flows through a stationary current channel in a manner that is sometimes referred to as deflagration-mode operation owing to the similarities to deflagration waves in gases. The plasma experiences electromagnetic acceleration as the plasma flows through the current channel towards the exit of the device. Quasi-steady MPD thrusters claim efficiencies up to 50% for certain propellants.4 There has been significant and sustained research over several decades on both gas-fed PPTs and quasi-steady MPD thrusters, however there have been pulsed thrusters that do not appear to exactly fit either classification, instead possessing a mixture of operational qualities characteristic of both thruster variants. The Coaxial High ENerGy (CHENG) thruster by Cheng, et al.5 operated on the short 10 μs timescales characteristic of PPTs, but claimed the high thrust densities, high efficiencies, and low electrode erosion rates that are more consistent with the MPD/deflagration mode of plasma acceleration. Gas-fed PPT research by Ziemer, et al. 3, 6 identified two separate regimes of performance in those thrusters. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency as a function of mass bit, while the second regime at very low mass bits (termed Mode II) exhibited an increase in efficiency with decreasing mass bit. Work by Poehlmann et al.7 and by Sitaraman and Raja8 sought to understand the performance of the CHENG thruster and the Mode I/Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes of plasma acceleration representing two distinct sets of solutions to the relevant conservation laws. In these works, it was proposed that the values of the various controllable parameters determined whether the accelerator would operate in detonation or deflagration mode. Our hypothesized view of the acceleration process in the CHENG thruster and in PPTs experiencing a transition from Mode I to Mode II is inspired by observations of the transition from the PPT mode of operation to the quasi-steady MPD mode. Specifically, the quasi-steady MPD was discovered by driving a PPT to extended pulse lengths. Above a certain pulse length threshold the transient plasma current sheet transitions into a stable plasma acceleration mode that 'replicates in every observable detail steady flow self-field magnetoplasmadynamic acceleration.' 9 In the present work, instead of treating the accelerator as if it were only operating in a single mode during a pulse, we consider the initial stage of the discharge in all cases as a current sheet forming at the breach of the accelerator and moving towards the exit as a detonation wave. If the current sheet reaches the exit of the accelerator before the discharge is completed, the view of the acceleration mode transitions to the deflagration mode-type found in quasi-steady MPD thrusters. In previous work10 we presented a modeling framework that first captured the time-evolution of the current sheet (detonation) mode of the thruster and then transitioned into the quasi-steady MPD (deflagration) mode of plasma acceleration. In the present work, variations of the controllable parameters - specifically the pulsed circuit properties, the amount of mass injected into the thruster, and the relative timing between the initial gas injection and the initiation of the plasma current sheet - will be used to explore the thruster performance. A range of parameters are explored to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗
DOE OSTI2022

Selectivity of Electrochemical Ion Insertion into Manganese Dioxide Polymorphs

The ion insertion redox chemistry of manganese dioxide has diverse applications in energy storage, catalysis, and chemical separations. Unique properties derive from the assembly of Mn–O octahedra into polymorphic structures that can host protons and nonprotonic cations in interstitial sites. Despite many reports on individual ion-polymorph couples, much less is known about the selectivity of electrochemical ion insertion in MnO 2 . In this work, we use density functional theory to holistically compare the electrochemistry of A x MnO 2 (where A = H + , Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ , Al 3+ ) in aqueous and nonaqueous electrolytes. We develop an efficient computational scheme demonstrating that Hubbard-U correction has a greater impact on calculating accurate redox energetics than choice of exchange-correlation functional. Using PBE+U, we find that for nonprotonic cations, ion selectivity depends on the oxygen coordination environments inside a polymorph. When H + is present, however, the driving force to form hydroxyl bonds is usually stronger. In aqueous electrolytes, only three ion-polymorph pairs are thermodynamically stable within water’s voltage stability window (Na + and K + in α-MnO 2 , and Li + in λ-MnO 2 ), with all other ion insertion being metastable. We find Al 3+ may insert into the δ, R, and λ polymorphs across the full 2-electron redox of MnO 2 at high voltage; however, electrolytes for multivalent ions must be designed to impede the formation of insoluble precipitates and facilitate cation desolvation. We also show that small ions coinsert with water in α-MnO 2 to achieve greater coordination by oxygen, while solvation energies and kinetic effects dictate water coinsertion in δ-MnO 2 . Finally, taken together, these findings explain reports of mixed ion insertion mechanisms in aqueous electrolytes and highlight promising design strategies for safe, high energy density electrochemical energy storage, desalination batteries, and electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗