Search NASA⌕ Search

SEARCH · Search NASA

Results for “Differential transport”

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

Production of Fischer-Tropsch Synfuels at Nuclear Plants

A case study analysis was performed to evaluate nuclear-powered synthetic fuel production in the midwestern United States (U.S.). A Fischer-Tropsch (FT) fuel synthesis plant design was used as the basis for the analysis. The FT plant design was configured to produce a product slate consisting of diesel fuel, jet fuel, and motor gasoline blend stocks from carbon dioxide (CO 2 ) and hydrogen (H 2 ) feedstocks. The CO 2 feedstock for the FT plant was assumed to be sourced from biorefineries in the region around a Midwest light water reactor (LWR) nuclear power plant (NPP). The analysis specifies that power from the LWR is used to produce H 2 via high-temperature steam electrolysis and to operate the FT synfuel production plant. Capital costs were estimated for the FT plant while capital costs for the electrolysis plant were based on previous Idaho National Laboratory (INL) studies. In addition to labor and maintenance costs for the FT and electrolysis plants, operating costs also include the costs for CO 2 feedstock transport. An analysis was performed to determine the cost of transporting CO 2 from the distributed biorefinery sources to the centralized fuel synthesis plant as a function of the synfuel plant capacity and corresponding CO 2 demand. The primary revenue streams are associated with sales of the synthetic fuel products. The synthetic fuel products will likely follow the same market trends as the conventional fuel products. The synfuel price data was thus based on projections made by the U.S. Energy Information Administration (EIA) 2021 Annual Energy Outlook (AEO) for conventional fuel products minus federal and state taxes, as well as marketing and distribution costs. The economic analysis also considered cases that included and excluded revenues from the 2022 Inflation Reduction Act (IRA) clean hydrogen production tax credit (PTC) of $\$ $3.00/kg for the first ten years of operation. The economic analysis calculated the net present value (NPV) for cases involving steady-state synfuel production for comparison with the NPV for a business-as-usual case in which NPP continues to sell only electric power to the grid. A synfuel production “Reference Case” was considered in addition to sensitivity cases in which the plant capacity, electricity price, and synthetic fuel product prices were perturbed. The synfuel production Reference Case considered a scenario in which the electrolysis and synfuel plants utilized a combined electrical load of 1000 megawatt electrical (MWe) from the LWR with the balance of the LWR power output being sold to the electric grid. The economic analysis suggests that the synfuel production Reference Case evaluated in this analysis would lead to considerable economic potential for near-term deployment of a nuclear-based synfuel production plant. Specifically, the economic analysis suggests that the deployment of a 1000 megawatt (MW) nuclear-powered synfuel plant could result in a NPV increase of approximately $\$ $1.7 billion for a case with no clean synfuel price premium relative to conventional petroleum fuels when accounting for the additional revenues from the 2022 IRA clean hydrogen PTCs of $\$ $3/kg. Sensitivity analysis was performed to evaluate the effect of perturbation of selected model input parameters on the NPV for the synfuel production Reference Case. The sensitivity analysis indicates that the plant capacity has the largest impact on the differential NPV, with a smaller synfuel production capacity resulting in a decrease in revenue when a larger fraction of the power from the NPP is sold to the grid and a smaller fraction of the power is used to produce synthetic fuel products. The synfuel product pricing has the next largest impact on the differential NPV, with lower synfuel prices resulting in decreased NPV from decreased synfuel sales revenue while higher synfuel prices result in increased NPV from increased synfuel sales revenue. Electricity pricing has a smaller effect on the NPV than the fuel sales price since, in the Reference Case, most of the energy from the NPP is used for synfuel production and a smaller amount of the system revenues are associated with electrical power sales. However, the electricity price sensitivity does indicate that the Synfuel Integrated Energy System (IES) would have a greater NPV than the business-as-usual case (e.g., grid power sales only) when electricity market prices are low, suggesting that synfuel production could provide a strategy for decreasing the economic risks to NPPs posed by a loss of revenues attributed to falling electricity market prices.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Validation of Jezebel Reactivity Coefficients and Sensitivity Analysis

Nuclear data validation is often performed today using criticality measurements. The gold standard for criticality measurements is the International Criticality Safety Benchmark Experiment Project (ICSBEP). The validation specifically focuses on the effective multiplication factor ($k_{eff}$). $K_{eff}$ is a relatively easy parameter to infer and has reduced uncertainty due to being at or above critical. However, while $k_{eff}$ is the most documented parameter and its uncertainties and sensitivities have been evaluated in great detail, it cannot be used as a standalone metric to determine inaccuracies in nuclear data (e.g., cross section data, PFNS, nu), which is based on theory, physics, and differential measurements. The Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) project aims to identify compensating errors in specific isotope nuclear data by optimally designing experiments that are, or are not sensitive to a suite of measurement parameters beyond $k_{eff}$. By identifying parameters that are sensitive to each other, oppositely sensitive, or have substantial magnitude differences in sensitivity, experiments can be designed to constrain questionable nuclear data. One sensitivity that is of particular interest to this project includes the sensitivity of reactivity coefficients. Reactivity coefficients compare reactivity, which is related to $k_{eff}$ at two different states therefore being sensitive to small changes in the system. The most common type of reactivity coefficient measurements is comparison to void for a small sample within the assembly. It is key that the sample sizes are small enough to not affect the flux of the full system. Reactivity coefficients were evaluated for many early experiments to better understand transport corrected cross sections. In fact, ICSBEP includes reactivity coefficient results as “Supplemental Measurements” in appendices for a handful of older benchmarks. One of those benchmarks is Jezebel, the bare Pu critical assembly. This paper compares new simulations of reactivity coefficients for Jezebel, and explores the sensitivity of reactivity coefficients to small changes in nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Using Co-Optimized Machine Learned Manifolds for Modeling Chemically Reacting Flows

Chemically reacting flows play a key role in a wide range of engineered systems, from chemical and polymer processing to combustion-based energy conversion technologies. Simulations of these flows involve solving a coupled set of partial differential equations for mass, momentum, energy, and all relevant chemical species in the system. Chemical reaction pathways may be extremely complex and involve hundreds or more intermediate species, with reactions that occur over timescales varying by several orders of magnitude - presenting a significant numerical stiffness challenge. The combination of these factors makes simulation of chemically reacting flows vastly more expensive than nonreactive simulations, and often makes direct solution of the governing equations intractable. It is necessary to apply lower-fidelity models in place of the detailed governing equations in order to reduce computational cost to enable reacting flow simulation tools to be used in the engineering design process. Many of the models employed for this purpose are based on reducing the dimension of the thermochemical state, motivated by the observation that the observed thermochemical states in a system lie on a low-dimensional manifold in thermochemical state space. This behavior occurs due to the fast equilibration of certain reactive and transport processes, and physics-based manifold models rely on idealized assumptions about the balance of timescales and the way in which chemistry and transport are coupled. In this work, we apply a novel method for data-driven manifold-based modeling that can leverage data from high-fidelity reacting flow simulations to improve model accuracy in cases where the physics-based modeling assumptions break down. The approach is designed to be broadly applicable across chemically reacting flow systems but is applied here to turbulent combustion modeling.

machine learning↗
DOE OSTI2023

Ice nucleation measurements from DRUM impactors during TRACER campaign in the Houston TX region from July to September 2022

During TRACER, three Davis Rotating Uniform size-cut Monitors (DRUM; DRUMAir 4-DRUM) were used to collect aerosols for ice nucleation measurements in the Brooks laboratory at Texas A&M University. The three instruments were located at AMF1 in La Porte, Ancillary site in Guy (ANC), and onboard the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V). ROAM-V was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and one of the Davis Rotating Uniform size-cut Monitors (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Each DRUM sampler was operated at a flow rate of 23 LPM. Each DRUM has four stages where the aerodynamic diameter size cuts are as follows: stage 1 larger than 3 μm, stage 2 from 3 to 1.2 μm, stage 3 from 1.2 to 0.34 μm, and stage 4 from 0.34 to 0.15 μm. Pretreated aluminum foil was used as a substrate on all stages. At AMF1 and ANC sites, the DRUMs were operated on the shared aerosol inlet, with generous support of the DOE ARM site staff. These instruments rotated 24 mm per day and only contained aluminum foil substrates. Substrates were changed weekly and transported to Texas A&M for storage in -80C freezer until analysis. The DRUM onboard ROAM-V was operated at a faster rotation rate of 150 mm per day to clearly separate the multiple deployment locations for ROAM-V. For the ROAM-V DRUM, substrates were changed every deployment and transported to Texas A&M for storage in -80C freezer until analysis. At Texas A&M, ice nucleation experiments were conducted to measure the ice nucleation temperature of the ambient aerosol samples collected from the three DRUMs using our previously established procedures (Alsante et al., 2023; Fornea et al., 2009; Matthews et al., 2023). For ice nucleation, only samples collected on stage 3 were analyzed, given that these are the most relevant size (1.2 to 0.34 μm diameter) for potential ice nucleating particles. For the AMF1 and ANC sites, we cut and analyzed 2 mm (2-hour) samples. We analyzed the time periods of the AMF1 site instrument when the ROAM-V was deployed. A 72-hour period from July 11th at 23:49 through July 15th at 1:49 was analyzed from the ANC site instrument. For the ROAM-V instrument, we separated the daily samples by site location. Between 1- and 6-hour independent samples were analyzed at each location. We also cut the ROAM-V samples in half to allow for compositional analysis of the aerosol on the other half of the substrate. All viable samples from the ROAM-V were analyzed. Analysis was done using a custom-built ice nucleation apparatus, recently updated to include an array of 16 individual samples (Matthews et al., 2023). On our experimental setup, we used 100 μL of Ultra-High-Performance Liquid Chromatography (UHPLC) water (Sigma Aldrich, >99.9% purity) to wash off the aerosol from the DRUM substrate. Then, we micropipetted 2 μL droplet samples into each of the 16 wells of the array, using hydrophobically coated microscope slides. Experiments would cycle 28 times from 10 C to -40C over a 20-hour period. During analysis of the TRACER samples, nine experiments were conducted with UHPLC water process blanks that followed an identical preparation procedure. The average freezing temperature and standard deviation for these process blanks is -27.7±2.4 C. This data was collected for ARM Field Campaign AFC07023 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Alsante, A. N., Thornton, D. C., & Brooks, S. D. (2023). Ice nucleation catalyzed by the photosynthesis enzyme RuBisCO and other abundant biomolecules. Communications Earth & Environment, 4(1), 51. Fornea, A. P., Brooks, S. D., Dooley, J. B., & Saha, A. (2009). Heterogeneous freezing of ice on atmospheric aerosols containing ash, soot, and soil. Journal of Geophysical Research: Atmospheres, 114(D13). DOI:10.1029/2009JD011958 Matthews, B. H., Alsante, A. N., & Brooks, S. D. (2023). Pollen Emissions of Subpollen Particles and Ice Nucleating Particles. ACS Earth and Space Chemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impacts of processing decisions on TNSL cross sections and their applications

Thermal neutron scattering law (TNSL) data describe low-energy neutrons scattering off of bound materials, and can have a significant impact on modeling any system with slow neutrons, including nuclear reactors. Previous work to introduce TNSL data to neutron transport codes at LLNL focused on COG and TART, with the limitation that these codes require highly specialized data processing and formatting. We have recently increased efforts to process TNSL data with the LLNL nuclear data processing code FUDGE. FUDGE reads and writes the evaluated and processed files using the generalized nuclear database structure (GNDS). This process uncovered some significant difficulties in processing TNSL data, and unearthed assumptions made in current TNSL data processing that we have found inadequate. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

General balance functions of identified charged hadron pairs of (π, Κ, ρ) in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

First measurements of balance functions (BFs) of all combinations of identified charged hadron (π, Κ, ρ) pairs in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV recorded by the ALICE detector are presented. The BF measurements are carried out as two-dimensional differential correlators versus the relative rapidity (Δy) and azimuthal angle (Δφ) of hadron pairs, and studied as a function of collision centrality. The Δφ dependence of BFs is expected to be sensitive to the light quark diffusivity in the quark–gluon plasma. While the BF azimuthal widths of all pairs substantially decrease from peripheral to central collisions, the longitudinal widths exhibit mixed behaviors: BFs of ππ and cross-species pairs narrow significantly in more central collisions, whereas those of KK and pp are found to be independent of collision centrality. This dichotomy is qualitatively consistent with the presence of strong radial flow effects and the existence of two stages of quark production in relativistic heavy-ion collisions. Finally, the first measurements of the collision centrality evolution of BF integrals are presented, with the observation that charge balancing fractions are nearly independent of collision centrality in Pb–Pb collisions. Overall, the results presented provide new and challenging constraints for theoretical models of hadron production and transport in relativistic heavy-ion collisions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Hygroscopic Growth of UO 2 F 2 Particles

Hygroscopicity is an important physicochemical property of aerosol that describes the ability of a particle to uptake water. The hygroscopic properties of uranyl fluoride (UO 2 F 2 ) aerosol generated from a UF 6 hydrolysis reactor was investigated for the first time using a custom-built Humidified Tandem Differential Mobility Analyzer (HTDMA). The HTDMA is capable of measuring UO 2 F 2 particle growth determined by mobility size over a wide range of atmospheric humidity from dry conditions at <10% relative humidity (RH) to 85% RH. The hygroscopic properties were determined for nanoparticles as small as 3.5 nm in this study. Anhydrous UO 2 F 2 particles with a mobility diameter of 3.5 nm were shown to be highly hygroscopic with a deliquescence relative humidity (DRH) of 10%. Hydrates with a larger mobility diameter of 10 to 80 nm were non-hygroscopic with no observable DRH and limited water uptake up to 85% RH. These results demonstrate the hygroscopic properties of UO 2 F 2 particles are highly variable and based on both the mobility size and hydration state. Hygroscopicity affects the physicochemical properties of UO 2 F 2 particles, including the aerosol phase state and viscosity, with impacts on aerosol growth, coagulation, and deposition that is critical for understanding the fate and transport of UO 2 F 2 particles in the atmosphere.

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

Permeability scaling relationships of volcanic tuff from core to field scale measurements

A recent chemical explosive test in P-Tunnel at the Nevada National Security Site, Nevada, USA, was conducted to better understand how signals propagate from explosions in the subsurface. A primary signal of interest is the migration of gases that can be used to differentiate chemical from nuclear explosions. Gas migration is highly dependent on the rock permeability which is notoriously difficult to determine experimentally in the field due to a potentially large dependence on the scale over which measurements are made. Here, we present pre-explosion permeability estimates to characterize the geologic units surrounding the recent test. Permeability measurements were made at three scales of increasing size: core samples (≈2 cm), borehole packer system tests (≈1 m), and a pre-shot cavity pressurization test (> 10 m) across ten tuff units. Permeability estimates based on core measurements showed little difference from borehole packer tests. However, permeability in most rock units calibrated from cavity pressurization tests resulted in higher permeability estimates by up to two orders of magnitude. Here, we demonstrate that the scale of the measurement significantly impacts the characterization efforts of hydraulic properties in volcanic tuff, and that local-scale measurements (< 10 m scale) do not incorporate enough heterogeneity to accurately predict field-scale flow and mass transport.

Environmental sciences↗
DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Analysis of an Arctic cold air outbreak during autumn and related air mass transformations forced by surface changes and advection in higher altitudes

During marine cold air outbreaks (MCAOs), cold and dry Arctic air masses are transported from the central Arctic southward across the closed sea ice and much warmer open oceans. They experience significant transformations including a rapid heating and moistening, often leading to cloud formation. While intense wintertime MCAOs have been analyzed widely, the air mass transformations during other seasons have been studied sparsely. We address this gap by investigating an MCAO case observed in September 2020. To study the transformation processes, we combine the fifth generation of atmospheric reanalyses of the global climate (ERA5), trajectory calculations, as well as shipborne and airborne measurements. In the central Arctic, observations acquired from aboard the research vessel (RV) Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition characterized the initial state of the air mass over closed sea ice. Trajectories indicated the pathway the air mass took from RV Polarstern southward to the Fram Strait. For the first 24 h of the southbound drift, the air masses remained quasi-stationary. Then, still 15 h ahead of the marginal sea ice zone, differential advection across the boundary layer flow introduced humidity and clouds at higher altitudes between 1.5 and 2.5 km. ERA5-derived temperature and humidity tendencies indicated complex vertical interactions. Radiative cloud-top cooling, entrainment, and turbulence were significantly reduced in the lower and enhanced in the upper advected cloud layer. Eventually, the lower cloud deck dissipated. After this confluence of 2 different air masses, observations gathered by Polar 5 in Fram Strait as part of the MOSAiC Airborne observations in the Central Arctic campaign revealed cloudy, moist layers throughout the lowest 3.5 km and an increasing boundary layer height. Comparing the initial with the final state 48 h later, the largest net heating of +8 K was found close to the surface, yet the largest net moistening of +2.5 g kg -1 at an altitude of 1 km, as the initial profile was exceptionally dry here. We conclude that the observed air mass transformations were driven by the surface changes from sea ice to open ocean but additionally strongly impacted by the differential advection of clouds and moisture across the near-surface MCAO flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Convergence Criteria for Multiphysics Simulations

The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Insights into negative differential resistance in MoS 2 Esaki diodes: A first-principles perspective

MoS 2 is a two-dimensional material with a band gap depending on the number of layers and tunable by an external electric field. The experimentally observed intralayer band-to-band tunneling and interlayer band-to-band tunneling in this material present an opportunity for new electronic applications in tunnel field-effect transistors. However, such a widely accepted concept has yet to be been supported by theoretical investigations based on first principles. In this paper, using density functional theory, in conjunction with nonequilibrium Green's function techniques and our electric field gating method, enabled by a large-scale computational approach, we study the relation between band alignment and transmission in planar and side-stack MoS 2 p–i–n junction configurations. Here, we demonstrate the presence of negative differential resistance for both in-plane and interlayer current, a staple characteristic of tunnel diode junctions, and analyze the physical origin of such an effect. Electrostatic potentials, the van der Waals barrier, and a complex band analysis are also examined for a thorough understanding of Esaki diodes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Horizontal gene transfer of the Mer operon is associated with large effects on the transcriptome and increased tolerance to mercury in nitrogen-fixing bacteria

Abstract Background Mercury (Hg) is highly toxic and has the potential to cause severe health problems for humans and foraging animals when transported into edible plant parts. Soil rhizobia that form symbiosis with legumes may possess mechanisms to prevent heavy metal translocation from roots to shoots in plants by exporting metals from nodules or compartmentalizing metal ions inside nodules. Horizontal gene transfer has potential to confer immediate de novo adaptations to stress. We used comparative genomics of high quality de novo assemblies to identify structural differences in the genomes of nitrogen-fixing rhizobia that were isolated from a mercury (Hg) mine site that show high variation in their tolerance to Hg. Results Our analyses identified multiple structurally conserved merA homologs in the genomes of Sinorhizobium medicae and Rhizobium leguminosarum but only the strains that possessed a Mer operon exhibited 10-fold increased tolerance to Hg. RNAseq analysis revealed nearly all genes in the Mer operon were significantly up-regulated in response to Hg stress in free-living conditions and in nodules. In both free-living and nodule environments, we found the Hg-tolerant strains with a Mer operon exhibited the fewest number of differentially expressed genes (DEGs) in the genome, indicating a rapid and efficient detoxification of Hg from the cells that reduced general stress responses to the Hg-treatment. Expression changes in S. medicae while in bacteroids showed that both rhizobia strain and host-plant tolerance affected the number of DEGs. Aside from Mer operon genes, nif genes which are involved in nitrogenase activity in S. medicae showed significant up-regulation in the most Hg-tolerant strain while inside the most Hg-accumulating host-plant. Transfer of a plasmid containing the Mer operon from the most tolerant strain to low-tolerant strains resulted in an immediate increase in Hg tolerance, indicating that the Mer operon is able to confer hyper tolerance to Hg. Conclusions Mer operons have not been previously reported in nitrogen-fixing rhizobia. This study demonstrates a pivotal role of the Mer operon in effective mercury detoxification and hypertolerance in nitrogen-fixing rhizobia. This finding has major implications not only for soil bioremediation, but also host plants growing in mercury contaminated soils.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2020

Pitch Angle Dependence of Electron and Ion Flux Changes During Local Magnetic Dipolarization Inside Geosynchronous Orbit

In this study we have statistically examined flux changes of 1–1,000-keV electrons, and hydrogen, helium, and oxygen ions during local magnetic dipolarization inside geosynchronous orbit, with a focus on their pitch angle dependence. Here, using 144 dipolarization events that were selected in a previous study with Van Allen Probes observations in 2012–2016, we have performed a superposed epoch analysis of differential flux changes after the dipolarization onset for each species. On average the electron flux increases primarily around pitch angle (α) = 90° at >80 keV and almost isotropically at 10–50 keV. The electron flux at <5 keV increases at α = 0° and 180° but slightly decreases (or remains unchanged) at α = 90°. On the other hand, the ion flux at >80 keV increases around α = 90°, while at <30 keV it decreases nearly independent of pitch angle. Only the low-energy (<5 keV) helium and oxygen ion flux changes indicate a strong field-aligned enhancement, which could be attributed to their outflows from the topside ionosphere. After the dipolarization onset, >5-keV electron and >30-keV ion fluxes exhibit a perpendicular anisotropy, which is most pronounced for 50–200-keV electrons. Both distributions become more isotropic with decreasing energy. A noticeable field-aligned anisotropy is seen for <5-keV ion fluxes. These statistical pitch angle-dependent electron and ion properties during dipolarizations may be explained by combining various processes, including adiabatic acceleration and/or transport, wave-particle interactions, and ion outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗