Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional expectations”

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 289 records · Page 16

Investigating parameter trainability in the SNAP-displacement protocol of a qudit system

Here, in this study, we explore the universality of Selective Number-dependent Arbitrary Phase (SNAP) and Displacement gates for quantum control in qudit-based systems. However, optimizing the parameters of these gates poses a challenging task. Our main focus is to investigate the sensitivity of training any of the SNAP parameters in the SNAP-Displacement protocol. We analyze conditions that could potentially lead to the Barren Plateau problem in a qudit system and draw comparisons with multi-qubit systems. The parameterized ansatz we consider consists of blocks, where each block is composed of hardware operations, namely SNAP and Displacement gates [Fösel et al 2020 Efficient cavity control with snap gates arXiv:2004.14256]. Applying Variational Quantum algorithm (VQA) with observable and gate cost functions, we utilize techniques similar to those in [McClean et al 2018 Barren plateaus in quantum neural network training landscapes Nat. Commun.9 1–6] and [Cerezo et al 2021 Cost function dependent barren plateaus in shallow parametrized quantum circuits. Nat. Commun.12 1–12] along with the concept of t − design. Through this analysis, we make the following key observations: (a) The trainability of a SNAP-parameter does not exhibit a preference for any particular direction within our cost function landscape, (b) By leveraging the first and second moments properties of Haar measures, we establish new lemmas concerning the expectation of certain polynomial functions, and (c) utilizing these new lemmas, we identify a general condition that indicates an expected trainability advantage in a qudit system when compared to multi-qubit systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interference of electrical double layers: Confinement effects on structure, dynamics, and screening of ionic liquids

Ionic liquids are widely used as electrolytes in electronic devices in which they are subject to nanoconfinement within nanopores or nanofilms. Because the intrinsic width of an electrical double layer is on the order of several nanometers, nanoconfinement is expected to fundamentally alter the double layer properties. Furthermore, in confined systems, a large portion of the ions are interfacial, e.g., at the electrode interface, leading to significant deviations of electrostatic screening and ion dynamics as compared to bulk properties. In this article, we systematically investigate the interference between electrical double layers for nanoconfined ionic liquids and the resulting influence on the structure, dynamics, and screening behavior. We perform molecular dynamics simulations for the ionic liquids [BMIm + ][BF$^-_4$] and [BMIm + ][PF$^-_6$] confined between two flat electrodes at systematic separation distances between 1.5 nm and 4.5 nm for both conducting and insulating boundary conditions. We find that while ion dynamics is expectedly slower than in the bulk (by ~2 orders of magnitude), there is an unexpected non-linear trend with the confinement length that leads to a local maximum in dynamic rates at ~3.5–4.5 nm confinement. We show that this nonlinear trend is due to the ion correlation that arises from the interference between opposite double layers. We further evaluate confinement effects on the ion structure and capacitance and investigate the influence of electronic polarization of the ionic liquid on the resulting properties. This systematic evaluation of the connection between electrostatic screening and structure and dynamics of ionic liquids in confined systems is important for the fundamental understanding of electrochemical supercapacitors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Efficiency-optimized relativistic plasma harmonics for extreme fields

Bright harmonic radiation from relativistically oscillating laser plasmas offers a direct route for generating extreme electromagnetic fields. Theory predicts that under optimized conditions, the plasma medium can support strong spatiotemporal compression of laser energy in a coherent harmonic focus (CHF), delivering intensity boosts many orders of magnitude greater than the incident driving laser pulse. Although diffraction-limited performance (spatial compression) and attosecond phase locking (temporal compression) have been demonstrated experimentally, efficient coupling of relativistically intense laser pulse energy into the emitted harmonic cone has not been realized so far. Here we demonstrate that this highly nonlinear interaction can be tailored to deliver the maximum conversion efficiencies predicted from simulations. By fine-tuning the temporal profile of the driving laser on sub-picosecond (<10 −12 s) timescales, energies >9 mJ between the 12th and 47th harmonics are observed. These results are in agreement with the theoretically expected efficiency dependence on harmonic order, verifying that optimal conditions have been achieved in the generation process. This is the important final element required to achieve the expected intensity boosts from a CHF in experiments. Although obtaining spatiotemporal compression and optimal efficiency simultaneously remains challenging, the path to realizing extreme optical field strengths approaching the critical field of quantum electrodynamics (the Schwinger limit at >10 16 V cm −1 or >10 29 W cm −2 ) is now open, permitting all-optical studies of the quantum vacuum and new frontiers for intense attosecond science.

High-harmonic generation↗

Improving Volatile Fatty Acid Productivity of Anaerobic Digestion

Typical anaerobic digestion (AD) focusses on complete conversion of waste to biogas (primarily carbon dioxide and methane). However, the intermediate metabolites of the AD process, which includes short- and long-chain volatile fatty acids (VFAs), that could serve as the precursors for useful industrial applications are typically ignored. The goal of this project is to eliminate production of biogas while enhancing the production of the intermediate VFAs. Using a mixed microbial consortium (from rumen sources and waste-water sludge), we have determined the optimal carbon loading (chemical oxygen demand, "COD") and optimal pH, that results in high VFA concentrations from food waste. The best VFA yields were obtained using 15 g COD/L and a pH of 9.0 for this substrate. pH 9 produced over 200% higher VFA titers than the controlled conditions (i.e., ph 7.0) after 35 d of digestion, in comparison to pH 5 that showed - 88% higher titers than the control digestion. As expected, the cumulative biogas production was highest in the pH 7.0 condition, in comparison to pH 5.0 or pH 9.0. We further observed that removal of VFAs using solid-liquid separation technique reduce the inhibitory effects of VFAs, thereby leading to overall improvement in conversion efficiency. 16s rRNA analysis is being carried out to explain and identify the biocatalysts that enable VFA production in these AD cultures. Additional efforts to improve VFA yields via increasing the total solid content, temperature optimizations, VFA removal via electrodialysis, and improving hydrolysis via microaeration will be presented.

anaerobic digestion↗

Assessment of Condition Monitoring Methods and Technologies for Inservice Inspection and Testing of Nuclear Power Plant Components

This report was prepared for the U.S. Nuclear Regulatory Commission (NRC) to explore the application of advanced technologies toward meeting the current and future regulatory requirements for maintenance and condition monitoring of structures, systems, and components. The advanced technologies considered in this work are advanced sensors and instrumentation, data analytics, machine learning and artificial intelligence (ML/AI), physics-based models, and digital twins (DT). The interest in the application of advanced technologies for condition monitoring in nuclear power plants continues to grow, and current and future licensees are expected to implement advanced technologies as part of their inservice inspection (ISI) and inservice testing (IST) programs. This report delineates the outcomes of an exploratory investigation into the implementation of advanced condition monitoring technologies to address ISI and IST requirements. A thorough review was conducted of the existing regulatory requirements for ISI and IST, along with an analysis of associated industry practices. Additionally, a state-of-the-art assessment was performed on advanced condition monitoring technologies frequently employed in non-nuclear sectors. This research incorporated two nuclear-specific case studies to illustrate the application of these technologies within the current nuclear fleet. The report provides an exhaustive discussion on the technical challenges, considerations, and opportunities associated with the deployment of advanced condition monitoring technologies. The following are key considerations in the application of advanced technologies for the ISI and IST of nuclear power plant components: • Developing adequate verification and validation procedures to confirm the functional and non-functional requirements, • Developing technical capabilities to conduct real-time asset condition monitoring, • Establishing guidance and protocol for modeling and simulation tools to continuously meet regulatory requirements, • Addressing trustworthiness, explainability, and interpretability of ML/AI methods, • Evaluating maintenance activities to maintain an adequate safety margin and avoid undesirable conditions, • Establishing cybersecure condition monitoring programs associated with a computer-based software system, and • Establishing standardized evaluation metrics for advanced condition monitoring programs. Interest in the use of advanced technologies for condition monitoring in ISI and IST programs continues to grow, and the technology is expected to experience rapid and wide industry adoption in the near future. Adoption of advanced technologies for condition monitoring could have novel and unique impacts on regulatory activities associated with ISI and IST programs. The NRC is continuing to explore the regulatory aspects of advanced technologies as part of ISI and IST programs by pursuing additional research in this technical area.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of Additively Manufactured Complex Tools for Autoclave Cure Composites

IACMI Project 4.9, Tooling for Composites with Washout Features Produced by Additive Manufacturing, assembled a team including the industry lead, Ability Composites, NREL and Colorado State University (CSU). Ability Composites had originally expressed interest in alternate methods of producing tooling for composite parts. In follow-up discussions, it became clear that one of the bigger tooling challenges revolved around small production volume composite parts that were tooled on washout material due to the complex geometry. To build an understanding of the potential, both from a technology and a cost perspective, for replacing conventional washout tooling with 3D printed thermoplastic tooling, a number of commercially available dissolvable FDM printing materials were evaluated, leading to tooling representative of commercial articles of interest to Ability Composites. Ultimately, Ability Composites was able to directly compare autoclave processed prepreg composite parts produced on conventional washout tooling to composite parts molded on 3D printed dissolvable tooling produced at CSU. Small, laboratory test specimens were developed to investigate the structural performance of the candidate materials under autoclave processing conditions, which were nominally 121 °C (250 °F) and 345 kPa (50 psi). In addition, several internal structural configurations (infills) were evaluated under autoclave conditions using model materials. The results of these tests indicated that two materials, Stratasys ST 130 and Infinite Materials Solutions Aquasys 180 (AQ 180), were the best candidates, given the specified autoclave processing conditions. ST-130 was slightly more robust than AQ-180; however, the AQ-180 was carried forward as it was dissolvable in water, not requiring the basic solution needed to dissolve ST-130. Based on the preliminary material and 3D printed structures evaluations, larger tools with a truncated square pyramid geometry were created to produce prepreg composite test articles for 3D printed dissolvable tool evaluation under standard autoclave fabrication conditions. Two tools were manufactured using ST-130 and one tool using traditional ceramic washout tooling media. The tools were evaluated for geometric fidelity and surface roughness changes before and after carbon fiber/epoxy prepreg composites were manufactured on the tooling. The autoclave processing did not impact the geometry significantly and was completed at 121 °C and 345 kPa, indicating satisfactory tool performance. The results from surface roughness testing of both the resulting composite and the associated tooling indicated that an adequate surface resulted without the need for a surface sealing step, as was required for the conventional washout tooling. Based on results of the truncated pyramid tests as a basis, ST-130, AQ-120 and AQ-180 materials were carried forward to the tool geometry of interest to Ability Composites. These hollow rectangular bent ducts, which were complex in nature and not extractable after cure, were used to understand the impacts of tool material and thickness. One ST-130 tool was produced as a partially solid part, with an enclosed 40% dense infill region to reduce weight and material use. This was the same approach evaluated in the truncated pyramid portion of the study. This tool was to be envelope vacuum bagged and directly compared to a monolithic tool of conventional washout material. The traditional monolithic ceramic tool was manufactured by Ability composites using CNC-based subtractive methods. An additional five dissolvable polymer tools, manufactured from ST-130, AQ-120, and AQ-180, using a hollow design were 3D printed and used to produce carbon fiber/epoxy composite evaluation articles. These hollow dissolvable tools were expected to be less influenced by the autoclave conditions as the wall was solid and vacuum bagging was inside and outside the tool. This alternative geometry was also evaluated as an option in techno-economic modeling. Print times were reduced from in excess of 3 days to under 30 hours, while surface quality and and tool integrity were substantially improved in the transition from the partially solid tool to the hollow tooling concept. Ability Composites produced autoclave-cured prepreg ducts on each of these tools. The autoclave conditions utilized were more severe than those of the initial trials, reaching temperatures of 160 °C and a pressure of 414 kPa. Under these conditions, the partially solid 3D printed tool with skin and 40% dense infill crushed significantly; however, the thicker ST-130 hollow tool showed good promise, deforming only slightly. The thinner hollow tool walls were unsuccessful as were the other materials. Overall, the hollow tool manufacturing process saved significant amounts of time and material in manufacturing as compared to the solid ducts and produced composite surface quality improvements compared to the traditional washout tooling. The TEM was developed to allow direct comparisons between conventional washout tool manufacturing processes and those developed at CSU. It also allowed for two separate 3D printed tool geometries to be analyzed and compared. In this case, the square bent duct tool geometry was determined to be representative of common washout tools. This geometry was compared with a scaled-up version of it to assess differences in the two manufacturing processes based on tool size. The model was developed to make use of user input in the form of geometry details, process steps, manufacturing parameters, bulk material costs, capital equipment costs, and general costs to calculate overall labor, material, capital equipment, and energy costs per manufactured tool for the conventional and additive manufacturing processes for the two representative geometries. It was also able to estimate step-by-step process times for the manufacturing process and geometries. Based on significant input from Ability Composites and CSU from their knowledge gained from hands-on manufacturing of the 3D printed bent duct tool geometry, costs and process times were calculated for the two manufacturing processes. Results showed that the additive manufacturing techniques developed at CSU can substantially reduce the costs of tool manufacturing by reducing labor times and material usage. This is because additive manufacturing is a relatively hands-off process and allows for the tool design to be optimized to reduce material usage. The disadvantage, however, is that process times for additive manufacturing are significantly longer. The three-dimensional (3D) printing process is slow if tight tolerances are required, but the analysis did show that print times could be reduced with the hollow tool geometry. Also, further advances in additive manufacturing could expedite the process. Costs and process times for the tool washout process were calculated separately. They showed that costs are relatively insignificant when compared to the overall tool manufacturing processes, but with increases in tool size, costs for the conventional manufacturing approach are larger than for additive manufacturing. Again, the washout process for conventional tools is very hands-on, whereas for additively manufactured tools the print medium is dissolved in an automated detergent bath at the sacrifice of process time. The analysis showed that optimizing the additively manufactured tools may also reduce washout times. Overall, Project 4.9 demonstrated that commercially available dissolvable 3D printing materials exist that can be used to produce dissolvable tooling capable of surviving prepreg composites fabrication under autoclave conditions of 121 °C (250 °F) and 345 kPa (50 psi). An alternative hollow dissolvable tool design was developed which was structurally superior to the initial concept and was cost and time effective versus conventional washout tooling. The 3D printed sacrificial tool required no added surface sealing steps prior to composite part layup and cure, offering a significant advantage over the porous conventional washout tooling.

36 MATERIALS SCIENCE↗

Integrated Sustainable Product Design With Warranty and End-of-Use Considerations

The concept of integrated sustainable product design has recently emerged, aiming to incorporate downstream life cycle performance into the initial product design to enhance sustainability. Various sustainable product design tools based on life cycle assessment or quality function deployment have been established while the impact of reliability on circular practices has received limited attention. Recognizing the critical role of product reliability in post-design performance, this paper develops a product design optimization model that considers the warranty performance and the effect of end-of-use options. The model takes into account the effect of uncertain operating conditions on product reliability. Two optimization goals including the minimization of expected unit life cycle cost and environmental impact are achieved by the model. To demonstrate the benefits of the integrated approach, the model is applied to an electric motor design problem. The results highlight that integrating end-of-use options in the early design phase leads to adjustments in component selection and reliability design. Moreover, the circular utilization of used products enables cost savings throughout the product’s life cycle and contributes to environmental impact reduction. Finally, the study analyzes the effects of operating conditions, warranty policies, and take-back prices for used products on design decisions, providing valuable insights for product designers.

Integrated design↗

Physics-based model of irradiation creep for ferritic materials under fusion energy operation conditions

Irradiation creep is known to be an important process for structural materials in nuclear environments, potentially leading to creep failure at temperatures where thermal creep is generally negligible. While there is a great deal of data for irradiation creep in steels and zirconium alloys in light water reactor conditions, much less is known for first wall materials under fusion energy conditions. Lacking suitable fusion neutron sources for detailed experimentation, modeling, and simulation can help bridge the dose-rate and spectral-effects gap and produce quantifiable expectations for creep deformation of first wall materials under standard fusion conditions. In this paper, we develop a comprehensive model for irradiation creep created from merging a crystal plasticity representation of the dislocation microstructure and a defect evolution simulator that accounts for the entire cluster dimensionality space. Both approaches are linked by way of a climb velocity that captures dislocation-biased defect absorption and a dislocation strengthening term that reflects the accumulation of defect clusters in the system. Here, we carry out our study in Fe under first wall fusion reactor conditions, characterized by a fusion neutron spectrum with average recoil energies of 20 keV and a damage dose rate of ≈3×10 –7 dpa/s at temperatures between 300 and 800 K.

36 MATERIALS SCIENCE↗

Irradiation creep model for α-Fe: Dislocation climb in the presence of stress fields (Final Report)

Irradiation creep is known to be an important process for structural materials in nuclear environments, potentially leading to creep failure at temperatures where thermal creep is generally negligible. While there is a great deal of data for irradiation creep in steels and zirconium alloys in light water reactor conditions, much less is known for first wall materials under fusion energy conditions. Lacking suitable fusion neutron sources for detailed experimentation, modeling, and simulation can help bridge the dose-rate and spectral-effects gap and produce quantifiable expectations for creep deformation of first wall materials under standard fusion conditions. In this work, we present a comprehensive model for irradiation creep created from merging a crystal plasticity representation of the dislocation microstructure and a defect evolution simulator that accounts for the entire cluster dimensionality space. Both approaches are linked by way of a climb velocity that captures dislocation-biased defect absorption and a dislocation strengthening term that reflects the accumulation of defect clusters in the system. We carry out our study in Fe under first wall fusion reactor conditions, characterized by a fusion neutron spectrum at temperatures between 300 and 800°K.

36 MATERIALS SCIENCE↗

Arctic shrub expansion generates regional variation in litter decomposition by altering litter quality and the decomposition environment

Abstract The expansion of deciduous shrubs into the graminoid‐dominated arctic tundra is expected to alter litter decomposition by changing litter quality and local abiotic and biotic conditions. However, it remains unclear how shrub expansion will affect litter decomposition at regional scales, where macroclimate is expected to be the dominant regulator of decomposition. To determine the relative influence of macroclimate and local controls on regional patterns of litter mass loss and nitrogen release, we conducted two hierarchical litter decomposition experiments across spatial scales. We decomposed leaf and root litter from a prominent graminoid ( Eriophorum vaginatum ) and three genera of deciduous shrubs ( Betula , Alnus and Salix ) for 1 year within replicated plots at five sites spanning a 160 km latitudinal gradient in northern Alaska. Using Eriophorum litter as a substrate, we found that macroclimate was the primary regulator of mass loss but had opposing effects on leaf and root litter. As summer temperature increased along the latitudinal gradient (11.9 to 13.9°C), leaf litter mass loss increased by 20% whereas root litter mass loss decreased by 33%. Leaf nitrogen release also increased with summer temperature. Conversely, root nitrogen release was controlled by the vegetation type of the decomposition environment. Using different shrub litters as substrates, we found that litter quality and its interaction with soil microclimate and macroclimate controlled decomposition. Overall, shrub root litters decomposed faster than Eriophorum root litter, losing 53% more mass and 190% more nitrogen across all sites and decomposition environments. For leaf litter, however, patterns varied by litter genus, with Salix losing more mass and Betula and Alnus losing less mass than Eriophorum . Our findings demonstrate that shrub expansion in the Arctic can regulate leaf and root litter decomposition at the regional scale through its effects on local controls, primarily litter quality. Ongoing increases in shrub cover are likely to accelerate the turnover of root litter carbon and nitrogen pools in tundra ecosystems. Therefore, including shrub‐related processes in Earth system models will improve our ability to predict regional‐scale litter decomposition and its effects on carbon and nutrient cycling in a warming Arctic. Read the free Plain Language Summary for this article on the Journal blog.

Vozzo, Justin T. [Department of Natural Resources ↗

FEA-aided investigation of the effective thermal conductivity in a medium with embedded spheres

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. Here, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing (TRISO) particles. A review on effective thermal conductivity (ETC) modeling is performed considering a random distribution of mono-sized spherical inclusions in a continuous matrix, with a primary focus on the analytical models. Finite element simulations are performed to evaluate each analytical model in varying material conditions. The model predictions are compared with the expected results obtained from the finite element predictions in addition to the Wiener and Hashin-Shtrikman bounds. Lastly, we included a practical discussion on the homogenization applied to a TRISO fuel pebble.

42 ENGINEERING↗

Neutron-proton pair transfer reactions and corresponding Weisskopf-type units

We introduce the concept of neutron-proton two-particle units (np-Weisskopf units) to be used in the analysis of the ( 3 He,p) and (p, 3 He) reactions on nuclei along the N=Z line. These are presented for the conditions relevant to the (n,j,ℓ) orbits expected from 16 O to 100 Sn. As is the case of the Weisskopf units for electromagnetic transitions, the np-WU's will provide a simple, yet robust, measure of isoscalar and isovector np pairing collective effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Worldwide benchmarking of cost-effective radiometers for direct and diffuse irradiance

Solar energy projects can benefit from direct normal irradiance (DNI) and diffuse horizontal irradiance (DHI) measurements during all project phases. Several commercial measurement systems for DNI and DHI are available. Sun trackers with pyranometers and pyrheliometers can provide highly accurate measurements but are often impractical in solar energy applications. For less expensive and more robust sensors, it is often unclear which accuracy can be expected under a project site's specific atmospheric conditions. We address this challenge through our dedicated experimental comparison of relevant sensor systems (rotating shadowband irradiometer [short RSI], Delta-T SPN1, EKO MS-90, PyranoCam, Sunto CaptPro, Kipp & Zonen CSD3) at up to six sites worldwide. The RSI systems (rRMSD 3 to 8.6%, DNI; 4.8 to 7.6%, DHI) and PyranoCam (rRMSD 2.6 to 5.2%, DNI; 4.4 to 5.8%, DHI) exhibit similar error metrics and are the most accurate systems in the test. Delta-T SPN1 and EKO MS-90 (rRMSD 6.8 to 15%, DNI; 10.6 to 20.1%, DHI) but especially Kipp & Zonen CSD3 and Sunto CaptPro show significant deviations (rRMSD 17.7 to 20%, DNI; 33 to 58%, DHI). We evaluate the influence of relevant atmospheric parameters on the sensors' accuracies by a rather unique measurement setup. MS-90's DNI errors depend on DNI itself, with overestimations for low reference DNI. The deviations of SPN1's DHI and DNI measurements increase sharply in situations with high circumsolar irradiance. Also CaptPro and CSD3's increased measurement errors are related to circumsolar irradiance. For RSI and PyranoCam, only moderate influences on the measurements are identified, indicating a general applicability of these instruments.

14 SOLAR ENERGY↗

Improving Time Step Convergence in an Atmosphere Model With Simplified Physics: The Impacts of Closure Assumption and Process Coupling

Convergence testing is a common practice in the development of dynamical cores of atmospheric models but is not as often exercised for the parameterization of sub-grid physics. An earlier study revealed that the stratiform cloud parameterizations in several predecessors of the Energy Exascale Earth System Model (E3SM) showed strong time-step sensitivity and slower-than-expected convergence when the model's time step was systematically refined. In this work, a simplified atmosphere model is configured that consists of the spectral-element dynamical core of the E3SM atmosphere model coupled with a large-scale condensation parameterization based on commonly used assumptions. This simplified model also resembles E3SM and its predecessors in the numerical implementation of process coupling and shows poor time-step convergence in short ensemble tests. We present a formal error analysis to reveal the expected time-step convergence rate and the conditions for obtaining such convergence. Numerical experiments are conducted to investigate the root causes of convergence problems. We show that revisions in the process coupling and closure assumption help to improve convergence in short simulations using the simplified model; the same revisions applied to a full atmosphere model lead to significant changes in the simulated long-term climate. This work demonstrates that causes of convergence issues in atmospheric simulations can be understood by combining analyses from physical and mathematical perspectives. Addressing convergence issues can help to obtain a discrete model that is more consistent with the intended representation of the physical phenomena.

54 ENVIRONMENTAL SCIENCES↗

Charge regulation mechanism in end-tethered weak polyampholytes

Weak polyampholytes, containing oppositely charged dissociable groups, are expected to be responsive to changes in ionic conditions. Here, we determine structural and thermodynamic properties, including the charged groups’ degrees of dissociation, of end-tethered weak polyampholyte layers as a function of salt concentration, pH, and the solvent quality. For diblock weak polyampholytes grafted by their acidic blocks, we find that the acidic monomers increase their charge while the basic monomers decrease their charge with decreasing salt concentration for pH values less than the pK a value of both monomers and vice versa when the pH > pK a . This complex charge regulation occurs because the electrostatic attraction between oppositely charged blocks is stronger than the repulsion between monomers with the same charge in both good and poor solvents when the screening by salt ions is weak. This is shown by the retraction of the top block into the bottom layer. In the case of poor solvent conditions to the basic block (the top block), we find lateral segregation of basic monomers into micelles, forming a two-dimensional hexagonal pattern on the surface at intermediate and high pH values for monovalent salt concentrations from 0.01 to 0.1 M. When the solvent is poor to both blocks, we find lateral segregation of the grafted acidic block into lamellae with longitudinal undulations of low and high acidic monomer density. By exploiting weak block polyampholytes, our work expands the parameter space for creating responsive surfaces stable over a wide range of pH and salt concentration.

36 MATERIALS SCIENCE↗

Model for photothermal ionization and molecular recombination during pulsed ablation of polyethylene

The surface erosion of polyethylene is studied under conditions characteristic of extremely high-rate heating expected in laser-driven and high energy density experiments. A coupled model for photoionization and volumetric photothermal degradation of polyethylene is developed that includes finite rate chemistry and molecular recombination of radicals. First, the model is used to explore the evolution of radicals and hydrocarbons during isochoric heating, for which it is found that polyethylene is not in chemical equilibrium for rates >108K/s. Then, the model is used to explore the cooperativity between photoionization and photothermal ionization during one-dimensional ablation from a pulsed heating source, accounting for coupled energy deposition, thermokinetics, thermochemistry, hydrodynamics, mechanics, and thermal conduction in the ALEGRA multi-physics code. The ablation depth per energy pulse is found to agree with an analytical model for instantaneous energy deposition and absent molecular recombination or photoionization. Otherwise, molecular recombination of radicals reduces the depth of ablation per pulse, owing to thermal stabilization on forming branched molecules. Radical formation through photoionization can be compensated for by recombination, but it accelerates ablation at high enough photoionization rates. Finally, the influence of molecular recombination on the ablation depth is found to decrease for long pulses, owing to decreased recombination reaction rates at lower vapor densities near the ablated surface. These and other findings highlight competition between the processes of scission, recombination, and photoionization during pulsed ablation of polyethylene.

Moore, N. W. (ORCID:0000000272259674)↗

Inferring demographic and selective histories from population genomic data using a 2-step approach in species with coding-sparse genomes: an application to human data

Abstract The demographic history of a population, and the distribution of fitness effects (DFE) of newly arising mutations in functional genomic regions, are fundamental factors dictating both genetic variation and evolutionary trajectories. Although both demographic and DFE inference has been performed extensively in humans, these approaches have generally either been limited to simple demographic models involving a single population, or, where a complex population history has been inferred, without accounting for the potentially confounding effects of selection at linked sites. Taking advantage of the coding-sparse nature of the genome, we propose a 2-step approach in which coalescent simulations are first used to infer a complex multi-population demographic model, utilizing large non-functional regions that are likely free from the effects of background selection. We then use forward-in-time simulations to perform DFE inference in functional regions, conditional on the complex demography inferred and utilizing expected background selection effects in the estimation procedure. Throughout, recombination and mutation rate maps were used to account for the underlying empirical rate heterogeneity across the human genome. Importantly, within this framework it is possible to utilize and fit multiple aspects of the data, and this inference scheme represents a generalized approach for such large-scale inference in species with coding-sparse genomes.

Soni, Vivak (ORCID:0000000294969562)↗