Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion problems”

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 1,153 records · Page 64

Multiscale approaches for optimizing the impact of strain on Na-ion battery cycle life

Abstract The high costs and geopolitical challenges inherent to the lithium-ion (Li-ion) battery supply chain have driven a rising interest in the development of sodium-ion (Na-ion) batteries as a potential alternative. Unfortunately, the larger ionic radius of Na limits the reversibility of cycling because of the extensive atomic rearrangements that accompany Na-ion insertion, which in turn limit diffusion and charging speed, and lead to rapid degradation of the electrodes. The Center for Strain Optimization for Renewable Energy (STORE) was established to address these challenges and develop new electrode materials for Na-ion cells. This article discusses the current state-of-the-art materials used in Na-ion cells and several directions that STORE believes are critical to understand and control the structural and volumetric changes during the reversible (de)insertion of large cations. Graphical abstract Highlights Understanding the fundamental way materials respond to localized strains at the atomic length-scale is a critical first step in the development of highly reversible, long cycle life, Na-ion insertion hosts. This perspective explores a variety of methods that can be employed to mitigate the detrimental effects of large strain. The insights gained from these investigations should help lay the foundation for the creation of more economical and sustainable batteries that could have immediate impact on global energy infrastructure. Discussion Although there is near universal agreement that electrochemical energy storage must be an integral part of a green-energy future, there is less agreement about how to reduce the cost of energy storage. Replacing high-cost lithium-ion cells with lower-cost sodium-ion batteries is one option frequently considered in future energy models, but the details of what can be achieve with optimized sodium cell performance remains unclear. Here we posit that developing methods to mitigating strain on the electrode particle length scale is a key factor for achieving long-cycle-life sodium-ion batteries. Mitigating strain on the atomic scale suppress electrode-level volume change. Allowing for fast cycling in materials without the problems of electrode cracking or delamination. We further posit that understanding volume change in sodium-ion electrodes at a fundamental level will lead to the designing new sodium-ion electrode materials that will allow for efficient, stable, lower-cost energy storage.

Brady, Michael J.↗

Large-Time Behavior of Solutions to Vlasov-Poisson-Fokker-Planck Equations: From Evanescent Collisions to Diffusive Limit

The present contribution investigates the dynamics generated by the two-dimensional Vlasov-Poisson-Fokker-Planck equation for charged particles in a steady inhomogeneous background of opposite charges. We provide global in time estimates that are uniform with respect to initial data taken in a bounded set of a weighted L{sup 2} space, and where dependencies on the mean-free path τ and the Debye length δ are made explicit. In our analysis the mean free path covers the full range of possible values: from the regime of evanescent collisions τ→∞ to the strongly collisional regime τ→0. As a counterpart, the largeness of the Debye length, that enforces a weakly nonlinear regime, is used to close our nonlinear estimates. Accordingly we pay a special attention to relax as much as possible the τ-dependent constraint on δ ensuring exponential decay with explicit τ-dependent rates towards the stationary solution. In the strongly collisional limit τ→0, we also examine all possible asymptotic regimes selected by a choice of observation time scale. Here also, our emphasis is on strong convergence, uniformity with respect to time and to initial data in bounded sets of a L{sup 2} space. Our proofs rely on a detailed study of the nonlinear elliptic equation defining stationary solutions and a careful tracking and optimization of parameter dependencies of hypocoercive/hypoelliptic estimates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coupled Neutron and Gamma Heating Calculation Based on VARIANT Nodal Transport Method

Accurate determination of temperature distributions in fuel pins and duct walls is essential for fuel-pin performance and assembly bowing analyses. Currently, the multi-assembly, sub-channel analysis code SE2-ANL is used in thermalhydraulic analysis of sodium cooled fast reactors (SFR). SE2- ANL is interfaced with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In this procedure, the neutron flux distribution is first calculated using the finite-difference diffusion theory option of the DIF3D code. Using the neutron flux solution, the GAMSOR code computes the gamma source distribution for the subsequent gamma flux calculation with DIF3D. However, in the current heating calculation scheme, the neutron and gamma flux calculations are limited to six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. In each assembly, the axial power profile is obtained by averaging the power densities of six triangular meshes at each axial mesh interval. The radial power profile is evaluated at a single axial mesh (e.g., at the core mid-plane) by approximating the neutron and gamma flux distributions by a linear or quadratic polynomial of two variables x and y, depending on the polynomial order specified by the user. The coefficients of the linear and quadratic polynomials are determined by the least squares fit of six and 24 triangular mesh fluxes, respectively. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. With the advancement of computer technology, the nodal transport theory code VARIANT is routinely used in core design analyses of SFR instead of DIF3D. It is therefore necessary to develop a new heating calculation procedure based on VARIANT neutron and gamma transport calculations. In this study, a new, coupled neutron and gamma heating calculation scheme has been developed to determine the pin power distributions more accurately by eliminating the aforementioned limitations of the current procedure. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. In the new heating calculation procedure, the neutron flux distribution is determined first by solving an eigenvalue transport problem using VARIANT. The obtained neutron flux is then used in GAMSOR to calculate the gamma source distribution. Using this gamma source distribution, the gamma flux distribution is determined by solving a fixed-source transport problem using VARIANT again. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at each fuel pin position and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The performance of the new heating calculation scheme based on VARIANT transport solutions has been tested against the existing method of SE2-ANL, using the ABR-1000 metal core design. The gamma source distribution and the assembly, fuel pin and duct wall power distributions were calculated and compared. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A GLOBAL GALACTIC DYNAMO WITH A CORONA CONSTRAINED BY RELATIVE HELICITY

We present a model for a global axisymmetric turbulent dynamo operating in a galaxy with a corona that treats the parameters of turbulence driven by supernovae and by magneto-rotational instability under a common formalism. The nonlinear quenching of the dynamo is alleviated by the inclusion of small-scale advective and diffusive magnetic helicity fluxes, which allow the gauge-invariant magnetic helicity to be transferred outside the disk and consequently to build up a corona during the course of dynamo action. The time-dependent dynamo equations are expressed in a separable form and solved through an eigenvector expansion constructed using the steady-state solutions of the dynamo equation. The parametric evolution of the dynamo solution allows us to estimate the final structure of the global magnetic field and the saturated value of the turbulence parameter α{sub m}, even before solving the dynamical equations for evolution of magnetic fields in the disk and the corona, along with α-quenching. We then solve these equations simultaneously to study the saturation of the large-scale magnetic field, its dependence on the small-scale magnetic helicity fluxes, and the corresponding evolution of the force-free field in the corona. The quadrupolar large-scale magnetic field in the disk is found to reach equipartition strength within a timescale of 1 Gyr. The large-scale magnetic field in the corona obtained is much weaker than the field inside the disk and has only a weak impact on the dynamo operation.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic and equilibrium reactions on natural and laboratory generation of thermogenic gases from Type II marine shale

The phenomenon that laboratory pyrolysis experiments produce much wetter gases than those in natural reservoirs is a long-recognized and debated problem in the investigation of natural gases in sedimentary basins. In this study, we explore the discrepancy by pyrolyzing a type II kerogen from the Woodford Shale in Oklahoma, compared with the previous results on the produced natural gases from the Arkoma Basin generated from the same source rock (Liu et al., 2019) with the discussion of gas and isotopic compositions at bulk and position-specific (PS) levels. An improved GC-pyrolysis-GC IRMS method is applied for the determination of PS δ 13 C of propane produced in the pyrolysis of the Woodford Shale at Easy %R o from 0.76 to 3.27. Kinetic and thermodynamic considerations of the chemical and isotopic compositions of the natural and laboratory pyrolysis gases suggest that the generation of light hydrocarbons involves uni-directional cracking reactions, exchange reactions with water, and likely reversible reactions among light hydrocarbons and other H-containing volatiles. After the gas generation in the unconventional Woodford Shale reservoirs, the C 1 -C 4 gases might have approached close to chemical equilibrium of C 1 -C 3 and isotope equilibrium of C 2 -C 1 and C 3 -C 1 pairs at their peak temperatures. The capping H for the generation of C 1 -C 4 in the Woodford Shale gases appears to have experienced at least partial exchange with the water, while that in the pyrolysis gases is only originated from organic-bound compounds with large kinetic isotope effects (KIE). Our findings indicate that elevated compound-specific and PS δ 13 C values of propane in the wet-gas cracking stage are significantly influenced by the breakdown of the thermally stable compounds (e.g., remaining kerogen, residues). A first synthesis of PS δ 13 C and δ 2 H isotopic compositions of propane from this study and the literature data suggests relatively similar isotopic structures of propane precursors in kerogens. Finally, this study demonstrates that PS isotope analysis of propane can contribute to identifying various geological (e.g., maturation, wet-gas cracking, H exchange, diffusion) and biodegradation processes.

58 GEOSCIENCES↗

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗

Corrosion of Zinc Cold Spray Coatings in a Wet Sweet and Sour Gas Environment

Internal corrosion is a problem for steel pipelines transporting natural gas or CO 2 containing water and partial pressures of H 2 S higher than 0.3 kPa (0.05 psi). This work aims to mitigate internal corrosion in steel pipelines transporting natural gas containing H 2 S using cold spray coatings. Two types of the cold spray binary metallic coatings (zinc chromium [ZnCr]. zinc niobium [ZnNb]) were studied using electrochemical techniques: potentiodynamic polarization, linear polarization resistance, and electrochemical impedance spectroscopy. The corrosion resistance of cold spray coatings (ZnCr, ZnNb) was evaluated in an environment containing 4 bar CO 2 pressure, simulating the partial pressures found in gas transmission lines over a solution of 3.5 wt% NaCl heated to 40°C. A concentration of 0.003 M Na 2 S 2 O 3 ·5H 2 O, corresponding to H 2 S partial pressures around 0.079 bar (1.146 psi), was used to simulate sour conditions. Postcorrosion surface characterization was performed using a scanning electron microscope (SEM) equipped with an energy-dispersive x-ray spectroscope (EDS) and x-ray diffraction analysis. The data showed that the presence of 0.003 M Na 2 S 2 O 3 ·5H 2 O shifted the corrosion potential to more anodic values and decreased the corrosion current density. Both coatings showed similar behavior after 1 h of exposure in the CO 2 /H 2 S environment, indicating that similar electrochemical reactions were occurring on ZnNb and ZnCr. SEM images and EDS surface analyses for specimens showed a significant change in the surface chemical composition of carbon steel coated with ZnNb and ZnCr after 24 h of immersion. In the presence of thiosulfate (under sour conditions), the formation of corrosion product layers (ZnCO 3 and ZnS) on top of ZnNb and ZnCr coatings increased their corrosion resistance, which helped to reduce their corrosion by a factor of 2. Under a sweet environment, the corrosion rates for steel coated with cold spray coatings after 14 d of exposure are lower than that for galvanized steel by a factor of 5 due to the ZnCO 3 layer formed on top of the coatings. The ZnCO 3 layer formed on the steel surface acts as a physical barrier against corrosion by blocking the diffusion of corrosive species to the surface. No localized attack was observed. ZnCr Cold spray coating with defect showed promising corrosion protection against CO 2 corrosion (sweet corrosion) after 14 d of exposure to a CO 2 environment. Here, the scratch on the coating simulated damage created in service, and it was deep enough to expose the substrate material (steel). The formation of zinc oxide (ZnO) and zinc carbonate (ZnCO 3 ) on the scratch confirmed the cathodic protection of the steel by ZnCr and ZnNb coatings.

36 MATERIALS SCIENCE↗

A Study of the von Neumann Stability Analysis of a Semi-implicit Numerical Method Applied Over the Radial Impurity Transport Equation in Tokamak Plasma

The radial impurity transport equation for tokamak plasma is a form of diffusion–convection–reaction equation. The impurity transport equation is solved to determine the distribution of impurity (non-fuel) ion species with different ionization states perpendicular to magnetic surfaces of tokamak plasma. The equation for each charge (ionization) state Z is a non-linear, second-order in space, first-order in time, parabolic partial differential equation coupled to the previous Z − 1 and the next Z + 1 charge states of the impurity species through its reaction term. The number of differential equations to be solved simultaneously is hence determined by the number of ionization states of the impurity species studied. The solution to the set of these coupled equations can be obtained using a semi-implicit numerical method applied on it. The present study describes the application of von Neumann stability analysis over the semi-implicit numerical method applied over the radial impurity transport equation and determines a generic stability criterion for the method. The stability analysis is further illustrated using the geometry of Aditya tokamak installed at the Institute for Plasma Research Gandhinagar, India as an example. The impurity species considered is oxygen (Atomic number = 8). This leads to a set of eight coupled equations for charge states Z = 1 to 8 over which von Neumann analysis is illustrated in present study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗

Novel Serpentine Robot Combinations for Inspection in Hard-to-Reach Areas of Damaged or Decommissioned Structures - 17335

U.S. government agencies, particularly the Dept. of Energy, need robotic tools to handle and inspect high-consequence materials - materials for which mishandling would result in great risk to both application operators and the general public. The remediation and documentation of nuclear waste has become a problem of national significance and it cannot be addressed without robotic tools. The use of robots bypasses the danger and decreases the preparation time to handle harmful materials. In the 1980's and 90's, the US Dept. of Energy (DOE) invested heavily in robotic technologies as part of a comprehensive program of research and development aimed at enabling safer and more efficient cleanup of nuclear waste sites. For the next generation of harder, more complex clean-up problems, there is a clear need for more R and D. There is a clear need for flexible, quickly deployable robotic technologies that can be configured to respond to unforeseen situations -- or can be purpose-built for planned activities -- while keeping human workers safe, yet effective. Modularity and inter-operability are core capabilities for inspection tasks in contaminated facilities with a large variety of inspection problems. With any tool, familiarity breeds ingenuity, and ingenuity leads to new innovations in efficiency, effectiveness, safety, and performance. For example, the simple screwdriver can be held fourteen different ways by an experienced machinist who knows exactly what he or she wants to accomplish at any given moment. As a tool, robots are no different. To maximize the efficiency, effectiveness, safety and performance, redeployment and reuse of familiar tools is essential. Modularity can assist with re-deployment and modularity manifests itself in both hardware and software. In some applications of inspection of high-consequence materials, a common concern has been the simplicity of the tool -- often correlated with cost (which is not always accurate) -- in case it becomes unrecoverable. In the case of robots, which have complex user interfaces and complex dynamics of interaction, the lack of familiarity can lead to clumsy operations and reduced accumulation of knowledge. As a result, while incremental successes may be achieved with sequential uses of different robots, there often is little knowledge gain that is generalizable across multiple uses. Total cost often goes down when a more sophisticated, but more familiar tool, can be used more frequently.An example is the inspection of the H-Canyon Exhaust Air Tunnel at the Savannah River Site nuclear waste facility in the United States [1]. This is a particularly nasty environment, due to the acidic gases, large amounts of erosion and debris, and wet, mucky puddles. Over the past twelve years, six inspections have been attempted with five different robotic vehicles. In fact, these inspections have been successful with their limited, but progressively increasing objectives, although most robots have been abandoned in the tunnel. While each vehicle has cost less than $75,000 (the sacrificial part), the total costs of development have been much higher. The robots used to investigate changed much from iteration to iteration and it is not clear a strong body of generalizable expertise was built up across trials. The authors of this study demonstrated the combination of two types of modular, serpentine, re-configurable robots for locomotion over complex terrains and the inspection of hard-to-reach areas. The MOTHERSHIP is a tread/limb/serpentine hybrid robot design for locomotion over complex terrains, while the CMU snake robot provides highly articulated inspection of highly confined spaces. The demonstration at the Portsmouth Gaseous Diffusion Plant involves a combination of the two whereby gross locomotion is achieved with the MOTHERSHIP and the CMU Snake deploys a sensor package to gather data from behind and under obstacles. Although statically configured in its present state, the combined robot is built of individual modules for which reconfiguration has been demonstrated in both labs at Purdue and CMU. In future work, much higher degrees of reconfiguration and re-purposing of both software and hardware are easily achievable with these unique platforms. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimating an executive summary of a time series: the tendency

In this paper, we revisit the problem of decomposing a signal into a tendency and a residual. The tendency describes an executive summary of a signal that encapsulates its notable characteristics while disregarding seemingly random, less interesting aspects. Building upon the Intrinsic Time Decomposition (ITD) and information-theoretical analysis, we introduce two alternative procedures for selecting the tendency from the ITD baselines. The first is based on the maximum extrema prominence, namely the maximum difference between extrema within each baseline. Specifically this method selects the tendency as the baseline from which an ITD step would produce the largest decline of the maximum prominence. The second method uses the rotations from the ITD and selects the tendency as the last baseline for which the associated rotation is statistically stationary. We delve into a comparative analysis of the information content and interpretability of the tendencies obtained by our proposed methods and those obtained through conventional low-pass filtering schemes, particularly the Hodrik–Prescott (HP) filter. Our findings underscore a fundamental distinction in the nature and interpretability of these tendencies, highlighting their context-dependent utility with emphasis in multi-scale signals. Through a series of real-world applications, we demonstrate the computational robustness and practical utility of our proposed tendencies, emphasizing their adaptability and relevance in diverse time series contexts.

Time series analysis↗

Planck 2018 results XI. Polarized dust foregrounds

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization in the quest for the curl-like B-mode polarization from primordial gravitational waves and the low-multipole E-mode polarization associated with the reionization of the Universe. We used the new Planck PR3 maps to characterize Galactic dust emission at high latitudes as a foreground to the CMB polarization and use end-to-end simulations to compute uncertainties and assess the statistical significance of our measurements. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for a set of six nested high-Galactic-latitude sky regions covering from 24 to 71% of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra, which for the largest sky region are α EE = -2.42 ± 0.02 and α BB = -2.54 ± 0.02, respectively. The spectra show that the TE correlation and E/B power asymmetry discovered by Planck extend to low multipoles that were not included in earlier Planck polarization papers due to residual data systematics. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we have determined the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component-separation procedure that is required for measuring the low-ℓ reionization CMB E-mode signal and detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified black-body emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean dust spectral index for dust polarization is β d P = 1.53±0.02. The difference between indices for polarization and total intensity is β d P -β d I = 0.05±0.03. By fitting multi-frequency cross-spectra between Planck data at 100, 143, 217, and 353 GHz, we examine the correlation of the dust polarization maps across frequency. We find no evidence for a loss of correlation and provide lower limits to the correlation ratio that are tighter than values we derive from the correlation of the 217- and 353 GHz maps alone. If the Planck limit on decorrelation for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then frequency decorrelation of dust polarization might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio r ≃ 0.01 at the recombination peak. However, the Planck sensitivity precludes identifying how difficult the component-separation problem will be for more ambitious experiments targeting lower limits on r.

79 ASTRONOMY AND ASTROPHYSICS↗

Rapid Patterning and Advanced Device Structures for Low Cost Manufacturable Crystalline Si IBC Cells

Recent record efficiency Si cells have had an interdigitated back contact (IBC) design where both positive and negative contacts are fabricated on the back side of the Si wafer. Many industry groups have investigated IBC cells but used photolithography to pattern the rear contacts although it is considered impractical for low cost, high volume manufacturing. Alternative patterning methods using lasers and mechanical masking have been utilized in fabrication of patterned regions in IBC Si solar cells to replace photolithography. Another exciting advance in Si solar cell technology has been replacing the high temperature diffusion of doped regions inside the crystalline Si wafer with low temperature deposition of very thin layers of doped amorphous hydrogenated Si (a-Si:H) on the Si wafer. This is called a Si heterojunction (HJ) device. The a-Si:H provides excellent passivation of surface defects and has produced the highest open circuit voltages (VOC) of any Si solar cell device structure. The Institute of Energy Conversion (IEC) demonstrated the first IBC-HJ solar cell combining these two strategies in 2007 and has fabricated IBC-HJ cells with 20% efficiency using three photolithography steps. Despite the demonstrated efficiency potential of this device by industry groups (>26%) there is no commercial production because of the challenges of patterning and processing the structure in an industrial environment. Our work sought to address that challenge. The objective of this 3 year project was to develop the processing for the IBC-HJ Si solar cell using lasers for patterning and contact formation instead of photolithography. Laser patterning enables rapid, contactless manufacturing of patterned regions. We intensively studied laser fired contacts, laser patterning of the a-Si multi-layer stacks and metal layers and the use of plasma shadow masks. After an exhaustive focus in the first year on the laser fired emitter (LFE) and contact (LFC), we were unable to achieve VOC greater than 660 mV, compared to the 720 mV required to meet our milestones. We discovered that an additional issue with our original IBC structure was the presence of an inversion layer at the back surface connecting the p and n regions which reduced the VOC and fill factor (FF). We developed innovative methods for characterizing the inversion layer. These two limitations in the original design lead to development of a new IBC-HJ process sequence and device structure which we called Plasma Masked Laser Processed (PMLP). It retained the original high efficiency features including manufacturability. The LFC was replaced with a standard n-type a-Si HJ contact. The inversion layer was eliminated by replacing the previous p-type stack in the gap with an n-type or SiN stack. The PMLP used laser ablation of a multi-layer a-Si stack followed by chemical etching to open the n-contact. It required deposition of a patterned stack through a mask in the plasma deposition chamber. This turned out to be a source of significant problems due to inevitable unwanted deposition ‘leakage’ under the mask. This formed a blocking contact on the emitter which significantly reduced Voc and FF. Several iterations in PMLP device structure and chemical etching steps resulted in a large improvements, e.g. the efficiency increasing from 3 to 15% and Voc from 450 to 660 mV but they were unable to completely eliminate it. The best IBC-HJ device we fabricated had only 15% efficiency. For perspective, our standard FHJ devices had 20% efficiency.

14 SOLAR ENERGY↗

Dust-trapping Vortices and a Potentially Planet-triggered Spiral Wake in the Pre-transitional Disk of V1247 Orionis

The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier. Here, we present ALMA 0.″04 resolution imaging of the pre-transitional disk of V1247 Orionis that reveals an asymmetric ring as well as a sharply confined crescent structure, resembling morphologies seen in theoretical models of vortex formation. The asymmetric ring (at 0.″17 = 54 au separation from the star) and the crescent (at 0.″38 = 120 au) seem smoothly connected through a one-armed spiral-arm structure that has been found previously in scattered light. We propose a physical scenario with a planet orbiting at ∼0.″3 ≈ 100 au, where the one-armed spiral arm detected in polarized light traces the accretion stream feeding the protoplanet. The dynamical influence of the planet clears the gap between the ring and the crescent and triggers two vortices that trap millimeter-sized particles, namely, the crescent and the bright asymmetry seen in the ring. We conducted dedicated hydrodynamics simulations of a disk with an embedded planet, which results in similar spiral-arm morphologies as seen in our scattered-light images. At the position of the spiral wake and the crescent we also observe {sup 12}CO(3-2) and H{sup 12}CO{sup +} (4-3) excess line emission, likely tracing the increased scale-height in these disk regions.

79 ASTRONOMY AND ASTROPHYSICS↗

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗

Technetium Stabilization in Low-Solubility Sulfide Phases: A Review

We report that technetium contamination remains a major environmental problem at nuclear reprocessing sites, e.g., the Hanford Site, Washington, USA. At these sites, Tc is present in liquid waste destined for immobilization in a waste form or has been released into the subsurface environment. The high environmental risk associated with Tc is due to its long half-life (214 000 years) and the mobility of the oxidized anionic species Tc(VII)O 4 - . Under reducing conditions, TcO 4 - is readily reduced to Tc(IV), which commonly exists as a relatively insoluble and therefore immobile, hydrous Tc-oxide (TcO 2 · n H 2 O). The stability of Tc(IV) sequestered as solid phases depends on the solubility of the solid and susceptibility to reoxidation to TcO 4 - , which in turn depend on the (biogeo)chemical conditions of the environment and/or nuclear waste streams. Unfortunately, the solubility of crystalline TcO 2 or amorphous TcO 2 ·H 2 O is still above the maximum contaminant level (MCL) established by the U.S. EPA (900 pCi/L), and the kinetics of TcO 2 oxidative dissolution can be on the order of days to years. In addition to oxygen, sulfur can form complexes that significantly affect the adsorption, solubility, and reoxidation potential of Tc, especially Tc(IV). The principal technetium sulfides are TcS 2 and Tc 2 S 7 , but much less is known about the mechanisms of formation, stabilization, and reoxidation of Tc-sulfides. A common assumption is that sulfides are less soluble than their oxyhydrous counterparts. Determination of the molecular structure of Tc 2 S 7 in particular has been hampered by the propensity of this phase to precipitate as an amorphous substance. Recent work indicates that the oxidation state of Tc in Tc 2 S 7 is Tc(IV), in apparent contradiction to its nominal stoichiometry. Technetium is relatively immobile in reduced sediments and soils, but in many cases the exact sink for Tc has not been identified. Experiments and modeling have demonstrated that both abiotic and biologic mechanisms can exert strong controls on Tc mobility and that Tc binding or uptake into sulfide phases can occur. These and similar investigations also show that extended exposure to oxidizing conditions results in transformation of sulfide-stabilized Tc(IV) to a Tc(IV)O 2 -like phase without formation of measurable dissolved TcO 4 - , suggesting a solid-state transformation in which Tc(IV)-associated sulfide is preferentially oxidized before the Tc(IV) cation. This transformation of Tc(IV)-sulfides to Tc(IV)-oxides may be the main process that limits remobilization of Tc as Tc(VII)O 4 - . The efficacy of the final waste form to retain Tc also strongly depends on the ability of oxidizing species to enter the waste and convert Tc(IV) to Tc(VII). Finally, many waste form designs are reducing (e.g., cementitious waste forms such as salt stone) and, therefore, attempt to restrict access of oxidizing species such that diffusion is the rate-limiting step in remobilization of Tc.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗