Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 397 records · Page 22

Classification of (2 + 1) D invertible fermionic topological phases with symmetry

Here we provide a classification of invertible topological phases of interacting fermions with symmetry in two spatial dimensions for general fermionic symmetry groups G f and general values of the chiral central charge c - . Here G f is a central extension of a bosonic symmetry group G b by fermion parity, (-1) F , specified by a second cohomology class [ω 2 ]∈ $\mathscr{H}^2$(Gb,$\mathbb{Z}_2$). Our approach proceeds by gauging fermion parity and classifying the resulting G b symmetry-enriched topological orders while keeping track of certain additional data and constraints. We perform this analysis through two perspectives, using G-crossed braided tensor categories and Spin(2c - ) 1 Chern-Simons theory coupled to a background G gauge field. These results give a way to characterize and classify invertible fermionic topological phases in terms of a concrete set of data and consistency equations, which is more physically transparent and computationally simpler than the more abstract methods using cobordism theory and spectral sequences. Our results also generalize and provide a different approach to the recent classification of fermionic symmetry-protected topological phases by Wang and Gu, which have chiral central charge c - = 0. We show how the tenfold way classification of topological insulators and superconductors fits into our scheme, along with general nonperturbative constraints due to certain choices of c - and G f . Mathematically, our results also suggest an explicit general parametrization of deformation classes of (2 + 1)D invertible topological quantum field theories with G f symmetry.

36 MATERIALS SCIENCE↗

Laser based 100 GeV electron acceleration scheme for muon production

High energy muons, due to their unique ability to penetrate deeply into matter, can enable radiography of structures that cannot be probed by other forms of radiation. Current terrestrial sources of muons require conventional GeV-TeV particle accelerators which are hundreds to thousands of meters in size. Laser wakefield acceleration (LWFA) can achieve acceleration gradients of two-to-three orders of magnitude greater than conventional accelerators, thus shrinking the accelerator to a number of meters. We propose a concept for a compact muon source based on the first self-consistent PIC simulations of an all optical LWFA that uses a guiding channel to achieve electron energies of 100 GeV in a distance of 6 m with a driving laser energy of 300 J in a single stage. From the resulting electron energy spectrum we estimate muon production for this source. We show that this accelerator, coupled with high average power laser driver technology, provides the basis for a high energy and high flux muon source.

laser-produced plasmas↗

MONACO: accurate biological network alignment through optimal neighborhood matching between focal nodes

Motivation: Alignment of protein–protein interaction networks can be used for the unsupervised prediction of functional modules, such as protein complexes and signaling pathways, that are conserved across different species. To date, various algorithms have been proposed for biological network alignment, many of which attempt to incorporate topological similarity between the networks into the alignment process with the goal of constructing accurate and biologically meaningful alignments. Especially, random walk models have been shown to be effective for quantifying the global topological relatedness between nodes that belong to different networks by diffusing node-level similarity along the interaction edges. However, these schemes are not ideal for capturing the local topological similarity between nodes. Results: Here, we propose MONACO, a novel and versatile network alignment algorithm that finds highly accurate pairwise and multiple network alignments through the iterative optimal matching of ‘local’ neighborhoods around focal nodes. Extensive performance assessment based on real networks as well as synthetic networks, for which the ground truth is known, demonstrates that MONACO clearly and consistently outperforms all other state-of-the-art network alignment algorithms that we have tested, in terms of accuracy, coherence and topological quality of the aligned network regions. Furthermore, despite the sharply enhanced alignment accuracy, MONACO remains computationally efficient and it scales well with increasing size and number of networks.

97 MATHEMATICS AND COMPUTING↗

𝛽−decay of neutron-rich 45 Cl located at the magic number 𝑁 = 28

Results from the β-decay study of 45 Cl, produced in the fragmentation of a 140-MeV/u 48 Ca beam, are presented. The half-life of 45 Cl is measured to be 513(36) ms. β – and β – 1n decay populate excited states in 45 Ar and 44 Ar, respectively, and decay schemes for both daughter nuclei are established based on γ-ray singles and γ-γ coincidence data. Comparison with shell-model calculations using the FSU interaction shows that the low-lying negative-parity states in 45 Ar are well described by single-neutron configurations near the Fermi surface, whereas positive-parity states require neutron excitations across the N = 20 shell gap, consistent with allowed Gamow–Teller decay from 45 Cl. The strong β-feeding to the 5/2 + state in 45 Ar supports a 3/2 + ground-state assignment for 45 Cl over 1/2 + . Furthermore, the high Q β – value enables population of 1p1h states above the neutron separation energy in 45 Ar, leading to positive-parity states in 44 Ar through neutron emission. Tentative spin-parity assignments for excited levels in 44 Ar are made for the first time based on comparison with shell-model predictions, and the 2978-keV level is identified as an excited 0 + state with a configuration distinct from the ground state.

39 ≤ A ≤ 58↗

CAS without SCF—Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree–Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Lastly, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semi-Automated, Object-Based Tomography of Dislocation Structures

The characterization of the three-dimensional arrangement of dislocations is important for many analyses in materials science. Dislocation tomography in transmission electron microscopy is conventionally accomplished through intensity-based reconstruction algorithms. Although such methods work successfully, a disadvantage is that they require many images to be collected over a large tilt range. Here, we present an alternative, semi-automated object-based approach that reduces the data collection requirements by drawing on the prior knowledge that dislocations are line objects. Our approach consists of three steps: (1) initial extraction of dislocation line objects from the individual frames, (2) alignment and matching of these objects across the frames in the tilt series, and (3) tomographic reconstruction to determine the full three-dimensional configuration of the dislocations. Drawing on innovations in graph theory, we employ a node-line segment representation for the dislocation lines and a novel arc-length mapping scheme to relate the dislocations to each other across the images in the tilt series. We demonstrate the method for a dataset collected from a dislocation network imaged by diffraction-contrast scanning transmission electron microscopy. Based on these results and a detailed uncertainty analysis for the algorithm, we discuss opportunities for optimizing data collection and further automating the method.

47 OTHER INSTRUMENTATION↗

Transverse momentum-dependent heavy-quark fragmentation at next-to-leading order

The transverse momentum-dependent fragmentation functions (TMD FFs) of heavy (bottom and charm) quarks, which we recently introduced, are universal building blocks that enter predictions for a large number of observables involving final-state heavy quarks or hadrons. They enable the extension of fixed-order subtraction schemes to quasi-collinear limits, and are of particular interest in their own right as probes of the nonperturbative dynamics of hadronization. In this paper we calculate all TMD FFs involving heavy quarks and the associated TMD matrix element in heavy-quark effective theory (HQET) to next-to-leading order in the strong interaction. Our results confirm the renormalization properties, large-mass, and small-mass consistency relations predicted in our earlier work. We also derive and confirm a prediction for the large-z behavior of the heavy-quark TMD FF by extending, for the first time, the formalism of joint resummation to capture quark mass effects in heavy-quark fragmentation. Our final results in position space agree with those of a recent calculation by another group that used a highly orthogonal organization of singularities in the intermediate momentum-space steps, providing a strong independent cross check. As an immediate application, we present the complete quark mass dependence of the energy-energy correlator (EEC) in the back-to-back limit at $\mathcal{O}\left({\alpha}_s\right)$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗

Coupling magnons to an opto-electronic parametric oscillator

Hybrid magnonic systems have emerged as versatile modular components for quantum signal transduction and sensing applications owing to their capability of connecting distinct quantum platforms. To date, the majority of the magnonic systems have been explored in a local, near-field scheme, due to the close proximity required for realizing a strong coupling between magnons and other excitations. This constraint greatly limits the applicability of magnons in developing remotely coupled, distributed quantum network systems. On the contrary, opto-electronic architectures hosting self-sustained oscillations have been a unique platform for long-haul signal transmission and processing. Here, we integrated an opto-electronic oscillator with a magnonic oscillator consisting of a microwave waveguide and a Y 3 Fe 5 O 12 (YIG) sphere, and demonstrated strong and coherent coupling between YIG’s magnon modes and the opto-electronic oscillator’s characteristic photon modes—revealing the hallmark anti-crossing gap in the measured spectrum. In particular, the photon mode is produced on-demand via a nonlinear, parametric process as stipulated by an external seed pump. Both the internal cavity phase and the external pump phase can be precisely tuned to stabilize either degenerate or nondegenerate auto-oscillations. Our result lays out a new, hybrid platform for investigating the long-distance coupling and nonlinearity in coherent magnonic phenomena, which may be found useful in constructing the future “distributed hybrid magnonic systems.”

36 MATERIALS SCIENCE↗

Two-dimensional shaping of Solov'ev equilibrium with vacuum using external coils

In this work, we demonstrate a method for constructing the Solov'ev equilibrium with any given 2D shape surrounded by a vacuum region using external poloidal field coils, based on the method previously developed by T. Xu and R. Fitzpatrick [Nuclear Fusion, 59, 064002 (2019)]. The computational domain consists of two parts: the plasma region, where the solution is the same as the Solov'ev solution, and the vacuum region, where the magnetic field generated by external coils as well as plasma current is determined using the Green function method through a matching condition near the separatrix. However, the method is not limited to the Solov'ev equilibrium in particular. In conclusion, the accuracy, efficiency, and robustness of such a scheme suggest that this method may be applied to the 2D shaping of a tokamak plasma with a vacuum region using external coils in general.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development and Performance Evaluation of of Active Insulation Systems using Solid-State Thermal Switches

Traditional building envelopes have passive insulation systems that cannot respond to dynamic changes in the environment. An Active Insulation System (AIS) consists of Active Insulation Materials (AIMs), which dynamically vary the thermal conductivity of the insulation system. Several researchers have evaluated the impact of AIS on building thermal and energy performance by using simulation tools. Up to 70% savings in annual heating and cooling energy and significant reductions in peak demand have been predicted for some climates with wall systems employing AIS. Sensitivity studies have provided requirements for the AIS properties and control schemes that can lead to those savings. However, materials and assembly development have not yet achieved a cost-effective product that achieves the required performance. In this study, we present the process to develop an AIS that we will install in a test hut for their performance evaluation. Minimum performance criteria of the AIS system are developed based on Rmin/Rmax ratio, required time and efficiency to switch states and cost estimate. The next step carried out during this study is creating the concept to meet the requirements, predicting the performance by simulations, developing the experimental setup for bench-scale testing, and finally, constructing a full-scale wall assembly and monitoring the performance when exposed to natural weather conditions. The selected approach uses off-the-shelf products to create a materials system that can switch R-value between R~=1 ft2·°F·h/BTU (0.18 m2·K/W) and R~=7 ft2·°F·h/BTU (1.23 m2·K/W and have a switching time of less than one minute between R-high and R-low.

Iffa, Emishaw↗

Fourteen Years of Energetic Neutral Atom Observations from IBEX

The Interstellar Boundary Explorer (IBEX) has been observing the outer heliosphere and its interactions with the very local interstellar medium (VLISM) via measurements of energetic neutral atoms (ENAs) for over 14 yr. We discovered the IBEX Ribbon—a structure completely unanticipated by any prior theory or model—that almost certainly resides beyond the heliopause in the VLISM. We also characterized the other major source of heliospheric ENAs, the globally distributed flux (GDF), produced largely in the heliosheath between the termination shock and heliopause. In this study, we make three major new contributions. First, we validate, provide, and analyze the most recent 3 yr of IBEX-Hi (0.5–6 keV FWHM) data (2020–2022) for the first time. Second, we link these observations to the prior 11 yr of observations, exploring long-term variations. Finally, we provide the first IBEX team-validated Ribbon/GDF separation scheme and separated maps. Because of the uncertainty in separating different line-of-sight integrated sources, we provide not just best guess (median) maps, but also maps with upper and lower reasonable values of Ribbon and GDF fluxes, along with bounding fluxes that add the uncertainties to the upper and lower values. This allows theories and models to be compared with a range of possible values that the IBEX team believes are consistent with data. These observations, along with the reanalysis of the prior 11 yr of IBEX-Hi data, provide new insights and even further develop our detailed understanding of the heliosphere’s interaction with the local interstellar medium unlocked by IBEX.

79 ASTRONOMY AND ASTROPHYSICS↗

PAMTRA 1.0: the Passive and Active Microwave radiative TRAnsfer tool for simulating radiometer and radar measurements of the cloudy atmosphere

Forward models are a key tool to generate synthetic observations given knowledge of the atmospheric state. In this way, they are an integral part of inversion algorithms that aim to retrieve geophysical variables from observations or in data assimilation. Their application for the exploitation of the full information content of remote sensing observations becomes increasingly important when these are used to evaluate the performance of cloud-resolving models (CRMs). Herein, CRM profiles or fields provide the input to the forward model whose simulation results are subsequently compared to the observations. This paper introduces the freely available comprehensive microwave forward model PAMTRA (Passive and Active Microwave TRAnsfer), demonstrates its capabilities to simulate passive and active measurements across the microwave spectral region for upward- and downward-looking geometries, and illustrates how the forward simulations can be used to evaluate CRMs and to interpret measurements to improve our understanding of cloud processes. PAMTRA is unique as it treats passive and active radiative transfer (RT) in a consistent way with the passive forward model providing upwelling and downwelling polarized brightness temperatures and radiances for arbitrary observation angles. The active part is capable of simulating the full radar Doppler spectrum and its moments. PAMTRA is designed to be flexible with respect to instrument specifications and interfaces to many different formats of input and output, especially CRMs, spanning the range from bin-resolved microphysical output to one- and two-moment schemes, and to in situ measured hydrometeor properties. A specific highlight is the incorporation of the self-similar Rayleigh–Gans approximation (SSRGA) for both active and passive applications, which becomes especially important for the investigation of frozen hydrometeors.

54 ENVIRONMENTAL SCIENCES↗

Graph-Based Representations and Applications to Process Simulation

Rapid and robust convergence of a process flowsheet is critical to enable large-scale simulations that address core scientific questions related to process design, optimization, and sustainability. However, due to the highly coupled and nonlinear nature of chemical processes, efficiently solving a flowsheet remains a challenge. In this work, we show that graph representations of the underlying physical phenomena in unit operations may help identify potential avenues to systematically reformulate the network of equations and enable more robust topology-based convergence of flowsheets. To this end, we developed graph abstractions of the governing equations of vapor-liquid and liquid-liquid equilibrium separation equipment. These graph abstractions consist of a mesh of interconnected variable nodes and equation nodes that are systematically generated through PhenomeNode, a new open-source library in Python developed in this study. We show that partitioning the graph into separate mass, energy, and equilibrium subgraphs can help decouple nonlinearities and guide decomposition algorithms. By employing the graph abstraction on an industrial separation process for separating glacial acetic acid from water, we implemented a new block decomposition scheme in BioSTEAM and demonstrated that this can accelerate convergence over a traditional sequential modular approach.

Distillation↗

Anomalous softening of 3D printed elastomeric foam irradiated under compressive strain

Elastomeric foam is an essential component in many industrial and technological settings, primarily as thermal insulators and as positional/mechanical support cushions. In particular, silicone foam is utilized in harsh environments due to exceptional thermal and chemical stability. Under service conditions within certain applications such material gets exposed to a high dosage of gamma radiation, which can permanently alter the material’s structural and mechanical response properties. Most studies on gamma-exposure under inert or oxidative atmosphere indicate hardening of silicone foam, which is attributed to radiation-induced enhancement in chemical cross-linking. Here we report two contrasting effects depending on whether (non-oxidative) radiation exposure is carried out with the foam under zero or finite compressive strain. While in the former case we observe radiation-hardening consistent with previous studies, in the latter case (50% porous foam under 30% uniaxial compression) we see a monotonic decrease in Young’s modulus with increasing dosage, although solvent swelling experiments on the constituent rubber indicate a net increase in cross-link density independent of the state of strain. We quantitatively model all dose-dependent data using the Ogden Hyperfoam strain-energy function within the framework of Tobolsky two-network scheme and attribute the above anomaly to a combined effect of radiation-induced thickness change (compression set) and inherent nonlinearity in the foam’s stress-strain response.

Coarse-grained models↗

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗

Recent advances in the quantification of uncertainties in reaction theory

Uncertainty quantification has become increasingly more prominent in nuclear physics over the past several years. In few-body reaction theory, there are four main sources that contribute to the uncertainties in the calculated observables: the effective potentials, approximations made to the few-body problem, structure functions, and degrees of freedom left out of the model space. In this work, we illustrate some of the features that can be obtained when modern statistical tools are applied in the context of nuclear reactions. This work consists of a summary of the progress that has been made in quantifying theoretical uncertainties in this domain, focusing primarily on those uncertainties coming from the effective optical potential as well as their propagation within various reaction theories. We use, as the central example, reactions on the doubly-magic stable nucleus 40 Ca, namely neutron and proton elastic scattering and single-nucleon transfer 40 Ca(d,p) 41 Ca. First, we show different optimization schemes used to constrain the optical potential from differential cross sections and other experimental constraints; we then discuss how these uncertainties propagate to the transfer cross section, comparing two reaction theories. Finally, we finish by laying out our future plans.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗