Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface coupling”

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 955 records · Page 53

Can electric fields drive chemistry for an aqueous microdroplet?

Abstract Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O–H bonds at the surface are ~16 MV/cm larger on average than that found for O–H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Joint Modeling/Experimental Approach to Characterize Ionization and Fragmentation of SOA Molecules with CIMS

The Molecular Dynamics technique Collective variable hyperdynamics (CVHD) interfaced with density functional tight-binding at the GFN1-xTB and GFN2-xTB levels of electronic structure theory has been applied to the H + GAG positively charged peptide system to explore the fragmentation of this simple tripeptide under relatively low temperature/energy conditions. Unlike conventional chemical dynamics simulations which can only be performed up to hundreds of picoseconds when coupled with semi-empirical Hamiltonians, CVHD is able to capture the long-time dynamics of multiple proton hopping, the formation of reversible ring structures, and ester rearrangement of the H + GAG system prior to fragmentation. The CVHD method applied to H + GAG is also able to uncover alternative fragmentation pathways not considered previously such as glycine fragmentation on the N-terminal side of the charged peptide and new cyclic cationic species. These findings in the H + GAG case have implications for other peptide systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and other new capabilities. The work described in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include implementing hash table matrix assembly for efficient sparsity pattern construction for contact in BISON, developing re-step testing infrastructure for ensuring the viability of overlapping domain coupling between SAM and Pronghorn, allowing unique preconditioners for single-input multi-system solves, supporting multi-system in MOOSE’s workhorse executioners, and many more smaller feature enhancements and bug fixes.

97 - MATHEMATICS AND COMPUTING↗

Nanoelectromechanical Control of Spin–Photon Interfaces in a Hybrid Quantum System on Chip

Color centers (CCs) in nanostructured diamond are promising for optically linked quantum technologies. Scaling to useful applications motivates architectures meeting the following criteria: C1 individual optical addressing of spin qubits; C2 frequency tuning of spin-dependent optical transitions; C3 coherent spin control; C4 active photon routing; C5 scalable manufacturability; and C6 low on-chip power dissipation for cryogenic operations. Here, we introduce an architecture that simultaneously achieves C1–C6. We realize piezoelectric strain control of diamond waveguide-coupled tin vacancy centers with ultralow power dissipation necessary. The DC response of our device allows emitter transition tuning by over 20 GHz, combined with low-power AC control. We show acoustic spin resonance of integrated tin vacancy spins and estimate single-phonon coupling rates over 1 kHz in the resolved sideband regime. Combined with high-speed optical routing, our work opens a path to scalable single-qubit control with optically mediated entangling gates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Nd, SbNd and Sb 3 Nd 4 and their interactions with the cladding alloy HT9

Lanthanide fission products, such as neodymium, formed during the irradiation of metallic fuels are known to cause deleterious effects from chemical interactions occurring at the fuel-cladding interface; a phenomenon known as fuel-cladding chemical interaction (FCCI). The use of fuel-based additives that bind with the lanthanide elements within the fuel meat, alleviating their interactions at the fuel-cladding interface, is one potential method proposed to mitigate the FCCI phenomenon and extend the burnup potential of such metallic fuel systems. In this study, antimony (Sb) is evaluated as one such additive, and neodymium (Nd) is used to represent the lanthanides. A Sb-Nd alloy is fabricated, which consists of two intermetallic phases, SbNd and Sb 3 Nd 4 . Isothermal diffusion couple experiments are carried out at 675 °C for 24 h between Nd and the Fe-12Cr based HT9 cladding alloy. The results are compared against similar diffusion couple experiments carried out between the Sb-Nd alloy and also HT9. Inter-diffusion between Nd/HT9 diffusion couple is characterized and the phases Fe 17 Nd 5 and Fe 17 Nd 2 were found to form whereas no interactions were observed in the Sb-Nd/HT9 diffusion couple. The lack of compound forming tendencies between Fe, the primary alloying constituent of HT9, and SbNd is elucidated through density functional theory (DFT) calculations on enthalpy of mixing, aligning well with the experimental observations. Lastly, the strong binding of Nd with Sb appears to be favorable, alleviating Nd interactions with HT9 constituent elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tardigrade, Version 1.1x

Tardigrade is a collection of functions, libraries, and stand-alone code which enables the simulation of materials using the higher-order micromorphic framework. The tools present are the interface to INL's MOOSE code which is used as the underlying FEA solver, several constitutive equations (linear elasticity, a pressure-sensitive elasto-plastic model, crystal plasticity, and additional models for geo-materials) useful for demonstration purposes and the simulation of some materials, the micromorphic filter which converts DNS information to micromorphic stress-deformation quantities, and the overlap coupling framework which enables multi-scale simulations of a coupled higher-order macroscale and a classically defined microscale.

Miller, Nathan↗

Stabilization of glucose-6-phosphate dehydrogenase oligomers enhances catalytic activity and stability of clinical variants

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a genetic trait that can cause hemolytic anemia. To date, over 150 nonsynonymous mutations have been identified in G6PD, with pathogenic mutations clustering near the dimer and/or tetramer interface and the allosteric NADP + -binding site. Recently, our lab identified a small molecule that activates G6PD variants by stabilizing the allosteric NADP + and dimer complex, suggesting therapeutics that target these regions may improve structural defects. Here, we elucidated the connection between allosteric NADP + binding, oligomerization, and pathogenicity to determine whether oligomer stabilization can be used as a therapeutic strategy for G6PD deficiency (G6PD def ). We first solved the crystal structure for G6PD K403Q , a mutant that mimics the physiological acetylation of wild-type G6PD in erythrocytes and demonstrated that loss of allosteric NADP + binding induces conformational changes in the dimer. These structural changes prevent tetramerization, are unique to Class I variants (the most severe form of G6PD def ), and cause the deactivation and destabilization of G6PD. We also introduced nonnative cysteines at the oligomer interfaces and found that the tetramer complex is more catalytically active and stable than the dimer. Furthermore, stabilizing the dimer and tetramer improved protein stability in clinical variants, regardless of clinical classification, with tetramerization also improving the activity of G6PD K403Q and Class I variants. These findings were validated using enzyme activity and thermostability assays, analytical size-exclusion chromatography (SEC), and SEC coupled with small-angle X-ray scattering (SEC-SAXS). Taken together, our findings suggest a potential therapeutic strategy for G6PD def and provide a foundation for future drug discovery efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin pinning effect to reconstructed oxyhydroxide layer on ferromagnetic oxides for enhanced water oxidation

Producing hydrogen by water electrolysis suffers from the kinetic barriers in the oxygen evolution reaction (OER) that limits the overall efficiency. With spin-dependent kinetics in OER, to manipulate the spin ordering of ferromagnetic OER catalysts (e.g., by magnetization) can reduce the kinetic barrier. However, most active OER catalysts are not ferromagnetic, which makes the spin manipulation challenging. In this work, we report a strategy with spin pinning effect to make the spins in paramagnetic oxyhydroxides more aligned for higher intrinsic OER activity. The spin pinning effect is established in oxideFM/oxyhydroxide interface which is realized by a controlled surface reconstruction of ferromagnetic oxides. Under spin pinning, simple magnetization further increases the spin alignment and thus the OER activity, which validates the spin effect in rate-limiting OER step. The spin polarization in OER highly relies on oxyl radicals (O∙) created by 1 st dehydrogenation to reduce the barrier for subsequent O-O coupling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A red-light-powered silicon nanowire biophotochemical diode for simultaneous CO 2 reduction and glycerol valorization

A bias-free photochemical diode, in which a p-type photocathode is connected to an n-type photoanode to harness light for driving photoelectrochemical reduction and oxidation pairs, serves as a platform for realizing light-driven fuel generation from CO 2 . However, the conventional design, in which cathodic CO 2 reduction is coupled with the anodic oxygen evolution reaction (OER), requires substantial energy input. Here we present a photochemical diode device that harnesses red light (740 nm) to simultaneously drive biophotocathodic CO 2 -to-multicarbon conversion and photoanodic glycerol oxidation as an alternative to the OER to overcome the above thermodynamic limitation. The device consists of an efficient CO 2 -fixing microorganism, Sporomusa ovata, interfaced with a silicon nanowire photocathode and a Pt–Au-loaded silicon nanowire photoanode. This photochemical diode operates bias-free under low-intensity (20 mW cm –2 ) red light irradiation with ~80% Faradaic efficiency for both the cathodic and anodic products. Finally, this work provides an alternative photosynthetic route to mitigate excessive CO 2 emissions and efficiently generate value-added chemicals from CO 2 and glycerol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HPC4Mfg with Carbon Inc.

Continuous Liquid Interface Production (CLIP) is poised to bring additive manufacturing to multiple American manufacturing sectors owing to its unique combination of rapid print speeds and material options that resemble injection molding thermoplastics. In spite of these benefits, CLIP is still a maturing process. To improve our understanding and control of the process, this project developed a multi-physics computational model that encompasses the coupled chemical-physical processes of photopolymerization and fluid flow to predict part outcomes. A physically predictive model can enable rapid optimization of CLIP and reduce the current cycle time and waste associated with optimization by four-fold. Ultimately, this effort will help open the U.S. $400B plastic manufacturing industry to the key benefits of additive manufacturing, namely mass customization, unlimited design space, and a cost- and energy-effective path to mainstream manufacturing. In the course of the project we produced a model which coupled fluid-flow, species transport, and photochemistry to model the Carbon CLIP process. As the development proceeded, we validated the code against theory and experiments. The resulting code was transferred to Carbon, and we expect it to be a useful part of their modeling capability. Future work includes improvements to the fluid solver’s robustness and performance, along with the addition of additional physics models.

42 ENGINEERING↗

Development of a Multi-physics Model to Optimize Continuous Liquid Interface Production (CLIP) for Additive Manufacturing (Carbon3D HPC4Mfg Final Report)

Continuous Liquid Interface Production (CLIP) is poised to bring additive manufacturing to multiple American manufacturing sectors owing to its unique combination of rapid print speeds and material options that resemble injection molding thermoplastics. In spite of these benefits, CLIP is still a maturing process. To improve our understanding and control of the process, this project developed a multi-physics computational model that encompasses the coupled chemical-physical processes of photopolymerization and fluid flow to predict part outcomes. A physically predictive model can enable rapid optimization of CLIP and reduce the current cycle time and waste associated with optimization by four-fold. Ultimately, this effort will help open the U.S. $400B plastic manufacturing industry to the key benefits of additive manufacturing, namely mass customization, unlimited design space, and a cost- and energy-effective path to mainstream manufacturing. In the course of the project we produced a model which coupled fluid-flow, species transport, and photochemistry to model the Carbon CLIP process. As the development proceeded, we validated the code against theory and experiments. The resulting code was transferred to Carbon, and we expect it to be a useful part of their modeling capability. Future work includes improvements to the fluid solver’s robustness and performance, along with the addition of additional physics models.

36 MATERIALS SCIENCE↗

An examination of the metal ion content in the active sites of human endonucleases CPSF73 and INTS11

Human cleavage and polyadenylation specificity factor (CPSF)73 (also known as CPSF3) is the endoribonuclease that catalyzes the cleavage reaction for the 3'-end processing of pre-mRNAs. The active site of CPSF73 is located at the interface between a metallo-β-lactamase domain and a β-CASP domain. Two metal ions are coordinated by conserved residues, five His and two Asp, in the active site, and they are critical for the nuclease reaction. The metal ions have long been thought to be zinc ions, but their exact identity has not been examined. Here we present evidence from inductively coupled plasma mass spectrometry and X-ray diffraction analyses that a mixture of metal ions, including Fe, Zn, and Mn, is present in the active site of CPSF73. The abundance of the various metal ions is different in samples prepared from different expression hosts. Zinc is present at less than 20% abundance in a sample expressed in insect cells, but the sample is active in cleaving a pre-mRNA substrate in a reconstituted canonical 3'-end processing machinery. Zinc is present at 75% abundance in a sample expressed in human cells, which has comparable endonuclease activity. We also observe a mixture of metal ions in the active site of the CPSF73 homolog INTS11, the endonuclease for Integrator. Taken together, our results provide further insights into the role of metal ions in the activity of CPSF73 and INTS11 for RNA 3'-end processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

User-oriented Improvements in the MOOSE Framework in Support of Multiphysics Simulation

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed or will be deployed in the coming months. The capabilities were implemented in the same order as they are covered in this report: chainable execution objects or executors, support for transfers between same-level applications in a coupling scheme, support for boundary/subdomain restricted transfers, support for transfers between applications with different coordinate or unit systems, support for MOOSE applications in the NEAMS Workbench, deployment of the MOOSE application of the Idaho National Laboratory (INL) high-performance computing (HPC) OnDemand platform, the addition of a triangular meshing library in libMesh, and increased support of face variables.

97 MATHEMATICS AND COMPUTING↗

Modeling Marangoni Convection using a Velocity-Time Relationship

The Waste Treatment and Immobilization Plant (WTP) at the U.S. Department of Energy’s (DOE’s) Handford site is designed to treat 56 million gallons of radioactive waste through vitrification. The melting vessels physical and chemical integrity is essential to maintain, as failure has a high safety and financial penalty. Monofrax® K-3 refractory, the material widely used for the melting vessel, undergoes significant corrosion over time, concentrated at the triple-point junction between the gaseous atmosphere, refractory, and glass melt, producing a characteristic neck-like profile that significantly limits refractory service life. Accurate prediction of the corrosion within the melter can be used to prevent this and extend the lifetime of the melters. A computational fluid dynamics (CFD) model was previously developed, and predicted subsurface refractory corrosion, but failed to predict the neck-line profile observed experimentally. To address this limitation, the glass-melt meniscus geometry at the glass-refractory interface was estimated from experimentally measured neck profiles and wettability experiments. A curve was fitted to the neck profile to identify the triple point location, and the meniscus shape was characterized by a contact angle of 1°. A velocity profile accounting for Marangoni convection-driven corrosion was derived from the estimated meniscus geometry. Surface tension measurements from LORPM14R1 glasses of varying composition were coupled with K-3 coupon dissolution rates to establish a surface tension–time relationship. This analytical velocity-time model was integrated into the existing CFD framework to improve prediction of the neck corrosion profile. Work was done in a computational setting with national collaboration, resulting in professional development through presentational and technical writing growth.

36 - MATERIALS SCIENCE↗