Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transport models”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Diagnosing the DARHT Electron Beam-Target Interaction and Hydrodynamic Expansion [Dissertation]

Quantitative electron beam-target interaction studies on the Dual-Axis Radiographic Hydrodynamic Test (DARHT) electron linear induction accelerators have only been performed recently. This includes characterization of the temperature, density, pressure, extent, and expansion velocity of the plasma plume. The results presented in this dissertation present a detailed and unified overview of the accelerator systems, target heating physics, beam transport, and diagnostic tools. Additional information includes calibration sources, radiation hydrodynamics, spectroscopic-quality radiation transport modeling, experimental measurements, and analyses of electron beam driven aluminum experiments. The first set of spatially and temporally resolved spectroscopic measurements of electron beam driven aluminum are presented. Contamination quantification analyses are used to understand the origin of the strong Na-I 3p-3s lines that are observed in absorption within the aluminum plasma continuum. These results inform the creation of the first spectroscopic-quality radiation transport model that links several atomic physics codes to interpret the conditions from which the Na-I lines originate. A good agreement is found between the surface analysis results and the model which confirms the concentration of the sodium present within the aluminum alloy foil material. It also demonstrates, for the first time on electron beam-driven target experiments, the ability to interpret plasma conditions from measured absorption lines. In a second experimental campaign which focuses on pure aluminum, the Al-I 3p-4s and 3p3d doublets are both measured in emission. A detailed analysis of the Al-I 3p-4s doublet reveals that the lines undergo moderate self-absorption. A simple model of the self-absorption effect is successfully used to match measured spectra at various temperature/density/plasma scale length combinations. These measurements led to the realization of the minimum density that can be resolved by the spectrometers for the Al-I 3p-4s lines due to the large slit width required to observe a signal on aluminum. These measurements also demonstrate the sensitivity of visible and long wave UV spectroscopy to minor changes in both temperature and density. The simple self-absorption model will be useful for analysis of other beam-target interaction experiments with spectra exhibiting either self-absorption or full self-reversal. Substantial headway has been made on the modeling front by linking together several codes needed to model both the energy deposition, hydrodynamic motion, and atomic kinetics to produce synthetic spectral calculations that are compared with experimental measurements. There exists ample space for improvement, especially with benchmarking the hydrodynamics codes and equation-of-state tables with experimental measurements. The X-ray diagnostics required to make these new measurements along with the simulation capabilities required to interpret the results are under development and will be the subject of future studies.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Electron Scattering from 1-Methyl-5-Nitroimidazole: Cross-Sections for Modeling Electron Transport through Potential Radiosensitizers

In this study, we present a complete set of electron scattering cross-sections from 1-Methyl-5-Nitroimidazole (1M5NI) molecules for impact energies ranging from 0.1 to 1000 eV. This information is relevant to evaluate the potential role of 1M5NI as a molecular radiosensitizers. The total electron scattering cross-sections (TCS) that we previously measured with a magnetically confined electron transmission apparatus were considered as the reference values for the present analysis. Elastic scattering cross-sections were calculated by means of two different schemes: The Schwinger multichannel (SMC) method for the lower energies (below 15 eV) and the independent atom model-based screening-corrected additivity rule with interferences (IAM-SCARI) for higher energies (above 15 eV). The latter was also applied to calculate the total ionization cross-sections, which were complemented with experimental values of the induced cationic fragmentation by electron impact. Double differential ionization cross-sections were measured with a reaction microscope multi-particle coincidence spectrometer. Using a momentum imaging spectrometer, direct measurements of the anion fragment yields and kinetic energies by the dissociative electron attachment are also presented. Cross-sections for the other inelastic channels were derived with a self-consistent procedure by sampling their values at a given energy to ensure that the sum of the cross-sections of all the scattering processes available at that energy coincides with the corresponding TCS. This cross-section data set is ready to be used for modelling electron-induced radiation damage at the molecular level to biologically relevant media containing 1M5NI as a potential radiosensitizer. Nonetheless, a proper evaluation of its radiosensitizing effects would require further radiobiological experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Phase Contrast Imaging Diagnostic on DIII-D (1994-2007); Phase Contrast Imaging Diagnostic on DIII-D and Alcator C-Mod (2007-2010); Experimental Study of Waves, Turbulence and Transport, and Model Validation with Phase Contrast Imaging on DIII-D and Alcator C-Mod (2010-2012); Phase Contrast Imaging Diagnostic Upgrade on DIII-D for Measuring High-k Turbulence and Alfven modes (2021-2016); Phase Contrast Imaging Diagnostic Upgrade on DIII-D for Full Multi-Scale Measurements (2016-2020) (Final Report)

This report describes the first twenty-five years of development and use of the Phase Contrast Imaging (PCI) diagnostic for measuring plasma fluctuations on the DIII–D tokamak. One of the key challenges in moving tokamak technology forward to the stage of a nuclear reactor remains understanding how energy and particles are transported across the magnetic field by the action of turbulence and waves, and PCI provides a highly flexible and capable measurement of these fluctuations. The PCI diagnostic uses a laser beam to probe the plasma and optical filtering to produce an amplitude variation on a linear detector array which provides an image of the density waves in the plasma. The PCI covers an exceptional range in frequency and wavelength of the measured fluctuations, allowing observation of plasma-sized magnetohydrodynamic (MHD) modes and ion- and electron-scale turbulence. Modulation of the beam amplitude allows measurement of radio frequency (RF) waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effective kinetics driven by dynamic concentration gradients under coupled transport and reaction

Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Quantification of Spore-forming Bacteria Carried by Dust Particles

In order to establish a biological contamination transport model for predicting the cross contamination risk during spacecraft assembly and upon landing on Mars, it is important to understand the relationship between spore-forming bacteria and their carrier particles. We conducted air and surface sampling in indoor, outdoor, and cleanroom environments to determine the ratio of spore forming bacteria to their dust particle carriers of different sizes. The number of spore forming bacteria was determined from various size groups of particles in a given environment. Our data also confirms the existence of multiple spores on a single particle and spore clumps. This study will help in developing a better bio-contamination transport model, which in turn will help in determining forward contamination risks for future missions.

particles↗
NASA NTRS2008

Impact of Surface Emissions to the Zonal Variability of Tropical Tropospheric Ozone and Carbon Monoxide for November 2004

The chemical and dynamical processes governing the zonal variability of tropical tropospheric ozone and carbon monoxide are investigated for November 2004 using satellite observations, in-situ measurements, and chemical transport models in conjunction with inverse-estimated surface emissions. Vertical ozone profile estimates from the Tropospheric Emission Spectrometer (TES) and ozone sonde measurements from the Southern Hemisphere Additional Ozonesondes (SHADOZ) network show the so called zonal 'wave-one' pattern, which is characterized by peak ozone concentrations (70-80 ppb) centered over the Atlantic, as well as elevated concentrations of ozone over Indonesia and Australia (60-70 ppb) in the lower troposphere. Observational evidence from TES CO vertical profiles and Ozone Monitoring Instrument (OMI) NO2 columns point to regional surface emissions as an important contributor to the elevated ozone over Indonesia. This contribution is investigated with the GEOS-Chem chemistry and transport model using surface emission estimates derived from an optimal inverse model, which was constrained by TES and Measurements Of Pollution In The Troposphere (MOPITT) CO profiles (Jones et al., 2007). These a posteriori estimates, which were over a factor of 2 greater than climatological emissions, reduced differences between GEOS-Chem and TES ozone observations by 30-40% and led to changes in GEOS-Chem upper tropospheric ozone of up to 40% over Indonesia. The remaining residual differences can be explained in part by upper tropospheric ozone produced from lightning NOx in the South Atlantic. Furthermore, model simulations from GEOS-Chem indicate that ozone over Indonesian/Australian is more sensitive to changes in surface emissions of NOx than ozone over the tropical Atlantic.

tropical tropospheric ozone↗
DOE OSTI2020

On the Impact of Income, Age, and Travel Distance on the Value of Time

The value of time (VOT) is a fundamental component used in transportation modeling, policy analysis, and economic appraisal. Decades of research and practice have empirically estimated the VOT across many factors (e.g., mode, purpose, time, comfort, etc.), yet little is known about its underlying form. Although it is well established that VOT can vary, it is still unclear whether patterns exist in this variation. The objective of this paper is not to merely estimate the VOT, but to model the VOT across multiple continuous and interacting variables. The purpose is to reveal its functional form with respect to mode, age, gender, purpose, income, and time of day to provide a generalizable understanding for future research and practice. Such an understanding can help develop simpler models and reduce the need for bespoke estimations for every conceivable variable perturbation. This research utilized a household travel survey containing 14,159 reported trips with imputed travel time and costs for the alternative mode choices. The average overall estimated VOT is 40.32 $/h, with results showing VOT varying log-linearly with income and trip distance, but following a Gaussian function (normal curve) with age. Overall, the results show that travel distance dominates VOT variation, which increases exponentially at a rate that is 3.61 times higher per mile of distance than per $10,000 of income, and that VOT by age peaks at age 54. This basic understanding of how the VOT varies sets the foundation for answering the subsequent question for why it might vary.

Engineering↗
DOE OSTI2021

In-medium $\Delta (1232)$ potential, pion production in heavy-ion collisions and the symmetry energy

Using the dcQMD transport model, the isoscalar and isovector in-medium potentials of the Δ (1232) baryon are studied and information regarding their effective strength is obtained from a comparison to experimental pion production data in heavy-ion collisions below 800 MeV/nucleon impact energy. The best description is achieved for an isoscalar potential moderately more attractive than the nucleon optical potential and a rather small isoscalar relative effective mass $m$$^{*}_{Δ}$ ≈ 0.45. For the isovector component only a constraint between the potential’s strength at saturation and the isovector effective mass difference can be extracted, which depends on quantities such as the slope of the symmetry energy and the neutron-proton effective mass difference. These results are incompatible with the usual assumption, in transport models, that the Δ (1232) and nucleon potentials are equal. The density dependence of symmetry energy can be studied using the high transverse momentum tail of pion multiplicity ratio spectra. Results are however correlated with the value of neutron-proton effective mass difference. Furthermore, this region of spectra is shown to be affected by uncertain model ingredients such as the pion potential or in-medium correction to inelastic scattering cross-sections at levels smaller than 10%. Extraction of precise constraints for the density dependence of symmetry energy above saturation will require experimental data for pion production in heavy-ion collisions below 800 MeV/nucleon impact energy and experimental values for the high transverse momentum tail of pion multiplicity ratio spectra accurate to better than 5%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Real-Time Atomic-Scale Structural Analysis Resolves the Amorphous to Crystalline CaCO 3 Mechanism Controversy

Amorphous calcium carbonate (ACC) occurs as a precursor to geological and biogenic calcium carbonate (CaCO 3 ), yet its transformation pathways and reaction mechanisms remain inconsistent and controversial. In this study, we investigated the transformation of ACC to calcite under both solution and dry conditions, in the presence and absence of impurity ions, utilizing operando time-resolved synchrotron X-ray diffraction (TRXRD) and reactive transport modeling. Results demonstrate that TRXRD techniques allow us to differentiate dissolution-reprecipitation versus solid-state transformation mechanisms for amorphous to crystalline phase transitions. Specifically, we observe that in environments with abundant water, ACC transforms to calcite through a dissolution-reprecipitation mechanism. This features an activation energy of 63 ± 2 kJ/mol and unit cell volume contraction during calcite crystal growth. Conversely, under water-limited conditions, ACC to calcite transformation proceeds via a solid-state transformation mechanism, with an activation energy of 210 ± 2 kJ/mol, three times greater than the dissolution-reprecipitation route, and a unit cell expansion during crystalline calcite growth. Further, to illustrate the magnitude of these effects, the rates of calcite growth were similar during dissolution-reprecipitation at 3 °C [0.00207(35) s –1 ] and solid-state transformation at 280 °C [0.00134(11) s –1 ]. Moreover, the incorporation of an impurity, strontium, significantly retards the rate of calcite growth while expanding its unit cell but whose incorporation is history dependent. Reactive transport modeling of the dissolution–precipitation kinetics suggests that ACC must be dissolving as compact aggregates. These various transformation mechanisms drive diverse geological and biological carbonate formations, impacting their use as paleoenvironmental markers and functional materials synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

How Well Will MODIS Measure Top of Atmosphere Aerosol Direct Radiative Forcing?

The new generation of satellite sensors such as the Moderate Resolution Imaging Spectroradiometer (MODIS) will be able to detect and characterize global aerosols with an unprecedented accuracy. The question remains whether this accuracy will be sufficient to narrow the uncertainties in our estimates of aerosol radiative forcing at the top of the atmosphere. Satellite remote sensing detects aerosol optical thickness with the least amount of relative error when aerosol loading is high. Satellites are less effective when aerosol loading is low. We use the monthly mean results of two global aerosol transport models to simulate the spatial distribution of smoke aerosol in the Southern Hemisphere during the tropical biomass burning season. This spatial distribution allows us to determine that 87-94% of the smoke aerosol forcing at the top of the atmosphere occurs in grid squares with sufficient signal to noise ratio to be detectable from space. The uncertainty of quantifying the smoke aerosol forcing in the Southern Hemisphere depends on the uncertainty introduced by errors in estimating the background aerosol, errors resulting from uncertainties in surface properties and errors resulting from uncertainties in assumptions of aerosol properties. These three errors combine to give overall uncertainties of 1.5 to 2.2 Wm-2 (21-56%) in determining the Southern Hemisphere smoke aerosol forcing at the top of the atmosphere. The range of values depend on which estimate of MODIS retrieval uncertainty is used, either the theoretical calculation (upper bound) or the empirical estimate (lower bound). Strategies that use the satellite data to derive flux directly or use the data in conjunction with ground-based remote sensing and aerosol transport models can reduce these uncertainties.

Remer, Lorraine A.↗
NASA NTRS2009

Use of NASA Satellite Data in Aiding Mississippi Barrier Island Restoration Projects

This presentation discusses a NASA Stennis Space Center project in which NASA-supported satellite and aerial data is being used to aid state and federal agencies in restoring the Mississippi barrier islands. Led by the Applied Science and Technology Project Office (ASTPO), this project will produce geospatial information products from multiple NASA-supported data sources, including Landsat, ASTER, and MODIS satellite data as well as ATLAS multispectral, CAMS multispectral, AVIRIS hyperspectral, EAARL, and other aerial data. Project objectives include the development and testing of a regional sediment transport model and the monitoring of barrier island restoration efforts through remote sensing. Barrier islands provide invaluable benefits to the State of Mississippi, including buffering the mainland from storm surge impacts, providing habitats for valuable wildlife and fisheries habitat, offering accessible recreational opportunities, and preserving natural environments for educating the public about coastal ecosystems and cultural resources. Unfortunately, these highly valued natural areas are prone to damage from hurricanes. For example, Hurricane Camille in 1969 split Ship Island into East and West Ship Island. Hurricane Georges in 1998 caused additional land loss for the two Ship Islands. More recently, Hurricanes Ivan, Katrina, Rita, Gustav, and Ike impacted the Mississippi barrier islands. In particular, Hurricane Katrina caused major damage to island vegetation and landforms, killing island forest overstories, overwashing entire islands, and causing widespread erosion. In response, multiple state and federal agencies are working to restore damaged components of these barrier islands. Much of this work is being implemented through federally funded Coastal Impact Assessment and Mississippi Coastal Improvement programs. One restoration component involves the reestablishment of the island footprints to that in 1969. Our project will employ NASA remote sensing data and products to support these federally funded efforts on multiple fronts. Landsat and ASTER data is being analyzed to assess changes in barrier island land cover over the last 35 years. ASTER, SRTM, and EAARL terrain products and other NASA airborne imagery are being applied in assessing changes in barrier island geomorphology and geospatial extent. MODIS data is being examined as a tool for sediment transport modeling by supplying geospatial data that quantifies in-water sediment concentrations. MODIS satellite data is being assessed for monitoring changes in the spatial extent of individual barrier islands. Results thus far indicate that NASA data products are useful in assessing barrier island conditions and changes. This value is enhanced with additional historical geospatial data, commercial high resolution satellite data, other non-NASA aerial imagery, and field survey data. The project s products are relevant to the Gulf of Mexico Alliance priority issues, including coastal habitat conservation, restoration and coastal community resilience. Such products will be available to state and federal agencies involved with coastal restoration. Potential end-users of these products include the National Park Service, U.S. Geological Survey, U.S. Army Corps of Engineers, Environmental Protection Agency, Mississippi Department of Environmental Quality, and Mississippi Department of Marine Resources.

Giardino, Marco↗
NASA NTRS2017

Multi-Decadal Records of Stratospheric Composition and Their Relationship to Stratospheric Circulation Change

Constituent evolution for 1990-2015 simulated using the Global Modeling Initiative chemistry and transport model driven by meteorological fields from the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) is compared with three sources of observations: ground-based column measurements of HNO3 and HCl from two stations in the Network for the Detection of Atmospheric Composition Change (NDACC, 1990- ongoing), profiles of CH4 from the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satellite (UARS, 1992-2005), and profiles of N2O from the Microwave Limb Sounder on the Earth Observing System satellite Aura (2005- ongoing). The differences between observed and simulated values are shown to be time dependent, with better agreement after 2000 compared with the prior decade. Furthermore, the differences between observed and simulated HNO3 and HCl columns are shown to be correlated with each other, suggesting that issues with the simulated transport and mixing cause the differences during the 1990s and that these issues are less important during the later years. Because the simulated fields are related to mean age in the lower stratosphere, we use these comparisons to evaluate the time dependence of mean age. The ongoing NDACC column observations provide critical information necessary to substantiate trends in mean age obtained using fields from MERRA-2 or any other reanalysis products.

Profiles of N2O from the Microwave Limb Sounde↗
NASA NTRS2019

Why Do Antarctic Ozone Recovery Trends Vary?

We use satellite ozone records and Global Modeling Initiative (GMI) chemistry transport model (CTM) simulations integrated with MERRA2 meteorology to identify a metric that accurately captures the trend in Antarctic ozone attributable to the decline in ozone depleting substances (ODSs). The GMI CTM Baseline simulation with realistically varying ODS levels closely matches observed interannual to decadal scale variations in Antarctic September ozone over the past 4 decades. The expected increase or recovery trend is obtained from the differences between the Baseline simulation and one with identical meteorology and fixed 1995 ODS levels. The differences show that vortex-averaged column O3 has the greatest sensitivity to ODS change from 1-20 Sep. The observed vortex-averaged column O3 during this period produces a trend consistent with the expected recovery attributable to ODS decline. Trends from dates after 20 Sep have smaller sensitivity to ODS decline and are more uncertain due to transport variability. Simulations show that the greatest decrease in O3 loss (i.e., recovery) occurs inside the vortex near the edge. The polar cap metrics have vortex size-dependent bias and do not consistently sample this region. Because the 60-90oS 220 DU O3 mass deficit (OMD) metric does not sample the edge region, its trend is lower than the expected trend; this is improved by area weighting. The 250 DU OMD metric samples more of the edge region, which increases its trend. Approximately 25% of the September Antarctic O3 increase is due to higher O3 levels in June prior to winter depletion.

Antarctic ozone↗
NASA NTRS2004

Sensitivity of Global Modeling Initiative Model Predictions of Antarctic Ozone Recovery to Input Meteorological Fields

We use the Global Modeling Initiative chemistry and transport model to simulate the evolution of stratospheric ozone between 1995 and 2030, using boundary conditions consistent with the recent World Meteorological Organization ozone assessment. We compare the Antarctic ozone recovery predictions of two simulations, one driven by an annually repeated year of meteorological data from a general circulation model (GCM), the other using a year of output from a data assimilation system (DAS), to examine the sensitivity of Antarctic ozone recovery predictions to the characteristic dynamical differences between GCM- and DAS-generated meteorological data. Although the age of air in the Antarctic lower stratosphere differs by a factor of 2 between the simulations, we find little sensitivity of the 1995-2030 Antarctic ozone recovery between 350 and 650 K to the differing meteorological fields, particularly when the recovery is specified in mixing ratio units. Percent changes are smaller in the DAS-driven simulation compared to the GCM-driven simulation because of a surplus of Antarctic ozone in the DAS-driven simulation which is not consistent with observations. The peak ozone change between 1995 and 2030 in both simulations is approx.20% lower than photochemical expectations, indicating that changes in ozone transport due to changing ozone gradients at 450 K between 1995 and 2030 constitute a small negative feedback. Total winter/spring ozone loss during the base year (1995) of both simulations and the rate of ozone loss during August and September is somewhat weaker than observed. This appears to be due to underestimates of Antarctic Cl(sub y) at the 450 K potential temperature level.

Considine, David B.↗
DOE OSTI2021

Modeling of transport in carrier-selective contacts in silicon heterojunction solar cells

High-performance silicon heterojunction (SHJ) solar cells use carrier-selective contact structures based on hydrogentated amorphous Si (a-Si:H) to maximize collection of photogenerated carriers. The high open circuit voltages observed experimentally in SHJ cells require that the carrier-selective contacts provide selectivity and passivation. However, a microscopic understanding of the dynamics of carrier transport through the a-Si layer is currently lacking. In this paper, we explicitly simulate the transport of holes across the a-Si:H(i) layer using a novel kinetic Monte Carlo approach. The hole-selective contact structure investigated in this paper uses p-type doped a-Si:H(p) and intrinsic a-Si:H(i) on an n-type crystalline silicon wafer, where the selectivity is provided by the a-Si:H(p) and the passivation is provided by the a-Si:H(i). However, in addition to the passivation provided by the a-Si:H(i), this layer also creates a potential barrier to the collection of photogenerated holes. There have been experimental studies in the literature that have suggested that multi-phonon processes are the main transport mechanism that assists in the transport of holes across the intrinsic a-Si:H barrier. Simulations presented here show that multi-phonon injection of holes into the a-Si:H(i) layer is the rate limiting step for transport across the a-Si:H(i) layer. Our results indicate that multi-phonon transport is strongly dependent on the electric field at the a-Si:H(i)/c-Si heterointerface as well. Transport simulations presented in this paper are consistent with experimental findings that multi-phonon processes limit transport across the a-Si:H(i) layer and are responsible for photocurrent suppression at the a-Si:H(i)/c-Si heterointerface when these processes are slower than the associated incident hole flux due to photo-excitation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae

The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.

79 ASTRONOMY AND ASTROPHYSICS↗