Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Oxidation kinetic modelling of Fe-based oxygen carriers for chemical looping applications: Impact of the topochemical effect

Chemical looping is a promising technology for fossil fuel utilization due to its high fuel conversion efficiency with in-situ CO 2 capture capability. Metal oxide are used as oxygen carriers (OCs) and circulate between a fuel reactor and an air reactor to perform reduction and oxidation reactions, respectively. In general, OC exiting the fuel reactor is not reduced fully to its metallic state due to many factors including carbon deposition and OC deactivation. Therefore, the effect of the initial reduction state on the OC oxidation in the air reactor is a significant parameter for consideration in developing the oxidation kinetic model. The objective of this work is to develop a physically significant kinetic model that applies to the oxidation of both initially fully and partially reduced OC with air. For this study, 1.5 mm Fe-based OC particle supported with TiO 2 was used as the model OC particle due to its complex multistep reaction nature. The oxidation kinetics were experimentally investigated in a thermogravimetric analyzer (TGA). Results indicate a significant difference in the oxidation rate profile for the OCs when oxidized from an initially fully reduced compared to an initially partially reduced state. Elemental mapping via energy-dispersive X-ray spectroscopy (EDS) reveals a shrinking-core type topochemical pattern across the OC particle, which was identified to be the cause of the dependency of kinetics on the initial reduction state. A generalized kinetic model was developed based on the observed shrinking-core behavior without presuming any rate-determining steps and experimentally validated over a broad range of temperatures (800–1000 °C) and oxygen concentrations (5, 7, 10, and 15 mol%). Impacts of particle porosity, size, and core-shell structure on the OC oxidation kinetics were analyzed in the developed oxidation kinetic model to suggest methods of improving the oxidation rate of the OC without modifying the chemical composition.

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

Approach to Evaluating Reorganization Energies of Interfacial Electrochemical Reactions

Reaction rate coefficients for electron-transfer processes at the electrode–electrolyte interface are commonly estimated by using the Butler–Volmer equation, but their values are inaccurate beyond a few tenths of volts of overpotential. The Marcus–Hush–Chidsey (MHC) formalism yields correct asymptotic behavior of the rate coefficients vs applied overpotential but has complex dependencies on the redox system’s intrinsic parameters, which can be difficult to model or measure. In this work, we bridge the two kinetics formalisms to estimate the reorganization energy, one of the important parameters for the MHC formalism, and investigate its dependence on other intrinsic parameters such as activation barriers, electronic coupling strength, and the density of states of the electrode surface. We examine the sensitivity of the reorganization energy to these parameters, establish some general relationships for accurately predicting rate coefficients using the MHC formalism over a wide range of applied overpotentials, and compare this approach to calculating MHC rate constants with other empirical approaches for the mechanisms of CO 2 reduction on different metal electrode surfaces.

Butler−Volmer↗
DOE OSTI2024

Direct Numerical Simulation of Involute Channel Turbulence

A direct numerical simulation (DNS) study was performed on turbulent flow in the high flux isotope reactor involute channel geometry to develop a numerical database and determine the differences compared with a flat parallel channel. The varying channel curvature along the walls was studied for differences in mean profiles. Parameters of interest include streamwise velocity, turbulent kinetic energy (TKE), and turbulence dissipation rate, as well as Reynolds stresses and turbulence transport terms. Profile sampling was carried out at 10 locations along the span of the involute. Additional DNS studies were performed on smaller domains of comparable curvature to the involute domain: a high curvature channel (high circular), a low curvature channel (low circular), and a flat channel (flat). Here, each of these four cases was compared against each other and to other DNS studies performed on parallel flows. The results indicate that the bulk involute channel flow does not differ significantly from a flat parallel channel flow and that the curvature of the walls does not significantly alter the mean flow parameters. However, the regions of the involute channel near the side walls exhibit relatively low magnitude twin recirculation structures driven toward the side walls from the centerline of the channel, which warrants further study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Energy-dependent optimization of the prompt fission neutron spectrum with CGMF

Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Oxidation kinetics of biochar from woody and herbaceous biomass

The goal of this study was to determine oxidation kinetics for biochar produced from fast pyrolysis of various biomass feedstocks. In particular, the role of inherent ash content on the oxidation rate was evaluated. Thermogravimetric analysis (TGA) and fluidized bed combustion experiments were used to explore oxidation kinetics of six fast pyrolysis produced biochars with diverse ash content. Here, the reaction rates varied by a factor of three under chemical kinetic-limited conditions, demonstrating inorganic content impacts oxidation rate. Chemical kinetic rate coefficients were proposed as a function of compositional parameters to determine overall fit and impact. Potassium content was found to have a positive correlation, best describing the differences in the oxidation kinetic rate coefficients. Additionally, feedstocks were subjected to a 1 M HCl acid wash mitigating the catalytic activity of the metals. Acid washed biochars had lower oxidation kinetic rates compared to their unwashed counterparts, indicating the removal of catalytically active metals reduced oxidation rate. Gas composition (CO/CO 2 ) was measured during fluidized bed experiments for both acid-washed and unwashed biochars, which varied between the six biochars. Formation of CO 2 was greatly affected by catalytic metals, finding potassium content to correlate well with a higher percentage of CO 2 formation as compared to CO. Comparison of oxidation rates were made between the two experimental apparatuses to measure the effect of attrition on biochar oxidation.

09 BIOMASS FUELS↗
DOE OSTI2020

A CALPHAD-informed approach to modeling constituent redistribution in Zr-based metallic fuels using BISON

Here, a CALPHAD-informed (Computer Coupling of Phase Diagrams and Thermochemistry) constituent redistribution model was developed for Zr-based metallic fuels and incorporated into the BISON fuel performance code. Three uncertain model parameters associated with β and γ phase kinetics were calibrated using integral test data from U-Zr fuel elements irradiated in Experimental Breeder Reactor II. The calibrated constituent redistribution model was shown to predict the behavior of U-Zr fuels with excellent accuracy. Model predictions for U-Pu-Zr fuels were physically reasonable but less accurate. Reduction of uncertainties in the ternary phase transition temperatures and collection of kinetic data for the ζ phase are expected to improve the model’s ternary predictions. Finally, the new model was coupled to existing thermomechanics models in BISON to simulate irradiation of an entire U-Zr fuel element, demonstrating its ability to accurately predict the behavior of U-Zr fuels with realistic geometries and mesh resolutions at the engineering scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

In Situ Visualization of Local Distortions in the High- T c Molecule-Intercalated Li x (C 5 H 5 N) y Fe 2– z Se 2 Superconductor

A time-resolved synchrotron X-ray total scattering study sheds light on the evolution of the different structural length scales involved during the intercalation of the layered iron–selenide host by organic molecular donors, aiming at the formation of the expanded-lattice Li x (C 5 H 5 N) y Fe 2– z Se 2 hybrid superconductor. The intercalates are found to crystallize in the tetragonal ThCr 2 Si 2 -type structure at the average level, however, with an enhanced interlayer iron–selenide spacing ( d = 16.2 Å) that accommodates the heterocyclic molecular spacers. Quantitative atomic pair distribution function (PDF) analysis at variable times suggests distorted FeSe 4 tetrahedral local environments that appear swollen with respect to those in the parent β-FeSe. Simultaneously acquired in situ synchrotron X-ray powder diffraction data disclose that secondary phases (α-Fe and Li 2 Se) grow significantly when a higher lithium concentration is used in the solvothermal reaction or when the solution is aged. These observations are in line with the strongly reducing character of the intercalation medium’s solvated electrons that mediate the defect chemistry of the expanded-lattice superconductor. In the latter, intralayer correlated local distortions indicate electron-donating aspects that reflect in somewhat enlarged Fe–Se bonds. Additionally, they also reveal a degree of relief of chemical pressure associated with a large distance between Fe and Se sheets (“taller” anion height) and a stretched Fe–Fe square planar topology. The elongation of the latter, derived from the in situ PDF study, speaks for a plausible increase in the Fe-site vacancy concentration. The evolution of the local structural parameters suggests an optimum reaction window where kinetically stabilized phases resemble the distortions of the edge-sharing Fe–Se tetrahedra, required for a high- T c in expanded-lattice iron-chalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hidden figures of photo-charging: a thermo-electrochemical approach for a solar-rechargeable redox flow cell system

Achieving high current densities without thermal performance degradation at high temperatures is one of the main challenges for enhancing the competitiveness of photo-electrochemical energy storage systems. We describe a system that overcomes this challenge by incorporating an integrated photoelectrode with a redox flow cell, which functions as a coolant for the excess heat from the photo-absorber. We perform quantitative analyses to theoretically validate and highlight the merit of the system. Practical operation parameters, including daily temperature and redox reaction kinetics, are modeled with respect to heat and charge transfer mechanisms. Our analyses show a profound impact on the resulting solar-to-chemical efficiencies and stored power, which are 21.8% higher than that of a conventional photovoltaic-assisted energy storage system. This paves the way for reassessing the merit of photovoltaic-integrated systems, which have hitherto been underrated as renewable energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Two-step spin-coating of vacancy-ordered double perovskites enables growth of thin films for electronic devices

Vacancy-ordered double perovskites (VODPs), such as Cs 2 TeX 6 (X = Cl, Br, I), are lead-free alternatives to conventional metal-halide perovskites (MHPs). One limitation of VODPs is the lack of processes to form thin films relevant for physical characterization and electronic devices. A two-step spin-coating method was developed for synthesizing high-quality films of Cs 2 TeBr 6 . Independently depositing CsBr and TeBr 4 enables high precursor concentrations and control over crystallization kinetics. By optimizing the spin-coating parameters, conversion of precursors to phase pure films was observed using structural and surface characterization methods. The growth of mixed-halide systems was investigated using alternative salts including CsCl and CsI. The formation of halide alloys was found to depend on the existence of routes to byproducts. Lastly, single carrier diodes of Cs2TeBr6 were designed following valence band characterization with photoelectron spectroscopy. Temperature-dependent space-charge-limited current measurements revealed that transport occurs by hopping and the hole mobility is 3.2 × 10 −5 cm 2 V −1 s −1 near room temperature. As a result, the insights from the 2-step procedure provide a pathway towards making semiconducting devices from VODPs.

Kuklinski, Owen [University of California, Santa B↗
DOE OSTI2024

Investigation of pedestal parameters and divertor heat fluxes in small ELM regimes in DIII-D

Abstract Divertor heat flux and its correlation with pedestal parameters within various small edge localized mode (ELM) regimes, including high beta poloidal, type-II and ELMs with negative triangularity H-modes were investigated in DIII-D. The parallel energy fluences of type-II and high beta poloidal small ELM regimes fall below the linear scaling with pedestal electron pressure for type-I ELMs put forward in Eich et al 2017 ( Nucl. Mater. Energy 12 84–90). The negative triangularity of H-mode ELMs follow the Eich scaling for type-I ELMs. The parallel heat flux and total heat loads to the divertor were determined using high-time resolution infrared thermography, while pedestal parameters were obtained through self-consistent kinetic equilibrium reconstructions. Linear regressions for the type-II and high beta poloidal regimes demonstrate that an equivalent 7.5 MA small ELM scenario in ITER would fall below the ~5 MJ m − 2 leading edge melting limit for tungsten (Gunn et al 2017 Nucl. Fusion 57 046025). Utilizing fast thermography, the scrape-off layer power fall-off length for both inter-ELM and intra-ELM was determined and compared to the Eich scaling with poloidal magnetic field in Eich et al (ASDEX Upgrade Team and JET EFDA Contributors 2013 Nucl. Fusion 53 093031). Except for the high beta poloidal scenario, all the small ELM regimes during both inter- and intra-ELM periods had power fall-off lengths ( λ q ) larger then would be expected from the B pol , MP − 1 scaling associated with type-I ELMs, signifying their potential in managing heat loads and offering a solution for core–edge integration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A forecast for large-scale structure constraints on Horndeski gravity with CO line intensity mapping

Here we consider the potential for line intensity mapping (LIM) of the rotational CO(1-0), CO(2-1), and CO(3-2) transitions to detect deviations from General Relativity from 0 < z < 3 within the framework of a very general class of modified gravity models, called Horndeski’s theories. Our forecast assumes a multitracer analysis separately obtaining information from the matter power spectrum and the first two multipoles of the redshift space distortion power spectrum. To achieve ±0.1 level constraints on the slope of the kinetic gravity braiding and Planck mass evolution parameters, a mm-wave LIM experiment would need to accumulate ≈10 8 –10 9 spectrometre-hours, feasible with instruments that could be deployed in the 2030s. Such a measurement would constrain the parameters of Horndeski’s theory at a level at worst competitive to and at best an order of magnitude tighter than existing constraints from the CMB and LSS. Our modelling code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Observation of ionization trends in a laboratory photoionized plasma experiment at Z

We report experimental and modeling results for the charge state distribution of laboratory photoionized neon plasmas in the first systematic study over nearly an order of magnitude range of ionization parameter ξ ∝ F/N e . The range of ξ is achieved by flexibility in the experimental platform to adjust either the x-ray drive flux F at the sample or the electron number density N e or both. Experimental measurements of photoionized plasma conditions over such a range of parameters enable a stringent test of atomic kinetics models used within codes that are applied to photoionized plasmas in the laboratory and astrophysics. From experimental transmission data, ion areal densities are extracted by spectroscopic analysis that is independent of atomic kinetics modeling. The measurements reveal the net result of the competition between photon-driven ionization and electron-driven recombination atomic processes as a function of ξ as it affects the charge state distribution. Results from radiation-hydrodynamics modeling calculations with detailed inline atomic kinetics modeling are compared with the experimental results. There is good agreement in the mean charge and overall qualitative similarities in the trends observed with ξ but significant quantitative differences in the fractional populations of individual ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Machine learning approaches for intentional materials engineering

In this article, the development of nanoporous metals and metallic composites through dealloying processes presents significant opportunities in materials engineering. However, designing multicomponent precursor alloys and establishing corresponding processing methods that yield predictable compositions and nanostructures remain a complex challenge. This article explores how machine learning (ML)-augmented computational and experimental methodologies can tackle these challenges by predicting precursor alloy compositions, final nanoporous structures, and mechanical properties, while integrating ML-enabled autonomous experimentation for material design and quantification. We highlight recent advancements in applying ML to nanostructured materials design via dealloying and discuss how techniques from other nanomaterial designs can be adapted for improved control over morphological and compositional outcomes in nanoporous and nanocomposite materials. Furthermore, we explore the role of ML in autonomous synchrotron x-ray experimentation, enabling real-time feedback between modeling and experimental setups. ML-driven approaches to microstructure characterization and mechanical property prediction are also examined, with a focus on modeling and advanced imaging techniques such as three-dimensional nanotomography. Finally, this article outlines future directions for ML-enhanced materials science, emphasizing the exploration of high-dimensional parameter spaces and the incorporation of materials kinetics into processing and property evaluation, ultimately advancing the design of nanoporous structures and materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗
DOE OSTI2020

Investigation of Mixing Law Efficacy for Gaseous Hydrodynamic Simulations

A computational simulation of various mixing laws for gaseous equations of state using planar traveling shocks for multiple mixtures in three dimensions is analyzed against nominal experimental data. Numerical simulations utilize the Sandia National Laboratories shock hydrodynamic code CTH and other codes including the thermochemical equilibrium code TIGER and the uncertainty qualification and sensitivity analysis code DAKOTA. The mixtures are a 1:1 and a 1:3 molar mixture of helium and sulfur hexafluoride. The mixing laws to be analyzed are the ideal gas law, Amagat's Law, Dalton's Law, the BKW EOS, the EXP6 EOS, and the JCZ3 EOS. Examination of the experimental data with TIGER revealed that the shock strength should not be strong enough to turn the mixture non-ideal as the compressibility factor z was essentially unity (z ≈ 1.02). Experimental results show that none of the equations of state are able to accurately predict the properties of the shocked mixture; similar discrepancies have been observed in previous works. Kinetic molecular theory appears to introduce a parameter which offers an explanation regarding the discrepancies. Implementation of the KMT pararneter into the EOS is left for future work.

42 ENGINEERING↗
DOE OSTI2025

Measurement of Energy Reduction of Inertial Alfvén Waves Propagating through Parallel Gradients in the Alfvén Speed

We have studied the propagation of inertial Alfvén waves through parallel gradients in the Alfvén speed using the Large Plasma Device at the University of California, Los Angeles. The reflection and transmission of Alfvén waves through inhomogeneities in the background plasma are important for understanding wave propagation, turbulence, and heating in space, laboratory, and astrophysical plasmas. Here we present inertial Alfvén waves under conditions relevant to solar flares and the solar corona. We find that the transmission of the inertial Alfvén waves is reduced as the sharpness of the gradient is increased. Any reflected waves were below the detection limit of our experiment, and reflection cannot account for all of the energy not transmitted through the gradient. Our findings indicate that, for both kinetic and inertial Alfvén waves, the controlling parameter for the transmission of the waves through an Alfvén speed gradient is the ratio of the Alfvén wavelength along the gradient divided by the scale length of the gradient. Furthermore, our results suggest that an as-yet-unidentified damping process occurs in the gradient.

79 ASTRONOMY AND ASTROPHYSICS↗