Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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 631 records · Page 35

Cellulose Structural Polymorphism in Plant Primary Cell Walls Investigated by High-Field 2D Solid-State NMR Spectroscopy and Density Functional Theory Calculations

The native cellulose of bacterial, algal, and animal origins has been well studied structurally using X-ray and neutron diffraction and solid-state NMR spectroscopy, and is known to consist of varying proportions of two allomorphs, Iα and Iβ, which differ in hydrogen bonding, chain packing, and local conformation. In comparison, cellulose structure in plant primary cell walls is much less understood because plant cellulose has lower crystallinity and extensive interactions with matrix polysaccharides. Here we have combined two-dimensional magic-angle-spinning (MAS) solid-state nuclear magnetic resonance (solid-state NMR) spectroscopy at high magnetic fields with density functional theory (DFT) calculations to obtain detailed information about the structural polymorphism and spatial distributions of plant primary-wall cellulose. 2D 13 C– 13 C correlation spectra of uniformly 13 C-labeled cell walls of several model plants resolved seven sets of cellulose chemical shifts. Among these, five sets (denoted a–e) belong to cellulose in the interior of the microfibril while two sets (f and g) can be assigned to surface cellulose. Importantly, most of the interior cellulose 13 C chemical shifts differ significantly from the 13 C chemical shifts of the Iα and Iβ allomorphs, indicating that plant primary-wall cellulose has different conformations, packing, and hydrogen bonding from celluloses of other organisms. 2D 13 C– 13 C correlation experiments with long mixing times and with water polarization transfer revealed the spatial distributions and matrix-polysaccharide interactions of these cellulose structures. Celluloses f and g are well mixed chains on the microfibril surface, celluloses a and b are interior chains that are in molecular contact with the surface chains, while cellulose c resides in the core of the microfibril, outside spin diffusion contact with the surface. Interestingly, cellulose d, whose chemical shifts differ most significantly from those of bacterial, algal, and animal cellulose, interacts with hemicellulose, is poorly hydrated, and is targeted by the protein expansin during wall loosening. To obtain information about the C6 hydroxymethyl conformation of these plant celluloses, we carried out DFT calculations of 13 C chemical shifts, using the Iα and Iβ crystal structures as templates and varying the C5–C6 torsion angle. Comparison with the experimental chemical shifts suggests that all interior cellulose favor the tg conformation, but cellulose d also has a similar propensity to adopt the gt conformation. In conclusion, these results indicate that cellulose in plant primary cell walls, due to their interactions with matrix polysaccharides, and has polymorphic structures that are not a simple superposition of the Iα and Iβ allomorphs, thus distinguishing them from bacterial and animal celluloses.

59 BASIC BIOLOGICAL SCIENCES↗

Overcoming element quality dependence of finite elements with adaptive extended stencil FEM (AES‐FEM)

Summary The finite element methods (FEMs) are important techniques in engineering for solving partial differential equations, but they depend heavily on element shape quality for stability and good performance. In this paper, we introduce the Adaptive Extended Stencil Finite Element Method (AES‐FEM) as a means for overcoming this dependence on element shape quality. Our method replaces the traditional basis functions with a set of generalized Lagrange polynomial basis functions , which we construct using local weighted least‐squares approximations. The method preserves the theoretical framework of FEM and allows imposing essential boundary conditions and integrating the stiffness matrix in the same way as the classical FEM. In addition, AES‐FEM can use higher‐degree polynomial basis functions than the classical FEM, while virtually preserving the sparsity pattern of the stiffness matrix. We describe the formulation and implementation of AES‐FEM and analyze its consistency and stability. We present numerical experiments in both 2D and 3D for the Poisson equation and a time‐independent convection–diffusion equation. The numerical results demonstrate that AES‐FEM is more accurate than linear FEM, is also more efficient than linear FEM in terms of error versus runtime, and enables much better stability and faster convergence of iterative solvers than linear FEM over poor‐quality meshes. Copyright © 2016 John Wiley & Sons, Ltd.

Conley, Rebecca↗

Characterization of shale using Helium and Argon at high pressures

In order to estimate the shale gas in place and the eventual recovery during shale gas operations, one of the key requirements is to accurately characterize the shale’s petrophysical and transport properties such as porosity, permeability, diffusivity, and storage capacity. Despite the many efforts reported in the technical literature aiming to provide an improved understanding of the complex pore structures and the associated fluid flow in gas shales, complete characterization of organic-rich shale samples still poses a big challenge. Here, we have characterized mass transfer and sorption in shale at different length scales using Helium (He) and Argon (Ar) as probe gases. Thermogravimetric analysis (TGA) with a shale cube of ~1 cm 3 in size and gas expansion experiments with a full-diameter core (3.5” in diameter) were used to measure sorption kinetics/isotherms and mass transfer, respectively. Both samples are from the same depth/location in the Marcellus shale formation. The TGA steady-state technique was initially used to generate excess sorption isotherms for Ar, while dynamic TGA experiments were used to study its sorption kinetics. The TGA experiments demonstrate that Ar, which has a similar sorption potential as Methane, but is generally assumed to be inert, adsorbs onto the surfaces of the mesoporous and microporous regions of the shale samples according to a Langmuir-type behavior. Helium expansion experiments, on the full-diameter core, were used to measure the overall porosity, on the basis that He is a non-sorbing and inert gas as compared to Ar. The He expansion experiments, furthermore, allow us to delineate the mass transfer of gas across the inherent hierarchy of pore sizes, including macropores (macro- and microcracks), mesopores and micropores. Similar expansion experiments were also performed with Ar to study the combined impact of mass transfer and sorption. A triple-porosity model (TPM) was utilized to interpret the He expansion experiments with the shale core and to extract (estimate) relevant transport parameters. We report and compare here the diffusivities and permeabilities of the whole core for both He and Ar, as calculated from the modeling and fitting of the experimental data. On the premise that the shale cube is representative of the matrix region of the core, the Ar sorption kinetics from the cube experiments were subsequently combined with the transport parameters extracted from the He experiments to predict the behavior of the Ar expansion test with the full-diameter core. An excellent agreement is observed between the model predictions and the experimental data. The experimental observations and their interpretation indicate that one must be cautious when using Ar to estimate the true porosity and permeability of shales. In addition, we demonstrate that He and Ar probe gases, when used in tandem, can be employed effectively as a tool to characterize shales in terms of mass transfer and sorption dynamics across scales.

04 OIL SHALES AND TAR SANDS↗

Limiting Retained Austenite Decomposition in Quenched and Tempered Steels: Influences of Rapid Tempering and Silicon

Tempering reactions are critical to microstructure and property control in martensitic steels. Here, retained austenite decomposition and cementite precipitation are monitored using Mössbauer spectroscopy in 4340 and 300-M steel under conventional and rapid tempering conditions. Tempering times are compared at a constant tempered hardness by increasing tempering temperatures associated with short time conditions to achieve equivalent matrix softening to that of longer tempering times. Time-temperature combinations that provide equivalent tempered hardness generated microstructures with similar dislocation densities and cementite precipitation fractions; these mechanisms are controlled by self-diffusion. However, systematic differences in retained austenite content were observed at a given degree of softening, where shorter tempering times exhibited higher levels of retained austenite compared to more conventional conditions. At low temperatures, the differences in retained austenite preservation between explored time-temperature conditions are attributed to corresponding differences in carbon diffusion distance (in austenite), the controlling diffusional process of retained austenite decomposition. At higher temperatures, retained austenite decomposition exhibits C-curve kinetic behavior in 4340. Thus, reduced thermodynamic driving force for cementite and ferrite formation at higher temperature is believed to play a role in restricting retained austenite decomposition within some short-time, high temperature tempering regimes. The addition of silicon pushes cementite precipitation and retained austenite decomposition to higher temperatures, although retained austenite decomposition is suppressed to a greater extent than cementite precipitation. Potential is illustrated for coupling rapid tempering with silicon alloying to produce appreciably tempered martensite (~490 HV) with relatively less retained austenite decomposition compared to conventional tempering conditions.

36 MATERIALS SCIENCE↗

TEMPORAL EVOLUTION OF SUB-NANOMETER COMPOSITIONAL PROFILES ACROSS THE GAMMA/GAMMA' INTERFACE IN A MODEL Ni-Al-Cr SUPERALLOY

Early-stage phase separation in a Ni-5.2 Al-14.2 Cr at.% superalloy, isothermally decomposing at 873 K, is investigated with atom-probe tomography. Sub-nanometer scale compositional profiles across the gamma/gamma'(L12) interfaces demonstrate that both the gamma-matrix and the gamma'-precipitate compositions evolve with time. Observed chemical gradients of Al depletion and Cr enrichment adjacent to the gamma'-precipitates are transient, consistent with well-established model predictions for diffusion-limited growth, and mark the first detailed observation of this phenomenon. Furthermore, it is shown that Cr atoms are kinetically trapped in the growing precipitates.

Sudbrack, Chantal K.↗

Evaluating the oxidative, photothermal and electrical stability of colloidal nanocrystal solids

IV-VI quantum dot (QD) solids are a novel class of granular electronic materials with great technological potential (e.g., in photodetectors, field-effect transistors (FETs), and solar cells), but their oxidative and thermal instability present a barrier to practical applications [1]. Poor stability is a fundamental issue facing many nanoscale materials due to high surface area and surface energy. Basic studies are needed to elucidate the most important mechanisms of degradation and develop robust countermeasures if QD materials are to become technologically important. This project determined the degradation mechanisms of IV-VI QD solids (primarily PbSe and PbS) and introduced new chemical strategies to drastically improve their performance, stability, and operating lifetimes [2-5]. Our approach was based on (1) detailed testing of QD thin film materials (principally FETs and solar cells) as a function of oxidative and thermal stress, and (2) the use of organic and inorganic approaches to link the QDs into strongly electronically coupled, high-mobility films, prevent their oxidation, and eliminate internal degrees of freedom that lead to film instability and degradation in response to electrical and thermal stress. Stability against oxidation and thermal degradation was the major focus of this project. We have evaluated the stability of QD thin films and interfaces at temperatures less than 100°C (the regime most relevant to solar and transistor applications). Low-temperature oxidation and sintering of QD films have been investigated using optical absorption spectroscopy, Fourier transform infrared spectroscopy (FTIR), current-voltage scanning of transistors and solar cells, X-ray photoelectron spectroscopy, and scanning Kelvin probe microscopy (SKPM). SKPM was used to map the potential profiles of operating QD FETs and solar cells as a function of bias and illumination, which provides detailed information on how the work functions, potential drops and electric field within these devices determine device operation, and set the stage for future studies targeted at understanding and preventing device failure. We pursued two strategies to fabricate QD films with stable electrical characteristics: (1) the use of robust molecular surface ligands [2], and (2) “matrix engineering,” i.e., infilling the QD solid with metal oxide or metal sulfide matrices by low-temperature atomic layer deposition (ALD) to passivate surface states, prevent oxidation, lock the QDs into position, inhibit diffusion, and tune the height and width of the inter-QD potential barriers that govern charge transport [5,6]. Poor stability is a common feature of nanoscale electronic materials, yet stability is all too rarely the focus of basic research. Fundamental studies are therefore needed to elucidate the most important mechanisms of degradation and develop simple yet effective countermeasures. By revealing both how QD solids degrade in response to environmental stresses (oxidative, photothermal, and electric) and how to prevent this degradation, the project has greatly improved our ability to develop stable, high-performance QD materials for real-world applications.

36 MATERIALS SCIENCE↗

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗

Study on hot deformation behavior and workability of squeeze-cast 20 vol%SiCw/6061Al composites using processing map

Highlights: • Prasad’s and Murty’s criteria were better than Gegel’s criterion in predicting the flow instability for the composites. • Partial melting and DRX at high temperature and high strain rate improved the workability of the composites. • The softening mechanism changed from DRV to DRX with the increase of deformation temperature and strain rate. • The safe deformation region of the composites was at higher temperature of 450-500°C and higher strain rate of 0.1-1s{sup -1}. • The safe deformation region predicted by processing maps was validated by uniaxial upsetting and two-directional forging. - Abstract: The hot deformation behavior and workability of squeeze-cast 20 vol%SiCw/6061Al composites were investigated by isothermal compression test and development of processing map, which was verified by uniaxial upsetting and two-directional forging of big cylindrical billets. Due to dispersed strengthening caused by SiC whiskers, the activation energy of 20 vol%SiCw/6061Al composites was much higher than the self-diffusion activation energy of pure Al. Both the DMM and MDMM were reliable in predicting the safe deformation regions. The instability regions predicted by the Prasad's and Murty's instability criterion were smaller than the actual ones, while the Gegel's instability criterion was not appropriate in predicting the instability regions of 20 vol%SiCw/6061Al composites. The optimum deformation region was at higher temperature of 450–500 °C and higher strain rate of 0.1–1 s{sup − 1}, in which the compressed samples exhibited perfect surface without cracks and microstructure with few of cavities due to the occurrence of partial melting and dynamic recrystallization. With the increase of deformation temperature and strain rate, the deformation mechanism changed from dynamic recovery to dynamic recrystallization. Squeezing casting technique with high filtration pressure may also contribute to the increase of strain rate of DRX, which was helpful for hot forging of Al matrix composites with high productivity.

36 MATERIALS SCIENCE↗

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗

Transient fuel performance analysis for the preliminary fuel concept of general atomics fast modular reactor

This study investigates the transient fuel performance of General Atomics Fast Modular Reactor (GA-FMR) during accident scenarios, focusing on the behavior of its innovative fuel system that combines high-assay low enriched uranium dioxide (HALEUO2) fuel with SiGA® ceramic matrix composite silicon carbide cladding. The preliminary fuel design’s response was analyzed during reactivity-initiated accidents (RIA) and loss of coolant accidents (LOCA) using BISON fuel performance analysis code, which included both the diffusion enhanced and BISON-FASTGRASS coupled UO 2 models. The RIA analysis demonstrated that effective reactivity control reduced fuel temperature, though with transient fission gas release resulting in additional tensile stress state on the cladding. LOCA simulations revealed differing predictions between the two models: the BISON UO 2 model showed more transient fission gas release but minimal pellet expansion, while the BISON-FASTGRASS UO 2 model predicted less pronounced fission gas release but more fuel swelling and thermal expansion, potentially leading to pellet-cladding mechanical interaction. Here, these findings highlight critical areas for fuel design optimization and identify knowledge gaps requiring further experimental and computational investigation to advance GA-FMR fuel development.

Lee, Soon K. [Argonne National Laboratory (ANL), A↗

H Atom Production and Reaction in the Gamma Radiolysis of Thermally Modified Boehmite

The correlation between the formation of H atoms and H2 was observed in the radiolysis of thermally modified boehmite (?-AlOOH). Boehmite was pretreated by heating to 300, 400, and 550°C, loaded with adsorbed water, and irradiated with gamma rays. Long-lived H atoms trapped in the interstitial matrix was observed by electron paramagnetic resonance techniques and the amount of gaseous H2 by gas chromatography. Heating the boehmite to 400°C before radiolysis results in the collapse of the matrix with loss of trapped H atoms and a substantial increase in the production of H2. Samples without added water post-heating showed no production of H2 suggesting that it is formed by the diffusion of H atoms to the surface followed by an abstraction reaction with water. Boehmite heated to 550°C has completely converted to alumina and the radiolytic responses of each are similar. Post-irradiation heating of boehmite to 300°C results in a doubling of the H2 yield at pretreated temperatures below 300°C due to the release of trapped H atoms.

La Verne, Jay A.↗

Assessment of a parabolic analysis for axisymmetric internal flows in rocket and turbomachinery ducts

The flow paths in gas turbine passages encompass a wide range of flow properties such as Reynolds number and Mach number as well as many other variable flow conditions such as swirl, free-stream turbulence, and laminar/turbulent transition. An existing computer program, the Axisymmetric Diffuser Duct (ADD) code, which calculates compressible, turbulent, swirling flow through axisymmetric ducts has been modified to include the effects of free-stream turbulence and laminar/turbulent transition. The program has been evaluated on a matrix of test cases to determine its accuracy, robustness, and limits of applicability. This improved version of the ADD code calculates solutions which compare well with available data and can now be applied to a wider range of problems than previously possible. In addition a concept called local enhancement was developed and tested on a simple two-dimensional geometry in order to demonstrate a method to reduce computer time. In this concept the pressure distribution is calculated on a coarse grid using the ADD code and the viscous layer is locally enhanced using a boundary layer analysis. By applying this concept, an order of magnitude reduction in computer time was possible without any loss in accuracy.

Power, G. D.↗

Zeolite-like MOF nanocrystals incorporated 6FDA-polyimide mixed-matrix membranes for CO 2 /CH 4 separation

MOF mixed-matrix membranes (MMMs) are regarded as potential candidates for energy-efficient natural gas purification. Our work reports the fabrication of high-performance 6FDA-polyimide MMMs, based on the incorporation of zeolite-like MOF (ZMOF) as fillers, and their associated permeation studies for CO 2 /CH 4 separation. To eliminate micron-sized crystals, a facile repeating sedimentation approach was used to harvest nanocrystals from the as-synthesized bulk ZMOF crystalline powder material. This enables the deployment of ZMOF nanocrystals with relatively uniform dimension and morphology in the polymer matrix. Typical 6FDA-polyimides encompassing distinct diamine moieties (6FDA-DAM, 6FDA-DETDA-DABA or PDMC) were explored as polymer matrices to disclose the transport property matching the hosted ZMOF filler. Mixed-gas permeation measurements revealed that the incorporation of the ZMOF filler affords a concurrent enhancement of the CO 2 permeability and the CO 2 /CH 4 selectivity for the three tested 6FDA-polyimides. Particularly, the highly permeable 6FDA-DAM showed a considerably enhanced performance for CO 2 /CH 4 that transcends the 2008 Robeson upper-bound. Detailed analysis of the sorption data and diffusion coefficients suggest that the enhanced transport property in the ZMOF-based MMM is plausibly attributed to the combination of the higher CO 2 sorption capacity and selectivity, and favorable gas diffusivity via the CO 2 -philic framework of ZMOF in moderately confined pores.

36 MATERIALS SCIENCE↗

The Identification of Complex Organic Molecules in the Interstellar Medium: Using Lasers and Matrix Isolation Spectroscopy to Simulate the Interstellar Environment

The Astrochemistry Group at NASA Ames Research Center is interested in the identification of large organic molecules in the interstellar medium Many smaller organic species (e.g. hydrocarbons, alcohols, etc.) have been previously identified by their radiofrequency signature due to molecular rotations. However, this becomes increasingly difficult to observe as the size of the molecule increases. Our group in interested in the identification of the carriers of the Diffuse Interstellar Bands (absorption features observed throughout the visible and near-infrared in the spectra of stars, due to species in the interstellar medium). Polycyclic Aromatic Hydrocarbons (PAHs) and related molecules are thought to be good candidates for these carriers. Laboratory experiments am performed at Ames to simulate the interstellar environment, and to compare spectra obtained from molecules in the laboratory to those derived astronomically. We are also interested in PAHs with respect to their possible connection to the UIR (Unidentified infrared) and ERE (Extended Red Emission) bands - emission features found to emanate from particular regions of our galaxy (e.g. Orion nebula, Red Rectangle, etc.). An old, "tried and proven spectroscopic technique, matrix isolation spectroscopy creates molecular conditions ideal for performing laboratory astrophysics.

Stone, Bradley M.↗

Microstructure, microtexture and precipitation in the ultrafine-grained surface layer of an Al-Zn-Mg-Cu alloy processed by sliding friction treatment

Precipitate redistribution and texture evolution are usually two concurrent aspects accompanying grain refinement induced by various surface treatment. However, the detailed precipitate redistribution characteristics and process, as well as crystallographic texture in the surface refined grain layer, are still far from full understanding. In this study, we focused on the microstructural and crystallographic features of the sliding friction treatment (SFT) induced surface deformation layer in a 7050 aluminum alloy. With the combination of transmission electron microscopy (TEM) and high angle angular dark field scanning TEM (HAADF-STEM) observations, a surface ultrafine grain (UFG) layer composed of both equiaxed and lamellar ultrafine grains and decorated by high density of coarse grain boundary precipitates (GBPs) were revealed. Further precession electron diffraction (PED) assisted orientation mapping unraveled that high angle grain boundaries rather than low angle grain boundaries are the most favorable nucleation sites for GBPs. The prominent precipitate redistribution can be divided into three successive and interrelated stages, i.e. the mechanically induced precipitate dissolution, solute diffusion and reprecipitation. The quantitative prediction based on pipe diffusion along dislocations and grain boundary diffusion proved the distribution feasibility of GBPs around UFGs. Based on PED and electron backscatter diffraction (EBSD) analyses, the crystallographic texture of the surface UFG layer was identified as a shear texture composed of major rotated cube texture (001) 〈110〉 and minor (111) 〈112〉, while that of the adjoining lamellar coarse grained matrix was pure brass. The SFT induced surface severe shear deformation is responsible for texture evolution. - Highlights: •The surface ultrafine grain layer in a 7050 aluminum alloy was focused. •Precipitate redistribution and texture evolution were discussed. •The quantitative prediction proved the distribution feasibility of GBPs. •Precession electron diffraction orientation mapping showed a shear texture.

36 MATERIALS SCIENCE↗

A Landau-fluid closure for arbitrary frequency response

The kinetic Landau-fluid (LF) closure which can be regarded as the exact closure is derived. For Maxwellian plasma, the kinetic closure is the same as Hammett-Perkins closure in static limit and totally the same with Chang-Callen closure. A new LF closure for arbitrary frequency response constructed with harmonic average technique is presented in this paper. This new LF closure bridges the existing LF closures in the low and high frequency limits: it recovers Hammett-Perkins closure when weight coefficient κ = 0 and converges to Chang-Callen closure at high frequency when weight coefficient κ = 1. By picking an appropriate κ, the harmonic average closure contains both nonlocal transport and local transport and the resulting fluid response function of a three moment fluid model well matches the exact kinetic response function within the entire frequency range. On the computational side, a sum of diffusion-convection solves (SDCS) method is developed to facilitate numerical implementation of the harmonic average LF closure. By using SDCS method, good agreement is achieved for the response functions between driven initial-value simulations and matrix eigenvalue calculations within the BOUT++ framework. The harmonic average closure of shifted-Maxwellian is also outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic wave propagation based upon spectral-element methodology in dispersive and attenuating media

SUMMARY We build on mathematical equivalences between Maxwell’s wave equations for an electromagnetic medium and elastic seismic wave equations. This allows us to readily model Maxwell’s wave propagation in the spectral-element codes SPECFEM2D and SPECFEM3D, written for acoustic, viscoelastic and poroelastic seismic wave propagation, providing the ability to handle complex geometries, inherent to finite-element methods and retaining the strength of exponential convergence and accuracy due to the use of high-degree polynomials to interpolate field functions on the elements, characteristic to spectral-element methods (SEMs). Attenuation and dispersion processes related to the frequency dependence of dielectric permittivity and conductivity are also included using a Zener model, similar to shear attenuation in viscoelastic media or viscous diffusion in poroelastic media, and a Kelvin–Voigt model, respectively. Ability to account for anisotropic media is also discussed. Here, we limit ourselves to certain dielectric permittivity tensor geometries, in order to conserve a diagonal mass matrix after discretization of the system of equations. Doing so, simulation of Maxwell’s wave equations in the radar frequency range based on SEM can be solved using explicit time integration schemes well suited for parallel computation. We validate our formulation with analytical solutions. In 2-D, our implementation allows for the modelling of both a transverse magnetic (TM) mode, suitable for surface based reflection ground penetration radar type of applications, and a transverse electric (TE) mode more suitable for crosshole and vertical radar profiling setups. Two 2-D examples are designed to demonstrated the use of the TM and TE modes. A 3-D example is also presented, which allows for the full TEM solution, different antenna orientations, and out-of-plane variations in material properties.

58 GEOSCIENCES↗

Aggregation of solutes in bosonic versus fermionic quantum fluids

Quantum fluid droplets made of helium-3 ( 3 He) or helium-4 ( 4 He) isotopes have long been considered as ideal cryogenic nanolabs, enabling unique ultracold chemistry and spectroscopy applications. The droplets were believed to provide a homogeneous environment in which dopant atoms and molecules could move and react almost as in free space but at temperatures close to absolute zero. Here, we report ultrafast x-ray diffraction experiments on xenon-doped 3 He and 4 He nanodroplets, demonstrating that the unavoidable rotational excitation of isolated droplets leads to highly anisotropic and inhomogeneous interactions between the host matrix and enclosed dopants. Superfluid 4 He droplets are laced with quantum vortices that trap the embedded particles, leading to the formation of filament-shaped clusters. In comparison, dopants in 3 He droplets gather in diffuse, ring-shaped structures along the equator. The shapes of droplets carrying filaments or rings are direct evidence that rotational excitation is the root cause for the inhomogeneous dopant distributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗