Search NASA⌕ Search

SEARCH · Search NASA

Results for “junction conditions”

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 19 records

Surface stress tensor and junction conditions on a rotating null horizon

The general form of the surface stress tensor of an infinitesimally thin shell located on a rotating null horizon is derived, when different interior and exterior geometries are joined there. Even though the induced metric on the surface must be the same approached from either side, the first derivatives of the metric need not be. Such discontinuities lead to a Dirac δ-distribution in the Einstein tensor localized on the horizon. For a general stationary axisymmetric geometry the surface stress tensor can be expressed in terms of two geometric invariants that characterize the surface, namely the discontinuities [κ] and [J] of the surface gravity κ and angular momentum density J. The Komar energy and angular momentum are given in coordinates adapted to the Killing symmetries, and the surface contributions to each determined in terms of [κ] and [J] . Guided by these, a simple modification of the original Israel junction conditions is verified directly from the Einstein tensor density to give the correct finite result for the surface stress, when the normal n to the surface is allowed to tend continuously to a null vector. The relation to Israel’s original junction conditions, which fail on null surfaces, is given. The modified junction conditions are suitable to the matching of a rotating “black hole” exterior to any interior geometry joined at the Kerr null horizon surface, even when the surface normal is itself discontinuous and the Barrabès-Israel formalism is also inapplicable. This joining on a rotating null horizon is purely of the matter shell type and does not contain a propagating gravitational shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonequilibrium effects on stability of hybrid stars with first-order phase transitions

The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely used Bardeen-Thorne-Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. Here, we apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that, like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.

79 ASTRONOMY AND ASTROPHYSICS↗

Two first-order phase transitions in hybrid compact stars: Higher-order multiplet stars, reaction modes, and intermediate conversion speeds

We study compact stars with hybrid equations of state consisting of a nuclear outer region and two nested quark phases, each separated from the lower density phase by a strong first-order phase transition. The stability of these models is determined by calculating their radial oscillation modes with different conversion rates between adjacent phases and hence junction conditions for the modes at the phase separation interface between them. In the case when the timescale of transition is faster than the period of oscillations, we recover the traditional stability criterion implying that ∂M/∂ρ c > 0 on the stable branch(es), where M is the mass and ρ c is the central density. In the opposite limit of slow conversion, we find stable stellar multiplets beyond triplets consisting of stars that are stable by the usual criterion plus slow-conversion (denoted by s) hybrid stars with ∂M/∂ρ c < 0 that are stabilized due to an alternative junction condition on the fluid displacement field at the interface reflecting the slow rate of conversion. We also study the properties of the reaction mode, the radial mode that only exists for stars with rapid (abbreviated by r) phase transitions, in stars with either two rapid phase transitions or alternating rapid and slow phase transitions for the first time. Here, the implications of alternative junction conditions are also examined, with these conditions generally being found to provide stability properties similar to those for a slow conversion rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Slowly rotating gravastars

We solve Einstein’s equations for slowly rotating gravitational condensate stars (gravastars) up to second order in the rotation by expanding about the spherically symmetric gravastar with de Sitter interior and Schwarzschild exterior matched at their common horizon. Requiring that the perturbations are finite on the null surface reduces the exterior geometry to that of a Kerr black hole, implying that a slowly rotating gravastar cannot be distinguished from a Kerr black hole by any measurement or observation restricted to the macroscopic spacetime exterior to the horizon. Here we determine the interior solution, the surface stress tensor, and the Komar mass and angular momentum localized on the slowly rotating horizon surface. With the interior equation of state fixed at p = –ρ, finite junction conditions on the null horizon surface necessarily lead to an interior solution with a singular core, where the perturbative expansion breaks down. Comparison to other models and implications for more rapidly rotating gravastars are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

The Effect of Grain Boundary Facet Junctions on Segregation and Embrittlement

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

36 MATERIALS SCIENCE↗

Graphite-Embedded High-Performance Insulated Metal Substrate for Wide-Bandgap Power Modules

Emerging wide-bandgap (WBG) semiconductor devices such as silicon carbide (SiC) metal–oxide semiconductor field-effect transistors (MOSFETs) and gallium nitride high-electron-mobility transistors can handle high power in reduced semiconductor areas better than conventional Si-based devices owing to superior material properties. With increased power loss density in a WBG-based converter and reduced die size in power modules, thermal management of power devices must be optimized for high performance. Here, we present a graphite-embedded insulated metal substrate (thermally-annealed-pyrolytic-graphite-embedded insulated metal substrate—IMSwTPG) designed for WBG power modules. Theoretical thermal performance analysis of graphite-embedded metal cores is presented, with design details for IMSwTPG with embedded graphite to replace a direct-bonded copper (DBC) substrate. The proposed IMSwTPG is compared with an aluminum nitride-based DBC substrate using finite-element thermal analysis for steady-state and transient thermal performance. The solutions’ thermal performances are compared under different coolant temperature and thermal loading conditions, and the proposed substrate's electrical performance is validated with static and dynamic characterization. Using graphite-embedded substrates, junction-to-case thermal resistance of SiC MOSFETs can be reduced up to 17%, and device current density can be increased by 10%, regardless of the thermal management strategy used to cool the substrate. Reduced transient thermal impedance of up to 40% of dies owing to increased heat capacity is validated in transient thermal simulations and experiments. The half-bridge power module's electrical performance is evaluated for on -state resistance, switching performance, and switching loss at three junction temperature conditions. The proposed substrate solution has minimal impact on conduction and switching performance of SiC MOSFETs.

42 ENGINEERING↗

Investigation of Off-Stoichiometry in Ternary Nitrides by EDS and HRTEM

We demonstrate an inverted metamorphic multijunction (IMM) photovoltaic cell comprising lattice-mismatched 1.2 eV AlGaInAs and 1.0 eV GaInAs junctions optimized for high-temperature thermophotovoltaic (TPV) applications. This device differs from traditional IMM solar cells because the mismatched junctions are grown at a single lattice constant. This architecture enables removal of the compositionally graded buffer which otherwise filters light from the junctions below and absorbs sub-bandgap light via free-carrier absorption. Sub-bandgap absorption dramatically reduces the efficiency of TPV systems using high reflectivity cells to enable bandedge spectrum filtering. Three components required development to enable this device: 1) a lattice-mismatched 1.2 eV AlGaInAs junction, 2) a metamorphic contact layer grown after the graded buffer, and 3) a transparent tunnel junction that sits in front of the 1.0 eV GaInAs junction. Growth conditions that minimize oxygen defect incorporation maximize AlGaInAs cell quality, enabling a 0.41 V bandgap open circuit voltage offset at 22 mA/cm2 under AM1.5D. A mismatched GaInAs:Se layer is developed as a low resistance contact. Lastly we develop a GaAsSb:C/GaInP:Se tunnel junction suitable for high-power densities with more transparency than the GaAsSb:C/GaInP:Se structure used in past IMM cells. We characterize the tandem device under a high-intensity spectrum that approximates the emission from a 2150 °C blackbody radiator and deduce a projected ideal TPV efficiency of 39.9% at ~30% of the blackbody irradiance, and 36% ideal TPV efficiency under the full 118 W/cm2 irradiance. Improvements to the back-surface reflectivity and series resistance are expected to increase the ideal TPV efficiency well above 40%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inverted Metamorphic AlGaInAs/GaInAs Tandem Thermophotovoltaic Cell Designed for Thermal Energy Grid Storage Application

We demonstrate an inverted metamorphic multijunction (IMM) photovoltaic cell comprising lattice-mismatched 1.2 eV AlGaInAs and 1.0 eV GaInAs junctions optimized for high-temperature thermophotovoltaic (TPV) applications. This device differs from traditional IMM solar cells because the mismatched junctions are grown at a single lattice constant. This architecture enables removal of the compositionally graded buffer that otherwise filters light from the junctions below and absorbs sub-bandgap light via free-carrier absorption. Sub-bandgap absorption dramatically reduces the efficiency of TPV systems using high reflectivity cells to enable band edge spectrum filtering. Three components required development to enable this device: (1) a lattice-mismatched 1.2 eV AlGaInAs junction, (2) a metamorphic contact layer grown after the graded buffer, and (3) a transparent tunnel junction that sits in front of the 1.0 eV GaInAs junction. Growth conditions that minimize oxygen defect incorporation maximize AlGaInAs cell quality, enabling a 0.41 V bandgap open circuit voltage offset at 22 mA/cm2 under AM1.5D. A mismatched GaInAs:Se layer is developed as a low resistance contact. Lastly, we develop a GaAsSb:C/GaInP:Se tunnel junction suitable for high-power densities with more transparency than the GaAsSb:C/GaInAs:Se structure used in past IMM cells. We characterize the tandem device under a high-intensity spectrum that approximates the emission from a 2150°C blackbody radiator and deduce a projected ideal TPV efficiency of 39.9% at ~30% of the blackbody irradiance and 36% ideal TPV efficiency under the full 118 W/cm2 irradiance. Improvements to the back-surface reflectivity and series resistance are expected to increase the ideal TPV efficiency well above 40%.

25 ENERGY STORAGE↗

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

Here we present a series of simple shear numerical simulations of dynamic recrystallization of two-phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (V p ) and fast S wave (V s ) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest V p quickly aligns with the maximum c-axis orientation. At the same time, the maximum c-axis development reduces V s in this orientation. When water is present, the developed maximum c-axis orientation is less intense, which results in lower V p and V s . At high percentage of water, V p does not align with the maximum c-axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of V s while V p remains close to the value for purely solid ice. These results suggest that the decrease in V s observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c-axis-dependent orientation and the magnitude of seismic anisotropy.

58 GEOSCIENCES↗

Diverse simulations of time-resolved photoluminescence in thin-film solar cells: A SnO 2 /CdSe y Te 1-y case study

Time-resolved photoluminescence (TRPL) is widely used to measure carrier lifetime in thin-film solar cell absorbers. However, the injection dependence of data and frequent non-exponential decay shapes complicate the interpretation. Here, we develop a numerical model to simulate injection-dependent TRPL measurements in a SnO 2 /CdSe y Te 1-y case study solar cell structure, considering parameters of interest to researchers in industry and academia. Previous simulations have shown that, in low injection, excess electrons and holes injected by the laser pulse are rapidly separated in the electric field formed by the pn junction. As a result, at early times, the PL signal can decay faster than the Shockley-Read-Hall lifetime in the absorber bulk (t_bulk). Prior simulations have shown that the charge stored in the junction can slowly leak out to affect decays at late times. However, it has not been clear if and to what degree charge storage can affect the slopes extracted from TRPL decays - t_2 - commonly cited as the TRPL-measured lifetime. Here, we show that charge storage can, in some cases, result in t_2 values that substantially overestimate t_bulk. Previous simulations indicate that high injection conditions can screen the junction field and minimize charge separation. Here, we show that continued injection increases can drive down t_2 below t_bulk as radiative recombination becomes dominant. We catalog charge storage and radiative recombination impacts for a diverse set of material parameters and compare results to double-heterostructure models.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,S↗

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗

Formation and Evolution of Metallocene Single-Molecule Circuits with Direct Gold-π Links

Single molecule circuits with group 8 metallocenes are formed without additional linker groups in Scanning Tunnelling Microscope-based break junction (STMBJ) measurements at cryogenic and room temperature conditions with gold (Au) electrodes. We investigate the nature of this direct gold-π binding motif and its affect on molecular conductance and persistence characteristics during junction evolution. The measurement technique under cryogenic conditions tracks molecular plateaus through the full cycle of extension and compression. Analysis reveals that junction persistence when the metal electrodes are pushed together correlates with whether electrodes are locally sharp or blunt, suggesting distinct scenarios for metallocene junction formation and evolution. The top and bottom surfaces of the "barrel" shaped metallocenes present the electron-rich π system of cyclopentadienyl rings which interacts with the gold electrodes in two distinct ways. An undercoordinated gold atom on a sharp tip forms a donor-acceptor bond to a specific carbon atom in the ring. However, a small, flat patch on a dull tip can bind more strongly to the ring as a whole through Van der Waals interactions. Density functional theory (DFT) based calculations of model electrode structures provide an atomic-scale picture of these scenarios, demonstrating the role of these bonding motifs during junction evolution and showing that the conductance is relatively independent of tip atomic-scale structure. Furthermore, the non-specific interaction of the cyclopentadienyl rings with the electrodes enables extended conductance plateaus, a mechanism distinct from that identified for the more commonly studied, rod-shaped organic molecular wires.

25 ENERGY STORAGE↗

Operando Temperature Measurements of Photovoltaic Laser Power Converter Devices Under Continuous High-Intensity Illumination

Photovoltaic devices that operate under extremely high irradiances, such as laser power converters (LPCs), may also operate at elevated temperatures, even under active cooling, as the result of large temperature gradients. We demonstrate the operation of GaAs LPC devices under orders of magnitude of irradiances up to 150 W/cm 2 in a monochromatic laser simulator with an active cooling stage. The steady-state open-circuit voltage (V OC ) as a function of irradiance is known to droop at high irradiance as the result of junction heating, but the junction temperature can be difficult to measure by conventional methods. Fast, transient V OC measurements under these extreme operating conditions are used here to determine the junction temperature. Empirical V OC temperature coefficients of the devices at each irradiance of interest are determined and used as an integral part of this technique. We show that the thermal design of different LPC devices strongly affects the operating temperature of the junctions. Knowledge of the operating temperature can be a strong tool for understanding the nature of loss mechanisms and improving the design for the operation of photovoltaic LPCs at high irradiances. Furthermore, this technique can be used for laboratory devices during initial design as well as to characterize mass-produced and packaged devices for quality control.

14 SOLAR ENERGY↗

Development of a Robust Reference Electrode in Aggressive Chemical and Radiation Environments in the Hanford Waste Tanks

The Hanford site stores more than 200 million liters of radioactive and chemically hazardous wastes from the production of weapons materials. The wastes are stored in 177 underground carbon-steel storage tanks, separated between 149 single shell tanks (SSTs) and 28 double shell tanks (DSTs). The DSTs provide critical retrieval and interim storage before the waste is vitrified in the Waste Treatment and Isolation Plant (WTP). The tanks will need to remain in-service far beyond the initial 40-year design life, and effective corrosion control practices must remain in force to extend the tanks’ lifespans. This effort includes direct measurements of corrosion rate (e.g., ultrasonic measurements and corrosion coupons) and electrochemical processes (e.g., linear polarization measurements and open circuit potential measurements). The Hanford site began monitoring the corrosion potential in select DSTs in 2008. Of the 45 reference electrodes that have been installed, 29 have failed and 6 others provided unreliable results. DOE-EM is supporting a 3-year program to develop a chemical and radiation resistant reference electrode for application in the Hanford tanks. The first year of the program focused on understanding the failure mechanism for the reference electrodes and identification of candidate construction materials that would mitigate degradation of the electrodes in the waste environment. During the second year of the program, the objectives were to: 1) test candidate materials under simulated waste conditions, 2) design components that will extend the service life of the electrode, 3) fabricate materials for prototype reference electrodes, and 4) assemble prototype reference electrodes for accelerated testing. The reference electrode is constructed of four principal parts: 1) junction, 2) casing, 3) inner chamber backfill materials, and 4) the sensing wire. Principally, improvements of the junction, casing, and inner chamber backfill materials are being pursued. The junction material at the interface between the waste and the inner chamber of the reference electrode was identified as a critical component in the failure of the reference electrodes. Nine candidate replacement junction materials were tested under simulated waste conditions to evaluate permeation rate. These materials included a variety of polymeric and ceramic materials, some of which were 3-D printed. Thus far, porous polyvinylidene fluoride materials have performed satisfactorily and are being considered for prototype development. The commercial electrode casing materials in general have performed well. Additionally, 3-D printing of chemically and mechanically stable materials is being investigated as a means for further improvement in fabrication consistency. SRNL has also investigated altering the reference electrode design to extend the service life. The new design of the interior of the reference electrode casing creates a longer, more tortuous path between the junction material and the electrode sensing wire. A finite element model was used to optimize the design without adversely impacting the circuit resistance of the electrode during the measurements, thus preserving the measurement accuracy while enhancing the service life. The inner chamber back fill materials are also critical to the performance of the reference electrode. Materials that are resistant to intruding tank waste and provide a conductive path to the sensing wire were investigated. Gel and powder materials that are interspersed with a conductive chloride bearing material were tested for their influence on diffusion and electrode resistance. All the investigated materials and components will be assembled, with collaboration from commercial vendors, to fabricate the initial prototypes. Accelerated testing of the prototypes will be initiated in Year 2 of the program and will be completed in Year 3. A recommendation on the materials of construction and the design of the new robust reference electrode will be presented to the Hanford tank farm facility.

Sykes, Kiana [Savannah River National Laboratory (↗

Measuring the Stress Factors for Photovoltaic (PV) Backsheet Degradation

Back sheet failure has resulted in power loss and large-scale recall of photovoltaic modules, resulting in billions of dollars in lost revenue. The light exposure on the backside of a photovoltaic module comes primarily from reflected light which alters the distribution of natural sunlight. Because of this, modelling the backside exposure and duplicating the exposure is much more difficult than modeling the frontside exposure. This project aims to study how various back sheets and junction box materials degrade under different conditions and to develop Python code to help model and predict degradation. The stress factors for back-sheet degradation must be quantified to extrapolate accelerated stress tests to the field. Test samples were placed in the A3, A4, and A5 conditions, as defined in IEC 62788-7-2, to assess the temperature and humidity dependence of ultraviolet (UV) induced degradation. We are utilizing a custom chamber with exposure from 0.5 UV-suns to 5 UV-suns to understand the dependence of degradation on light intensity. A group of samples put in the A3 condition had glass filters with 50% UV cut-offs of 320 nm, 335 nm, and 360 nm to assess the wavelength dependence of UV degradation. All this data is necessary to assess the impact of non-standard UV light exposure. The material evaluation tests include gloss measurements, attenuated total internal reflectance Fourier transform infrared spectroscopy (ATR-FTIR), UV-visible reflectance/transmittance utilizing a Cary Ci7000 spectrophotometer, and a nano-indenter for surface hardness and modulus measurements. Alongside the experimental work, there is a computational effort using raytracing and Python open-source tools in PVDeg , PVLib, and Bifacial_Radiance. This code will create specific exposure scenarios and enable the evaluation of chamber degradation relative to field degradation. Equation 1 is a strawman equation used to model degradation on the backside of a PV module. We will create simplified code, based on the results of ray-tracing calculations, which uses a view factor approach to provide fast calculations for the most common exposure scenarios.

14 SOLAR ENERGY↗

Optimizing energy yield of monolithic perovskite/silicon tandem solar cells in real-world Conditions: The impact of luminescent coupling

Efficient light management is key to maximizing power conversion efficiency (PCE) in monolithic perovskite/silicon tandem solar cells. Achieving peak efficiency requires closely matched current generation in all junctions, especially in integrated configurations. However, real-world conditions vary significantly due to factors such as sunlight spectrum, diffuse-to-direct sunlight ratio, angular distribution of light, subcell temperature coefficients, and ground reflection. This study introduces a comprehensive optical and device simulation to optimize perovskite/silicon tandem cells, considering experimental luminescent coupling (LC) efficiency and its dependence on working conditions, alongside variations in radiative recombination, effect of temperature on absorptivity spectra, and cloud cover. Our results show potential energy yield improvements of up to 1.4 % with LC, based on current perovskite radiative recombination records, and up to 4 % with direct bandgap materials. Although radiative recombination's dependence on excitation intensity reduces output power and requires thicker absorbers, LC compensates for these losses. LC also lowers the optimized bandgap for the perovskite top cell from 1.72 eV to 1.64–1.68 eV, or even lower in regions with redshifted irradiance. Additionally, optimization revealed that thinner silicon bottom cells require a higher perovskite top cell bandgap, impacting the balance between fabrication cost and cell stability.

14 SOLAR ENERGY↗

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗