Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Friction-mediated phase transition in confined active nematics

Using a minimal continuum model, we investigate the interplay between circular confinement and substrate friction in active nematics. Upon increasing the friction from low to high, we observe a dynamical phase transition from a circulating flow phase to an anisotropic flow phase in which the flow tends to align perpendicular to the nematic director at the boundary. We demonstrate that both the flow structure and dynamic correlations in the latter phase differ from those of an unconfined, active turbulent system and may be controlled by the prescribed nematic boundary conditions. Here, our results show that substrate friction and geometric confinement act as valuable control parameters in active nematics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Interactions Remove the Quantization of the Chiral Photocurrent at Weyl Points

The chiral photocurrent or circular photogalvanic effect (CPGE) is a photocurrent that depends on the sense of circular polarization. In a disorder-free, non-interacting chiral Weyl semimetal, the magnitude of the effect is approximately quantized with a material-independent quantum $e^3/h^2$ for reasons of band topology. We study the first-order corrections due to the Coulomb and Hubbard interactions in a continuum model of a Weyl semimetal in which known corrections from other bands are absent. We find that the inclusion of interactions generically breaks the quantization. The corrections have a weaker dependence on the form of the cutoff than previously studied interaction corrections to the (non-topological) linear optical conductivity of graphene, and a potentially observable frequency dependence. Here, we conclude that, unlike the quantum Hall effect in gapped phases or the chiral anomaly in field theories, the quantization of the CPGE in Weyl semimetals is not protected but has perturbative corrections in interaction strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Twist Angle-Dependent Interlayer Exciton Lifetimes in van der Waals Heterostructures

In van der Waals (vdW) heterostructures formed by stacking two monolayers of transition metal dichalcogenides, multiple exciton resonances with highly tunable properties are formed and subject to both vertical and lateral confinement. We investigate how a unique control knob, the twist angle between the two monolayers, can be used to control the exciton dynamics. We observe that the interlayer exciton lifetimes in MoSe 2 /WSe 2 twisted bilayers (TBLs) change by one order of magnitude when the twist angle is varied from 1° to 3.5°. Additionally, using a low-energy continuum model, we theoretically separate two leading mechanisms that influence interlayer exciton radiative lifetimes. The shift to indirect transitions in the momentum space with an increasing twist angle and the energy modulation from the moiré potential both have a significant impact on interlayer exciton lifetimes. We further predict distinct temperature dependence of interlayer exciton lifetimes in TBLs with different twist angles, which is partially validated by experiments. While many recent studies have highlighted how the twist angle in a vdW TBL can be used to engineer the ground states and quantum phases due to many-body interaction, our studies explore its role in controlling the dynamics of optically excited states, thus, expanding the conceptual applications of “twistronics”.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Exact Diagonalization for Magic-Angle Twisted Bilayer Graphene

We report on finite-size exact-diagonalization calculations in a Hilbert space defined by the continuum-model flat moiré bands of magic angle twisted bilayer graphene. For moiré band filling 3>|ν|>2, where superconductivity is strongest, we obtain evidence that the ground state is a spin ferromagnet. Near |ν|=3, we find Chern insulator ground states that have spontaneous spin, valley, and sublattice polarization, and demonstrate that the anisotropy energy in this order-parameter space is strongly band-filling-factor dependent. We emphasize that inclusion of the remote band self-energy is necessary for a reliable description of magic angle twisted bilayer graphene flat band correlations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Reentrant Correlated Insulators in Twisted Bilayer Graphene at 25 T ( 2 π Flux)

Twisted bilayer graphene (TBG) is remarkable for its topological flat bands, which drive strongly interacting physics at integer fillings, and its simple theoretical description facilitated by the Bistritzer-MacDonald Hamiltonian, a continuum model coupling two Dirac fermions. Because of the large moiré unit cell, TBG offers the unprecedented opportunity to observe reentrant Hofstadter phases in laboratory-strength magnetic fields near 25 T. This Letter is devoted to magic angle TBG at 2π flux where the magnetic translation group commutes. Here we use a newly developed gauge-invariant formalism to determine the exact single-particle band structure and topology. We find that the characteristic TBG flat bands reemerge at 2π flux, but, due to the magnetic field breaking C 2z T, they split and acquire Chern number ±1. We show that reentrant correlated insulating states appear at 2π flux driven by the Coulomb interaction at integer fillings, and we predict the characteristic Landau fans from their excitation spectrum.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magic Angles and Fractional Chern Insulators in Twisted Homobilayer Transition Metal Dichalcogenides

We explain the appearance of magic angles and fractional Chern insulators in twisted K-valley homobilayer transition metal dichalcogenides by mapping their continuum model to a Landau level problem. Here, our approach relies on an adiabatic approximation for the quantum mechanics of valence band holes in a layer-pseudospin field that is valid for sufficiently small twist angles and on a lowest Landau level approximation that is valid for sufficiently large twist angles. It provides a simple qualitative explanation for the nearly ideal quantum geometry of the lowest moir´e miniband at particular twist angles, predicts that topological flat bands occur only when the valley-dependent moir´e potential is sufficiently strong compared to the interlayer tunneling amplitude, and provides a convenient starting point for the study of interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Moiré Exchange Effect in Twisted WSe 2 /WS 2 Heterobilayer

Moiré superlattices of layered transition metal dichalcogenides are proven to host periodic electron crystals due to strong correlation effects. These electron crystals can also be intertwined with intricate magnetic phenomena. Here, in this Letter, we present our findings on the moiré exchange effect, resulting from the modulation of local magnetic moments by electron crystals within well-aligned WSe 2 /WS 2 heterobilayers. Employing polarization-resolved magneto-optical spectroscopy, we unveil a high-energy excitonic resonance near one hole per moiré unit cell (v=–1), which possesses a giant g factor several times greater than the already very large g factor of the WSe 2 A exciton in this heterostructure. Supported by continuum model calculations, these high-energy states are found to be dark excitons brightened through Umklapp scattering from the moiré mini-Brillouin zone. When the carriers form a Mott insulating state near v=–1, the Coulomb exchange between doped carriers and excitons forms an effective magnetic field with moiré periodicity. This moiré exchange effect gives rise to the observed giant g factor for the excitonic Umklapp state.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonmonotonic Band Flattening near the Magic Angle of Twisted Bilayer MoTe 2

Twisted bilayer MoTe 2 (tMoTe 2 ) is an emergent platform for exploring exotic quantum phases driven by the interplay between nontrivial band topology and strong electron correlations. Direct experimental access to its momentum-resolved electronic structure is essential for uncovering the microscopic origins of the correlated topological phases therein. Here, we report angle-resolved photoemission spectroscopy measurements of tMoTe 2 , revealing pronounced twist-angle-dependent band reconstruction shaped by orbital character, interlayer coupling, and moiré potential modulation. Density functional theory captures the qualitative evolution, yet underestimates key energy scales across twist angles, highlighting the importance of electronic correlations. Notably, the hole effective mass at the 𝐾 point exhibits a nonmonotonic dependence on twist angle, peaking near 2°, consistent with band flattening at the magic angle predicted by continuum models. Via electrostatic gating and surface dosing, we further visualize the evolution of electronic structure versus doping, enabling direct observation of the conduction band minimum and confirm tMoTe 2 as a direct band gap semiconductor. These results establish a spectroscopic foundation for modeling and engineering emergent quantum phases in this moiré platform.

2D materials↗
DOE OSTI2024

Higgs-confinement transitions in QCD from symmetry protected topological phases

In gauge theories with fundamental matter there is typically no sharp way to distinguish confining and Higgs regimes, e.g. using generalized global symmetries acting on loop order parameters. It is standard lore that these two regimes are continuously connected, as has been explicitly demonstrated in certain lattice and continuum models. We point out that Higgsing and confinement sometimes lead to distinct symmetry protected topological (SPT) phases -- necessarily separated by a phase transition -- for ordinary global symmetries. We present explicit examples in 3+1 dimensions, obtained by adding elementary Higgs fields and Yukawa couplings to QCD while preserving parity P and time reversal T. In a suitable scheme, the confining phases of these theories are trivial SPTs, while their Higgs phases are characterized by non-trivial P- and T-invariant theta-angles θ f , θ g = π for flavor or gravity background gauge fields, i.e. they are topological insulators or superconductors. Finally, we consider conventional three-flavor QCD (without elementary Higgs fields) at finite U(1) B baryon-number chemical potential μ B , which preserves P and T. At very large μ B , three-flavor QCD is known to be a completely Higgsed color superconductor that also spontaneously breaks U(1) B . We argue that this high-density phase is in fact a gapless SPT, with a gravitational theta-angle θ g = π that safely co-exists with the U(1) B Nambu-Goldstone boson. We explain why this SPT motivates unexpected transitions in the QCD phase diagram, as well as anomalous surface modes at the boundary of quark-matter cores inside neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Helium Bubble Growth and Retention Characteristics in Aging Palladium Tritide

The helium-related characteristics of aged palladium tritide, including tritium loss, the formation of He bubbles from tritium decay-induced He atoms, and He thermal desorption behavior are computed using a continuum model of bubble evolution that incorporates the log-normal bubble spacing distribution deduced from 3 He nuclear magnetic resonance (NMR) measurements. This spacing distribution produces significant differences between mean (expectation) values of the bubble size and pressure and the average (median) values of these aging characteristics obtained by a previous calculation. The new calculations find much of the He is retained in a large quantity of smaller bubbles with higher bubble pressures, improving the overall stability of the bubbles and extending the age of He retention. By contrast, the integrated tritide material characteristics of swelling and tritium pressure are found relatively insensitive to details of the bubble spacing distribution. Inter-bubble fracture is predicted to begin with small, closely-spaced bubbles and progress to include larger bubbles near the onset of rapid He release. The critical concentration for this onset increases almost linearly with material tensile strength and decreases with shear modulus. Release of He from the aging solid is modeled as a growing network of fractured inter-bubble ligaments. For small particle or small grain material, the physical size of the fracture cluster, relative to the material dimensions, becomes important. In contrast with normal aging, inter-bubble fracture occurring during thermal desorption appears to begin with large, widely-spaced bubbles, indicating elevated temperature techniques may be of limited use for evaluating development of the bubble fracture network.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Multimode Metastructures: Novel Hybrid 3D Lattice Topologies

With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Magnetism in Moiré Materials

This project focused on the theory of the electronic properties of magnetic moire materials. Moire materials can be formed by overlaying two or more two-dimensional crystals that are either semiconductors or semimetals and establishing a moire pattern. Low-energy electronic properties are then described by a continuum model with crystalline periodicity. Moire materials are therefore artificial two-dimensional crys- tals whose periodicity is at the tens of nanometers scale, increasing the unit cell area by a factor of about 10,000 compared to real crystals. They are important because the main tuning knob of chemistry, the number of electrons per effective atom or molecule (the analog of the number of electrons per unit cell in real crystals), can be shifted by more than one simply by using electrical gates. The research covered two diffeerent classes of moire materials, graphene multilayers twisted to a rotation angle at which correlations are strong, and transition metal dichalcogenide semiconductors. The proposal identifed one project in each area both related very broadly to magnetism: i) Advancing understanding of broken spin-valley and sublattice symmetries in Magic Angle Twisted Bilayer Graphene and ii) Developing the theory of itinerant electron magnetism in Gamma-valley and K-valley transition metal dichalcogenide moire materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Scaling Turbulent Combustion Fields in Explosions

We considered the topic of explosions from spherical high-explosive (HE) charges. We studied how the turbulent combustion fields scale. On the basis of theories of dimensional analysis by Bridgman and similarity theories of Sedov and Barenblatt, we found that all fields scaled with the explosion length scale r0. This included the blast wave, the mean and root mean squared (RMS) profiles of thermodynamic variables, combustion variables, velocities, vorticity, and turbulent Reynolds stresses. This was a consequence of the formulation of the problem and our numerical method, which both satisfied the similarity conditions of Sedov. We performed numerical simulations of 1 g charges and 1 kg charges; the solutions were identical (within roundoff error) when plotted in scaled variables. We also explored scaling laws related to three-phase pyrotechnic explosions. We show that although the scaling formally broke down, the fireball still essentially scaled with the explosion length scale r0. However, the discrete Lagrange particles (DLP) (phase 2) and the heterogeneous continuum model (HCM) of the DLP wakes (phase 3) did not scale with r0, and mean and RMS profiles could differ by a factor of 10 in some regions. This was because the DLP particles and wakes introduced an additional scale that broke the similarity conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

65M75, 70H05, 70G65↗
NASA NTRS1973

Contact problem for an elastic reinforcement bonded to an elastic plate

The stiffening layer is treated as an elastic membrane and the base plate is assumed to be an elastic continuum. The bonding between the two materials is assumed to be either one of direct adhesion ro through a thin adhesive layer which is treated as a shear spring. The solution for the simple case in which both the stiffener and the base plate are treated as membranes is also given. The contact stress is obtained for a series of numerical examples. In the direct adhesion case the contact stress becomes infinite at the stiffener ends with a typical square root singularity for the continuum model, and behaving as a delta function for the membrane model. In the case of bonding through an adhesive layer the contact stress becomes finite and continuous along the entire contact area.

Erdogan, F.↗