Search NASA⌕ Search
Home/Archive/DOE OSTI

DOE OSTI · 3023365

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Leng, Yu [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:000000026128409X)·Buechler, Miles Allen [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000156354594)·Manzini, Gianmarco [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000336263112)
View original source ↗Read source document ↗Export metadata ↓

Abstract

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

Explore related subjects

97 MATHEMATICS AND COMPUTINGFinite volume methodVirtual element methodhourglass controlstaggered-grid hydrodynamics

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Leng, Yu [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:000000026128409X), Buechler, Miles Allen [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000156354594), Manzini, Gianmarco [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000336263112). 2026-03-12. Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques. https://doi.org/10.1016/j.jcp.2026.114842

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

DOE OSTI2026

TANTE: Time-adaptive operator learning via neural Taylor expansion

Operator learning for time-dependent partial differential equations (PDEs) has seen rapid progress in recent years, enabling efficient approximation of complex spatiotemporal dynamics. However, most existing methods rely on fixed time step sizes during rollout, which limits their ability to adapt to varying temporal complexity and often leads to error accumulation. In this work, we propose the Time-Adaptive Transformer with Neural Taylor Expansion (TANTE), a novel operator-learning framework that produces continuous-time predictions with adaptive step sizes. TANTE predicts future states by performing a Taylor expansion at the current state, where neural networks learn both the higher-order temporal derivatives and the local radius of convergence. This allows the model to dynamically adjust its rollout based on the local behavior of the solution, thereby reducing cumulative error and improving computational efficiency. We demonstrate the effectiveness of TANTE across a wide range of PDE benchmarks, achieving superior accuracy and adaptability compared to fixed-step baselines, delivering accuracy gains of 60-80 % and speed-ups of 30-40 % at inference time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Structured illumination for surface-resolved grazing-incidence X-ray scattering

Grazing-incidence (GI) scattering techniques are widely used to characterize thin films, offering high surface sensitivity and insight into morphology and structure. However, these approaches typically provide statistical averaged information due to elongated footprint or limited spatial resolution due to beam size. Here we introduce a method that combines structured illumination with GI X-ray scattering and leverages our computational imaging approach to resolve local structural details. We demonstrate that our method captures local features of an organic semiconductor thin film without the need for sample rotation as in tomography. The method expands GI techniques from statistical averaging to high-resolution imaging, thereby providing the capability for detailed analysis of local material properties, such as domain shape, orientation and polymorphism, which are critical for advancing material design towards more efficient and tailored materials.

97 MATHEMATICS AND COMPUTING↗