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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Sunspots and the physics of magnetic flux tubes. I - The general nature of the sunspot. II - Aerodynamic drag

Analysis of the dynamical stability of a large flux tube suggests that the field of a sunspot must divide into many separate tubes within the first 1000 km below the surface. Buoyancy of the Wilson depression at the visible surface and probably also a downdraft beneath the sunspot hold the separate tubes in a loose cluster. Convective generation of Alfven waves, which are emitted preferentially downward, cools the tubes. Aerodynamic drag on a slender flux tube stretched vertically across a convective cell is also studied. Since the drag is approximately proportional to the local kinetic energy density, the density stratification weights the drag in favor of the upper layers. Horizontal motions concentrated in the bottom of the convective cell may reverse this density effect. A downdraft of about two km/sec through the flux tubes beneath the sunspot is hypothesized.

Parker, E. N.↗
NASA NTRS1984

The far reaches of the solar wind - Pioneer 10 and Pioneer 11 plasma results

Selected plasma parameters observed by Pioneer 10 and Pioneer 11 between launch (1972 and 1973) and the end of 1979 are used to find the large-scale radial structure of the solar wind. Comparison of data from the two spacecraft is used to separate temporal from spatial variations. The average bulk speed is found to remain constant at about 430 km/s, with stream structure still evident, though of diminished amplitude, at 20.5 AU (Pioneer 10's distance by the end of 1979). Proton density, flux, pressure, and kinetic energy flux are found to have radial profiles consistent with 1/R-squared. Proton temperatures decrease as R to the -0.6 power, too slowly for an adiabatic expansion.

Kayser, S. E.↗
NASA NTRS1995

A turbulence module for the NPARC code

A turbulence module is developed for the 2D version of the NPARC code which is currently restricted to planar or axisymmetric flows without swirling. Four turbulence models have been built into the module: Baldwin-Lomax, Chien, Shih-Lumley and CMOTT models. The first is a mixing-length eddy-viscosity model which is mainly used for initialization of computational fields and the last three are the low Reynolds number two-equation models. Unlike chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +), an advantage for separated flow calculations. Contrary to the NPARC and most other compressible codes, the non-delta form of transport equations is used which leads to a simpler linearization and is more effective than using the delta form in ensuring the positiveness of the turbulent kinetic energy and its dissipation rate. To reduce numerical diffusion while maintaining necessary stability, a second-order accurate and bounded scheme is used for the convective terms of the turbulent transport equations. This scheme is implemented in a deferred correction manner so that the main coefficients of the resulting difference equations are always positive, thus making the numerical solutions process unconditionally stable. The system of equations is solved via a decoupled method and by the alternating direction TDMA of Thomas. The module can be easily linked to the NPARC code for turbulent flow calculations.

Zhu, J.↗
DOE OSTI2026

Effect of NO on DME-Methanol HCCI Experimental Observations

Methanol is an attractive fuel for the maritime sector due to its wide availability. Its direct use as a fuel, however, is accompanied by challenges such as high latent heat of vaporization and low cetane number. A potential solution to overcome the ignition properties of methanol could be through on-board generation of dimethyl ether (DME) via catalytic dehydration of methanol. The resulting mixture from dehydration can be mixed in with the intake air to generate a homogenous charge compression ignition (HCCI) preburn for subsequent direct injection (DI) and successful ignition methanol at diesel–like timescales. However, if the preburn species are treated separately from the complete methanol MCCI approach the preburn heat release rate (HRR) phasing and behavior do not replicate the preburn behavior of the complete approach. Thus, the presence of the main methanol mixing controlled compression ignition (MCCI) combustion event influences the DME/methanol preburn kinetics. Specifically, it was found that trapped residual temperature alone was insufficient alone to be responsible for the observed differences, and that trace species concentrations of NO in the trapped residual gas also influenced DME/methanol kinetics increasing low temperature heat release (LTHR) magnitude and advancing high temperature heat release (HTHR) phasing. NO is not present in HCCI combustion of neat DME or DME/methanol blends nor is elevated gas temperature in the trapped residuals; both of which result in failure to accurately predict the HCCI combustion of DME/methanol blends when coupled with subsequent DI methanol MCCI. This work experimentally explores the effect of trapped residual temperature and NO on DME and DME/methanol HCCI combustion.

Jatana, Gurneesh [ORNL] (ORCID:0000000288903225)↗
DOE OSTI2020

Kinetic Isotope Effects: Interpretation and Prediction Using Degrees of Rate Control

Kinetic isotope effects (KIEs) have been used for decades in catalysis research as a tool for clarifying reaction mechanisms. Significant primary kinetic isotope effects have usually been interpreted as being a result of isotope substitution at a site of bond breaking (or forming) in the rate-determining step in the reaction mechanism. However, quantitative analysis of the magnitude of the KIE in complex multistep reaction mechanisms is seldom reported. We prove here that the logarithm of the rate ratio for two isotopes is the weighted average over all species in the mechanism of their standard-state free energy difference between the two isotopes, divided by RT. The weighting factor is the degree of rate control (DRC) for that species (e.g., transition state, intermediate, reactant) when the rate is measured separately for each isotope. It is instead the degree of selectivity ratio control (DSRC) when the KIE is measured as the product selectivity in a parallel competition between two isotopes within the same reactant molecule. Since only a few species have non-zero DRCs (or DSRCs) for most reactions, only these few contribute to this weighted average and the KIE. We show that this provides a simple way to interpret and quantitatively predict kinetic isotope effects that is powerful in the insights it provides, allowing one to evaluate directly which species contribute most to the KIE. By applying it to H/D KIEs in several example mechanisms, we further show that the traditional way of interpreting KIEs that focuses only on the rate-determining step can easily lead to misunderstanding of KIE and the reaction mechanism. Finally, this highlights the importance to consider the effect of isotope substitution on the energies of all species with large DRCs (i.e, those whose energies are kinetically relevant). This method also offers opportunities for quantitative validation of mechanism-based microkinetic models.

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

A slip-spring framework to study relaxation dynamics of entangled wormlike micelles with kinetic Monte Carlo algorithm

Hypothesis: Wormlike micelles (WLMs) formed due to the self-assembly of amphiphiles in aqueous solution have similar viscoelastic properties as polymers. Owing to this similarity, in this work, it is postulated that kinetic Monte Carlo (kMC) sampling of slip-springs dynamics, which is able to model the rheology of polymers, can also be extended to capture the relaxation dynamics of WLMs. Theory: The proposed modeling framework considers the following relaxation mechanisms: reptation, union-scission, and constraint release. Specifically, each of these relaxation mechanisms is simulated as separate kMC events that capture the relaxation dynamics while considering the living nature of WLMs within the slip-spring framework. As a case study, the model is implemented to a system of sodium oleate and sodium chloride to predict the linear rheology and the characteristic relaxation times associated with the individual relaxation mechanisms at different pH and salt concentrations. Findings: Linear rheology predictions were found to be in good agreement with experimental data. Furthermore, the calculated relaxation times highlighted that reptation contributed to a continuous increase in viscosity while union-scission contributed to the decrease in viscosity of WLM solutions at a higher salinity and pH. Finally, this manifests the proposed model’s capability to provide insights into the key processes governing WLM’s rheology.

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

Bridging the thermodynamics and kinetics of temperature-induced morphology evolution in polymer/fullerene organic solar cell bulk heterojunction

The performance of organic solar cells (OSC) critically depends on the morphology of the active layer. After deposition, the active layer is in a metastable state and prone to changes that lead to cell degradation. Here, a high efficiency fullerene:polymer blend is used as a model system to follow the temperature-induced morphology evolution through a series of thermal annealing treatments. Electron microscopy analysis of the nano-scale phase evolution during the early stages of thermal annealing revealed that spinodal decomposition, i.e. spontaneous phase separation with no nucleation stage, is possibly responsible for the formation of a fine scale bicontinuous structure. In the later evolution stages, large polycrystalline fullerene aggregates are formed. Optical microscopy and scattering revealed that aggregate-growth follows the Johnson–Mehl–Avrami–Kolmogorov equation indicating a heterogeneous transformation process, i.e., through nucleation and growth. These two mechanisms, spinodal decomposition vs. nucleation and growth, are mutually exclusive and their co-existence is surprising. This unexpected observation is resolved by introducing a metastable monotectic phase diagram and showing that the morphology evolution goes through two distinct and consecutive transformation processes where spinodal decomposition of the amorphous donor:acceptor blend is followed by nucleation and growth of crystalline acceptor aggregates. Lastly, this unified thermodynamic and kinetic mechanism allows us to correlate the morphology evolution with OSC degradation during thermal annealing.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Observations of circumstellar carbon monoxide and evidence for multiple ejections in red giants

Observations of the fundamental 4.6 micron band of CO in nine red giants are presented. A common feature is multiple absorption lines which are identified as products of separate components or shells. Column densities are derived; the relative values should be free of the uncertainties inherent in determining the absolute scale. These column densities are well fitted by single excitation temperatures for each absorption component; these excitation temperatures are identified with the local kinetic temperatures. There is no correlation of CO column density with either gas or dust column density nor of the expansion velocity of the component with its distance from the star. The evidence is reviewed, and it is concluded that mass loss from red giants is most likely episodic in nature.

Bernat, A. P.↗
NASA NTRS2004

Demonstration of Single Axis Combined Attitude Control and Energy Storage Using Two Flywheels

The energy storage and attitude control subsystems of the typical satellite are presently distinct and separate. Energy storage is conventionally provided by batteries, either NiCd or NiH, and active attitude control is accomplished with control moment gyros (CMGs) or reaction wheels. An overall system mass savings can be realized if these two subsystems are combined using multiple flywheels for simultaneous kinetic energy storage and momentum transfer. Several authors have studied the control of the flywheels to accomplish this and have published simulation results showing the feasibility and performance. This paper presents the first experimental results showing combined energy storage and momentum control about a single axis using two flywheels.

Kenny, Barbara H.↗
NASA NTRS2004

Application of Simulated Annealing and Related Algorithms to TWTA Design

Simulated Annealing (SA) is a stochastic optimization algorithm used to search for global minima in complex design surfaces where exhaustive searches are not computationally feasible. The algorithm is derived by simulating the annealing process, whereby a solid is heated to a liquid state and then cooled slowly to reach thermodynamic equilibrium at each temperature. The idea is that atoms in the solid continually bond and re-bond at various quantum energy levels, and with sufficient cooling time they will rearrange at the minimum energy state to form a perfect crystal. The distribution of energy levels is given by the Boltzmann distribution: as temperature drops, the probability of the presence of high-energy bonds decreases. In searching for an optimal design, local minima and discontinuities are often present in a design surface. SA presents a distinct advantage over other optimization algorithms in its ability to escape from these local minima. Just as high-energy atomic configurations are visited in the actual annealing process in order to eventually reach the minimum energy state, in SA highly non-optimal configurations are visited in order to find otherwise inaccessible global minima. The SA algorithm produces a Markov chain of points in the design space at each temperature, with a monotonically decreasing temperature. A random point is started upon, and the objective function is evaluated at that point. A stochastic perturbation is then made to the parameters of the point to arrive at a proposed new point in the design space, at which the objection function is evaluated as well. If the change in objective function values (Delta)E is negative, the proposed new point is accepted. If (Delta)E is positive, the proposed new point is accepted according to the Metropolis criterion: rho((Delta)f) = exp((-Delta)E/T), where T is the temperature for the current Markov chain. The process then repeats for the remainder of the Markov chain, after which the temperature is decremented and the process repeats. Eventually (and hopefully), a near-globally optimal solution is attained as T approaches zero. Several exciting variants of SA have recently emerged, including Discrete-State Simulated Annealing (DSSA) and Simulated Tempering (ST). The DSSA algorithm takes the thermodynamic analogy one step further by categorizing objective function evaluations into discrete states. In doing so, many of the case-specific problems associated with fine-tuning the SA algorithm can be avoided; for example, theoretical approximations for the initial and final temperature can be derived independently of the case. In this manner, DSSA provides a scheme that is more robust with respect to widely differing design surfaces. ST differs from SA in that the temperature T becomes an additional random variable in the optimization. The system is also kept in equilibrium as the temperature changes, as opposed to the system being driven out of equilibrium as temperature changes in SA. ST is designed to overcome obstacles in design surfaces where numerous local minima are separated by high barriers. These algorithms are incorporated into the optimal design of the traveling-wave tube amplifier (TWTA). The area under scrutiny is the collector, in which it would be ideal to use negative potential to decelerate the spent electron beam to zero kinetic energy just as it reaches the collector surface. In reality this is not plausible due to a number of physical limitations, including repulsion and differing levels of kinetic energy among individual electrons. Instead, the collector is designed with multiple stages depressed below ground potential. The design of this multiple-stage collector is the optimization problem of interest. One remaining problem in SA and DSSA is the difficulty in determining when equilibrium has been reached so that the current Markov chain can be terminated. It has been suggested in recent literature that simulating the thermodynamic properties opecific heat, entropy, and internal energy from the Boltzmann distribution can provide good indicators of having reached equilibrium at a certain temperature. These properties are tested for their efficacy and implemented in SA and DSSA code with respect to TWTA collector optimization.

Radke, Eric M.↗
NASA NTRS2004

Modeling of Isotope Fractionation in Stratospheric CO2, N2O, CH4, and O3: Investigations of Stratospheric Chemistry and Transport, Stratosphere-Troposphere Exchange, and Their Influence on Global Isotope Budgets

We investigated the isotopic fractionation of CH4 and hydrogen (H2) in the stratosphere by incorporating isotope-specific rate coefficients into the Lawrence Livermore National Laboratory (LLNL) 2D model and comparing the model results with new observations from the NASA ER-2 aircraft (funded through a separate task under the Upper Atmosphere Research Program). The model results reveal that fractionation which occurs in the stratosphere has a significant influence on isotope compositions in the free troposphere, an important point which had previously been ignored, unrecognized or unquantified for many long-lived trace gases, including CH4 and H2 which we have focused our efforts on to date. Our analyses of the model results and new isotope observations have also been used to test how well the kinetic isotope effects are known, at least to within the uncertainties in model chemistry and transport. Overall, these results represent an important step forward in our understanding of isotope fractionation in the atmosphere and demonstrate that stratospheric isotope fractionation cannot be ignored in modeling studies which use isotope observations in the troposphere to infer the global budgets of CH4 (an important greenhouse gas) and of H2 (a gas whose atmospheric budget must be better quantified, particularly before a large human perturbation from fuel cell use is realized). Our analyses of model results and observations from the NASA ER-2 aircraft are briefly summarized separately below for CH4, H2, and H2O and for the contribution of these modeling studies to date to our understanding of isotope fractionation for N2O, CO2, and O3 as well.

Boering, Kristie A.↗
DOE OSTI2023

3D-Printed Graded Electrode with Ultrahigh MnO 2 Loading for Non-Aqueous Electrochemical Energy Storage

Electrolytic manganese dioxide is one of the promising cathode candidates for electrochemical energy storage devices due to its high redox capacity and ease of synthesis. Yet, high-loading MnO 2 often suffers from sluggish reaction kinetics, especially in non-aqueous electrolytes. The non-uniform deposition of MnO 2 on a porous current collectors also makes it difficult to fully utilize the active materials at high mass loading. Here, a 3D printed graded graphene aerogel (3D GA) that contains sparsely separated exterior ligaments is developed to create large open channels for mass transport as well as densely arranged interior ligaments providing large ion-accessible active surface. The unique structural design homogenizes the thickness of electro deposited MnO 2 even at an ultrahigh mass loading of ≈70 mg cm -2 . The electrode achieves a remarkable volumetric capacity of 29.1 mA h cm -3 in the non-aqueous electrolyte. A Li-ion hybrid capacitor device assembled with a graded 3D GA/MnO2 cathode and graded 3D GA/VO x anode exhibits a wide voltage window of 0–4 V and a superior volumetric energy density of 20.2 W h L -1 . The findings offer guidance on 3D printed electrode design for supporting ultrahigh loading of active materials and developments of high energy density energy storage devices.

electrochemical energy storage↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Weak double layers

The characteristics of weak double layers with potential steps approximately equal to five are investigated in a triple plasma device and compared with a kinetic model. The double layers considered here differ from those previously investigated in two significant ways: (1) the double layers are much weaker; and (2) trapped ions play a much less important role. In this study, two plasmas at different potentials but with similar characteristics are separated from a target chamber by grids so that ions enter only at the high potential side and electron beams enter at the low potential side. What is more, thermal electrons enter at the high potential side and are trapped by the double layer. The model calculations illustrate the role of the trapped ions in determining the position of the double layer, and they demonstrate the dependence of the length of the double layer on the entering ion and electron fluxes.

Hershkowitz, N.↗
NASA NTRS1984

Basic governing equations for the flight regimes of aeroassisted orbital transfer vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are given explicitly. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, is derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, J.-H.↗
NASA NTRS1990

An analysis of the contact sintering process in III-V solar cells

The possibility of eliminating the gross mechanical changes that take place during sintering while retaining the processes that reduce contact resistivity is investigated. Before considering the electrical aspects, however, it is necessary to understand the relevant metallurgy. As the sintering process proceeds in the InP-Au system, dramatic color changes occur that enable a precise measurement of the kinetics of each step. A detailed insight into the processes that control each stage in the InP-Au interaction has resulted. It is shown that contact sintering in the InP-Au system consists of three consecutive stages, each of which is controlled by a separate mechanism.

Weizer, Victor G.↗
NASA NTRS1985

Basic Governing Equations for the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are explicitly given. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, are derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, Jong-Hun↗