Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subdivision method”

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.

9 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗
DOE OSTI2024

At-power subcritical multiplication in the Advanced Test Reactor during nuclear requalification testing

Power division information during nuclear requalification of the Advanced Test Reactor (ATR) is of considerable interest as an importance function for observed changes to core reactivity. The degree to which a given physical subdivision of a critical reactor acts as a neutron source for other lobes is not analytically characterized for general application. When ATR operates at power, individual power-producing lobes rely on each other as neutron sources in order to maintain constant power, which in general requires either exactly critical multiplication within a reactor or an external neutron source. Here, this work shows that fuel element and lobe powers in ATR can be related with subcritical multiplication theory. Subcritical multiplication factors are computed with a physically validated analytical method based on actual at-power operation, quantifying for each lobe its dependence on other lobes as an external neutron source. This explanation is significant for ATR due to the desire to irradiate a large variety of experiments simultaneously, each having its impact on the core neutron population. For any physical subdivision of any other critical reactor, it is likewise true that the subdivision undergoes only subcritical multiplication.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Normalized Solar-Induced Fluorescence Responds Earlier Than Vegetation Indices to the 2019 North China Plain Drought

Recently, solar-induced chlorophyll fluorescence (SIF) from satellites has shown potential for evaluating vegetation status and stress responses. Fluorescence quantum yield ($Φ_F$) is essentially linked to vegetation stress. However, the complex physiological and structural responses of SIF and $Φ_F$ to drought need further study. This study normalized SIF as SIFn to account for angular variations and fluctuations in photosynthetically active radiation (PAR), aiming for more accurate drought monitoring. SIFn anomalies were compared to historical baselines (2019–2021 averages) of vegetation indices (VIs), raw SIF, and $Φ_F$ during a 2019 drought in the North China Plain (NCP). Here, the results show SIFn provides an effective method for drought monitoring, showing the earliest decline compared to raw SIF, VIs, and $Φ_F$. In the first two weeks of drought, SIFn decreased by 8.2%, 7.0%, 12.5%, and 8.2% across the four NCP subdivisions. SIFn outperformed other indicators, proving sensitive to early drought detection. SIFn was also examined for tracking drought alleviation by rainfall. The uncertainty under different viewing geometries was quantified. SIFn anomalies showed a strong correlation with rainfall anomalies (R: 0.45 ~ 0.52) and meteorological factors like PAR (R: 0.80 ~ 0.84) and relative humidity (R:0.52 ~ 0.54). The correlation of near-infrared reflectance (NIRv) and $Φ_F$ anomalies with SIF was weak during drought onset (R: 0.16 ~ 0.32) but strong at the end (R: 0.83 ~ 0.87). These suggest both canopy structure (mainly characterized by NIRv) and vegetation chlorophyll ($Φ_F$) are impacted by drought and influence SIF at different stages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Informed unsupervised machine learning analysis of dislocation microstructure from high-resolution differential aperture X-ray structural microscopy data

This study leverages high-resolution differential-aperture X-ray structural microscopy (DAXM) to probe the local dislocation structure in deformed 304L-stainless steel at small strain, by measuring the lattice rotation and deviatoric elastic strain with a sub-micron resolution. For a single grain in a polycrystalline specimen, the measured lattice rotation field over the measured volume exhibited a multimodal distribution while the deviatoric elastic strain showed a single-mode distribution. An unsupervised Cauchy mixture machine learning model was developed to resolve the multimodal distribution of the lattice rotation. By mapping the lattice rotation data associated with each Cauchy peak in the model back onto the measured volume, we identify contiguous regions of the crystal rotated near the average values corresponding to the peaks of the overall rotation distribution. These regions represent the grain subdivision in the microstructure. Finally, the dislocation density tensor was also computed and its norm was laid over the rotation field to detect the subgrain boundaries. This step provided a validation of the Cauchy mixture model for the analysis of the lattice rotation distribution. The current study highlights the integration of advanced X-ray microscopy techniques with data-driven analysis methods to uncover detailed microstructure scales in deformed crystals.

Machine learning; Lattice rotation; High-energy X-↗
DOE OSTI2025

Defect formation and microstructure tuning via proton irradiation to control electrochemical and phase reversibility in layered battery materials

The reversibility of phase transformation influences the functionality of electrode materials in batteries. In many battery materials, nanosized grains favor phase reversibility but at the cost of cyclability due to aggravated side reactions with the electrolyte. Here, in this study, we present a novel approach to enhance the phase transformation reversibility of layered oxide cathodes, exemplified by Na 2/3 Fe 1/2 Mn 1/2 O 2 through proton irradiation. In addition to forming defects, proton irradiation at sufficiently high doses can subdivide single grains into multiple nanodomains without physically rupturing them. Hence, the single grains of the material assume a pseudo-secondary particle nature without reducing the overall grain size. Preserving the grain size is advantageous, as it reduces side reactions, which is not possible with conventional grain size reduction methods. While chemical transformations and defect formation induced through proton irradiation can influence the stability of battery materials, it is expected that structural reorganization due to cycling-induced phase transformation will be contained within these nanodomains. Such confinement of phase transformation is potentially responsible for enhancing the reversibility of layered oxide materials in our study. Thus, our study suggests that grain subdivision could become an effective microstructure tuning strategy for managing electrochemical cycling-induced phase changes in battery electrodes.

25 ENERGY STORAGE↗
DOE OSTI2025

Scale-up study of enhancing algal growth through bioaugmentation with the indole-3-acetic acid producing bacteria Azospirillum brasilense

Algae-based wastewater treatment technologies can recover nutrients while generating algal biomass with diverse application potentials. While numerous bacteria have demonstrated the ability to enhance algal growth through symbiotic relationships in laboratory settings, scaling this benefit to larger applications remains challenging due to complex wastewater conditions and microbe-algae dynamics. This study explores the effectiveness of bioaugmentation with Azospirillum brasilense, a plant growth-promoting bacterium known for producing the phytohormone indole-3-acetic acid (IAA), in boosting algal productivity within an Algaewheel wastewater treatment system. Through a series of experiments conducted at lab-, pilot-, and full-scale levels, we evaluated the impact of bioaugmentation on algal growth and developed a bioaugmentation strategy for an Algaewheel system. The results indicate that a biweekly dosing, achieving a final bioaugmenting cell density of 0.2–4 × 10 9 cells/L, effectively doubles the biomass productivity in a residential subdivision wastewater treatment plant in Northern Illinois. The bacterial production of IAA was the primary mechanism driving this enhancement, as demonstrated by the similar growth and yield improvements observed with chemically supplemented IAA and bioaugmentation across all experimental scales. High-throughput sequencing of 18S rRNA genes revealed that the algal community in the bioaugmented tank exhibited more stable biodiversity than the control. Additionally, bioaugmentation improved nutrient removal efficiency during winter, highlighting its potential to enhance the overall performance and sustainability of algae-based wastewater systems. The developed bioaugmentation and chemical treatment methods offer further operational solutions for managing yield, optimizing biochemical profiles, and enhancing biofuel production potential.

Algae-bacterial symbiotic relationship↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗