Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “qualitative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

Role of effective mass and long-range interactions in the band-gap renormalization of photoexcited semiconductors

Understanding how to control changes in the electronic structure and related dynamical renormalizations by external driving fields is the key for understanding ultrafast spectroscopy and applications in electronics. Here, we focus on the band gap's modulation by external electric fields and uncover the effect of band dispersion on the gap renormalization. We employ the Green's function formalism using the real-time Dyson expansion to account for dynamical correlations induced by photodoping. The many-body formalism captures the dynamics of systems with long-range interactions, carrier mobility, and variable electron and hole effective mass. We also demonstrate that mean-field simulations based on the Hartree-Fock Hamiltonian, which lacks dynamical correlations, yields a qualitatively incorrect picture of band-gap renormalization. We find the trend that increasing effective mass, thus decreasing mobility, leads to as much as a 6% enhancement in band-gap renormalization. Further, the renormalization is strongly dependent on the degree of photodoping. As the screening induced by free electrons and holes effectively reduces any long-range and interband interactions for highly excited systems, we show that there is a specific turnover point with a minimal band gap. Here, we further demonstrate that the optical gap renormalization follows the same trend though its magnitude is altered by the Moss-Burstein effect.

Approximation methods for many-body systems↗

Pyrogenic Organic Matter Laboratory Experiment: Aerobic Respiration and Geochemistry from Variably Inundated Stream Sediments (v3)

This dataset supports a broader study examining the effects of variable inundation and pyrogenic organic matter on ecosystem respiration. The dataset provides data generated from a laboratory batch experiment investigating the interaction between variable inundation conditions (wet and dry sediment) and pyrogenic organic matter (burned and unburned treatments). The contents include time series dissolved oxygen, sediment geochemistry data, and field metadata (including qualitative information on instream and river corridor characteristics). This data package was originally published in November 2025. It was updated in April 2026 (v2; new and modified files) and May 2026 (v3; modified files). See the change history section in the readme for more details For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) international generic sample number (IGSN) mapping file; (5) readme; (6) field protocol; (7) sample name metadata; (8) an environmental context picture for the dry and inundated sampling locations; and (9) a subfolder with sample data from the sediment incubation experiment. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC); (2) total nitrogen (TN); (3) gravimetric moisture; (4) partial pressure and production rates of carbon dioxide, methane, and nitrous oxide; (5) field wet sediment mass, dry sediment mass, water mass, and field wet sediment volume in incubation and sediment NPOC/TN vials; (6) methods codes; (7) respiration rates, pH, and temperature from after the incubation, raw time series dissolved oxygen and temperature, and a subfolder containing associated plots and scripts; (8) ions; (9) FTICR-MS methods; and (10) a subfolder of 12 Tesla (12T) FTICR-MS data. This folder contains the CoreMS processed data and three subfolders, one containing the .xml files, one containing the CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .xml, .html, .Rmd, .py, .cal, .json, or .jpg.

54 ENVIRONMENTAL SCIENCES↗

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗

Direct reactions with the AT-TPC

Direct reactions are crucial tools for accessing properties of the atomic nucleus. Fundamental and exotic phenomena such as collective modes, pairing, weakbinding effects and evolution of single-particles energies can be investigated in peripheral collisions between a heavy nucleus and a light target. The necessity of using inverse kinematics to reveal how these structural properties change with isospin imbalance renders direct reactions a challenging technique when using the missing mass method. In this scenario, Active Target Time Projection Chambers (AT-TPC) have demonstrated an outstanding performance in enabling these types of reactions even under conditions of very low beam intensities. The AT-TPC of the Facility for Rare Isotope Beams (FRIB) is a next generation multipurpose Active Target. When operated inside a solenoidal magnet, direct reactions benefit from the measurement of the magnetic rigidity that enables particle identification and the determination of the excitation energy with high resolution without the need of auxiliary detectors. Additionally, the AT-TPC can be coupled to a magnetic spectrometer improving even further its spectroscopic investigation capability. In this contribution, we discuss inelastic scattering and transfer reaction data obtained via the AT-TPC and compare them to theory. In particular, we present the results for the 14 C(p,p′) and 12 Be (p,d) 11 Be reactions. For 14 C, we compare the experimental excitation energy of the first 1 – excited state with coupled-cluster calculationsbased on nuclear interactions from chiral effective field theory and with available shell-model predictions. For 12 Be, we determine the theoretical spectroscopic factors of the 12 Be (p,d) 11 Be transfer reaction in the shell modeland compare them to the experimental excitation spectrum from a qualitative standpoint.

active target↗

Thermally-Conductive Metallic Coatings and Applications for Heat Removal on In-Space Cryogenic Vehicles

For large in-space cryogenic upper stages, substantial axial heat removal from a forward skirt by vapor-based heat interception may not be achieved by simple attachment methods unless sufficient thermal conductance from the skirt to the cooling fluid can be achieved. Preferable methods would allow for the addition of the cooling system to existing structure with minimal impact on the structure. Otherwise, significant modification to the basic structural design andor novel and complex attachment mechanisms with high effective thermal conductance are likely to be required. The approach being pursued by evolvable Cryogenics (eCryo) is to increase the thermal performance of a relatively simple attachment system by applying metallic or other thermally conductive material coatings to the mating surface area of the fluid channel where it is attached the skirt wall. The expectation of candidate materials is that the dramatic increase in conductivity of pure metals at temperatures close to liquid hydrogen vapor temperature will compensate for the reduced actual contact area typical of mechanical joints. Basic contact conductance data at low temperatures for candidate interface materials is required to enable the test approach. A test rig was designed at NASA Glenn Research Center to provide thermal contact resistance testing between small sample coupons coated with conductive material via electron beam evaporation, a low-temperature option that will not affect physical properties of base materials. Average coating thicknesses were 10 k. The test fixture was designed to mount directly to a cryocooler cold head within a vacuum test chamber. The purpose of this test was to determine qualitative contact conductance between various test samples. Results from this effort will be implemented in a sub-scale vapor-based heat interception test, where the applicability for increased heat removal on large structural skirts will be considered.

small scale testing↗

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗

Evaluating the Importance of Conformers for Understanding the Vacuum-Ultraviolet Spectra of Oxiranes: Experiment and Theory

Vacuum-ultraviolet (VUV) absorption spectroscopy enables electronic transitions that offer the unambiguous identification of molecules. As target molecules become more complex, multifunctional species present a great challenge to both experimental and computational spectroscopy. This research reports both experimental and theoretical studies of oxiranes. Computationally, the nuclear ensemble approach has been used to accurately predict experimental spectra for a variety of molecules. However, this approach incurs great computational cost, as ensembles generally consist of thousands of geometries. The present study aims to drastically reduce the ensemble by evaluating the significance of the conformers to the predicted spectra. This approach was applied to 11 substituted oxiranes using the Conformer Rotamer Ensemble Sampling Tool (CREST) of Grimme to generate an ensemble of unique conformers determined by their Boltzmann populations. Five TD-DFT functionals (BMK, CAM-B3LYP, M06-2X, MN15, ωB97X-D) and EOM-CCSD were used to simulate the spectrum of each substituted oxirane ensemble. Computed spectra were then compared to the experiment using both qualitative and quantitative metrics. Based on these metrics, it was observed that certain conformers may not be necessary to characterize this set of oxiranes despite the temperature (323 K) of the experiment. A single conformer can then be used with TD-DFT and EOM-CCSD to replicate the experimental spectra of these medium-sized combustion species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Planetary and satellite x ray spectroscopy: A new window on solid-body composition by remote sensing

The rings and most of the satellites of the outer planets orbit within the radiation belts of their parent bodies. This is an environment with intense fluxes of energetic electrons. As a result, these objects are strong emitters of X-rays. The characteristic X-ray lines from these bodies depend on atomic composition, but they are not sensitive to how the material is arranged in compounds or mixtures. X-ray fluorescence spectral analysis has demonstrated its unique value in the laboratory as a qualitative and quantitative analysis tool. This technique has yet to be fully exploited in a planetary instrument for remote sensing. The characteristic X-ray emissions provide atomic relative abundances. These results are complementary to the molecular composition information obtained from IR, visible, and UV emission spectra. The atomic relative abundances are crucial to understanding the formation and evolution of these bodies. They are also crucial to the proper interpretation of the molecular composition results from the other sensors. The intensities of the characteristic X-ray emissions are sufficiently strong to be measured with an instrument of modest size. Recent developments in X-ray detector technologies and electronic miniaturization have made possible space-flight X-ray imaging and nonimaging spectrometers of high sensitivity and excellent energy resolution that are rugged enough to survive long-duration space missions. Depending on the application, such instruments are capable of resolving elemental abundances of elements from carbon through iron. At the same time, by measuring the bremsstrahlung intensity and energy spectrum, the characteristics of the source electron flux can be determined. We will discuss these concepts, including estimated source strengths, and will describe a small instrument capable of providing this unique channel of information for future planetary missions. We propose to build this instrument using innovative electronics packaging methods to minimize size and weight.

Chenette, D. L.↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗

Extracting Oxygen from Lunar Simulant Using a Transparent Furnace Pulsed Fluidized Bed

In the event that humans return to the moon, utilizing the local materials will be beneficial for extended stays. Rather than transporting resources, such as oxygen, from Earth, NASA is investigating methods of extracting it from lunar regolith. One promising process is hydrogen reduction. In the hydrogen reduction process, lunar regolith is heated to 1000 C in the presence of hydrogen. The iron oxide (Fe-O) bonds, found in lunar material, are broken and the hydrogen attracts the oxygen to produce water vapor [Allen et al., 1996]. FeO + H2 (right arrow) Fe +H2O. The water vapor is then captured, cleaned, and electrolyzed. The hydrogen is recycled back to the reduction process and the oxygen is stored until consumed by an end user (propulsion, life support, etc.). To obtain a good oxygen yield, the majority of lunar regolith must be exposed to the hydrogen gas and have a high rate of heat transfer from heat source to particle. This is achieved with good solids mixing via fluidization or mechanical agitation. In Generation II of the ROxygen program, the ROxygen Team at Johnson Space Center (JSC) investigated the feasibility of gas only pulsed fluidization as the only means to mix synthetic lunar regolith (simulant) at high temperatures. Fluidized beds have been used in industry to effectively process powders for decades. They consist of gas flowing upward through a bed of particles. The stirring action continuously moves the grains around to achieve uniform mixing of gas, solids, and heat [Geldart, 1986]. A transparent furnace unit was developed by Thoughventions Unlimited LLC (TvU) to aid in the qualitative observation of the fluidization behavior at high temperatures. Multipoint thermocouples and pressure sensors provided quantitative information regarding the quality of mixing. The water produced was measured using humidity sensors and captured using a NASA designed and built condenser. Once the simulant was processed, pneumatically transporting the 'hot' simulant out of the furnace was investigated.

Oryshchyn, L.↗

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗

Skeletal muscle responses to unloading in humans

This study examined the effects of unloading on skeletal muscle structure. Method: Eight subjects walked on crutches for six weeks with a 110 cm elevated sole on the right shoe. This removed weight bearing by the left lower limb. Magnetic resonance imaging of both lower limbs and biopsies of the left m. vastus laterallis (VL) were used to study muscle structure. Results: Unloading decreased (P less than 0.05) muscle cross-sectional areas (CSA) of the knee extensors 16 percent. The knee flexors showed about 1/2 of this response (-7 percent, P less than 0.05). The three vasti muscles each showed decreases (P less than 0.05) of about 15 percent. M. rectus femoris did not change. Mean fiber CSA in VL decreased (P less than 0.05) 14 percent with type 2 and type 1 fibers showing reductions of 15 and 11 percent respectively. The ankle extensors showed a 20 percent decrease (P less than 0.05) in CSA. The reduction for the 'fast' m. gastrocnemius was 27 percent compared to the 18 percent decrease for the 'slow' soleus. Summary: The results suggest that decreases in muscle CSA are determined by the relative change in impact loading history because atrophy was (1) greater in extensor than flexor muscles, (2) at least as great in fast as compared to slow muscles or fibers, and (3) not dependent on single or multi-joint function. They also suggest that the atrophic responses to unloading reported for lower mammals are quantitatively but not qualitatively similar to those of humans.

Dudley, G.↗

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Biofouling and Corrosion Study for a Novel Linear Guided Wave Energy Converter (Post Access Report)

The overall objective of this project was to examine the reliability and performance of antibiofouling coatings used for a wave energy converter (WEC) developed by E-Wave Technologies. The particular coatings were selected for their low toxicity and potential compatibility with aquaculture. The aim of this work was to 1) test coating solutions to prevent biofouling growth and saltwater corrosion on the static (paddle and attachment frame surface) components of the WEC that are submerged, 2) determine adhesion of the coatings to system components, and 3) assess the ease and effectiveness of biofouling cleaning to insure long term performance of the system. An analysis of commercial coatings was performed using methods to examine the prevention of biofouling and coating adhesion properties on two key materials of the WEC, which were 316L low carbon marine grade stainless steel (SS) and Ultra High Molecular Weight Polyethylene (PE). Three marine antifouling paints were selected based on their unique properties to test how different paint styles perform on different materials. The selected paints were ePaint Ecominder self-polishing paint with Zinc Omadine for slime control, Pettit ECO HRT Copper-Free ablative antifouling with Econea biocide, and Intersleek 1100SR foul release. Pacific Northwest National Laboratory (PNNL) prepared PE and SS substrates coated with the three paints and compared the performance against uncoated substrates when submerged in raw seawater for 3-, 6-, and 9-month (m) time periods. Results in adhesion testing indicated that Pettit and ePaint materials clearly bonded strongly to SS, but did not bond comparably well to PE. It was noted during adhesion testing that the Intersleek surfaces were especially difficult to test as the paint highly resists bonding to the epoxy adhesives used with the adherence testing platform. The wear rate of the coatings was not measured under this study; however, based on adhesion testing, coatings in the sliding regions of the Ewave device are expected to wear rapidly. Sandia National Laboratories (SNL) evaluated the adhesion of three different paints to PE and SS substrates which were exposed to a marine environment for time intervals of 0, 3, 6, and 9 months. From qualitative visual analysis of the 3 in2 coupons when pulled from the tank, the 3 in2 coupons generally only appeared to have biofouling consisting of filamentous algae or diatoms, which all have relatively low mass and can be easily wiped from the surface of coupons. Qualitative visual analysis indicated that ECO HRT and Unpainted were consistently worse than Intersleek and Ecominder at all time points. Results provide insight to aid with down-selection of commercial coatings under static conditions to support reliability of the WEC and potential maintenance schedules. This investigation was conducted using small coupon samples suspended in seawater and the development of testing rigs for dynamic component level testing is needed for future work. In addition, the potential toxicity of these commercial coatings on aquaculture has not been determined by this study. One recommendation is to conduct toxicity investigations at the Environmental Toxicity Laboratory at Oak Ridge National Laboratory.

16 TIDAL AND WAVE POWER↗

A Gridded Solar Irradiance Ensemble Prediction System Based on WRF-Solar EPS and the Analog Ensemble

The WRF-Solar Ensemble Prediction System (WRF-Solar EPS) and a calibration method, the analog ensemble (AnEn), are used to generate calibrated gridded ensemble forecasts of solar irradiance over the contiguous United States (CONUS). Global horizontal irradiance (GHI) and direct normal irradiance (DNI) retrievals, based on geostationary satellites from the National Solar Radiation Database (NSRDB) are used for both calibrating and verifying the day-ahead GHI and DNI predictions (GDIP). A 10-member ensemble of WRF-Solar EPS is run in a re-forecast mode to generate day-ahead GDIP for three years. The AnEn is used to calibrate GDIP at each grid point independently using the NSRDB as the “ground truth”. Performance evaluations of deterministic and probabilistic attributes are carried out over the whole CONUS. The results demonstrate that using the AnEn calibrated ensemble forecast from WRF-Solar EPS contributes to improving the overall quality of the GHI predictions with respect to an AnEn calibrated system based only on the deterministic run of WRF-Solar. In fact, the calibrated WRF-Solar EPS’s mean exhibits a lower bias and RMSE than the calibrated deterministic WRF-Solar. Moreover, using the ensemble mean and spread as predictors for the AnEn allows a more effective calibration than using variables only from the deterministic runs. Finally, it has been shown that the recently introduced algorithm of correction for rare events is of paramount importance to obtain the lowest values of GHI from the calibrated ensemble (WRF-Solar EPS AnEn), qualitatively consistent with those observed from the NSRDB.

14 SOLAR ENERGY↗