Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Evaluation of exchange-correlation functionals with multiple-shock conductivity measurements in hydrogen and deuterium at the molecular-to-atomic transition

The temperature ( T ) and density (ρ) conditions at which hydrogen undergoes a molecular-to-atomic (MA) transition is crucial to our understanding of the gas-giant planets such as Jupiter and Saturn. First-principles (FP) calculations suggest that this transition is coincident with metallization and acts as a catalyst for hydrogen-helium demixing, which has significant consequences for models of planetary interiors. Prediction of this transition boundary has proven to be difficult using FP methods. In particular, detailed comparisons of finite temperature density functional theory (FT-DFT) calculations of the MA transition in both the high- T , low-ρ regime, where the transition is largely T driven, and the low- T , high-ρ regime, where the transition is largely ρ driven, suggest that the transition is very sensitive to the exchange-correlation (xc) functional used in the calculation. Here we present a detailed comparison of previous multiple-shock electrical conductivity measurements with FT-DFT calculations employing various xc functionals to probe a regime where both T and ρ play an important role in the transition. The measurement results are found to be inconsistent with the semilocal xc functional PBE and are in much better agreement with the nonlocal xc functionals vdW-DF1 and vdW-DF2. Furthermore, we show that the inconsistency with PBE likely stems from pressure errors associated with the PBE xc functional, resulting in calculated pressures that are too low at these T and ρ conditions. Together with previous comparisons at high- T , low-ρ and low- T , high-ρ these results provide a consistent picture for the MA transition over a wide T and ρ range. Here, this picture may also provide insight into differences in experimental observations of the metallization of liquid hydrogen and deuterium in the low- T regime.

08 HYDROGEN↗
DOE OSTI2015

A Discontinuous Petrov-Galerkin Methodology for Adaptive Solutions to the Incompressible Navier-Stokes Equations

The discontinuous Petrov-Galerkin methodology with optimal test functions (DPG) of Demkowicz and Gopalakrishnan [18, 20] guarantees the optimality of the solution in an energy norm, and provides several features facilitating adaptive schemes. Whereas Bubnov-Galerkin methods use identical trial and test spaces, Petrov-Galerkin methods allow these function spaces to differ. In DPG, test functions are computed on the fly and are chosen to realize the supremum in the inf-sup condition; the method is equivalent to a minimum residual method. For well-posed problems with sufficiently regular solutions, DPG can be shown to converge at optimal rates—the inf-sup constants governing the convergence are mesh-independent, and of the same order as those governing the continuous problem [48]. DPG also provides an accurate mechanism for measuring the error, and this can be used to drive adaptive mesh refinements. We employ DPG to solve the steady incompressible Navier-Stokes equations in two dimensions, building on previous work on the Stokes equations, and focusing particularly on the usefulness of the approach for automatic adaptivity starting from a coarse mesh. We apply our approach to a manufactured solution due to Kovasznay as well as the lid-driven cavity flow, backward-facing step, and flow past a cylinder problems.

adaptive finite elements↗
DOE OSTI2014

Optically tunable acoustic wave band-pass filter

The acoustic properties of a hybrid composite that exhibits both photonic and phononic behavior are investigated numerically with finite-element and finite-difference time-domain simulations. The structure is constituted of a periodic array of photonic resonant cavities embedded in a background superlattice. The resonant cavities contain a photo-elastic chalcogenide glass that undergoes atomic-scale structural reorganization when irradiated with light having energy close to its band-gap. Photo-excitation of the chalcogenide glass changes its elastic properties and, consequently, augments the acoustic transmission spectrum of the composite. By modulating the intensity of light irradiating the hybrid photonic/phononic structure, the position and spectral width of phonon passing-bands can be controlled. This demonstration offers the technological platform for optically-tunable acoustic wave band-pass filters.

36 MATERIALS SCIENCE↗
DOE OSTI2019

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

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

58 GEOSCIENCES↗
DOE OSTI2016

SU-F-I-50: Finite Element-Based Deformable Image Registration of Lung and Heart

Purpose: Photodynamic therapy (PDT) is used after surgical resection to treat the microscopic disease for malignant pleural mesothelioma and to increase survival rates. Although accurate light delivery is imperative to PDT efficacy, the deformation of the pleural volume during the surgery impacts the delivered light dose. To facilitate treatment planning, we use a finite-element-based (FEM) deformable image registration to quantify the anatomical variation of lung and heart volumes between CT pre-(or post-) surgery and surface contours obtained during PDT using an infrared camera-based navigation system (NDI). Methods: NDI is used during PDT to obtain the information of the cumulative light fluence on every cavity surface point that is being treated. A wand, comprised of a modified endotrachial tube filled with Intralipid and an optical fiber inside the tube, is used to deliver the light during PDT. The position of the treatment is tracked using an attachment with nine reflective passive markers that are seen by the NDI system. Then, the position points are plotted as three-dimensional volume of the pleural cavity using Matlab and Meshlab. A series of computed tomography (CT) scans of the lungs and heart, in the same patient, are also acquired before and after the surgery. The NDI and CT contours are imported into COMSOL Multiphysics, where the FEM-based deformable image registration is obtained. The NDI and CT contours acquired during and post-PDT are considered as the reference, and the Pre-PDT CT contours are used as the target, which will be deformed. Results: Anatomical variation of the lung and heart volumes, taken at different times from different imaging devices, was determined by using our model. The resulting three-dimensional deformation map along x, y and z-axes was obtained. Conclusion: Our model fuses images acquired by different modalities and provides insights into the variation in anatomical structures over time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗
DOE OSTI2018

Study of Dissipative Losses in AC-Biased Mo/Au Bilayer Transition-Edge Sensors

We are developing kilo-pixel arrays of transition-edge sensors (TESs) for the X-ray Integral Field Unit on ESA’s Athena observatory. Previous measurements of AC-biased Mo/Au TESs have highlighted a frequency-dependent loss mechanism that results in broader transitions and worse spectral performance compared to the same devices measured under DC bias. In order to better understand the nature of this loss, we are now studying TES pixels in different geometric configurations. We present measurements on devices of different sizes and with different metal features used for noise mitigation and X-ray absorption. Our results show how the loss mechanism is strongly dependent upon the amount of metal in close proximity to the sensor and can be attributed to induced eddy current coupling to these features. We present a finite element model that successfully reproduces the magnitude and geometry dependence of the losses. Using this model, we present mitigation strategies that should reduce the losses to an acceptable level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

A Modified Johnson–Cook Constitutive Model for the Compressive Flow Behaviors of the SnSbCu Alloy at High Strain Rates

The SnSbCu alloy is widely used as the material for the main bearing in low-speed marine engines, and an accurate constitutive model is the foundation for studying the frictional behaviors of bearings. In this work, the compressive flow behaviors of the SnSbCu alloy were considered under the different strain rates (1000-5000 s{sup −1}) and temperatures (20-110 °C) by quasi-static and split-Hopkinson pressure bar dynamic compression tests. First, the original Johnson–Cook model was used to describe the constitutive relation of the SnSbCu alloy at high strain rates, and the results predicted by the original model showed relatively large errors compared with the experimental results since the coupled effect of temperature and strain rate was omitted. Then, a modified Johnson–Cook constitutive model was developed to describe the compressive flow behaviors of the SnSbCu alloy, and the results predicted by this modified model agreed well with the experimental data. Moreover, a finite element analysis was also conducted to verify the accuracy of the modified Johnson–Cook model.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SAMBA: Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos

A new arbitrary Polynomial Chaos (aPC) method is presented for moderately high-dimensional problems characterised by limited input data availability. The proposed methodology improves the algorithm of aPC and extends the method, that was previously only introduced as tensor product expansion, to moderately high-dimensional stochastic problems. The fundamental idea of aPC is to use the statistical moments of the input random variables to develop the polynomial chaos expansion. This approach provides the possibility to propagate continuous or discrete probability density functions and also histograms (data sets) as long as their moments exist, are finite and the determinant of the moment matrix is strictly positive. For cases with limited data availability, this approach avoids bias and fitting errors caused by wrong assumptions. In this work, an alternative way to calculate the aPC is suggested, which provides the optimal polynomials, Gaussian quadrature collocation points and weights from the moments using only a handful of matrix operations on the Hankel matrix of moments. It can therefore be implemented without requiring prior knowledge about statistical data analysis or a detailed understanding of the mathematics of polynomial chaos expansions. The extension to more input variables suggested in this work, is an anisotropic and adaptive version of Smolyak's algorithm that is solely based on the moments of the input probability distributions. It is referred to as SAMBA (PC), which is short for Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos. It is illustrated that for moderately high-dimensional problems (up to 20 different input variables or histograms) SAMBA can significantly simplify the calculation of sparse Gaussian quadrature rules. SAMBA's efficiency for multivariate functions with regard to data availability is further demonstrated by analysing higher order convergence and accuracy for a set of nonlinear test functions with 2, 5 and 10 different input distributions or histograms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Using ARM Measurements to Understand and Reduce the Double ITCZ Biases in the Community Atmospheric Model

1. Understanding of the observed variability of ITCZ in the equatorial eastern Pacific. The annual mean precipitation in the eastern Pacific has a maximum zonal band north of the equator in the ITCZ where the maximum SST is located. During the boreal spring (referring to February, March, and April throughout the present paper), because of the accumulated solar radiation heating and oceanic heat transport, a secondary maximum of SST exists in the southeastern equatorial Pacific. Associated with this warm SST is also a seasonal transitional maximum of precipitation in the same region in boreal spring, exhibited as a weak double ITCZ pattern in the equatorial eastern Pacific. This climatological seasonal variation, however, varies greatly from year to year: double ITCZ in the boreal spring occurs in some years but not in other years; when there a single ITCZ, it can appear either north, south or at the equator. Understanding this observed variability is critical to find the ultimate cause of the double ITCZ in climate models. Seasonal variation of ITCZ south of the eastern equatorial Pacific: By analyzing data from satellites, field measurements and atmospheric reanalysis, we have found that in the region where spurious ITCZ in models occurs, there is a “seasonal cloud transition” — from stratocumulus to shallow cumulus and eventually to deep convection —in the South Equatorial Pacific (SEP) from September to April that is similar to the spatial cloud transition from the California coast to the equator. This seasonal transition is associated with increasing sea surface temperature (SST), decreasing lower tropospheric stability and large-scale subsidence. This finding of seasonal cloud transition points to the same source of model errors in the ITCZ simulations as in simulation of stratocumulus-cumulus-deep convection transition. It provides a test for climate models to simulate the relationships between clouds and large-scale atmospheric fields in a region that features a spurious double Inter-tropical Convergence Zone (ITCZ) in most models. This work is recently published in Yu et al. (2016). Interannual variation of ITCZ south of the eastern equatorial Pacific: By analyzing data from satellites, field measurements and atmospheric reanalysis, we have characterized the interannual variation of boreal spring precipitation in the eastern tropical Pacific and found the cause of the observed interannual variability. We have shown that ITCZ in this region can occur as a single ITCZ along the Equator, single ITCZ north of the Equator, single ITCZ south of the Equator, and double ITCZ on both sides of the Equator. We have found that convective instability only plays a secondary role in the ITCZ interannual variability. Instead, the remote impact of the Pacific basin-wide SST on the horizontal gradient of surface pressure and wind convergence is the primary driver of this interannual variability. Results point to the need to include moisture convergence in convection schemes to improve the simulation of precipitation in the eastern tropical Pacific. This result has been recently submitted for publication (Yu and Zhang 2016). 2. Improvement of model parameterizations to reduce the double ITCZ bias We analyzed the current status of climate model performance in simulating precipitation in the equatorial Pacific. We have found that the double ITCZ bias has not been reduced in CMIP5 models relative to CMIP4 models. We have characterized the dynamic structure of the common bias by using precipitation, sea surface temperature, surface winds and sea-level. Results are published in Zhang et al. (2015): Since cumulus convection plays a significant role in the double ITCZ behavior in models, we have used measurements from ARM and other sources to carry out a systematic analysis of the roles of shallow and deep convection in the CAM. We found that in both CAM4 and CAM5, when the intensity of deep convection decreases as a result of parameterization change, the intensity of shallow convection increases, leading to very different changes in precipitation partitions but little change in the total precipitation. The different precipitation partition however can manifest themselves in other measures of model performances including temperature and humidity. This study points to the need to treat model physical parameterizations as integrated system rather than individual components. Results from this study are published in Wang and Zhang (2013). Since shallow convection interacts with the deep convection scheme and surface turbulence to trigger the double ITCZ, we studied methods to improve the shallow convection scheme in climate models. We investigated the bulk budgets of the vertical velocity and its parameterization in convective cores, convective updrafts, and clouds by using large-eddy simulation (LES) of four shallow convection cases including one from ARM. We proposed optimal forms of the Simpson and Wiggert equation to calculate the vertical velocity in bulk mass flux convection schemes for convective cores, convective updrafts, and convective clouds as parameterization schemes. The new scheme is published in Wang and Zhang (2014). By using long-term radar-based ground measurements from ARM, we derived a scale-aware inhomogeneity parameterization of cloud liquid water in climate models. We found a relationship between the inhomogeneity parameter and the model grid size as well as atmospheric stability. This relationship is implemented in the CESM to describe the subgrid-scale cloud inhomogeneity. Relative to the default CESM with the finite-volume dynamic core at 2-degree resolution, the new parameterization leads to smaller cloud inhomogeneity and larger cloud liquid-water path in high latitudes, and the opposite effect in low latitudes, with the regional impact on shortwave cloud radiation effect of up to 10 W/m 2 . This is due to both the smaller (larger) grid size in high (low) latitudes in the longitude-latitude grid setting of CESM and the more stable (unstable) atmosphere. This parameterization is expected lead to more realistic simulation of tropical precipitation in high-resolution models. Results from this study are reported in Xie and Zhang (2015).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Quarkonium in-medium properties from realistic lattice NRQCD

Here, we present the final results of our high statistics study on the properties of bottomonium and charmonium at finite temperature. We focus on the temperature range around the crossover transition 150 ≤ T ≤ 410MeV, relevant for current heavy ion collision experiments. The QCD medium degrees of freedom, which consist of dynamical u,d, and s quarks and gluons are captured by realistic state-of-the art (m π ≈ 161MeV) lattice QCD simulations of the HotQCD collaboration. For the heavy quarks we deploy the non-relativistic effective field theory of QCD, NRQCD. The in-medium properties of quarkonium are deduced from their spectral functions, which are reconstructed using improved and novel Bayesian approaches. Through a systematic analysis we shed light on the origin of the discrepancies in melting temperatures previously reported in the literature, showing that they are owed to underestimated methods uncertainties of the deployed spectral reconstructions. Our simulations corroborate a picture of sequential in-medium modification, ordered according to the vacuum binding energy of the states. As a central quantitative result, our study reveals how the mass of the heavy quarkonium ground state reduces as temperature increases. The observed spectral modifications are interpreted in the light of, and compared to previous studies based on the complex lattice potential for heavy quarkonium. Furthermore for the first time we provide a robust picture of in-medium heavy quarkonium modification in the quark-gluon plasma consistent among different non-relativistic methods. We also critically discuss the perspectives for improving on these results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Antisites in III-V semiconductors: Density functional theory calculations

Density functional based simulation, corrected for finite size effects, is used to investigate systematically the formation of antisite defects in III-V semiconductors (III = Al, Ga, and In and V = P, As, and Sb). Different charge states are modelled as a function of the Fermi level and under different growth conditions. The formation energies of group III antisites (III{sub V}{sup q}) decrease with increasing covalent radius of the group V atom though not group III radius, whereas group V antisites (V{sub III}{sup q}) show a consistent decrease in formation energies with increase in group III and group V covalent radii. In general, III{sub V}{sup q} defects dominate under III-rich conditions and V{sub III}{sup q} under V-rich conditions. Comparison with equivalent vacancy formation energy simulations shows that while antisite concentrations are always dominant under stoichiometric conditions, modest variation in growth or doping conditions can lead to a significantly higher concentration of vacancies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Manufacture of gradient micro-structures of magnesium alloys using two stage extrusion dies

This paper aims to manufacture magnesium alloy metals with gradient micro-structures using hot extrusion process. The extrusion die was designed to have a straight channel part combined with a conical part. Materials pushed through this specially-designed die generate a non-uniform velocity distribution at cross sections inside the die and result in different strain and strain rate distributions. Accordingly, a gradient microstructure product can be obtained. Using the finite element analysis, the forming temperature, effective strain, and effective strain rate distributions at the die exit were firstly discussed for various inclination angles in the conical die. Then, hot extrusion experiments with a two stage die were conducted to obtain magnesium alloy products with gradient micro-structures. The effects of the inclination angle on the grain size distribution at cross sections of the products were also discussed. Using a die of an inclination angle of 15°, gradient micro-structures of the grain size decreasing gradually from 17 μm at the center to 4 μm at the edge of product were achieved.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effects of cross-section on mechanical properties of Au nanowire

The aim of this paper is study of the effects of multiple cross-section of Au nanowire on mechanical properties. Different cross-section models of Au nanowires including circular, hexagonal, pentagonal and rectangular were simulated by finite element modeling using ABAQUS. In this study, the bending technique was applied so that both ends of the model were clamped with mid-span under loading condition. The cross-sections had the length of 400 nm and the diameter of 40 nm, except the circular cross-section while the rest of the cross-sections had an equivalent diameter. Von Misses stresses distribution were used to define the stress distribution in the cross-section under loading condition, and elastic deformation was analyzed by the beam theory. The results disclosed that the circular and the rectangular models had highest and lowest strengths against plastic deformation, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

A Generalized Backward Euler algorithm for the numerical integration of a viscous breakage model

This paper discusses the formulation and the numerical performance of a fully implicit algorithm used to integrate a rate-dependent model defined within a breakage mechanics framework. For this purpose, a Generalized Backward Euler (GBE) algorithm has been implemented according to two different linearization strategies: The former is derived by a direct linearization of the constitutive equations, while the latter introduces rate effects through a consistency parameter. Additionally, the accuracy and efficiency of the GBE algorithm have been investigated by (1) performing material point analyses and (2) solving initial boundary value problems. In both cases, the overall performance of the underlying algorithm is inspected for a range of loading rates, thus simulating comminution from slow to fast dynamic problems. As the viscous response of the breakage model can be recast through a viscous nucleus function, the presented algorithm can be considered as a general framework to integrate constitutive equations relying on the overstress approach typical of Perzyna-like viscoplastic models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prediction of softening kinetics and recrystallization texture in non-isothermally annealed bulged tubes using CPFEM and CA models

A hierarchically coupled cellular automata (CA) model, crystal plasticity finite element method (CPFEM), and thermal finite element (FE) model is developed to predict the softening kinetics of the bulged steel tube during non-isothermal annealing. Through the developed model, the kinetics of softening mechanisms including static recovery (SRV) and static recrystallization (SRX), as well as the recrystallization texture are predicted. Later, the Johnson-Mehl-Avrami-Kohnogorov (JMAK) model based on the predicted SRX data is developed to interpret the recrystallization behavior of the material. To perform this study, diverse experimental tests including tube hydroforming (THF), annealing, uniaxial tensile test, hardness test, as well as microstructure observations through optical microscopy and Electron Backscatter Diffraction (EBSD) tests on steel tube are performed. The obtained experimental data are utilized to calibrate and verify the implemented CPFEM model for simulation of THF process, thermal FE model for prediction of the local temperature over annealing time, and CA algorithm for modeling of the softening kinetics and texture evolution throughout the annealing process. The study shows that the predicted deformation characteristics, softening kinetics, recrystallization texture and temperature profile during non-isothermal annealing are in good agreement with experimental data. During the annealing process, a total of four stages for the kinetics of softening mechanisms is observed: No softening; SRV only; SRV dominant; and SRX dominant. During the progress of SRX, the behavior of recrystallization is abruptly changed, confirming that two different mechanisms are controlling the kinetics of transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Topological supercurrents interaction and fluctuations in the multiterminal Josephson effect

We study the Josephson effect in the multiterminal junction of topological superconductors. We use the symmetry-constrained scattering matrix approach to derive band dispersions of emergent subgap Andreev bound states in a multidimensional parameter space of superconducting phase differences. We find distinct topologically protected band crossings that serve as monopoles of finite Berry curvature. Particularly, in a four-terminal junction the admixture of 2π and 4π periodic levels leads to the appearance of finite-energy Majorana-Weyl nodes. This topological regime in the junction can be characterized by a quantized nonlocal conductance that measures the Chern number of the corresponding bands. In addition, we calculate current-phase relations, variance, and cross correlations of topological supercurrents in multiterminal contacts and discuss the universality of these transport characteristics. At the technical level these results are obtained by integrating over the group of a circular ensemble that describes the scattering matrix of the junction. Furthermore, we briefly discuss our results in the context of observed fluctuations of the gate dependence of the critical current in topological planar Josephson junctions and comment on the possibility of parity measurements from the switching current distributions in multiterminal Majorana junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

MONTE CARLO SIMULATIONS OF THE PHOTOSPHERIC EMISSION IN GAMMA-RAY BURSTS

We studied the decoupling of photons from ultra-relativistic spherically symmetric outflows expanding with constant velocity by means of Monte Carlo simulations. For outflows with finite widths we confirm the existence of two regimes: photon-thick and photon-thin, introduced recently by Ruffini et al. (RSV). The probability density function of the last scattering of photons is shown to be very different in these two cases. We also obtained spectra as well as light curves. In the photon-thick case, the time-integrated spectrum is much broader than the Planck function and its shape is well described by the fuzzy photosphere approximation introduced by RSV. In the photon-thin case, we confirm the crucial role of photon diffusion, hence the probability density of decoupling has a maximum near the diffusion radius well below the photosphere. The time-integrated spectrum of the photon-thin case has a Band shape that is produced when the outflow is optically thick and its peak is formed at the diffusion radius.

79 ASTRONOMY AND ASTROPHYSICS↗