SEARCH · Search NASA
Results for “Diffuse interface methods”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles
In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.
NASA Tech Briefs, September 2005
Topivs include: Diamond-Coated Carbon Nanotubes for Efficient Field Emission; Improved Anode Coatings for Direct Methanol Fuel Cells; Advanced Ablative Insulators and Methods of Making Them; PETIs as High-Temperature Resin-Transfer-Molding Materials; Stable Polyimides for Terrestrial and Space Uses; Low-Density, Aerogel-Filled Thermal-Insulation Tiles; High-Performance Polymers Having Low Melt Viscosities; Nonflammable, Hydrophobic Aerogel Composites for Insulation; Front-Side Microstrip Line Feeding a Raised Antenna Patch; Medium-Frequency Pseudonoise Georadar; Facilitating Navigation Through Large Archives; Program for Weibull Analysis of Fatigue Data; Comprehensive Micromechanics-Analysis Code - Version 4.0; Component-Based Visualization System; Software for Engineering Simulations of a Spacecraft; LabVIEW Interface for PCI-SpaceWire Interface Card; Path Following with Slip Compensation for a Mars Rover; International Space Station Electric Power System Performance Code-SPACE; Software for Automation of Real-Time Agents, Version 2; Software for Optimizing Plans Involving Interdependent Goals; Computing Gravitational Fields of Finite-Sized Bodies; Custom Sky-Image Mosaics from NASA's Information Power Grid; ANTLR Tree Grammar Generator and Extensions; Generic Kalman Filter Software; Alignment Stage for a Cryogenic Dilatometer; Rugged Iris Mechanism; Treatments To Produce Stabilized Aluminum Mirrors for Cryogenic Uses; Making AlNx Tunnel Barriers Using a Low-Energy Nitrogen-Ion Beam; Making Wide-IF SIS Mixers with Suspended Metal-Beam Leads; Sol-Gel Glass Holographic Light-Shaping Diffusers; Automated Counting of Particles To Quantify Cleanliness; Phase Correction for GPS Antenna with Nonunique Phase Center; Compact Infrasonic Windscreen; Broadband External-Cavity Diode Laser; High-Efficiency Solar Cells Using Photonic-Bandgap Materials; Generating Solid Models from Topographical Data; Computationally Lightweight Air-Traffic-Control Simulation; Spool Valve for Switching Air Flows Between Two Beds; Partial Model of Insulator/ Insulator Contact Charging; Asymmetric Electrostatic Radiation Shielding for Spacecraft; and Reusable Hybrid Propellant Modules for Outer-Space Transport.
Process-dependent anisotropic thermal conductivity of laser powder bed fusion AlSi10Mg: impact of microstructure and aluminum-silicon interfaces
Purpose AlSi10Mg alloy is commonly used in laser powder bed fusion due to its printability, relatively high thermal conductivity, low density and good mechanical properties. However, the thermal conductivity of as-built materials as a function of processing (energy density, laser power, laser scanning speed, support structure) and build orientation, are not well explored in the literature. This study aims to elucidate the relationship between processing, microstructure, and thermal conductivity. Design/methodology/approach The thermal conductivity of laser powder bed fusion (L-PBF) AlSi10Mg samples are investigated by the flash diffusivity and frequency domain thermoreflectance (FDTR) techniques. Thermal conductivities are linked to the microstructure of L-PBF AlSi10Mg, which changes with processing conditions. The through-plane exceeded the in-plane thermal conductivity for all energy densities. A co-located thermal conductivity map by frequency domain thermoreflectance (FDTR) and crystallographic grain orientation map by electron backscattered diffraction (EBSD) was used to investigate the effect of microstructure on thermal conductivity. Findings The highest through-plane thermal conductivity (136 ± 2 W/m-K) was achieved at 59 J/mm 3 and exceeded the values reported previously. The in-plane thermal conductivity peaked at 117 ± 2 W/m-K at 50 J/mm 3 . The trend of thermal conductivity reducing with energy density at similar porosity was primarily due to the reduced grain size producing more Al-Si interfaces that pose thermal resistance. At these interfaces, thermal energy must convert from electrons in the aluminum to phonons in the silicon. The co-located thermal conductivity and crystallographic grain orientation maps confirmed that larger colonies of columnar grains have higher thermal conductivity compared to smaller columnar grains. Practical implications The thermal properties of AlSi10Mg are crucial to heat transfer applications including additively manufactured heatsinks, cold plates, vapor chambers, heat pipes, enclosures and heat exchangers. Additionally, thermal-based nondestructive testing methods require these properties for applications such as defect detection and simulation of L-PBF processes. Industrial standards for L-PBF processes and components can use the data for thermal applications. Originality/value To the best of the authors’ knowledge, this paper is the first to make coupled thermal conductivity maps that were matched to microstructure for L-PBF AlSi10Mg aluminum alloy. This was achieved by a unique in-house thermal conductivity mapping setup and relating the data to local SEM EBSD maps. This provides the first conclusive proof that larger grain sizes can achieve higher thermal conductivity for this processing method and material system. This study also shows that control of the solidification can result in higher thermal conductivity. It was also the first to find that the build substrate (with or without support) has a large effect on thermal conductivity.
Global Energy and Aviation Concerns
Renewable energy sources are usually diffuse and require large facilities. Biofuels work better, are more economical to produce for ground transportation, but sharply increase competition for food croplands. Noble laureate Richard Smalley (deceased-2005) conceptual 20 TWe power generation covers hundreds x hundreds of miles. Combined with Fuller s superconducting power grid system would enable renewable planetary energy. A solar-wind project in Australia will have a 7km diameter collector interfacing with a 1 km tower to extract 200 MW from wind turbines mounted at the base. GE Energy s 3.5MW Wind Turbine is large and placing this in perspective, it is as if one were rotating a Boeing 747-200; the blade diameter is that large. Wind turbines are rapidly gaining popularity in Europe and photovoltaic (PV) is expected to also expand rapidly. It becomes clear that we need (and still have time) to develop new sources of energy. Hf 178 bombarded by X-rays produces Gamma-rays for heating. The reaction stops when the X-rays stop; the half life is about 30 years and seems manageable vs 30 000 years. Water splitting needs to be perused as do ultra fast ultra intense laser applications in terms of fusion and new materials developments including new ways to strip and re-bind hydrogen into fuels. New methods and tools for development are being found in quantum mechanical applications to macro-systems and need to be developed into a set of new tool boxes for development of these new energy sources.
Optimization of the compositions of polyanionic sodium-ion battery cathode NaFe 2-x V x (PO 4 )(SO 4 ) 2
Sodium (Na) super ionic conductor (NASICON) polyanionic compounds have recently attracted much attention from the battery community because of their electroactive properties and reasonably high ionic conductivities, leading to their use as a cathode in sodium-ion batteries. This article describes the compositional optimizations, crystallographic evaluations, and electrochemical behavior of a new mixed NASICON polyanionic compound, NaFe 2-x V x (PO 4 )(SO 4 ) 2 . By doping the characteristic Fe 3+ sites of the FeO 6 octahedrons with varying amounts of V 3+ , the electrochemical stability and charge transport in NaFe 2 (PO 4 )(SO 4 ) 2 were enhanced. The resulting best composition, with crystal structure NaFe 1.4 V 0.6 (PO 4 )(SO 4 ) 2 resolved through the Rietveld method, exhibited a stable capacity compared with the other synthesized compositions. In situ powder x-ray diffraction measurements, a single-phase intercalation/deintercalation mechanism of the NASICON structure in the measured sodium concentration window was observed with no impurity phase formation. Further electrochemical assessments revealed the interfacial charge transfer kinetics to be the rate-limiting step in the sodium concentration window. Also, the measured sodium-ion diffusivity values in the range of 6 × 10 -11 to 7 × 10 -11 cm 2 /s in the measured sodium concentration range. The results reported here highlight the potential of compositionally and morphologically optimized NaFe 1.4 V 0.6 (PO 4 )(SO 4 ) 2 with higher particle surface areas as a cathode material for high-performance sodium-ion batteries.
Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)
This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.
Realization of a high luminosity muon beam at the Jefferson Lab
The proton is the major component of the visible matter of our universe, so being able to determine its fundamental properties, such as its charge radius, is of vital importance. In recent years, several experiments have been performed for the accurate measurement of the proton radius with different experimental techniques. The main and most common are based on two different techniques: ep scattering and Lamb shift spectroscopy of regular hydrogen atoms. The latter can also be performed on muonic hydrogen atoms with a significant improvement in accuracy. Recent measurements that use this method have highlighted a discrepancy between the results obtained with the diffusion of electrons and those obtained with the spectroscopic technique leading to the so-called ?proton radius puzzle?. The explanation of why this difference between the measured values exists is still a matter of discussion. The last attempt to address the puzzle was the PRad Experiment. It was performed with a unique experimental setup that gave the scientists an excellent control on the systematic errors that are usually present in electron-proton scattering experiments. This has significantly increased the precision of the measurement, bringing it in the same area as the spectroscopic ones done on muonic hydrogen atoms, which are considered the most precise. However, to solve the puzzle once and for all, even more precise electron scattering measurements are needed. Alternatively, new methods can be used, such as eg. muon-proton scattering at low transferred momenta in order to lower the systematic errors that comes with the measurement. This thesis tries to address this problem by studying the set up of an intense muon beam at the Jefferson Laboratory (JLab) that can be used to carry out a high-precision m ? p scattering experiment. Using Montecarlo simulations I studied the muon production resulting from the interaction of the primary electron beam with the beam-dump at JLab Hall-A. The simulations were performed using the FLUKA package via the FLAIR graphical interface. In this study, I modified an existing geometry (built for the ?Beam Dump eXperiment?) and performed several studies in order to optimize the muon beam parameters. In order to estimate the beam intensity, I studied the muon generation and transportation through the concrete bunker surrounding the Hall-A beam-dump. Detailed studies on muon attenuation and dispersion at different depth allowed me to track the beam profile in space and momentum as a function of the distance from the beam-dump. Last, but not least, I studied a possible tagging system for the muon beam. The tagging system will allow us to know with good precision (O(1=1000)) the energy of each muon of the beam. Based on the results obtained, I created an event generator that incorporates the detailed description of the muon beam (energy, position, emittance) useful to define and optimize the experimental setup for a future measurement of the proton radius at the JLab.
Effect of Sulfonation Level on the Percolated Morphology and Proton Conductivity of Hydrated Fluorine-Free Copolymers: Experiments and Simulations
Using all-atom molecular dynamics simulations and a variety of experimental methods, we previously reported on a linear polyethylene with pendant phenyl sulfonated groups precisely on every fifth carbon along the backbone. With increasing relative humidity this fluorine-free polymer self-assembled to form nanoscale water channels and exhibited exceptional proton conductivity. Expanding upon those findings, here we explore partially sulfonated random copolymers, referred to as p 5PhSH-Y. Using either acetyl sulfate or sulfuric acid, a wide range of sulfonation levels were prepared ( Y = 34−98%) corresponding to ion-exchange capacities (IEC) of 2.0−4.4 mmol/g. Combining experimental techniques and all-atom molecular dynamics simulations, we study the effect of Y on water uptake, nanoscale morphology, and the proton/water transport properties of p5PhSH- Y . The proton conductivity of p 5PhSH- Y increases with relative humidity and with Y and achieves values in excess of 0.1 S/cm. These high conductivities are attributed to high IEC and welldeveloped nanoscale percolated hydrophilic domains made possible by the flexible backbone. We quantitatively describe the nature of the water channels using the characteristic distance, channel width distribution, the area per sulfonate group at the hydrophilic/ hydrophobic interface, and the fractal dimension. Notably, the channel widths and the areas per sulfonate group are nominally independent of the level of sulfonation, while depending significantly on the level of hydration. The fractal dimension of the water channels correlates strongly with the water diffusion coefficients calculated from the molecular dynamics (MD) simulations. These findings demonstrate that the p 5PhSH- Y hydrocarbon copolymers can be modified to tune properties, particularly proton conductivity.
Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers
We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.
Probing phonon focusing, thermomechanical behavior, and moiré patterns in van der Waals architectures using surface acoustic waves
Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS 2 . Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS 2 bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS 2 monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS 2 by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS 2 layers.
Magnetic damping in ferromagnetic/heavy-metal systems: The role of interfaces and the relation to proximity-induced magnetism
Damping and spin transport in spintronic multilayered systems continues to be a topic of active research. The enhancement of damping in ferromagnet (FM)/spacer layer (SL)/heavy-metal (HM) thin-film systems was studied for Co 25 Fe 75 /SL/Pt with a nonmagnetic (NM) SL of either Au or Cu with variable thickness, in order to understand the correlation with proximity-induced magnetism (PIM) in the HM. Structural, PIM and magnetic damping measurements were undertaken on the same samples. Specifically, secondary ion mass spectroscopy, element specific x-ray magnetic reflectivity and x-ray magnetic circular dichroism at the Pt and Au L 3 edges, and ferromagnetic resonance methods were used. With increasing thickness of a Cu or Au SL directly between the FM and the Pt layer, the Pt PIM and the damping both fall rapidly, with a relationship between damping and PIM that depends on the SL material. The PIM observed in the Au layer showed a complex dependence on the layer thickness, suggesting some hybridization with the Pt. The role of the number and location of interfaces on the damping was demonstrated with the addition of a SL within the Pt layer, which showed that the specific details of the NM/HM interface also affects the damping. The insertion of a Cu SL within the Pt showed a measurable increase in the overall enhancement of the damping while the insertion of a Au SL into Pt had almost no effect on the damping. Finally, together these results demonstrate the role of both PIM and of additional interfaces in the enhancement of damping in FM/HM systems, which is not fully accounted for by existing theory
Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses
In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.
The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas
The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.
Novel and effective surface passivation for high efficiency n- and p-type Silicon solar cell
The project objective was to develop a novel Si surface passivation method using chalcogens, sulfur (S) and/or selenium (Se), as passivating elements, to withstand industry-standard high temperature contacting and metallization schemes for p-type Si based passivated emitter and rear contact (p-PERC) solar cells. The back surface passivation of PERC cells has been improved drastically with the invention and successful application of an Al 2 O 3 passivation layer. However, the front n + diffused junction surface is still poorly passivated by the standard amorphous silicon nitride (SiNx) anti-reflection coating (ARC) layer. This project sought to address the passivation challenges of both front n+ emitter and undiffused p-Si back surface. Improved p-PERC solar cell performance with open circuit voltage (V OC ) > 680 mV and efficiency of 22% were targeted to validate superior defect passivation properties as compared to standard SiO 2 / Al 2 O 3 passivation. During this project, we systematically investigated process-structure-properties-performance relationships of this novel advanced defect passivation approach. The S/Se passivation was carried out by reacting industrial Czochralski (Cz) Si wafers in H 2 S and H 2 Se gases in a chemical vapor deposition (CVD) reactor at temperatures up to 700°C. After an exhaustive optimization of the process parameters (temperature, time, and gas concentration), we established an optimized process and demonstrated extremely low surface recombination velocities (SRVs) of 1.5 cm/s and 8 cm/s on n-type and p-type Si, respectively, by S-passivation. In-depth surface and interface characterization were performed using soft x-ray and photoelectron spectroscopies (XPS, UPS, XES), combined with capacitance-voltage-frequency (C-V-f) measurements, to decipher the surface chemical/electronic structure and interface defect state densities. These measurements provided critical understanding of the defect passivation mechanism and elucidated the presence of surface S-Si bonds, a reduction of surface dipoles, and low interface state densities (D it ) < 10 11 cm -2 ev -1 . We also found that the Se-passivation is inferior to the S-passivation (by at least one order of magnitude in SRV). Application of the optimized S-passivation to the n+ diffused emitter surface led to a low surface recombination current density, J0 ≈ 40 fA/cm 2 (~ 1/4 of the industry-standard SiNx-passivation), and high implied V OC (686 mV) in p-PERC solar cell structures. The S-passivation process also was found to improve the bulk quality of the p-type Si, better than the SiO 2 or Al 2 O 3 passivation processes. After successful demonstration of efficient passivation of Si surface defects by S, we extensively studied the air, thermal, and illumination stability of the passivation structure. S-passivation itself degrades in air due to competing reactions with moisture and oxygen to form oxides, which can be eliminated by a SiNx capping layer (also acting as a anti-reflective coating). After SiNx process optimization, we demonstrated illumination and thermally stable S-passivation with SRV < 5 cm/s and J 0 < 80 fA/cm 2 . These enhancements in Si passivation, incorporated into p-PERC cells, achieved an efficiency of 19.93% with V OC = 649 mV, using manufacturing metallization and contacting schemes. The low cell performance (cell V OC is much less than the implied V OC = 686 mV, anticipated from surface passivation) was identified due to degradation of S-passivation during the metal firing step (out-diffusion of S from the Si interface to the SiNx surface). The S-passivation of Si surfaces shows significant promise with excellent passivation quality, essential for high performance (high V OC , high efficiency) solar cells. Integration of this innovative defect passivation into devices, however, demands further development of the capping layer, low temperature (<700°C) metallization process, and/or engineering of advanced device structures. Surface passivation-dominated advanced Si solar cells, such as tunnel oxide passivated contacts and Si heterojunctions, are increasingly of interest due to their high-performance potential and will have a growing photovoltaic market share in the near future.
Weak scaling of the contact distance between two fluctuating interfaces with system size
A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.
MEUMAPPS (C++ Version)
Many materials, metal alloys in particular, have features on the on micrometer or nanometer scale that have a large impact on the properties of the material. These features are known as the microstructure of the material. Understanding why and how the microstructure forms in a material is of fundamental scientific interest as well as of significant technological interest. The capability to predict microstructure evolution in a material allows the intentional design of microstructures and hence the intentional design of material properties. The phase-field method is one of the leading methods for predicting microstructure evolution. One of the most significant problems for phase-field models is their computational expense. Even limited phase-field simulations can easily require thousands of CPU core-hours to complete, which significantly limits their use. This code provides both a general framework for creating scalable, GPU-accelerated phase-field model applications as well as several applications themselves. The code is capable of using hundreds of GPUs efficiently, which greatly reduces the time required to perform simulations. The code is written with an emphasis on performance portability, that is the ability for the code to run efficiently on a number of different computing architectures without modification of the source code. The performance portability of this code is primarily enabled through the use of two libraries, Kokkos (performance portable data structures and execution patterns) and heFFTe (performance portable distributed 3D fast Fourier transforms). The code consists of a core library, applications, and tests. The core library includes shared functionality between applications. This includes interfaces with fast Fourier transform (FFT) libraries such as heFFTe, data structures based on Kokkos, file input and output capabilities, and a solver for infinitesimal strain mechanical equilibrium problems. Five applications are included in the code. The flagship application is the MEUMAPPS-SS application, which implements the Kim-Kim-Suzuki phase-field model for precipitation for an arbitrary number of phases and components in a metal alloy. Five simpler applications are also included that solve the Eshelby inclusion problem, Allen-Cahn equation, the coupled Allen-Cahn and diffusion equations, and the Cahn-Hilliard equation. The code includes two applications to solve the Cahn-Hilliard equation, one with constant-step-size first-order time integration and the second with adaptive high-order time integration.
Multi-technique characterization of spray coated and roll-to-roll coated gas diffusion fuel cell electrodes
Here this work reports the characterization of catalyst layer (CL) structure and composition for a set of gas-diffusion electrodes (GDEs) fabricated by different coating methods with the goal of developing a better understanding of CL processing-structure-property-performance relationships. The CL coating techniques used in this study were ultrasonic spray coating, and two roll-to-roll (R2R) methods - gravure and slot die. Scanning transmission electron microscopy (STEM) coupled with X-ray energy dispersive spectroscopy (EDS) analysis of GDE cross-sections provided bulk information, while scanning electron microscopy (SEM) with EDS and X-ray photoelectron spectroscopy (XPS) were used to investigate the near surface and surface composition of the CLs that will be in direct contact with the membrane once the membrane electrode assembly (MEA) is constructed. This study demonstrates that common quantification parameters, namely the F:Pt ratio extracted from cross-sectional STEM-EDS maps and top-down CL measurements with SEM and XPS, along with roughness parameters quantified from SEM micrographs are useful parameters in explaining and potentially predicting performance trends. Additionally, it highlights the importance of understanding key CL surface properties for advancing the manufacturability of high-performance components for polymer electrolyte membrane (PEM) technologies.