Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗

Artificial Gravity as a Bone Loss Countermeasure in Simulated Weightlessness

The impact of microgravity on the human body is a significant concern for space travelers. We report here initial results from a pilot study designed to explore the utility of artificial gravity (AG) as a countermeasure to the effects of microgravity, specifically to bone loss. After an initial phase of adaptation and testing, 15 male subjects underwent 21 days of 6 head-down bed rest to simulate the deconditioning associated with space flight. Eight of the subjects underwent 1 h of centrifugation (AG, 1 gz at the heart, 2.5 gz at the feet) each day for 21 days, while 7 of the subjects served as untreated controls (CN). Blood and urine were collected before, during, and after bed rest for bone marker determinations. At this point, preliminary data are available on the first 8 subjects (6 AG, and 2 CN). Comparing the last week of bed rest to before bed rest, urinary excretion of the bone resorption marker n-telopeptide increased 95 plus or minus 59% (mean plus or minus SD) in CN but only 32 plus or minus 26% in the AG group. Similar results were found for another resorption marker, helical peptide (increased 57 plus or minus 0% and 35 plus or minus 13% in CN and AG respectively). Bone-specific alkaline phosphatase, a bone formation marker, did not change during bed rest. At this point, sample analyses are continuing, including calcium tracer kinetic studies. These initial data demonstrate the potential effectiveness of short-radius, intermittent AG as a countermeasure to the bone deconditioning that occurs during bed rest.

Smith, S. M.↗

Small Scale Response and Modeling of Periodically Forced Turbulence

The response of the small scales of isotropic turbulence to periodic large scale forcing is studied using two-point closures. The frequency response of the turbulent kinetic energy and dissipation rate, and the phase shifts between production, energy and dissipation are determined as functions of Reynolds number. It is observed that the amplitude and phase of the dissipation exhibit nontrivial frequency and Reynolds number dependence that reveals a filtering effect of the energy cascade. Perturbation analysis is applied to understand this behavior which is shown to depend on distant interactions between widely separated scales of motion. Finally, the extent to which finite dimensional models (standard two-equation models and various generalizations) can reproduce the observed behavior is discussed.

Bos, Wouter↗

Capture efficiency and bias from the defect dynamics near grain boundaries in BCC Fe using mesoscale simulations

The capture efficiency describes the capability of a sink, such as a grain boundary (GB), dislocation, and void, to absorb point defects (PDs). The bias defines the difference in capture efficiency between the absorption of a vacancy and dumbbell at a sink. Complete kinetic information on PDs, including diffusion barriers and diffusion orientations, as well as accurate saddle points, are needed to determine the capture efficiency and bias at a sink accurately, which is computationally demanding. Here, in the present study, the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method was used to investigate the defect dynamics of PDs near different types of grain boundaries (GBs) (with both $\langle$100$\rangle$ and $\langle$110$\rangle$ families) accurately in body-centered cubic (BCC) iron (Fe). The capture efficiency, sink strength, and bias factor of different types of GBs were determined in Fe, which, different from traditional rate theory estimation, showed a distinct capture efficiency, sink strength, and bias in different GBs. The results demonstrate a strong positive correlation between the capture efficiency and the GB strain width, instead of the GB misorientation, GB energy, or GB-PD binding energy, which have been investigated previously. This work provides valuable insight into the radiation-induced microstructural evolution of GBs.

36 MATERIALS SCIENCE↗

Differential roles of kinetic on- and off-rates in T-cell receptor signal integration revealed with a modified Fab’-DNA ligand

Antibody-derived T-cell receptor (TCR) agonists are commonly used to activate T cells. While antibodies can trigger TCRs regardless of clonotype, they bypass native T cell signal integration mechanisms that rely on monovalent, membrane-associated, and relatively weakly binding ligand in the context of cellular adhesion. Commonly used antibodies and their derivatives bind much more strongly than native peptide major histocompatibility complex (pMHC) ligands bind their cognate TCRs. Because ligand dwell time is a critical parameter that tightly correlates with physiological function of the TCR signaling system, there is a general need, both in research and therapeutics, for universal TCR ligands with controlled kinetic binding parameters. To this end, we have introduced point mutations into recombinantly expressed α-TCRβ H57 Fab to modulate the dwell time of monovalent Fab binding to TCR. When tethered to a supported lipid bilayer via DNA complementation, these monovalent Fab’-DNA ligands activate T cells with potencies well-correlated with their TCR binding dwell time. Single-molecule tracking studies in live T cells reveal that individual binding events between Fab'-DNA ligands and TCRs elicit local signaling responses closely resembling native pMHC. The unique combination of high on- and off-rates of the H57 R97L mutant enables direct observations of cooperative interplay between ligand binding and TCR-proximal condensation of the linker for activation of T cells, which is not readily visualized with pMHC. This work provides insights into how T cells integrate kinetic information from TCR ligands and introduces a method to develop affinity panels for polyclonal T cells, such as cells from a human patient.

Science & Technology - Other Topics↗

Asymmetric ether solvents for high-rate lithium metal batteries

Recent electrolyte solvent design based on weakening lithium-ion solvation have shown promise in enhancing cycling performance of Li-metal batteries. However, they often face slow redox kinetics and poor cycling reversibility at high rate. Here we report using asymmetric solvent molecules substantially accelerates Li redox kinetics. Asymmetric ethers (1-ethoxy-2-methoxyethane, 1-methoxy-2-propoxyethane) showed higher exchange current densities and enhanced high-rate Li 0 plating/stripping reversibility compared to symmetric ethers. Adjusting fluorination levels further improved oxidative stability and Li 0 reversibility. The asymmetric 1-(2,2,2-trifluoro)-ethoxy-2-methoxyethane, with 2 M lithium bis(fluorosulfonyl)imide, exhibited high exchange current density, oxidative stability, compact solid–electrolyte interphase (~10 nm). This electrolyte exhibited superior performance among state-of-the-art electrolytes, enabling over 220 cycles in high-rate Li (50 μm)||LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811, 4.9 mAh cm −2 ) cells and for the first time over 600 cycles in anode-free Cu | |Ni95 pouch cells (200 mAh) under electric vertical take-off and landing cycling protocols. Our findings on asymmetric molecular design strategy points to a new pathway towards achieving fast redox kinetics for high-power Li-metal batteries.

batteries↗

Boundary layer description of directional polymer crystallisation

Nearly fifty years ago Lovinger and Gryte suggested that the directional crystallization of a polymer was analogous to the quiescent isothermal crystallization experiment but at a supercooling where the crystal growth velocity was equal to the velocity of the moving front. Our experiments showed that this equivalence holds in a detailed manner at low directional velocities. To understand the underlying physics of these situations, we modeled the motion of a crystallization front in a liquid where the left side boundary is suddenly lowered below the melting point (Stefan's problem) but with the modification that the crystallization kinetics follow a version of the Avrami model. Our numerical results surprisingly showed that the results of the polymer analog track with the Stefan results which were derived for a simple liquid that crystallizes completely at its melting point; in particular, the position of the crystal growth-front evolved with time exactly as in the Stefan problem. The numerical solution also showed that the temperature in the immediate vicinity of the growth-front decreased with increasing front velocity, which is in line with Lovinger and Gryte's ansatz. To provide a clear theoretical understanding of these numerical results we derive a boundary layer solution to the governing coupled differential equations of the polymer problem. The analytical results are in agreement with our observations from experiments and numerical computations but show that this equivalence between the small molecule and polymer analog only holds in the limit where the crystallization enthalpy is much larger than the rate at which heat is conducted away in the polymer. In particular, in the context of the temperature profile, the enthalpy generated by the crystallisation process which is spread out over a narrow spatial region can be approximated as a point source whose location and temperature correspond to the Lovinger–Gryte ansatz.

36 MATERIALS SCIENCE↗

Tuning Point Defects by Elastic Strain Modulates Nanoparticle Exsolution on Perovskite Oxides

Exsolution generates stable and catalytically active metal nanoparticles via phase precipitation out of a host oxide. An ability to control the size and dispersion of the exsolution particles is desirable for design of nanostructured (electro)catalysts. In this work, we demonstrate that tuning point defects by lattice strain affects both the thermodynamics and the kinetics of iron (Fe 0 ) exsolution on La 0.6 Sr 0.4 FeO 3 (LSF) thin film model. By combining in situ surface characterization and ab initio defect modeling, we show oxygen vacancy and Schottky defects to be the primary point defects formed upon Fe0 exsolution. Lattice strain tunes the formation energy, and thus the abundance of these defects, and alters the amount and size of the resulting exsolution particles. In addition, we find that the density of exsolved nanoparticles matches the concentration of oxygen vacancy pairs, thus pointing to the surface oxygen vacancy pairs as preferential nucleation sites for exsolution. The tensile-strained LSF with a facile formation of these critical point defects results in a higher Fe 0 metal concentration, a larger density of nanoparticles, and a reduced particle size at its surfaces. These results provide important mechanistic insights and highlight the role of point-defect engineering in designing nanostructured catalysts in energy and fuel conversion technologies.

36 MATERIALS SCIENCE↗

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope effects on the structural transformation and relaxation of deeply supercooled water

Herein, we have examined the structure of supercooled liquid D 2 O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D 2 O agree remarkably well with those for H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗

Detailed Chemical Kinetic Modeling of Hydrazine Decomposition

The purpose of this research project is to develop and validate a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. Hydrazine is used extensively in aerospace propulsion, and although liquid hydrazine is not considered detonable, many fuel handling systems create multiphase mixtures of fuels and fuel vapors during their operation. Therefore, a thorough knowledge of the decomposition chemistry of hydrazine under a variety of conditions can be of value in assessing potential operational hazards in hydrazine fuel systems. To gain such knowledge, a reasonable starting point is the development and validation of a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. A reasonably complete mechanism was published in 1996, however, many of the elementary steps included had outdated rate expressions and a thorough investigation of the behavior of the mechanism under a variety of conditions was not presented. The current work has included substantial revision of the previously published mechanism, along with a more extensive examination of the decomposition behavior of hydrazine. An attempt to validate the mechanism against the limited experimental data available has been made and was moderately successful. Further computational and experimental research into the chemistry of this fuel needs to be completed.

Meagher, Nancy E.↗

Effects of Gas-Phase Radiation and Detailed Kinetics on the Burning and Extinction of a Solid Fuel

This is the first attempt to analyze both radiation and detailed kinetics on the burning and extinction of a solid fuel in a stagnation-point diffusion flame. We present a detailed and comparatively accurate computational model of a solid fuel flame along with a quantitative study of the kinetics mechanism, radiation interactions, and the extinction limits of the flame. A detailed kinetics model for the burning of solid trioxane (a trimer of formaldehyde) is coupled with a narrowband radiation model, with carbon dioxide, carbon monoxide, and water vapor as the gas-phase participating media. The solution of the solid trioxane diffusion flame over the flammable regime is presented in some detail, as this is the first solution of a heterogeneous trioxane flame. We identify high-temperature and low-temperature reaction paths for the heterogeneous trioxane flame. We then compare the adiabatic solution to solutions that include Surface radiation only and gas-phase and surface radiation using a black surface model. The analysis includes discussion of detailed flame chemistry over the flammable regime and, in particular, at the low stretch extinction limit. We emphasize the low stretch regime of the radiatively participating flame, since this is the region representative of microgravity flames. When only surface radiation is included, two extinction limits exist (the blow-off limit, and the low stretch radiative limit), and the burning rate and maximum flame temperatures are lower, as expected. With the inclusion of surface and gas-phase radiation, results show that, while flame temperatures are lower, the burning rate of the trioxane diffusion flame may actually increase at low stretch rate due to radiative feedback from the flame to the surface.

Rhatigan, Jennifer L.↗

Critical analysis of nitramine decomposition data: Activation energies and frequency factors for HMX and RDX decomposition

A summary of a literature review on thermal decomposition of HMX and RDX is presented. The decomposition apparently fits first order kinetics. Recommended values for Arrhenius parameters for HMX and RDX decomposition in the gaseous and liquid phases and for decomposition of RDX in solution in TNT are given. The apparent importance of autocatalysis is pointed out, as are some possible complications that may be encountered in interpreting extending or extrapolating kinetic data for these compounds from measurements carried out below their melting points to the higher temperatures and pressure characteristic of combustion.

Schroeder, M. A.↗

1D Simulation of Avalanche to Streamer to Spark Transition of Plasma Discharge in Ammonia-Air Combustion

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of ..alpha.., the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek's criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗

Analysis of two-equation turbulence models for recirculating flows

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Thangam, S.↗

A parametric analysis of radiative structure in aerobrake shock layers

A broad-spectrum version of the NEQAIR code was modified to account for self-absorption and applied to AFE flowfields calculated by the LAURA code with a variety of kinetic models. The resulting radiative fluxes were obtained in a decoupled fashion from the flowfield solver along the vehicle's stagnation streamline. The radiative flux obtained was broken down by causative process to study the radiative structure of the AFE's flowfield for the various kinetic models. In addition, the radiative fluxes for several points on a typical AFE trajectory were analyzed to examine how the radiative structure changes as the vehicle completes its aeropass. Only two radiative processes dominated the stagnation radiative flux, and the flow field conditions near the wal were found to exert considerable influence over the radiative flux to the wall.

Greendyke, Robert B.↗

Termolecular ion-molecule reactions in Titan's atmosphere. IV. A search made at up to 1 micron in pure hydrocarbons

The results of a study of ion-molecule reactions occurring in pure methane, acetylene, ethylene, ethane, propyne, propene, propane, and diacetylene at pressures up to 40 microns of pressure are reported. A variety of experimental methods are used: The standard double resonance in an ICR, for determination of the precursor ions and the modulated double resonance ejection in an ICR, for the determination of the daughter ions. The FA-SIFT technique was used for validation and examination of termolecular reactions with rate coefficients that are less than 10(-26) cm(6) s(-1). An extensive database of reaction kinetics already exists for many of these reactions. The main point of this study was the determination of the accuracy of this database and to search for any missing reactions and reaction channels that may have been omitted from earlier investigations. A specific objective of this work was to extend the study to the highest pressures possible to find out if there were any important termolecular reaction channels occurring. A new approach was used here. In the pure hydrocarbon gases the mass spectra were followed as a function of the pressure changes of the gas. An initial guess was first made using the current literature as a source of the reaction kinetics that were expected. A model of the ion abundances was produced from the solution of the partial differential equations in terms of reaction rate coefficients and initial abundances. The experimental data was fitted to the model for all of the pressures by a least squares minimization to the reaction rate coefficients and initial abundances. The reaction rate coefficients obtained from the model were then compared to the literature values. Several new channels and reactions were discovered when the modeled fits were compared to the actual data. This is all explained in the text and the implications of these results are discussed for the Titan atmosphere.

Hydrocarbons, Acyclic/analysis↗