Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “warping”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Restoration of weak localization in bilayer graphene by a molecular thin film

Quantum coherent effects can be probed in multilayer graphene through electronic transport measurements at low temperatures. In particular, bilayer graphene (BLG) is known to be susceptible to quantum interference corrections of the conductivity, presenting weak localization at all electronic densities, and dependent on different scattering mechanisms such as those related to the trigonal warping of the electron dispersion near the K and K′ valleys. Proximity effects with a molecular thin film influence these scattering mechanisms, which can be quantified through the known theory of magnetoconductance for BLG. Here, we present electronic transport measurements in a copper-phthalocyanine (CuPc) / BLG / hexagonal boron nitride (h-BN) heterostructure that suggest the restoration of weak localization in BLG, associated to a reduction of trigonal warping effects, that are known to suppress weak localization in BLG. Additionally, we observe a charge transfer of 3.6×10 12 cm −2 from the BLG to the molecules, as well as a very small degradation of the mobility of the BLG/h-BN heterostructure upon the deposition of CuPc. The molecular arrangement of the CuPc thin film is characterized in a control sample through transmission electron microscopy, that we relate to the electronic transport results.

bilayer graphene↗
DOE OSTI2020

Topological band structure transitions and goniopolar transport in honeycomb antimonene as a function of buckling

The electronic band topology of monolayer β-Sb (antimonene) is studied from the flat honeycomb to the equilibrium buckled structure using first-principles calculations and analyzed using a tight-binding model and low-energy Hamiltonians. In flat monolayer Sb, the Fermi level occurs near the intersection of two warped Dirac cones, one associated with the pz orbitals, and one with the {p x , p y } orbitals. The differently oriented threefold warping of these two cones leads to an unusually shaped nodal line, which leads to anisotropic in-plane transport properties and goniopolarity. A slight buckling opens a gap along the nodal line except at six remaining Dirac points, protected by symmetry. Under increasing buckling, pairs of Dirac points of opposite winding number annihilate at a critical buckling angle. At a second critical angle, the remaining Dirac points disappear when the band structure opens a gap. Spin-orbit coupling and edge states are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonperturbative resolution of strong coupling singularities in 4D 𝒩 = 1 heterotic M-theory

We investigate the interior of the moduli space of four-dimensional 𝒩 = 1 theories of gravity arising from compactifications of the 𝐸 8 × 𝐸 8 heterotic string on Calabi-Yau threefolds. By studying the threshold corrections to the coupling of the heterotic gauge groups, we infer the existence of a strong coupling singularity for one of the perturbative heterotic gauge groups, which effectively yields an additional finite distance boundary of the classical scalar field space. In heterotic M-theory, this boundary maps to a domain wall solution for which the gauge coupling and the warp factor on one of the Horava-Witten 9-branes diverge, thus highlighting the gravitational origin of the classical strong coupling singularity. The divergence of the warp factor is, however, regulated once nonperturbative effects are taken into account, as we demonstrate by studying the instanton corrections to the 5D BPS domain wall equations. This regularization implies that the classical strong coupling boundary of the scalar field 4D 𝒩 = 1 heterotic M-theory is resolved, indicating that, at the quantum level, the field space can be extended beyond this classical boundary.

Cvetič, Mirjam↗
DOE OSTI2020

A High-Throughput Solver for Marginalized Graph Kernels on GPU

Here, we present the design and optimization of a solver for efficient and high-throughput computation of the marginalized graph kernel on General Purpose GPUs. The graph kernel is computed using the conjugate gradient method to solve a generalized Laplacian of the tensor product between a pair of graphs. To cope with the large gap between the instruction throughput and the memory bandwidth of the GPUs, our solver forms the graph tensor product on-the-fly without storing it in memory. This is achieved by using threads in a warp cooperatively to stream the adjacency and edge label matrices of individual graphs by small square matrix blocks called tiles, which are then staged in registers and the shared memory for later reuse. Warps across a thread block can further share tiles via the shared memory to increase data reuse. We exploit the sparsity of the graphs hierarchically by storing only non-empty tiles using a coordinate format and nonzero elements within each tile using bitmaps. We propose a new partition-based reordering algorithm for aggregating nonzero elements of the graphs into fewer but denser tiles to further exploit sparsity. We carry out extensive theoretical analyses on the graph tensor product primitives for tiles of various density and evaluate their performance on synthetic and real-world datasets. Our solver delivers three to four orders of magnitude speedup over existing CPU-based solvers such as GraKeL and GraphKernels. The capability of the solver enables kernel-based learning tasks at unprecedented scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Smaller feature size impairs densification in digital light processing 3D printed ceramic samples

In this study, we investigate how geometric miniaturization affects the sintering densification behavior of alumina cylinders with various wall thicknesses (down to ∼500 µm in green body) fabricated via digital light processing (DLP). As the wall thickness decreased, the density of the sintered samples consistently declined. Microstructural characterization revealed the presence of a crack-rich surface skin that is ∼500 µm thick. As the wall thickness decreases, the skin's volume fraction rises, and its poor local densification increasingly dominates the part average, leading to an overall decrease in the density of the sintered samples. For the thinnest wall (500 µm) cylinder, macroscopic warping occurred during sintering. We suggest that warping further disrupts particle packing and induces through-thickness nonuniformity, further impairing densification in addition to the skin-fraction effect. These results highlight how geometric miniaturization may adversely influence the sintering behavior and final density of DLP-printed ceramics through distinct densification-limiting mechanisms.

Materials science↗
DOE OSTI2024

Divergence Reduction in Monte Carlo Neutron Transport with On-GPU Asynchronous Scheduling

While Monte Carlo Neutron Transport (MCNT) is near-embarrasingly parallel, the effectively unpredictable lifetime of neutrons can lead to divergence when MCNT is evaluated on GPUs. Divergence is the phenomenon of adjacent threads in a warp executing different control flow paths; on GPUS, it reduces performance because each work group may only execute one path at a time. The process of Thread Data Remapping (TDR) resolves these discrepancies by moving data across hardware such that data in the same warp will be processed through similar paths. A common issue among prior implementations of TDR is the synchronous nature of its remapping and processing cycles, which exhaustively sort data produced by prior processing passes and exhaustively evaluate the sorted data. In another work, we defined a method of remapping data through an asynchronous scheduler which allows for work to be stored in shared memory and deferred arbitrarily until that work is a viable option for low-divergence evaluation. This article surveys a wider set of cases, with the goal of characterizing performance trends across a more comprehensive set of parameters. These parameters include cross sections of scattering/capturing/fission, use of implicit capture, source neutron counts, simulation time spans, and tuned memory allocations. Across these cases, we have recorded minimum and average execution times, as well as a heuristically tuned near-optimal memory allocation size for both synchronous and asynchronous scheduling. Across the collected data, it is shown that the asynchronous method is faster and more memory efficient in the majority of cases, and that it requires less tuning to achieve competitive performance.

Computer Science↗
DOE OSTI2023

The combined effect of temperature and seawater on the compression properties of carbon fiber vinyl ester composites for sandwich structures

For naval applications, composite sandwich structures are of significant interest, and are often manufactured using thin (2 to 4 mm) composite facings made from carbon or glass fiber reinforcement, attached to a thick (25 to 50 mm) section of PVC cellular foam or balsa wood-based core materials using a suitable polymeric resin. In the present study, we focus on the hygrothermal effect on the fiber-dominated compression properties of carbon fiber reinforced vinyl ester resin based polymeric composite (CF/VE), used as “skin” for a polymeric composite sandwich material. Hygrothermal conditioning is achieved by saturating samples in simulated seawater at 40°C. Compression properties are evaluated for coupons extracted along warp and fill undergoing- no conditioning, conditioning till saturation (up to 6 months), and long-term conditioning (2 years). Sea-water saturation yields in up to 12% drop in compression strength with a further 3–4% drop resulting from long-term conditioning. No statistically significant modulus degradation is noticed due to short or long-term hygrothermal exposure. The failure mechanism of the warp extracted coupon, which fails in a splitting failure mode originating due to the delamination between the 0/90 interface, or the fill extracted coupon, which fails due to the instability caused by tow micro-buckling, remains unchanged due to combined exposure (short or long-term) of seawater and temperature. The loss in strength is attributed to the degradation of the fiber-matrix interface, which is validated via conducting single fiber push-in tests with a nominal diameter of 7 micron for conditioned and unconditioned coupons.

Engineering↗
DOE OSTI2020

Simulation of Nanoscale High-Efficiency Thermionic Converters

LBNL provided development and feedback for Warp's mesh refinement solver and cutcell capabilities, authoring adjusted implementations for use with the dielectric solver in 3D. This work has been documented in a separate branch of the RadiaSoft Warp fork, and is currently under review. LBNL has also provided feedback on the new tests developed to study dielectric charging and shaped emitters. The LBNL contributions enabled the Radiasoft team to progress steadily and meet the goal of the projects for the duration of the collaboration. Results have been reported in one refereed paper, three proceeding papers, four oral and six poster presentions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Graphics processor unit with opportunistic inter-path reconvergence

A graphics processing unit and methods for comping and executing instructions with opportunistic inter-path reconvergence are provided. A graphics processing unit may access computer executable instructions mapped to code blocks of a control flow for a warp. The code blocks may include an immediate dominator block and an intermediate post dominator block. The graphics processing unit may store a first thread mask associated with the first code block. The first thread mask may include a plurality of bits indicative of the active or non-active status for the threads of the warp, respectively. The graphics processing unit may a second thread mask corresponding to an intermediate code block between the immediate dominator block and intermediate post dominator block. The graphics processing unit may execute, with threads indicated as active by the first thread mask, instructions of the intermediate code block with a first operand or a second operand depending on the second thread mask.

Kulkarni, Milind↗
DOE OSTI2020

Time Series Classification for Locating Forced Oscillation Sources

Here, this article presents a machine learning based time-series classification method for using synchrophasor measurements to locate the source of forced oscillation (FO) for fast disturbance removal. First, multivariate time series (MTS) matrices are constructed by the most informative measurements selected by sequential feature selection from each power plant. Then, the Mahalanobis matrix is trained such that the Mahalanobis distance between the MTSs from the same class (i.e., with the same FO source location) are minimized and from different classes (i.e., with different FO source locations) are maximized. This allows MTSs to be classified by classifiers with class membership corresponding to the location of each FO source. To meet the runtime requirements of online matching, class templates are constructed to reduce data size and improve matching efficiency. To account for uncertainty in identifying the exact beginning of an FO event, dynamic time warping is used to align the out-of-sync MTSs. IEEE 39bus and WECC 179bus systems are used for algorithm development and validation. Simulation results demonstrate that the algorithm meets online operation runtime requirement with high accuracy using misaligned data sets.

42 ENGINEERING↗
DOE OSTI2022

Large Magnetic Gap in a Designer Ferromagnet–Topological Insulator–Ferromagnet Heterostructure

Abstract Combining magnetism and nontrivial band topology gives rise to quantum anomalous Hall (QAH) insulators and exotic quantum phases such as the QAH effect where current flows without dissipation along quantized edge states. Inducing magnetic order in topological insulators via proximity to a magnetic material offers a promising pathway toward achieving the QAH effect at a high temperature for lossless transport applications. One promising architecture involves a sandwich structure comprising two single‐septuple layers (1SL) of MnBi 2 Te 4 (a 2D ferromagnetic insulator) with ultrathin few quintuple layer (QL) Bi 2 Te 3 in the middle, and it is predicted to yield a robust QAH insulator phase with a large bandgap greater than 50 meV. Here, the growth of a 1SL MnBi 2 Te 4 /4QL Bi 2 Te 3 /1SL MnBi 2 Te 4 heterostructure via molecular beam epitaxy is demonstrated and the electronic structure probed using angle‐resolved photoelectron spectroscopy. Strong hexagonally warped massive Dirac fermions and a bandgap of 75 ± 15 meV are observed. The magnetic origin of the gap is confirmed by the observation of the exchange‐Rashba effect, as well as the vanishing bandgap above the Curie temperature, in agreement with density functional theory calculations. These findings provide insights into magnetic proximity effects in topological insulators and reveal a promising platform for realizing the QAH effect at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Joint statistics of cosmological constant and SUSY breaking in flux vacua with nilpotent Goldstino

We obtain the joint distribution of the gravitino mass and the cosmological constant in KKLT and LVS models with anti-D3 brane uplifting described via the nilpotent goldstino formalism. Moduli stabilisation (of both complex structure and Kähler moduli) is incorporated so that we sample only over points corresponding to vacua. Our key inputs are the distributions of the flux superpotential, the string coupling and the hierarchies of warped throats. In the limit of zero cosmological constant, we find that both in KKLT and LVS the distributions are tilted favourably towards lower scales of supersymmetry breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Local SU(2) × U(1) quark flavor symmetry in the RS bulk

We propose a model of quark flavor based on an additional SU(2) × U(1) local symmetry in a warped extra dimensional bulk. In contrast to other works, we break the additional gauge symmetry in the bulk via two complex scalars which acquire bulk vevs, rather than relying on brane-localized symmetry breaking. A gauge-covariant Kaluza-Klein decomposition of a theory with a bulk spontaneously broken gauge symmetry is performed, and exact expressions for the bulk profiles of all physical particles in such systems are given. The SM quark masses and mixings are then recreated using gauge-covariant bulk quark mass terms and Yukawa-like couplings to the new bulk scalars. A numerical sampling of points in the model parameter space that recreate the quark masses and mixings is performed at a KK scale of M KK = 5 TeV. We then compute the ΔF = 2 4-quark operators arising from our new flavor gauge bosons and scalars, and those arising from Kaluza-Klein modes of SM gauge bosons. By decoupling one of our bulk scalar fields to all quark fields except the right-handed up-like sector, we find that it is possible to greatly suppress tree-level contributions to the highly constrained Kaon mixing parameters. Instead, the dominant constraints on the model emerge from neutral Bd and D meson mixing. These constraints are explored with our numerical sampling of the model parameter space, and the specific contribution of the new flavor gauge bosons and scalars is discussed. We find that for a significant range of realistic flavor gauge couplings, the new gauge bosons compete with the normally dominant gluon flavor-changing currents, but flavor-changing operators emerging from the bulk scalar fields are highly suppressed. Finally, we briefly comment on flavor constraints that are independent of the flavor gauge sector arising from the Z$\overline{b}$ L b L coupling and rare top decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Search for resonant pair production of Higgs bosons in the $\textrm{b}\overline{\textrm{b}}\textrm{b}\overline{\textrm{b}}$ final state using large-area jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search is presented for the resonant production of a pair of standard model-like Higgs bosons using data from proton-proton collisions at a centre-of-mass energy of 13 TeV, collected by the CMS experiment at the CERN LHC in 2016–2018, corresponding to an integrated luminosity of 138 fb −1 . The final state consists of two b quark-antiquark pairs. The search is conducted in the region of phase space where at least one of the pairs is highly Lorentz-boosted and is reconstructed as a single large-area jet. The other pair may be either similarly merged or resolved, the latter reconstructed using two b-tagged jets. The data are found to be consistent with standard model processes and are interpreted as 95% confidence level upper limits on the product of the cross sections and the branching fractions of the spin-0 radion and the spin-2 bulk graviton that arise in warped extradimensional models. The limits set are in the range 9.74–0.29 fb and 4.94–0.19 fb for a narrow radion and a graviton, respectively, with masses between 1 and 3 TeV. For a radion and for a bulk graviton with widths 10% of their masses, the limits are in the range 12.5–0.35 fb and 8.23–0.23 fb, respectively, for the same masses. These limits result in the exclusion of a narrow-width graviton with a mass below 1.2 TeV, and of narrow and 10%-width radions with masses below 2.6, and 2.9 TeV, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Search for high-mass resonances in a final state comprising a gluon and two hadronically decaying W bosons in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for high-mass resonances decaying into a gluon, g, and two W bosons is presented. A Kaluza-Klein gluon, g$_{KK}$, decaying in cascade via a scalar radion R, g$_{KK}$ → gR → gWW, is considered. The final state studied consists of three large-radius jets, two of which contain the products of hadronically decaying W bosons, and the third one the hadronization products of the gluon. The analysis is performed using proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the CMS experiment at the CERN LHC during 2016–2018, corresponding to an integrated luminosity of 138 fb$^{−1}$. The masses of the g$_{KK}$ and R candidates are reconstructed as trijet and dijet masses, respectively. These are used for event categorization and signal extraction. No excess of data events above the standard model background expectation is observed. Upper limits are set on the product of the g$_{KK}$ production cross section and its branching fraction via a radion R to gWW. This is the first analysis examining the resonant WW+jet signature and setting limits on the two resonance masses in an extended warped extra-dimensional model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Kinetic mixing, dark photons and extra dimensions. Part III. Brane localized dark matter

Extra dimensions have proven to be a very useful tool in constructing new physics models. In earlier work, we began investigating toy models for the 5-D analog of the kinetic mixing/vector portal scenario where the interactions of dark matter, taken to be, e.g., a complex scalar, with the brane-localized fields of the Standard Model (SM) are mediated by a massive U(1) D dark photon living in the bulk. These models were shown to have many novel features differentiating them from their 4-D analogs and which, in several cases, avoided some well-known 4-D model building constraints. However, these gains were obtained at the cost of the introduction of a fair amount of model complexity, e.g., dark matter Kaluza-Klein excitations. In the present paper, we consider an alternative setup wherein the dark matter and the dark Higgs, responsible for U(1) D breaking, are both localized to the ‘dark’ brane at the opposite end of the 5-D interval from where the SM fields are located with only the dark photon now being a 5-D field. The phenomenology of such a setup is explored for both flat and warped extra dimensions and compared to the previous more complex models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Mapping out the internal space in AdS/BCFT with Wilson loops

We study Wilson loops in string theory realizations of AdS/BCFT and wedge holography. The field theories are based on 3d N = 4 long quiver gauge theories engineered by D3, D5 and NS5 branes, and on BCFTs involving 4d N = 4 SYM coupled to such 3d theories. The holographic duals have geometry AdS 4 × S 2 × S 2 warped over a strip. We identify the holographic representation of antisymmetric Wilson loops associated with individual 3d gauge nodes in terms of probe D5-branes. The expectation values obtained holographically are matched to supersymmetric localization computations. Our results yield an identification of regions in the internal space with individual 3d gauge nodes. Connecting to bottom-up braneworld models, this gives a concrete notion of which parts of the 10d solutions correspond to the end-of-the-world brane and which to the remaining bulk. We also construct supersymmetric Janus on the brane embeddings in AdS 5 ×S 5 , which describe surface defects with boundaries and interfaces in N = 4 SYM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗