INTERFACE MECHANICS FOR EXPLICIT EULERIAN TO EXPLICIT LAGRANGIAN COUPLING USING THE MODIFIED IMMERSED FINITE ELEMENT METHOD
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A system includes a first optical sensor sensitive to both a parameter of interest, Parameter1, and at least one confounding parameter, Parameter2 and a second optical sensor sensitive only to the confounding parameter. Measurement circuitry measures M1 in response to light scattered by the first optical sensor, where M1=value of Parameter1+K*value of Parameter2. The measurement circuitry also measures M2 in response to light scattered by the second optical sensor, where M2=value of Parameter2. Compensation circuitry determines a compensation factor, K, for the confounding parameter based on measurements of M1 and M2 taken over multiple load/unload cycles or over one or more thermal cycles. The compensation factor is used to determine the parameter of interest.
A quantum dot structure having a split-gate geometry is provided. The quantum dot is configured for incorporation into a quantum dot array of a quantum processing unit. A gap between a reservoir accumulation gate and a quantum dot accumulation gate provides a tunnel barrier between an electric charge reservoir and a quantum dot well. An electrical potential applied to the gates defines a tunnel barrier height, width and charge tunneling rate between the well and the reservoir without relying on any barrier gate to control the charge tunneling rate.
Explore the source record for details and available documents.
Motionally averaged dipolar couplings are an important tool for understanding the complex dynamics of catalysts, polymers, and biomolecules. While there is a plethora of solid-state NMR pulse sequences available for their measurement, in can be difficult to gauge the methods’ strengths and weaknesses. In particular, there has not been a comprehensive comparison of their performance in natural abundance samples, where 1H homonuclear dipolar couplings are important and the use of large MAS rotors may be required for sensitivity reasons. In this work, we directly compared some of the more common methods for measuring C–H dipolar couplings in natural abundance samples using L-alanine (L-Ala) and the N-formyl-L-methionyl-L-leucyl-L-phenylalanine (fMLF) tripeptide as model systems. Here, we evaluated their performance in terms of accuracy, resolution, sensitivity, and ease of implementation. We found that, despite the presence of 1H homonuclear dipolar interactions, all methods, with the exception of REDOR, were able to yield the reasonable dipolar coupling strengths for both mobile and static moieties. Of these methods, PDLF provides the most convenient workflow and precision at the expense of low sensitivity. In low-sensitivity cases, MAS-PISEMA and DIPSHIFT appear to be the better options.
Here a benchmark computational study of K-edge core-ionization energies of third-row elements using relativistic delta-coupled-cluster (ΔCC) methods and a revised core-valence separation (CVS) scheme is reported. High-level relativistic (HLR) corrections beyond the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), including the contributions from two-electron picture-change effects, spin–orbit coupling, the Breit term, and quantum electrodynamics effects, have been taken into account and demonstrated to play an important role. Relativistic ΔCC calculations are shown to provide accurate results for core-ionization energies of third-row elements. The SFX2C-1e-CVS-ΔCC results augmented with HLR corrections show a maximum deviation of less than 0.5 eV with respect to experimental values.
CONSPECTUS: Electrochemical organic oxidation reactions are highly appealing because protons are often effective terminal electron acceptors, thereby avoiding undesirable stoichiometric oxidants. These reactions are often plagued by high overpotentials, however, that greatly limit their utility. Single-electron transfer (SET) from organic molecules generates high energy radical-cations. Formation of such intermediates often requires electrode potentials far above the thermodynamic potentials of the reaction and frequently causes decomposition and/or side reactions of ancillary functional groups. In this Account, we show how electrocatalytic electron-proton transfer mediators (EPTMs) address this challenge. EPTMs bypass the formation of radical-cation intermediates by supporting mechanisms that operate at electrode potentials much lower (=1 V) than analogous direct electrolysis reactions. The stable aminoxyl radical TEMPO (2,2,6,6-tetramethylpiperidine N-oxyl) is an effective mediator for electrochemical alcohol oxidation, and we have employed such processes for applications ranging from pharmaceutical synthesis to biomass conversion. A complementary electrochemical alcohol oxidation method employs a cooperative Cu/TEMPO mediator system that operates at 0.5 V lower electrode potential than the TEMPO-only mediated process. This difference, which arises from a different catalytic mechanism, rationalizes the broad functional group tolerance of Cu/TEMPO-based aerobic alcohol oxidation catalysts. Aminoxyl mediators address long-standing challenges in the "Shono oxidation", an important method for a-C–H oxidation of tertiary amides and carbamates. Shono oxidations are initiated by a high-potential SET step that limits their utility. Aminoxyl-mediated Shono-type oxidations have been developed that operate at much lower potentials and tolerate diverse functional groups. Analogous reactivity underlies a-C–H cyanation of secondary cyclic amines, a new method that enables efficient diversification of piperidine-based pharmaceutical building blocks and preparation of non-natural amino acids. Electrochemical oxidations of benzylic C–H bonds are commonly initiated by SET to generate radical cations, but such methods are again plagued by large overpotentials. Mediated electrolysis methods that promote hydrogen-atom-transfer (HAT) from benzylic C–H bonds to Fe-oxo species and phthalimide N-oxyl (PINO) support C–H oxygenation, iodination, and oxidative-coupling reactions. A complementary method merges photochemistry with electrochemistry to achieve amidation of C(sp3)–H bonds. This unique process operates at much lower overpotentials compatible with diverse functional groups. These results have broad implications for organic electrochemistry, highlighting the importance of "overpotential" considerations and the prospects for expanding synthetic utility by using mediators to bypass high-energy outer-sphere electron-transfer mechanisms. Principles demonstrated here for oxidation are equally relevant to electrochemical reductions.
The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.
Quantum transduction is the process of upconverting microwave quantum signals into optical signals to develop quantum networks through optical fibers outside the dilution refrigerator, enabling connections between quantum technologies on the quantum internet. In this project, upconversion is achieved by directing an optical laser and the microwave quantum signal into an electro-optic bulk crystal. The crystal is housed within a Superconducting Radio Frequency (SRF) cavity designed to maximize the overlap between the microwave field and the crystal volume. In the presence of microwaves, the refractive index of the crystal changes through the Pockels effect. This change in refractive index modifies the propagation of the optical field within the crystal, allowing the quantum information carried by the microwave field to be transferred to the optical field. This study focuses on coupling laser light from suspended waveguide chips into the whispering-gallery modes (WGMs) of the crystal to enable transduction. The coupling efficiency between the optical field in the suspended waveguide and the WGM depends on the position of the laser spot on the crystal. To address this challenge, a feedback-based algorithm is being developed to fine-tune the waveguide position so that the optical signal is coupled efficiently into the crystal. After passing through the crystal, the optical signal is detected by a photodiode connected to an oscilloscope. The algorithm evaluates the coupling quality and iteratively adjusts the waveguide position to maximize coupling efficiency. The expected outcome is an automated waveguide alignment method that improves optical coupling and enables more efficient microwave-to-optical quantum transduction.
The complete and incomplete fusion cross sections for 6 Li+ 209 Bi were measured using the in-beam γ-ray method around the Coulomb barrier. The cross sections of (deuteron captured) incomplete fusion (ICF) products were re-quantified experimentally for this reaction system. Here, the results reveal that the ICF cross section is equivalent to that of complete fusion (CF) above the Coulomb barrier and dominant near or below the Coulomb barrier. A theoretical calculation based on the continuum discretized coupled channel (CDCC) method was performed for the aforementioned CF and ICF cross sections; the result is consistent with the experiments. The universal fusion function (UFF) was also compared with the measured CF cross section for different barrier parameters, demonstrating that the CF suppression factor is significantly influenced by the choice of potential, which can reflect both dynamic and static effects of breakup on the fusion process.
In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.
This dataset represents results from an experimental study designed to compare how the chemical composition of organic matter changes across different burn conditions and feedstock materials. The dataset provides both solid and dissolved phase bulk concentration and organic matter characterization data from experimentally generated chars. Chars were created in a closed muffle furnace or on an open burn table from four different feedstock species representing vegetation commonly impacted by fire regimes across the Pacific Northwest, USA. This data can be used to compare how different burn conditions may influence resultant organic matter chemistry and help further our understanding of potential biogeochemical impacts on river corridors post-fire. This dataset is comprised of one data package readme, one data dictionary (dd), one file level metadata (flmd), fourteen burn table videos, burn table video metadata and three folders containing (A) data; (B) metadata and protocols; and (C) photos. The folder names and the file name of the data package readme include a version number which will be updated with future iterations of this data package. The data folder includes (1) solid carbon and solid nitrogen; (2) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) and total dissolved nitrogen (TN); (3) pH; (4) thermocouple time series temperature; (5) methods codes; (6) installation methods; (7) excitation emissions matrix (EEM) methods information; (8) a folder of excitation emissions matrix (EEM) fluorescence and absorbance spectra in dissolved organic matter and EEMs processing instructions; (9) solid state carbon-13 and solution state phosphorus nuclear magnetic resonance (13-C NMR and 31-P NMR) data and methods; (10) benzene polycarboxylic acid (BPCA) concentration and stable isotope data; (11) FTICR-MS methods; (12) Inductively coupled plasma (ICP) data for total calcium, magnesium, iron, aluminum, potassium, phosphorus, sodium, and sulfur along with sodium hydroxide-ethylenediaminetetraacetic acid (sodium hydroxide-EDTA) extractable calcium, magnesium, iron, aluminum, potassium, phosphorus, and sulfur; (13) a folder of phosphorus, carbon, and nitrogen X-ray absorption near edge structure (P-XANES, N-XANES, C-XANES) data for samples and standards; (14) P-XANES, N-XANES, C-XANES methods; (15) molybdate reactive phosphorus; and (16) folder of high resolution characterization of organic matter via 21 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) generated through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). The FTICR folder contains .txt data files and a subfolder containing instruments for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. The metadata and protocols folder includes (1) international geo-sample number (IGSN) mapping file (2) burn and laboratory metadata; (3) burn protocol; (4) laboratory protocol; (5) vegetation collection metadata; and (6) vegetation collection protocol. The folder contains photos of the solid chars. All files are .csv, .txt, .pdf, .jpg, .jpeg, .R, .ref, or .mp4. The data package was originally published October 2022 (v1). It was updated April 2023 (v2; new data files), September 2023 (v3; new and corrected data files), and September 2024 (v4; new and added/updated files). Metadata files were also updated to reflect these changes. See the change history section in the readme for more details.
Advection-dispersion equations (ADEs) are commonly used to describe transport phenomena in porous media. Even though mature discretization-based numerical methods for ADEs exist, some challenges still remain, especially when it comes to solving advection-dominated forward ADEs and diffusion-dominated backward ADEs. The latter problem usually arises in the source identification context and leads to numerically unstable grid-based solutions that require a form of regularization or should be treated as an inverse problem that is computationally more expensive because it requires solving the forward problem multiple times. In this study, we propose a discretization-free approach based on the physics-informed neural network (PINN) method for solving coupled ADE and Darcy flow equations with space-dependent hydraulic conductivity. In this approach, the hydraulic conductivity, hydraulic head, and concentration fields are approximated with deep neural networks (DNNs). We assume that the conductivity field is given by its values on a grid, and we use these values to train the conductivity DNN. The head and concentration DNNs are trained by minimizing the residuals of the flow equation and ADE and using the initial and boundary conditions as additional constraints. The PINN method is applied to one- and two-dimensional forward ADE problems, where its performance for various P\'{e}clet numbers ($Pe$) is compared with the analytical and numerical solutions. We find that the PINN method is accurate with errors of less than 1\% and outperforms some conventional discretization-based methods for $Pe$ larger that 100. Next, we demonstrate that the PINN method remains accurate for the backward ADEs, with the relative errors in most cases staying under 5\% compared to the reference concentration field. Finally, we show that when available, the concentration measurements can be easily incorporated in the PINN method and significantly improve (by more than 50\% in the considered cases) the accuracy of the PINN solution of the backward ADE.
We present a sixth order finite difference scheme for Poisson’s equation when discretized with the compact 27-point stencil based on Mehrstellen corrections of the forcing function term f. Our approach results in a sixth order accurate solution error as opposed to a fourth-order error imposed by the classical Mehrstellen correction for the 19-point and 27-point stencils. The present study is a continuation of former work of Spotz and Carey (1996) on compact finite difference schemes for Poisson’s equation where sixth order convergence may be obtained under the assumption that the fourth order derivatives of f are determined analytically. Specifically, we show that sixth order convergence can still be attained when only values of f at grid points are available. The sixth order Mehrstellen scheme is further coupled with a Method of Local Corrections (MLC) 3D Poisson solver improving to sixth order accuracy the results reported in Kavouklis and Colella (2019). The MLC test case considered involves an adaptive grid that comprises 7.5 billion cells.
In principle, a patterned dielectric back contact structure in a GaAs thermophotovoltaic device boosts the sub-bandgap reflectance. Rigorous coupled wave analysis method is used to study the three dimensional (3D) periodic grating effects of metal point contact diameter and spacing on sub-bandgap reflectance of the device.
Electronic structure methods typically benefit from symmetry breaking and restoration, specially in the strong correlation regime. Here, the same goes for ansätze on a quantum computer. We develop a unitary coupled cluster method based on the antisymmetrized geminal power (AGP)—a state formally equivalent to the number-projected Bardeen–Cooper–Schrieffer wavefunction. We demonstrate our method for the single-band Fermi–Hubbard Hamiltonian in one and two dimensions. We also explore post-selection as a state preparation step to obtain correlated AGP and prove that it scales no worse than O(√M) in the number of measurements, thereby making it a less expensive alternative to gauge integration to restore particle number symmetry.
Multi-scale coupling has increasingly drawn attention as a promising approach for modeling thermal systems. Traditional system codes provide fast and robust predictions at the plant scale, while high-fidelity computational fluid dynamics (CFD)-based tools resolve localized flow and heat transfer phenomena with greater accuracy. By combining these complementary strengths, co-simulations enable multi-scale analysis that would otherwise be computationally prohibitive for a standalone CFD code. Here, this work introduces a robust and problem-agnostic domain overlapping (DO) coupling between the system thermal-hydraulic (STH) code System Analysis Module (SAM) and the coarse-mesh CFD code Pronghorn. Both applications belong to the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite, a code suite in active development at the Idaho National Laboratory (INL), tailored for multi-physics analysis of advanced reactors. Unlike previous approaches, BlueCRAB supports an agnostic interface between codes based on different fidelity, while its coupling formulation can address arbitrary flow geometries with multiple inlets and outlets in coupled components. The implemented method leads to consistent pressure drops, enthalpies, and scalar concentrations between coupled SAM and Pronghorn simulations. The methodology is demonstrated through two verification tests, which ensure the numerical consistency and conservation across the codes, and through one validation test against experimental data. The proposed problems explore different physical aspects inherent to thermal systems, with particular attention given to nuclear reactor analysis. These include buoyancy-driven flows, complex flow patterns, and setups with multiple inlets and outlets, representing challenges in advanced reactor applications.