SEARCH · Engineering Papers

Results for “Unstructured”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗

DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

DOE OSTI2026

Preliminary neutronic and thermal analysis of Inertial Fusion Energy systems with a focus on thick liquid wall concepts

This paper provides an assessment of the neutronic behavior of Inertial Fusion Energy systems using the Serpent 2 Monte Carlo code. The study explores different types of reactor concepts, including dry-wall and thick liquid wall systems. In particular, it characterizes the tritium breeding and energy multiplication performance, as well as the effect of long-term irradiation on structural materials. For this purpose, the displacement per atom (DPA) and the gas production rate in structural materials are calculated. A simplified geometry with varying geometrical parameters is used as a representation for the different kinds of devices. Additionally, the case of the HYLIFE-II concept is used to assess modeling inaccuracies arising from certain approximations. A comparison between the simplified geometry and a realistic one is performed, thanks to Serpent’s capability to run transport calculations on unstructured meshes. The study then addresses the impact of these discrepancies on the temperature distribution in the first wall using the GeN-Foam code for conjugate heat transfer analysis. The results highlight the necessity to consider the heterogeneous nature of the liquid wall, as simplified models tend to underestimate radiation damage metrics and heat loads.

Conjugate heat transfer↗

DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗

DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗

DOE OSTI2025

Graph neural networks for mechanical property prediction of 2D fiber composites

This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.

Brittle strength↗

DOE OSTI2025

Coupled beam-target-moderator optimization for the Second Target Station

To support the design of the Second Target Station, that aims to provide the world’s highest peak brightness of cold neutrons, studies that optimize the dimensions of the target and moderators are invaluable. In this work, we investigate the influence of the target dimensions and beam profile on the performance and optimal size of the moderators. We perform optimization runs with detailed MCNP6.2 simulations using high-fidelity unstructured mesh geometries generated from parametrized CAD models. We demonstrate that small changes in target height do not influence the moderator performance if the beam dimensions are chosen adequately. We quantify the effect of beam footprint and target width on the moderator performance and show that the optimal moderator dimensions are insensitive to limited changes in target and beam profile.

Dakota↗

DOE OSTI2024

Debunking common myths in coastal circulation modeling

Despite tremendous progress in algorithm development, computational efficiency and transition into operations over the past two decades, coastal modeling still lacks scientific rigor due to proliferation of many ‘gray’ areas related to various modeling choices made by modelers. Here, in this paper, we propose some guiding principles for the modeling community to improve performance, and we also debunk commonly held myths that make the coastal modeling lack rigor. Using our own experience in developing seamless cross-scale unstructured-grid based models for the past two decades, we describe in unprecedented detail the end-to-end modeling process (i.e., from digital elevation models (DEMs) to mesh generation to post analysis), and demonstrate that defensible modeling is within reach for any end user by following three guiding principles: (1) Bathymetry is a first order forcing in coastal domains and thus should be respected in all aspects of modeling; (2) Oceanographic processes are driven across multiple spatial scales and so models should enable appropriate resolution as needed; and (3) Model assessment should focus on physical processes. Through qualitative and quantitative model assessments, we demonstrate the fundamental role played by bathymetry/topography as embedded in DEMs in making the results defensible, which is unfortunately glossed over in many modeling studies. Focusing on process-based assessment simplifies the calibration process. A major conclusion of this work is that model developers and operators should maximize the scientific rigor for in silico oceanography by avoiding some common pitfalls that rely on error compensation at the expense of representation of physical system processes. We present some best practice procedures for defensive and trustworthy numerical modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Total Water Level Prediction at Continental Scale: Coastal Ocean

We demonstrate recent progress made in the simulation of total water level (TWL) at continental scale, using the coastal ocean of US East Coast/Gulf of Mexico coast as an example. A key difference between the continental-scale and small-scale modeling is that the former requires a more accurate vertical datum. Using a geoid-based datum (xGEOID20b), a satellite altimetry product, and a state-of-the-art 3D unstructured-grid model, we significantly improve the accuracy for TWL both near- and off-shore. The average root-mean-square error at all NOAA stations is 14 cm. The non-tidal signals are found to be sensitive to the representation of a large-scale current system near the boundary and extending the domain extent to accommodate this system improves these signals.

Cui, Linlin↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

3D Poisson equation↗

DOE OSTI2025

Characterization of Gulf Stream energy resources: A 30-year high-resolution simulation study

The strength and proximity of the Gulf Stream (GS) along the U.S. coast offers opportunities for marine energy extraction. Leveraging high-resolution (up to 400 m) long-term (30-year) hindcasts based on an unstructured grid model, this study characterizes the energy potential of the GS in two viable locations: the Florida Straits (FS) and Cape Hatteras (CH). Results reveal distinct characteristics in the kinetic energy flux (KEF) and kinetic energy density (KED) for the two regions. The 30-year mean KEF at CH is 27.50 GW, surpassing the 19.74 GW estimated for the FS, which is attributed to the wider cross-section of the GS in CH. Additionally, CH exhibits higher KEF variability due to GS’s larger path shifts and meanders. Estimated KEDs at ADCP site locations reach up to 2,908 W/m 2 in the FS and 1,512 W/m 2 in the CH at a depth of 20 m. High-energy zones are closer to shore in the FS, whereas energy-rich areas at CH are located further offshore. A depth-based KED analysis shows significant decreases across typical energy extraction depths (20 to 100 m), with reductions up to 44%. Finally, this study provides high-resolution resource characterization, forming a foundation for optimized site selection and marine hydrokinetic system design.

16 TIDAL AND WAVE POWER↗

DOE OSTI2024

Street-level temperature estimation using graph neural networks: Performance, feature embedding and interpretability

Estimating street-level air temperature is a challenging task due to the highly heterogeneous urban surfaces, canyon-like street morphology, and the diverse physical processes in the built environment. Though pioneering studies have embarked on investigations via data-driven approaches, many questions remain to be answered. Here, in this study, we leveraged an innovative framework and redefined the street-level temperature estimation problem using Graph Neural Networks (GNN) with spatial embedding techniques. The results showed that GNN models are more capable and consistent of estimating street-level temperature among tested locations, benefiting from its unique strength in handling extensive data over unstructured graph topology. In addition, we conducted in-depth analysis of feature importance to enhance the model interpretability. Among the urban features analyzed in this study, the time-variant canopy density and meter-level land use data emerge as crucial factors. Our findings highlight GNN 's high potential in capturing the complex dynamics between urban elements and their impacts on microclimate, thus offering valuable insights for comprehensive urban data collection and urban climate modeling in general. Collectively, this study also contributes to urban planning and policy by providing avenues to enhance city resilience against climate change, thereby advancing the agenda for environmental stewardship and urban sustainability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

MechBERT: Language Models for Extracting Chemical and Property Relationships about Mechanical Stress and Strain

Language models are transforming materials-aware naturallanguage processing by enabling the extraction of dynamic, context-rich information from unstructured text, thus, moving beyond the limitations of traditional information-extraction methods. Moreover, small language models are on the rise because some of them can perform better than large language models (LLMs) when given domain-specific questionanswer tasks, especially about an application area that relies on a highly specialized vernacular, such as materials science. We therefore present a new class of MechBERT language models for understanding mechanical stress and strain in materials. These employ Bidirectional Encoder Representations for transformer (BERT) architectures. We showcase four MechBERT models, all of which were pretrained on a corpus of documents that are textually rich in chemicals and their stress–strain properties and were fine-tuned on question-answering tasks. We evaluated the level of performance of our models on domain-specific as well as general English-language question-answer tasks and also explored the influence of the size and type of BERT architectures on model performance. We find that our MechBERT models outperform BERT-based models of the same size and maintain relevancy better than much larger BERT-based models when tasked with domain-specific question-answering tasks within the stress–strain engineering sector. These small language models also enable much faster processing and require a much smaller fraction of data to pretrain them, affording them greater operational efficiency and energy sustainability than LLMs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Static and Dynamic Thermomechanical Properties of Phase-Separated Epoxy Networks with Tuned Microstructures

Here, polymerization-induced phase separation is a useful method for the construction of heterogeneous epoxy networks with properties exceeding their homogeneous counterparts. In this work, we examine the static and dynamic thermomechanical properties of phase-separated epoxy networks salient to their application as encapsulants. Three heterogeneous epoxy-amine networks with nano-, meso-, and macro-phase-separated morphologies comprised of hard and soft domains are compared to a rigid, unstructured network. The glass transition profiles of the heterogeneous networks are complex, spanning many decades in the frequency domain. The nanophase-separated morphology leads to higher coefficient of thermal expansion, yet surprisingly is characterized by reduced residual stress. Under both quasi-static and dynamic compression (strain rates of order 10 –3 and 10 3 s –1 , respectively), the nanophase-separated network also exhibits higher modulus and strength. In split-Hopkinson bar experiments, the energy dissipation characteristics of the epoxy networks were nearly identical. Curiously, however, the Hugoniot response of the macro-phase-separated network determined by ballistic shockwave analysis indicates a remarkable ability of this material to mitigate shockwave propagation in comparison to many homogeneous and heterogeneous polymer materials. Collectively, this work reveals several previously unreported phenomena with respect to structure–property relationships in phase-separated epoxy networks, illustrating the potential value of systematically tuned microstructures for optimization of application-specific physical properties.

36 MATERIALS SCIENCE↗

DOE OSTI2026

The Energy Exascale Earth System Model Version 3: 2. Overview of the Coupled System

The Energy Exascale Earth System Model version 3 (E3SMv3) represents the latest advancement in Earth system modeling developed by the U.S. Department of Energy (DOE). Building upon previous versions, E3SMv3 introduces significant updates across its coupled components to enhance capability and improve fidelity. The atmosphere component incorporates advancements in chemistry, aerosol-cloud interactions, convection, and microphysics. The ocean features a new time-stepping scheme and a higher-resolution unstructured mesh with sub-ice-shelf cavities, while the sea ice model integrates advanced snow and ice physics for more realistic cryospheric simulations. The land model introduces prognostic vegetation dynamics and a new sub-grid topographic treatment of solar radiation. A new tri-grid configuration harmonizes the horizontal grids of the land and river components for improved process coupling. It is enabled by a new non-linear remapping between the atmosphere and land. E3SMv3 underwent extensive testing through a comprehensive simulation campaign, including pre-industrial control, idealized CO 2 experiments, and historical simulations spanning 1850–2024. The model demonstrates significant improvements in simulating the evolution of the historical surface temperature, particularly addressing the “pothole cooling” bias in earlier versions. Reduced aerosol-related forcing contributes to more realistic radiative forcing and better alignment with the observational record. Ocean heat content (OHC) and sea ice trends are also improved as a result.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Systematic benchmarking demonstrates large language models have not reached the diagnostic accuracy of traditional rare-disease decision support tools

Large language models (LLMs) show promise in supporting differential diagnosis, but their performance is challenging to evaluate due to the unstructured nature of their responses, and their accuracy compared to existing diagnostic tools is not well characterized. To assess the current capabilities of LLMs to diagnose genetic diseases, we benchmarked these models on 5213 previously published case reports using the Phenopacket Schema, the Human Phenotype Ontology and Mondo disease ontology. Prompts generated from each phenopacket were sent to seven LLMs, including four generalist models and three LLMs specialized for medical applications. The same phenopackets were used as input to a widely used diagnostic tool, Exomiser, in phenotype-only mode. The best LLM ranked the correct diagnosis first in 23.6% of cases, whereas Exomiser did so in 35.5% of cases. While the performance of LLMs for supporting differential diagnosis has been improving, it has not reached the level of commonly used traditional bioinformatics tools. Future research is needed to determine the best approach to incorporate LLMs into diagnostic pipelines.

Reese, Justin T. [Lawrence Berkeley National Labor↗

DOE OSTI2025

Space-time crystals from particle-like topological solitons

Time crystals are unexpected states of matter that spontaneously break time-translation symmetry either in a discrete or continuous manner. However, spatially mesoscale space-time crystals that break both space and time symmetries have not been reported. Here we report a continuous space-time crystal in a nematic liquid crystal driven by ambient-power, constant-intensity unstructured light. Our numerically constructed four-dimensional configurations exhibit good agreement with these experimental findings. Although meeting the established criteria to identify time-crystalline order, both experiments and computer simulations reveal a space-time crystallization phase formed by particle-like topological solitons. The robustness against temporal perturbations and spatiotemporal dislocations shows the stability and rigidity of the studied space-time crystals, which relates to their locally topological nature and many-body interactions between emergent spontaneously twisted, particle-like solitonic building blocks. Their potential technological utility includes optical devices, photonic space-time crystal generators, telecommunications and anti-counterfeiting designs, among others.

Liquid crystals↗