Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic methods”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Validity of approximate methods in molecular scattering. III - Effective potential and coupled states approximations for differential and gas kinetic cross sections

Two dimensionality-reducing approximations, the j sub z-conserving coupled states (sometimes called the centrifugal decoupling) method and the effective potential method, were applied to collision calculations of He with CO and with HCl. The coupled states method was found to be sensitive to the interpretation of the centrifugal angular momentum quantum number in the body-fixed frame, but the choice leading to the original McGuire-Kouri expression for the scattering amplitude - and to the simplest formulas - proved to be quite successful in reproducing differential and gas kinetic cross sections. The computationally cheaper effective potential method was much less accurate.

Monchick, L.↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

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

Ground Vibration Test Planning and Pre-Test Analysis for the X-33 Vehicle

This paper describes the results of the modal test planning and the pre-test analysis for the X-33 vehicle. The pre-test analysis included the selection of the target modes, selection of the sensor and shaker locations and the development of an accurate Test Analysis Model (TAM). For target mode selection, four techniques were considered, one based on the Modal Cost technique, one based on Balanced Singular Value technique, a technique known as the Root Sum Squared (RSS) method, and a Modal Kinetic Energy (MKE) approach. For selecting sensor locations, four techniques were also considered; one based on the Weighted Average Kinetic Energy (WAKE), one based on Guyan Reduction (GR), one emphasizing engineering judgment, and one based on an optimum sensor selection technique using Genetic Algorithm (GA) search technique combined with a criteria based on Hankel Singular Values (HSV's). For selecting shaker locations, four techniques were also considered; one based on the Weighted Average Driving Point Residue (WADPR), one based on engineering judgment and accessibility considerations, a frequency response method, and an optimum shaker location selection based on a GA search technique combined with a criteria based on HSV's. To evaluate the effectiveness of the proposed sensor and shaker locations for exciting the target modes, extensive numerical simulations were performed. Multivariate Mode Indicator Function (MMIF) was used to evaluate the effectiveness of each sensor & shaker set with respect to modal parameter identification. Several TAM reduction techniques were considered including, Guyan, IRS, Modal, and Hybrid. Based on a pre-test cross-orthogonality checks using various reduction techniques, a Hybrid TAM reduction technique was selected and was used for all three vehicle fuel level configurations.

Bedrossian, Herand↗
NASA NTRS1987

Wall relaxation in growing stems: comparison of four species and assessment of measurement techniques

This study was carried out to develop improved methods for measuring in-vivo stress relaxation of growing tissues and to compare relaxation in the stems of four different species. When water uptake by growing tissue is prevented, in-vivo stress relaxation occurs because continued wall loosening reduces wall stress and cell turgor pressure. With this procedure one may measure the yield threshold for growth (Y), the turgor pressure in excess of the yield threshold (P-Y), and the physiological wall extensibility (phi). Three relaxation techniques proved useful: "turgor-relaxation", "balance-pressure" and "pressure-block". In the turgor-relaxation method, water is withheld from growing tissue and the reduction in turgor is measured directly with the pressure probe. This technique gives absolute values for P and Y, but requires tissue excision. In the balance-pressure technique, the excised growing region is sealed in a pressure chamber, and the subsequent reduction in water potential is measured as the applied pressure needed to return xylem sap to the cut surface. This method is simple, but only measures (P-Y), not the individual values of P and Y. In the pressure-block technique, the growing tissue is sealed into a pressure chamber, growth is monitored continuously, and just sufficient pressure is applied to the chamber to block growth. The method gives high-resolution kinetics of relaxation and does not require tissue excision, but only measures (P-Y). The three methods gave similar results when applied to the growing stems of pea (Pisum sativum L.), cucumber (Cucumis sativus L.), soybean (Glycine max (L.) Merr.) and zucchini (Curcubita pepo L.) seedlings. Values for (P-Y) averaged between 1.4 and 2.7 bar, depending on species. Yield thresholds averaged between 1.3 and 3.0 bar. Compared with the other methods, relaxation by pressure-block was faster and exhibited dynamic changes in wall-yielding properties. The two pressure-chamber methods were also used to measure the internal water-potential gradient (between the xylem and the epidermis) which drives water uptake for growth. For the four species it was small, between 0.3 and 0.6 bar, and so did not limit growth substantially.

NASA Program Space Biology↗
NASA NTRS1984

Fast algorithms for combustion kinetics calculations: A comparison

To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.

Radhakrishnan, K.↗
DOE OSTI2024

The Effect of Air Separations on Fast Pyrolysis Products for Forest Residue Feedstocks

This study investigates the intricate relationship between biomass preprocessing and pyrolysis product yields, employing the air classification technique for the treatment of loblolly pine residues with varying moisture content. A comprehensive exploration of the physicochemical properties of air-classified loblolly pine informs a sophisticated pyrolysis simulation model. Given the complex and multifaceted nature of biomass pyrolysis, operating across diverse temporal and spatial scales, a pyrolysis kinetics-based CFD–DEM simulation method is employed to predict product yields. Results showed that the elevated moisture content amplifies particle adhesiveness, necessitating augmented air velocities for effective separation, thereby influencing the efficiency of the separation process. While carbon and hydrogen contents exhibit relative stability across diverse moisture contents and blower frequencies, the oxygen content undergoes noticeable changes. For example, the oxygen contents were measured as 29.2 and 38.6 wt% in the light fraction of 30% moisture content sample at blower frequencies of 10 and 20 Hz, respectively. An intriguing finding emerges from pyrolysis simulation, indicating that a lower blower frequency in air classification moderately enhances bio-oil yield and significantly improves its quality, particularly in terms of water content. For instance, the water content in the bio-oil was about 1.5% and 10% in the heavy and light fractions, respectively from 10% moisture sample under 15 Hz blower frequency. In summary, a detailed understanding and strategic manipulation of critical material attributes in biomass through efficient fractionation techniques are imperative for advancing fast pyrolysis as a sustainable avenue for renewable energy and chemical production.

09 BIOMASS FUELS↗
NASA NTRS1998

Experiments on Nucleation in Different Flow Regimes

The vast majority of metallic engineering materials are solidified from the liquid phase. Understanding the solidification process is essential to control microstructure, which in turn, determines the properties of materials. The genesis of solidification is nucleation, where the first stable solid forms from the liquid phase. Nucleation kinetics determine the degree of undercooling and phase selection. As such, it is important to understand nucleation phenomena in order to control solidification or glass formation in metals and alloys. Early experiments in nucleation kinetics were accomplished by droplet dispersion methods. Dilitometry was used by Turnbull and others, and more recently differential thermal analysis and differential scanning calorimetry have been used for kinetic studies. These techniques have enjoyed success; however, there are difficulties with these experiments. Since materials are dispersed in a medium, the character of the emulsion/metal interface affects the nucleation behavior. Statistics are derived from the large number of particles observed in a single experiment, but dispersions have a finite size distribution which adds to the uncertainty of the kinetic determinations. Even though temperature can be controlled quite well before the onset of nucleation, the release of the latent heat of fusion during nucleation of particles complicates the assumption of isothermality during these experiments. Containerless processing has enabled another approach to the study of nucleation kinetics. With levitation techniques it is possible to undercool one sample to nucleation repeatedly in a controlled manner, such that the statistics of the nucleation process can be derived from multiple experiments on a single sample. The authors have fully developed the analysis of nucleation experiments on single samples following the suggestions of Skripov. The advantage of these experiments is that the samples are directly observable. The nucleation temperature can be measured by noncontact optical pyrometry, the mass of the sample is known, and post-processing analysis can be conducted on the sample. The disadvantages are that temperature measurement must have exceptionally high precision, and it is not possible to isolate specific heterogeneous sites as in droplet dispersions. Levitation processing of refractory materials in ultra high vacuum provides an avenue to conduct these kinetic studies on single samples. Two experimental methods have been identified where ultra high vacuum experiments are possible; electrostatic levitation in ground-based experiments and electromagnetic processing in low earth orbit on TEMPUS. Such experiments, reported here, were conducted on zirconium. Liquid zirconium is an excellent solvent and has a high solubility for contaminants contained in the bulk material as well as those contaminants found in the vacuum environment. Oxides, nitrides, and carbides do not exist in the melt, and do not form on the surface of molten zirconium, for the materials and vacuum levels used in this study. Ground-based experiments with electrostatic levitation have shown that the statistical nucleation kinetic experiments are viable and yield results which are consistent with classical nucleation theory. The advantage of low earth orbit experiments is the ability to vary the flow conditions in the liquid prior to nucleation. The put-pose of nucleation experiments in TEMPUS was to examine.

Bayuzick, R. J.↗
DOE OSTI2026

Steric Mapping, Ligand Dynamics, and Cycloisomerization Catalysis with Redox Robust Mn I/0/-I Dicarbenes

Manganese is perhaps the most electronically versatile element, yet the redox properties, reactivity, and catalytic applications of low-valent Mn 0 /Mn −I complexes remain underexplored due to the propensity for Mn 0 to dimerize, quenching highenergy metalloradicals. We report a series of redox-active monometallic Mn I , Mn 0 and Mn −I complexes containing a BH 2 - bridged dicarbene, characterized using a suite of experimental and cutting-edge computational (DFT) methods. Slow electron transfer kinetics at Mn I/0 are observed, with computations and electrochemical simulations in excellent agreement with experimental values. Despite the lack of steric bulk at the BH 2 -bridged Mn 0 , the t Bu groups at the dicarbene provide adequate steric protection to prevent dimerization, with percent buried volume (% V bur ) serving as a valuable steric ranking tool. We also show that a %V bur > 83% prevents dimerization for a diverse array of Mn 0 complexes from the literature. Ligand sterics of BPh 2 -and BH 2 -bridged complexes dictate reaction outcomes when Mn I and Mn −I are exposed to nucleophiles and electrophiles, respectively, while Mn0 facilitates the radical cycloisomerization catalysis of 6-iodo-1- hexene at room temperature. Furthermore, this work underscores the importance of ligand sterics in rationalizing reactivity patterns at Mn and provides valuable insights for designing chelating ligands that can selectively leverage Mn I/0/‑I states in redox-mediated catalytic reactions.

Ligands↗
DOE OSTI2024

Formation of hierarchically structured martensites in pure iron with ultrahigh strength and stiffness

Strong steels are primarily fabricated by introducing spatial obstacles (e.g., stacking faults and precipitates) that inhibit dislocation slips under stress to achieve high strength. However, for most low-carbon steels, such obstacles are difficult to form mainly because the martensitic transition is kinetically unfavorable by conventional methods, which precludes the attainment of high-strength materials in these steels with low solute contents. Here, we report an innovative high-pressure preparation of martensitic pure Fe with involving nano-effect, which leads to the formation of ultrastrong bulk iron with exceptionally high yield strength, ultimate strength, and hardness of 2.9 GPa, 3.7 GPa, and 9.0 GPa, respectively, exceeding those of high-speed steels. Such extraordinary mechanical properties are closely attributed to its high-density martensites with unique multiscale hierarchical structures formed due to complex phase transitions under pressure.

Science & Technology - Other Topics↗
DOE OSTI2025

Lilith: An Enduring Plutonium Critical Assembly [Slides]

A design of an enduring Pu critical assembly has been produced. This is very important to DOE/NNSA missions. This assembly would result in advances/improvements in nuclear data validation, analytical methods validation, dosimetry, reactor kinetics, materials, and more.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1977

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 2: The two-dimensional model

The two-dimensional model of stratospheric constituents is presented in detail. The derivation of pertinent transport parameters and the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations that represent the chemical kinetic terms, and appropriate methods for simulating the diurnal variations of the solar zenith angle and species concentrations are discussed. Predicted distributions of tracer constituents (ozone, carbon 14, nitric acid) are compared with observed distributions.

Whitten, R. C.↗
NASA NTRS1986

Noise-driven sidebranching in the boundary-layer model of dendritic solidification

In the context of local models of solidification, it is suggested that the sidebranching which appears to be a characteristic feature of dendritic processes may be generated by selective amplification of noise near the tip and that, under some circumstances, thermal fluctuations may be strong enough to explain the observed effects. This proposal is discussed based on the two-dimensional boundary-layer model with kinetic crystalline anisotropy. Experimental methods to distinguish this stochastic process from a deterministic, dynamical sidebranching mechanism are suggested, including the possibility of space-based experiments.

Pieters, R.↗
NASA NTRS1991

Molecules and chains in a strong magnetic field - Statistical treatment

A Thomas-Fermi-Dirac-Weizsaecker statistical model is developed and employed to investigate diatomic molecules and infinite molecular chains in strong magnetic fields. The standard magnetic Thomas-Fermi-Dirac kinetic, potential, and exchange energy functionals are supplemented by a gradient correction to the kinetic energy. The numerical method used for solving this system in two spatial dimensions is detailed. Numerical solutions for a wide range of magnetic strengths and elements are presented to demonstrate the robustness, as well as the limitations, of the statistical approach. These calculations qualitatively reproduce many of the results of detailed quantum mechanical treatments. For example, the fractional binding energy is greatest for low atomic numbers and for strong magnetic fields.

Abrahams, Andrew M.↗
NASA NTRS1994

Direct measurement of free-energy barrier to nucleation of crystallites in amorphous silicon thin films

A method is introduced to measure the free-energy barrier W(sup *), the activation energy, and activation entropy to nucleation of crystallites in amorphous solids, independent of the energy barrier to growth. The method allows one to determine the temperature dependence of W(sup *), and the effect of the preparation conditions of the initial amorphous phase, the dopants, and the crystallization methds on W(sup *). The method is applied to determine the free-energy barrier to nucleation of crystallites in amorphous silicon (a-Si) thin films. For thermally induced nucleation in a-Si thin films with annealing temperatures in the range of from 824 to 983 K, the free-energy barrier W(sup *) to nucleation of silicon crystals is about 2.0 - 2.1 eV regardless of the preparation conditions of the films. The observation supports the idea that a-Si transforms into an intermediate amorphous state through the structural relaxation prior to the onset of nucleation of crystallites in a-Si. The observation also indicates that the activation entropy may be an insignificant part of the free-energy barrier for the nucleation of crystallites in a-Si. Compared with the free-energy barrier to nucleation of crystallites in undoped a-Si films, a significant reduction is observed in the free-energy barrier to nucleation in Cu-doped a-Si films. For a-Si under irradiation of Xe(2+) at 10(exp 5) eV, the free-energy barrier to ion-induced nucleation of crystallites is shown to be about half of the value associated with thermal-induced nucleation of crystallites in a-Si under the otherwise same conditions, which is much more significant than previously expected. The present method has a general kinetic basis; it thus should be equally applicable to nucleation of crystallites in any amorphous elemental semiconductors and semiconductor alloys, metallic and polymeric glasses, and to nucleation of crystallites in melts and solutions.

Shi, Frank G.↗
NASA NTRS1995

Calculation of kinetic rate constants from thermodynamic data

A new scheme for relating the absolute value for the kinetic rate constant k to the thermodynamic constant Kp is developed for gases. In this report the forward and reverse rate constants are individually related to the thermodynamic data. The kinetic rate constants computed from thermodynamics compare well with the current kinetic rate constants. This method is self consistent and does not have extensive rules. It is first demonstrated and calibrated by computing the HBr reaction from H2 and Br2. This method then is used on other reactions.

Marek, C. John↗