Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling coefficient”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Interdiffusion in Ternary Magnesium Solid Solutions of Aluminum and Zinc

Al and Zn are two of the most common alloying elements in commercial Mg alloys, which can improve the physical properties through solid solution strengthening and precipitation hardening. Diffusion plays a key role in the kinetics of these and other microstructural design relevant to Mg-alloy development. However, there is a lack of multicomponent diffusion data available for Mg alloys. Through solid-to-solid diffusion couples, diffusional interactions of Al and Zn in ternary Mg solid-solution at 400° and 450 °C were examined by an extension of the Boltzmann-Matano analysis based on Onsager s formalism. Concentration profiles of Mg-Al-Zn ternary alloys were determined by electron probe microanalysis, and analyzed to determine the ternary interdiffusion coefficients as a function of composition. Zn was determined to interdiffuse the fastest, followed by Mg and Al. Appreciable diffusional interactions among Mg, Al, and Zn were observed by variations in sign and magnitude of cross interdiffusion coefficients. In particular, Zn was found to significantly influence the interdiffusion of Mg and Al significantly: the and ternary cross interdiffusion coefficients were both negative, and large in magnitude, in comparison to and , respectively. Al and Mg were observed influence the interdiffusion of Mg and Al, respectively, with positive and interdiffusion coefficients, but their influence on the Zn interdiffusion was negligible.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗
DOE OSTI2015

Collisional excitation of NH(X{sup 3}Σ{sup −}) by Ne: Potential energy surface, scattering calculations, and comparison with experiments

We present a new three-dimensional potential energy surface (PES) for the NH(X{sup 3}Σ{sup −})–Ne van der Waals system, which explicitly takes into account the NH vibrational motion. Ab initio calculations of the NH–Ne PES were carried out using the open-shell single- and double-excitation coupled cluster approach with non-iterative perturbational treatment of triple excitations [RCCSD(T)]. The augmented correlation-consistent quadruple zeta (aug-cc-pVQZ) basis set was employed. Mid-bond functions were also included in order to improve the accuracy in the van der Waals well. Using this new PES, we have studied the collisional excitation of NH(X{sup 3}Σ{sup −}) by Ne. Close-coupling calculations of the collisional excitation cross sections of the fine-structure levels of NH by Ne are performed for energies up to 3000 cm{sup −1}, which yield, after thermal average, rate coefficients up to 350 K. The propensity rules between fine-structure levels are reported, and it is found that F-conserving cross sections are larger than F-changing cross sections even if the propensity rules are not as strong as for the NH–He system. The calculated rate coefficients are compared with available experimental measurements at room temperature and a fairly good agreement is found between experimental and theoretical data, confirming the good quality of the scattering calculations and also the accuracy of the potential energy surface used in this work.

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

Effective Field Theory islands from perturbative and nonperturbative four-graviton amplitudes

Theoretical data obtained from physically sensible field and string theory models suggest that gravitational Effective Field Theories (EFTs) live on islands that are tiny compared to current general bounds determined from unitarity, causality, crossing symmetry, and a good high-energy behavior. In this work, we present explicit perturbative and nonperturbative 2 → 2 graviton scattering amplitudes and their associated low-energy expansion in spacetime dimensions D ≥ 4 to support this notion. Our new results include a first example of gravity weakly coupled to a nonperturbative effective action. We show that, at energies below the mass of its nonperturbative matter, the D = 4, $\mathcal{N}$ = 1 supersymmetric field theory in the confined phase lies on the same islands identified using four-dimensional perturbative models based on string theory and minimally-coupled matter circulating a loop. Furthermore, we generalize the previous four-dimensional perturbative models based on string theory and minimally-coupled massive spin-0 and spin-1 states circulating in the loop to D dimensions. Remarkably, we again find that the low-energy EFT coefficients lie on small islands. These results offer a useful guide towards constraining possible extensions of Einstein gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Global anomalies and effective field theory

Here, we show that matching anomalies under large gauge transformations and large diffeomorphisms can explain the appearance and non-renormalization of couplings in effective field theory. We focus on thermal effective field theory, where we argue that the appearance of certain unusual Chern-Simons couplings is a consequence of global anomalies. As an example, we show that a mixed global anomaly in four dimensions fixes the chiral vortical effect coefficient (up to an overall additive factor). This is an experimentally measurable prediction from a global anomaly. For certain situations, we propose a simpler method for calculating global anomalies which uses correlation functions rather than eta invariants.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

4-Clique network minor embedding for quantum annealers

Quantum annealing is a quantum algorithm for computing solutions to combinatorial optimization problems. This study proposes a method for minor embedding optimization problems onto sparse quantum annealing hardware graphs called 4-clique network minor embedding. This method is in contrast to the standard minor embedding technique of using a path of linearly connected qubits in order to represent a logical variable state. The 4-clique minor embedding is possible on Pegasus graph connectivity, which is the native hardware graph for some of the current D-Wave quantum annealers. The Pegasus hardware graph contains many cliques of size 4, making it possible to form a graph composed entirely of paths of connected 4-cliques on which a problem can be minor-embedded. The 4-clique chains come at the cost of additional qubit usage on the hardware graph, but they allow for stronger coupling within each chain, thereby increasing chain integrity, reducing chain breaks, and allow for greater usage of the available energy scale for programming logical problem coefficients on current quantum annealers. The 4-clique minor embedding technique is compared with the standard linear path minor embedding with experiments on two D-Wave quantum annealing processors with Pegasus hardware graphs. We show proof-of-concept experiments where the 4-clique minor embeddings can use weak chain strengths while successfully carrying out the computation of minimizing random all-to-all spin glass problem instances. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-principles studies of the unconventional spin-phonon coupling mediated by spin-orbit coupling in pyrochlore Cd 2 Os 2 O 7

Here, we perform systematic first-principles calculations to investigate the spin-phonon coupling (SPC) mediated by spin-orbit coupling in pyrochlore Cd 2 Os 2 O 7 which has the peculiar all-in-all-out noncollinear magnetic ground state. A large SPC coefficient up to 34 cm –1 is found, showing the potential of using this material in thermal spin devices.We further split the SPC coefficients to contributions from symmetric exchange interaction, Dzyaloshinskii-Moriya interaction (DMI), and the single ion anisotropy. Unlike what was assumed in previous studies, the DMI contributions are rather important for most modes. These results provide insights for the understanding of the SPC in 5d transition metal oxides and useful guidance for the search of excellent spin caloritronic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Development of an Integrated Multidiagnostic to Assess the High-Z Impurity Fluxes in the Metallic Environment of WEST Using IMAS

WEST is an actively cooled, long-pulse tokamak with nearly all plasma-facing components (PFC) made of tungsten. One of the aims of WEST is to study plasma operations with tungsten PFCs in preparation for long-pulse operations on high-Z divertor devices, such as ITER. For long-pulse operation, the high-Z impurity content and transport to the core plasma are critical concerns that require further measurement and interpretation to improve plasma performance and PFC durability. This work details the impurity influxes in WEST during a series of discharges in which the lower hybrid (LH) injected power was incrementally increased. An analysis has been performed of measurements collected from an array of edge diagnostics. Visible spectroscopy was utilized to measure the spectral radiances generated by fuel particles (D) and impurities (W, O, and C) at the divertor and at the antennas with a newly developed spectral peak-fitting tool used to analyze the data in WEST. The scrape-off layer (SOL) plasma conditions (n e and T e ) measured at the divertor target with flush Langmuir probes and near the outer mid-plane (OMP) using reciprocating Langmuir probes (RCPs) are used to evaluate the number of ionizations per photons (S/XB) coefficients required to estimate the impurity fluxes obtained with the collisional-radiative model ColRadPy. Here, the array of edge diagnostics discussed in this work, coupled with SOL plasma modeling tools, represents a multidiagnostic interpretative modeling workflow that will continue to be applied to upcoming experimental campaigns on the WEST experiment to assess high-Z impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A higher order approximate static condensation method for multi-material diffusion problems

The paper studies an approximate static condensation method for the diffusion problem with discontinuous diffusion coefficients. The method allows for a general polygonal mesh which is unfitted to the material interfaces. Moreover, the interfaces can be discontinuous across the mesh edges as typical for numerical reconstructions using the volume or moment-of-fluid methods. We apply a mimetic finite difference method to solve local diffusion problems and use P 1 (mortar) edge elements to couple local problems into the global system. The condensation process and the properties of the resulting algebraic system are discussed. It is demonstrated that the method is second order accurate on smooth solutions and performs well for problems with high contrast in diffusion coefficients. Experiments also show the robustness with respect to position of the material interfaces against the underlying mesh.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Development of a water source heat pump hardware-in-the-loop (HIL) testing facility for smart building applications

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while maintaining the thermal comfort of building occupants. As the pivotalpoint between building electricity consumption and indoor thermal comfort, high-efficiency electrical heatpumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on the heatpump operation andperformance under different weather, building load, and grid requests calls for investigation and verification via experimental tests with actual heat pumps integrated with real-time building and grid responses. This study presents a Water-Source Heat Pump Hardware-in-The-Loop (HIL) Test Facility developed with the capability to perform such tests. The hardware configuration for this testfacility introduces a hydronic system that emulates the conditions for the heat pump water-side, and a duct system that emulates conditions for the heat pump airside. Both data acquisition and emulator control are implemented through the National Instruments (NI) LabVIEW software running on an NI PXIplatform. The HIL mechanism based on the hardware-software integration that allows the testbed to communicate with a generic simulation environment is also discussed. Currently, the test facility setup includes a single heatpump and virtual building model in EnergyPlus coupled with an occupant behavioral model in MATLAB. Preliminary test results of the current setup demonstrate the building load emulator's ability to track the simulated gone temperature with a Root Mean Square Deviation (RSME) below 0.12°C (0.216°F). An uncertainty analysis based on sensor accuracies shows that the heat pump coefficient of performance (COP) can be measured with a relative uncertainty of 10.4% in cooling and 3.7% in heating. Apartfrom the current testing on a single heat pump, the test facility also provides the flexibility to include additional heat pumps to form a heat pump cluster, as well as coupling the heat pump with active thermal storage to provide enhanced demandflexibility.

Calfa, Caleb↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

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

The influence of respiratory motion on CT image volume definition

Purpose: Radiotherapy treatments are based on geometric and density information acquired from patient CT scans. It is well established that breathing motion during scan acquisition induces motion artifacts in CT images, which can alter the size, shape, and density of a patient's anatomy. The aim of this work is to examine and evaluate the impact of breathing motion on multislice CT imaging with respiratory synchronization (4DCT) and without it (3DCT). Methods: A specific phantom with a movable insert was used. Static and dynamic phantom acquisitions were obtained with a multislice CT. Four sinusoidal breath patterns were simulated to move known geometric structures longitudinally. Respiratory synchronized acquisitions (4DCT) were performed to generate images during inhale, intermediate, and exhale phases using prospective and retrospective techniques. Static phantom data were acquired in helical and sequential mode to define a baseline for each type of respiratory 4DCT technique. Taking into account the fact that respiratory 4DCT is not always available, 3DCT helical image studies were also acquired for several CT rotation periods. To study breath and acquisition coupling when respiratory 4DCT was not performed, the beginning of the CT image acquisition was matched with inhale, intermediate, or exhale respiratory phases, for each breath pattern. Other coupling scenarios were evaluated by simulating different phantom and CT acquisition parameters. Motion induced variations in shape and density were quantified by automatic threshold volume generation and Dice similarity coefficient calculation. The structure mass center positions were also determined to make a comparison with their theoretical expected position. Results: 4DCT acquisitions provided volume and position accuracies within ±3% and ±2 mm for structure dimensions >2 cm, breath amplitude ≤15 mm, and breath period ≥3 s. The smallest object (1 cm diameter) exceeded 5% volume variation for the breath patterns of higher frequency and amplitude motion. Larger volume differences (>10%) and inconsistencies between the relative positions of objects were detected in image studies acquired without respiratory control. Increasing the 3DCT rotation period caused a higher distortion in structures without obtaining their envelope. Simulated data showed that the slice acquisition time should be at least twice the breath period to average object movement. Conclusions: Respiratory 4DCT images provide accurate volume and position of organs affected by breath motion detecting higher volume discrepancies as amplitude length or breath frequency are increased. For 3DCT acquisitions, a CT should be considered slow enough to include lesion envelope as long as the slice acquisition time exceeds twice the breathing period. If this requirement cannot be satisfied, a fast CT (along with breath-hold inhale and exhale CTs to estimate roughly the ITV) is recommended in order to minimize structure distortion. Even with an awareness of a patient's respiratory cycle, its coupling with 3DCT acquisition cannot be predicted since patient anatomy is not accurately known.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Electronic and optical properties of mixed perovskites CsSn{sub x}Pb{sub (1−x)}I{sub 3}

The electronic structure and optical properties of the new solar cells absorber material: mixed perovskites CsSn{sub x}Pb{sub (1−x)}I{sub 3} are studied by the first-principle calculations with mBJ + SOC (modified Beak Johnson approximation plus spin-orbit coupling) method. The band gap of the serial of compounds almost quasi-linearly reduces with increasing Sn content from 0.96 eV (x = 0) to 0.16 eV (x = 1). Optical absorption coefficient revealed a progressive red shift with the increment of the Sn content, accompanying with the absorption edge broadening. The absorption coefficient and Ideal Power Absorption Coefficient (IPAC) increase greatly with the Pb atoms being partially substituted by Sn atoms. The pure CsSnI{sub 3} has the highest IPAC, but it is unstable in the air because the Sn{sup 2+} will be oxidized to Sn{sup 4+}. So our results indicate that partially substituted CsSn{sub x}Pb{sub (1−x)}I{sub 3} might be the good solar cell absorption material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, Part II: A liquid LBE loop under a high vacuum

To provide US domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. Part I of this work discusses the design of a liquid lead–bismuth eutectic (LBE) windowless target for an accelerator-based neutron source development. Part II discusses how to couple this windowless target to an accelerator operating at an ultra-high vacuum and the subcritical UTA cooled by water at room temperature. Due to the windowless design of the target, the liquid LBE flow shares an ultra-high vacuum (<1.3 × 10 -7 Pa or 10 -9 Torr) space with the accelerator. As a result of this shared vacuum space, the LBE system must operate at a high vacuum (10 -3 ~10 -6 Pa or 10 -5 ~10 -8 Torr). A magnetic rotary motion feedthrough unit utilizes magnetic fluid to allow rotation of the pump while maintaining a high vacuum environment. Prior to testing the LBE system under vacuum, a pump curve measurement is performed to estimate flowrate in the system. This measurement also generates data on orifice loss coefficients, which are compared to correlations in literature. The second experiment investigates vacuum level in the LBE system during operation. High vacuum is maintained (10 -3 ~10 -5 Pa or 10 -5 ~10 -7 Torr) during system operation, and a residual gas analyzer (RGA) scan shows that partial pressures of residual gases in the LBE system lower over the duration of LBE system operations. The third experiment investigates the gravity driven liquid LBE flowing out of the target chamber in the return line, which is partially full. If the liquid LBE is not drained quickly enough, flooding in the target chamber could occur. The coefficient n in the Manning equation is found to be around 0.008 s/m 1/3 . The last experiment performed is a demonstration that a vacuum jacket could provide sufficient thermal insulation to allow coupling between 300 °C LBE loop and a water tank at room temperature. In conclusion, the results from these experiments have influenced the development of the neutron source for the future commercial scale Mo-99 production system.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

High-Temperature Thermodynamics of Cerium Silicates, A-Ce 2 Si 2 O 7 , and Ce 4.67 (SiO 4 ) 3 O

Lanthanide disilicates and oxyapatites have potential roles in high temperature applications as thermal (TBC) and environmental barrier coatings (EBC), or possible alteration phases in geological nuclear waste repositories. However, those Ce 3+ -bearing silicates have only been limitedly studied. In this work, we performed detailed structural and thermodynamic investigations on A-Ce 2 Si 2 O 7 (tetragonal, P4 1 ) and Ce 4.67 (SiO 4 ) 3 O (hexagonal, P6 3 /m). The high temperature structural behaviors and coefficients of thermal expansion were determined by in situ high temperature synchrotron X-ray diffraction (HT-XRD) implemented with Rietveld analysis and thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC). A-Ce 2 Si 2 O 7 was found to be stable in N 2 and air up to ~1483 K with anisotropic thermal expansion along the a and c axes (α a = 12.3 × 10 -6 K -1 and α c = 12.4 × 10 -6 K -1 ). Ce 4.67 (SiO 4 ) 3 O had a slow partial oxidation between 533 K and 873 K to a new nonstoichiometric phase Ce 3+ 1.67- xCe4 + x Ce 3+ 3 (SiO 4 ) 3 O 1+0.5x , followed by a thermal decomposition to CeO 2 and SiO 2 at ~1000 K in air. By using high temperature oxide melt solution calorimetry at 973 K with lead borate as the solvent, the standard enthalpy of formation was determined for A-Ce 2 Si 2 O 7 (-3825.1 ± 6.0 kJ/mol) and Ce 4.67 (SiO 4 ) 3 O (-7391.3 ± 9.5 kJ/mol). Finally, these thermodynamic parameters were compared with those of CeO 2 , CeSiO 4 , and other silicate oxyapatites for examining their chemical stability in high temperature environments relevant for aeronautical applications, mineral formation, and nuclear fuel cycle.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Chemical kinetics and relaxation of non-equilibrium air plasma generated by energetic photon and electron beams

The comprehension of electromagnetic perturbations of electronic devices, due to air plasma-induced electromagnetic field, requires a thorough study on air plasma. In the aim to understand the phenomena at the origin of the formation of non-equilibrium air plasma, we simulate, using a volume average chemical kinetics model (0D model), the time evolution of a non-equilibrium air plasma generated by an energetic X-ray flash. The simulation is undertaken in synthetic air (80% N{sub 2} and 20% O{sub 2}) at ambient temperature and atmospheric pressure. When the X-ray flash crosses the gas, non-relativistic Compton electrons (low energy) and a relativistic Compton electron beam (high energy) are simultaneously generated and interact with the gas. The considered chemical kinetics scheme involves 26 influent species (electrons, positive ions, negative ions, and neutral atoms and molecules in their ground or metastable excited states) reacting following 164 selected reactions. The kinetics model describing the plasma chemistry was coupled to the conservation equation of the electron mean energy, in order to calculate at each time step of the non-equilibrium plasma evolution, the coefficients of reactions involving electrons while the energy of the heavy species (positive and negative ions and neutral atoms and molecules) is assumed remaining close to ambient temperature. It has been shown that it is the relativistic Compton electron beam directly created by the X-ray flash which is mainly responsible for the non-equilibrium plasma formation. Indeed, the low energy electrons (i.e., the non-relativistic ones) directly ejected from molecules by Compton collisions contribute to less than 1% on the creation of electrons in the plasma. In our simulation conditions, a non-equilibrium plasma with a low electron mean energy close to 1 eV and a concentration of charged species close to 10{sup 13 }cm{sup −3} is formed a few nanoseconds after the peak of X-ray flash intensity. 200 ns after the flash application, the electrons are thermalized and their concentration has decreased from about 10{sup 13 }cm{sup −3} down to about 10{sup 12 }cm{sup −3} leaving positive and negative ionised species and atomic radicals whose recombination characteristic times are much longer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗