SEARCH · Engineering Papers
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.
A large number of model particles are shown necessary for quantitatively correct simulations of the kinetic beam-plasma instability with the clouds-in-cells method. The required number of particles scales inversely with the expected growth rate, as only a narrow interval of beam velocities is resonant with the wave in the kinetic regime.
The low-frequency limit of Maxwell equations is considered in the Maxwell-Vlasov system. This limit produces a neutral Vlasov system that captures essential features of plasma dynamics, while neglecting radiation effects. Euler-Poincaré reduction theory is used to show that the neutral Vlasov kinetic theory possesses a variational formulation in both Lagrangian and Eulerian coordinates. By construction, the new model recovers all collisionless neutral models employed in plasma simulations. Then, comparisons between the neutral Vlasov system and hybrid kinetic-fluid models are presented in the linear regime.
Linear gyro-kinetic simulations of the classical tearing mode in three-dimensional toroidal geometry were performed using the global gyro-kinetic turbulence code, GKW. The results were benchmarked against a cylindrical ideal MHD and analytical theory calculations. The stability, growth rate, and frequency of the mode were investigated by varying the current profile, collisionality, and the pressure gradients. Both collisionless and semi-collisional tearing modes were found with a smooth transition between the two. A residual, finite, rotation frequency of the mode even in the absence of a pressure gradient is observed, which is attributed to toroidal finite Larmor-radius effects. When a pressure gradient is present at low collisionality, the mode rotates at the expected electron diamagnetic frequency. However, the island rotation reverses direction at high collisionality. The growth rate is found to follow a η{sup 1∕7} scaling with collisional resistivity in the semi-collisional regime, closely following the semi-collisional scaling found by Fitzpatrick. The stability of the mode closely follows the stability analysis as performed by Hastie et al. using the same current and safety factor profiles but for cylindrical geometry, however, here a modification due to toroidal coupling and pressure effects is seen.
A response to “Comment on ‘The ion-kinetic D'Angelo mode’” by Aman-ur-Rehman, Shaukat Ali Shan, and Hamid Saleem is given. We find that all Comments are grounded on the misinterpretations of the results of our paper “The ion-kinetic D'Angelo mode” and are erroneous.
In this paper, the effect of an external inhomogeneous magnetic field on the high intensity laser absorption rate in a sub-critical plasma has been investigated by employing a relativistic electromagnetic 1.5 dimensional particle-in-cell code. Relying on the effective nonlinear phenomena such as phase-mixing and scattering, this study shows that in a finite-size plasma the laser absorption increases with inhomogeneity of the magnetic field (i.e., reduction of characteristic length of inhomogeneous magnetic field, λ{sub p}) before exiting a considerable amount of laser energy from the plasma due to scattering process. On the other hand, the presence of the external inhomogeneous magnetic field causes the maximum absorption of laser to occur at a shorter time. Moreover, study of the kinetic results associated with the distribution function of plasma particles shows that, in a special range of the plasma density and the characteristic length of inhomogeneous magnetic field, a considerable amount of laser energy is transferred to the particles producing a population of electrons with kinetic energy along the laser direction.
We discuss the complete theory of spin-1/2 electron-positron quantum plasmas, when electrons and positrons move with velocities mach smaller than the speed of light. We derive a set of two fluid quantum hydrodynamic equations consisting of the continuity, Euler, spin (magnetic moment) evolution equations for each species. We explicitly include the Coulomb, spin-spin, Darwin and annihilation interactions. The annihilation interaction is the main topic of the paper. We consider the contribution of the annihilation interaction in the quantum hydrodynamic equations and in the spectrum of waves in magnetized electron-positron plasmas. We consider the propagation of waves parallel and perpendicular to an external magnetic field. We also consider the oblique propagation of longitudinal waves. We derive the set of quantum kinetic equations for electron-positron plasmas with the Darwin and annihilation interactions. We apply the kinetic theory to the linear wave behavior in absence of external fields. We calculate the contribution of the Darwin and annihilation interactions in the Landau damping of the Langmuir waves. We should mention that the annihilation interaction does not change number of particles in the system. It does not related to annihilation itself, but it exists as a result of interaction of an electron-positron pair via conversion of the pair into virtual photon. A pair of the non-linear Schrodinger equations for the electron-positron plasmas including the Darwin and annihilation interactions is derived. Existence of the conserving helicity in electron-positron quantum plasmas of spinning particles with the Darwin and annihilation interactions is demonstrated. We show that the annihilation interaction plays an important role in the quantum electron-positron plasmas giving the contribution of the same magnitude as the spin-spin interaction.
Localized vibrational mode spectroscopy measurements on Czochralski silicon (Cz-Si) samples subjected to isothermal annealing at 450 °C are reported. First, we studied the effect of carbon (C) and tin (Sn) isovalent dopants on the aggregation kinetics of oxygen. It is determined that the reduction rate of oxygen is described by the Johnson-Mehl-Avrami equation in accordance with previous reports. The activation energy related with the reaction rate constant of the process is calculated to increase from Cz-Si, to C-doped Cz-Si (CCz-Si), to Sn-doped Cz-Si contained C (SnCz-Si). This is attributed to the presence of the isovalent dopants that may impact both the kinetics of the oxygen atoms and also may lead to the formation of other oxygen-related clusters. Second, we studied the effect of Sn on the formation and evolution of carbon-oxygen (C-O) defects. It was determined that the presence of Sn suppresses the formation of the C-O defects as indicated by the reduction in the strength of the 683, 626, and 586 cm{sup −1} well-known bands of C{sub s}O{sub i} defect. The phenomenon is attributed to the association of Sn with C atoms that may prevent the pairing of O with C. Third, we investigated the effect of C and Sn on the formation of thermal donors (TDs). Regarding carbon our results verified previous reports that carbon suppresses the formation of TDs. Interestingly, when both C and Sn are present in Si, very weak bands of TDs were observed, although it is known that Sn alone suppress their formation. This may be attributed to the competing strains of C and Sn in the Si lattice.
A quantum kinetic theory of molecule formation is presented which includes three-body recombination and radiative association for a thermodynamically closed system which may or may not exchange energy with its surrounding at a constant temperature. The theory uses a Sturmian representation of a two-body continuum to achieve a steady-state solution of a governing master equation which is self-consistent in the sense that detailed balance between all bound and unbound states is rigorously enforced. The role of quasibound states in catalyzing the molecule formation is analyzed in complete detail. The theory is used to make three predictions which differ from conventional kinetic models. These predictions suggest significant modifications may be needed to phenomenological rate constants which are currently in wide use. Implications for models of low and high density systems are discussed.
In a recent paper, Gaelzer et al. [Phys. Plasmas 22, 032310 (2015)] revisited the second-order nonlinear kinetic theory for turbulence propagating in directions parallel/anti-parallel to the ambient magnetic field. The original work was according to Yoon and Fang [Phys. Plasmas 15, 122312 (2008)], but Gaelzer et al. noted that the terms pertaining to discrete-particle effects in Yoon and Fang's theory did not enjoy proper dimensionality. The purpose of Gaelzer et al. was to restore the dimensional consistency associated with such terms. However, Gaelzer et al. was concerned only with linear wave-particle interaction terms. The present paper completes the analysis by considering the dimensional correction to nonlinear wave-particle interaction terms in the wave kinetic equation.
The correlation energy of interaction is an elusive and sought-after interaction between molecular systems. By partitioning the response function of the system into subsystem contributions, the Frozen Density Embedding (FDE)-vdW method provides a computationally amenable nonlocal correlation functional based on the adiabatic connection fluctuation dissipation theorem applied to subsystem density functional theory. In reproducing potential energy surfaces of weakly interacting dimers, we show that FDE-vdW, either employing semilocal or exact nonadditive kinetic energy functionals, is in quantitative agreement with high-accuracy coupled cluster calculations (overall mean unsigned error of 0.5 kcal/mol). When employing the exact kinetic energy (which we term the Kohn-Sham (KS)-vdW method), the binding energies are generally closer to the benchmark, and the energy surfaces are also smoother.
Reflection electron energy loss spectra from some insulating materials (CaCO{sub 3}, Li{sub 2}CO{sub 3}, and SiO{sub 2}) taken at relatively high incoming electron energies (5–40 keV) are analyzed. Here, one is bulk sensitive and a well-defined onset of inelastic excitations is observed from which one can infer the value of the band gap. An estimate of the band gap was obtained by fitting the spectra with a procedure that includes the recoil shift and recoil broadening affecting these measurements. The width of the elastic peak is directly connected to the mean kinetic energy of the atom in the material (Doppler broadening). The experimentally obtained mean kinetic energies of the O, C, Li, Ca, and Si atoms are compared with the calculated ones, and good agreement is found, especially if the effect of multiple scattering is taken into account. It is demonstrated experimentally that the onset of the inelastic excitation is also affected by Doppler broadening. Aided by this understanding, we can obtain a good fit of the elastic peak and the onset of inelastic excitations. For SiO{sub 2}, good agreement is obtained with the well-established value of the band gap (8.9 eV) only if it is assumed that the intensity near the edge scales as (E − E{sub gap}){sup 1.5}. For CaCO{sub 3}, the band gap obtained here (7 eV) is about 1 eV larger than the previous experimental value, whereas the value for Li{sub 2}CO{sub 3} (7.5 eV) is the first experimental estimate.
Photo induced birefringent materials can be used to develop optical and conversion energy devices, and consequently, the study of the variables that influences the creation and relaxation of birefringence should be carefully analyzed. In this work, the parameters of birefringence creation and relaxation kinetics curves obtained on layer-by-layer (LBL) films, prepared from azo-polyectrolyte poly[1-[4-(3-carboxy-4 hydroxyphenylazo) benzene sulfonamido]-1,2-ethanediyl, sodium salt] (PAZO) and poly(allylamine hydrochloride)(PAH), are related with the presence of counterions and the degree of ionization of the polyelectrolytes. Those kinetics curves obtained on PAH/PAZO LBL films, prepared from PAH solutions with different pHs and maintaining the pH of PAZO solution constant at pH = 9, were analyzed taking into account the films composition which was characterized by X-ray photoelectron spectroscopy. The creation and relaxation birefringence curves are justified by two processes: one associated to local mobility of the azobenzene with a characteristic time 30 s and intensity constant and other associated with polymeric chains mobility with the characteristic time and intensity decreasing with pH. These results allow us to conclude that the birefringence creation process, associated to local mobility of azobenzenes is independent of the degree of ionization and of number of counterions or co-ions present while the birefringence creation process associated to mobility of chains have its characteristic time and intensity dependent of both degree of ionization and number of counterions. The birefringence relaxation processes are dependent of the degree of ionization. The analysis of the films composition revealed, in addition, the presence of a protonated secondary or tertiary amine revealing that PAZO may have positive charges and consequently a zwitterionic behavior.
The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.
In this work, we use GaN with different deposited Pt nanostructures as a controllable model system to investigate the kinetics of photo-generated charge carriers in hybrid photocatalysts. We combine conductance and contact potential difference measurements to investigate the influence of Pt on the processes involved in the capture and decay of photo-generated charge carriers at and close to the GaN surface. We found that in the presence of Pt nanostructures the photo-excitation processes are similar to those found in Pt free GaN. However, in GaN with Pt nanostructures, photo-generated holes are preferentially trapped in surface states of the GaN covered with Pt and/or in electronic states of the Pt and lead to an accumulation of positive charge there, whereas negative charge is accumulated in localized states in a shallow defect band of the GaN covered with Pt. This preferential accumulation of photo-generated electrons close to the surface is responsible for a dramatic acceleration of the turn-off charge transfer kinetics and a stronger dependence of the surface photovoltage on light intensity when compared to a Pt free GaN surface. Our study shows that in hybrid photocatalysts, the metal nanostructures induce a spatially inhomogeneous surface band bending of the semiconductor that promotes a lateral drift of photogenerated charges towards the catalytic nanostructures.
In this work, we deal with general reactive systems involving N species and M elementary reactions under applicability of the mass-action law. Starting from the dynamic variables introduced in two previous works [P. Nicolini and D. Frezzato, J. Chem. Phys. 138(23), 234101 (2013); 138(23), 234102 (2013)], we turn to a new representation in which the system state is specified in a (N × M){sup 2}-dimensional space by a point whose coordinates have physical dimension of inverse-of-time. By adopting hyper-spherical coordinates (a set of dimensionless “angular” variables and a single “radial” one with physical dimension of inverse-of-time) and by examining the properties of their evolution law both formally and numerically on model kinetic schemes, we show that the system evolves towards the equilibrium as being attracted by a sequence of fixed subspaces (one at a time) each associated with a compact domain of the concentration space. Thus, we point out that also for general non-linear kinetics there exist fixed “objects” on the global scale, although they are conceived in such an abstract and extended space. Moreover, we propose a link between the persistence of the belonging of a trajectory to such subspaces and the closeness to the slow manifold which would be perceived by looking at the bundling of the trajectories in the concentration space.
We have found that thermodynamic state and kinetic process co-determine the dual ferromagnetic (FM) orders in high-Si content FeMnP{sub 1−x}Si{sub x} (0.25 < x < 0.5). Alloys undergoing high temperature annealing and quenching process prefer a high magnetic moment FM state in a chemically partial disordered structure with low c/a ratio. This mechanism is suggested to be responsible for the often discussed virgin effect as well. A chemically ordered structure obtained by a slow cooling process from a relatively low annealing temperature and the increase in Si content stabilize a metastable lattice with high c/a ratio and FM order with low magnetic moment. The non-simultaneity of the magnetic and structural transitions can be responsible for the occurrence of FM state in the high c/a range. Thus, a c/a ratio that changes from high to low is physically plausible to stabilize the metastable FM order at low temperature. Our theoretical observations indicate that suitable thermodynamic state and kinetic diffusion process is crucial for optimizing magnetocaloric properties and exploring feasible magnetocaloric materials.
A new high-temperature NMR technique for measurements of the phase transition kinetics in solids has been developed. The technique allows measuring the time evolution of the volume of the appearing phase at controlled cooling rates. Developed method was applied to study the phase transition kinetics in the superionic conductor LiNaSO{sub 4}. It was revealed that the phase transition in LiNaSO{sub 4} is governed by the diffusion-controlled growth of nuclei (“germs”). An effect of the crystallite rearrangement in the LiNaSO{sub 4} powder after cooling through the phase transition was also revealed. This effect was studied by means of high-temperature XRD and NMR.
The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766–1793 (1996); ibid. 56, 1794–1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence.