The Schwarz Alternating Method for ROM-FOM and ROM-ROM Coupling.
Abstract not provided.
Tezaur, Irina↗
SEARCH · Engineering Papers
Results for “coupling method”
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.
Turbine airfoil with leading edge cooling passage(s) coupled via plenum to film cooling holes, and related method
A turbine airfoil includes a body including a wall defining pressure and suction sides, and a leading edge extending between the pressure and suction sides. A cooling circuit inside the wall of the body includes at least one of: a) a suction side to pressure side cooling sub-circuit including a first cooling passage(s) extending from the suction side to the pressure side around the leading edge to a first plenum, and a plurality of first film cooling holes communicating with the first plenum and extending through the wall on the pressure side; and b) a pressure side to suction side cooling sub-circuit including second cooling passage(s) extending from the pressure side to the suction side around the leading edge to a second plenum, and a plurality of second film cooling holes communicating with the second plenum and extending through the wall on the suction side.
Lacy, Benjamin Paul↗
Advances in an r-adaptive high-order method for high-speed flows: multi-domain coupling and uncertainty quantification
Explore the source record for details and available documents.
Van Heyningen, Robert Loek [Sandia National Labora↗
Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective
The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.
42 ENGINEERING↗
Model Reduction by Generalized Falk Method for Efficient Field-Circuit Simulations
The Generalized Falk Method (GFM) for coordinate transformation, together with two model-reduction strategies based on this method, are presented for efficient coupled field-circuit simulations. Each model-reduction strategy is based on a decision to retain specific linearly-independent vectors, called trial vectors, to construct a vector basis for coordinate transformation. The reduced-order models are guaranteed to be stable and passive since the GFM is a congruence transformation of originally symmetric positive definite systems. We also show that, unlike the Pad´e-via-Lanczos (PVL) method, the GFM does not generate unstable positive poles while reducing the order of circuit problems. Further, the proposed GFM is also faster when compared to methods of the type Lanczos (or Krylov) that are already widely used in circuit simulations for electrothermal and electromagnetic problems. The concept of response participation factors is introduced for the selection of the trial vectors in the proposed model-reduction methods. Further, we present methods to develop simple equivalent circuit networks for the field component of the overall field-circuit system. The implementation of these equivalent circuit networks in circuit simulators is discussed. With the proposed model-reduction strategies, significant improvement on the efficiency of the generalized Falk method is illustrated for coupled field-circuit problems.
42 ENGINEERING↗
Interface learning in fluid dynamics: Statistical inference of closures within micro–macro-coupling models
Many complex multiphysics systems in fluid dynamics involve using solvers with varied levels of approximations in different regions of the computational domain to resolve multiple spatiotemporal scales present in the flow. The accuracy of the solution is governed by how the information is exchanged between these solvers at the interface and several methods have been devised for such coupling problems. In this article, we construct a data-driven model by spatially coupling a microscale lattice Boltzmann method (LBM) solver and macroscale finite difference method (FDM) solver for reaction-diffusion systems. The coupling between the micro-macro solvers has one to many mapping at the interface leading to the interface closure problem, and we propose a statistical inference method based on neural networks to learn this closure relation. The performance of the proposed framework in a bifidelity setting partitioned between the FDM and LBM domain shows its promise for complex systems where analytical relations between micro-macro solvers are not available.
42 ENGINEERING↗
Saddle point search with dynamic active volume
Sampling potential energy surface (PES) is critical for many problems in materials science, chemistry, physics, and biology and requires highly efficient saddle point searches (SPS). In the study, we introduce the concept of dynamic active volume (DAV) in addition to the active volume in self-evolving atomistic kinetic Monte Carlo (SEAKMC). The DAV method has further reduced the dimensionality of the PES at the elevation stage of a SPS. At the subsequent converging stage, the dynamic boundary in DAV is lifted to allow the system to converge to the right location of a saddle point. Coupled with the dimer method, the DAV method not only significantly reduces the time cost for a given search attempt, but also dramatically increases the probability of finding relevant saddle points for PES sampling. Further, a Python software package within the framework SEAKMC (SEAKMC_py) with the DAV method has been developed.
36 MATERIALS SCIENCE↗
Coherent electrical control of a single high-spin nucleus in silicon
Nuclear spins are highly coherent quantum objects. In large ensembles, their control and detection via magnetic resonance is widely exploited, for example, in chemistry, medicine, materials science and mining. Nuclear spins also featured in early proposals for solid-state quantum computers and demonstrations of quantum search and factoring algorithms. Scaling up such concepts requires controlling individual nuclei, which can be detected when coupled to an electron. However, the need to address the nuclei via oscillating magnetic fields complicates their integration in multi-spin nanoscale devices, because the field cannot be localized or screened. Control via electric fields would resolve this problem, but previous methods relied on transducing electric signals into magnetic fields via the electron–nuclear hyperfine interaction, which severely affects nuclear coherence. In this work, we demonstrate the coherent quantum control of a single 123 Sb (spin-7/2) nucleus using localized electric fields produced within a silicon nanoelectronic device. The method exploits an idea proposed in 1961 but not previously realized experimentally with a single nucleus. Our results are quantitatively supported by a microscopic theoretical model that reveals how the purely electrical modulation of the nuclear electric quadrupole interaction results in coherent nuclear spin transitions that are uniquely addressable owing to lattice strain. The spin dephasing time, 0.1 seconds, is orders of magnitude longer than those obtained by methods that require a coupled electron spin to achieve electrical driving. These results show that high-spin quadrupolar nuclei could be deployed as chaotic models, strain sensors and hybrid spin-mechanical quantum systems using all-electrical controls. Integrating electrically controllable nuclei with quantum dots could pave the way to scalable, nuclear- and electron-spin-based quantum computers in silicon that operate without the need for oscillating magnetic fields.
42 ENGINEERING↗
Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems
Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.
Basis sets↗
Cascade of vestigial orders in two-component superconductors: Nematic, ferromagnetic, s -wave charge- 4 e , and d -wave charge- 4 e states
Electronically ordered states that break multiple symmetries can melt in multiple stages, similarly to liquid crystals. In a partially melted phase, known as vestigial phase, a bilinear made out of combinations of the multiple components of the primary order parameter condenses. Multicomponent superconductors are thus natural candidates for vestigial order since they break both the U(1)-gauge and also time-reversal or lattice symmetries. Here, we use group theory to classify all possible real-valued and complex-valued bilinears of a generic two-component superconductor on a tetragonal or hexagonal lattice. While the more widely investigated real-valued bilinears correspond to vestigial nematic or ferromagnetic order, the little explored complex-valued bilinears correspond to a vestigial charge-4e condensate, which itself can have an underlying s-wave, d x 2 –y 2 -wave, or d xy -wave symmetry. To properly describe the fluctuating regime of the superconducting Ginzburg-Landau action and thus access these competing vestigial phases, we employ both a large-N and a variational method. We show that while vestigial order can be understood as a weak-coupling effect in the large-N approach, it is akin to a moderate-coupling effect in the variational method. Despite these distinctions, both methods yield similar results in wide regions of the parameter space spanned by the quartic Landau coefficients. Specifically, we find that the nematic and ferromagnetic phases are the leading vestigial instabilities, whereas the various types of charge-4e order are attractive albeit subleading vestigial channels. The only exception is for the hexagonal case, in which the nematic and s-wave charge-4e vestigial states are degenerate. Finally, we discuss the limitations of our approach, as well as the implications of our results for the realization of exotic charge-4e states in material candidates.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework
We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.
Computer science↗
Bootstrap embedding for interacting electrons in phonon coherent-state mean field
Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.
Islam, Shariful [North Carolina State University, ↗
Complex paths around the sign problem
The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of physics, from QCD and nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half filling, and dynamical and nonequilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called sign problem. Here, a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it is reviewed. The mathematical ideas underpinning this approach, as well as the algorithms developed thus far, are described together with nontrivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks
Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation timeframes and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford.
Lines, Amanda M.↗
Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks - 20380
Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation time frames and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Simulation scaling studies of reactor core two-phase flow using direct numerical simulation
Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Hexagonal Geometries in MPACT
The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Islet: interpolation semi-Lagrangian element-based transport
Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.
58 GEOSCIENCES↗