Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape function”

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

Transverse momentum dependent shape function for $J/ψ$ production in SIDIS

It has been shown previously that the transverse momentum dependent (TMD) factorization of heavy quarkonium production requires a TMD shape function. Its perturbative tail can be extracted by matching the cross sections valid at low and high transverse momenta. In this article we compare the order-$α$ s TMD expressions with the order-${α}_{s}^{2}$ collinear ones for $J/ψ$ production in semi-inclusive deep inelastic scattering (SIDIS), employing nonrelativistic QCD in both cases. In contrast to previous studies, we find that the small transverse momentum limit of the collinear expressions contains discontinuities. We demonstrate how to properly deal with them and include their finite contributions to the TMD shape functions. Moreover, we show that soft gluon emission from the low transverse momentum Born diagrams provide the same leading order TMD shape functions as required for the matching. Their revised perturbative tails have a less divergent behavior as compared to the TMD fragmentation functions of light hadrons. Finally, we investigate the universality of TMD shape functions in heavy quarkonium production, identify the need for process dependent factorization and discuss the phenomenological implications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Improved C 1 shape functions for simplex meshes

Using continuously-differentiable (C 1 ) interpolation methods previously developed for computer visualization and other applications, a new class of shape functions is constructed and applied to electrostatic (ES) Particle-In-Cell (PIC) simulations on simplex (triangle or tetrahedral) meshes. The resulting shape functions are non-negative and partition unity,and lead to a reduction of shot noise and aliasing associated with particle interpolation to and from the mesh. For an energy-conserving implementation as chosen, C 1 interpolation of the potential leads to a continuous (C 0 ) electric field, with associated improved fidelity of particle motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Establishing Models for Digital Twin of Hydropower Systems Using Probability Density Function Shaping

This paper introduces a digital twin modeling method for hydropower systems with Kaplan turbines using probability density function (PDF) shaping. We first use multilayer perceptron (MLP) model to build the discretized openloop Kaplan unit, where the MLP is trained by historical data. Then we use a proportional integral double derivative (PIDD) controller and a lead-lag exciter to test the obtained digital twin model in a closed-loop fashion. Simulation results show that the proposed digital twin modeling method can accurately capture the dynamics of the Kaplan hydropower unit. Finally, we show that the obtained digital twin can help to optimize the PIDD parameters. Compared with the original PIDD controller, the optimized one can achieve an over 90% improvement on the mean square tracking error.

Yin, Zhun [New York University]↗
DOE OSTI2022

Coral and Seawater Metagenomes Reveal Key Microbial Functions to Coral Health and Ecosystem Functioning Shaped at Reef Scale

The coral holobiont is comprised of a highly diverse microbial community that provides key services to corals such as protection against pathogens and nutrient cycling. The coral surface mucus layer (SML) microbiome is very sensitive to external changes, as it constitutes the direct interface between the coral host and the environment. Here, we investigate whether the bacterial taxonomic and functional profiles in the coral SML are shaped by the local reef zone and explore their role in coral health and ecosystem functioning. The analysis was conducted using metagenomes and metagenome-assembled genomes (MAGs) associated with the coral Pseudodiploria strigosa and the water column from two naturally distinct reef environments in Bermuda: inner patch reefs exposed to a fluctuating thermal regime and the more stable outer reefs. The microbial community structure in the coral SML varied according to the local environment, both at taxonomic and functional levels. The coral SML microbiome from inner reefs provides more gene functions that are involved in nutrient cycling (e.g., photosynthesis, phosphorus metabolism, sulfur assimilation) and those that are related to higher levels of microbial activity, competition, and stress response. In contrast, the coral SML microbiome from outer reefs contained genes indicative of a carbohydrate-rich mucus composition found in corals exposed to less stressful temperatures and showed high proportions of microbial gene functions that play a potential role in coral disease, such as degradation of lignin-derived compounds and sulfur oxidation. The fluctuating environment in the inner patch reefs of Bermuda could be driving a more beneficial coral SML microbiome, potentially increasing holobiont resilience to environmental changes and disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Voltage Probability Density Function Shaping Control Strategy Considering Grid Operational Uncertainties

It is well-known that power systems operation always affected by various uncertainties which make the bus voltage a random process that can be characterized by its probability density function (PDF) at any time instant. In this context, this paper presents a novel PDF-based voltage control framework for power systems. By modeling voltage as a stochastic process, we formulate a stochastic differential equationthat captures grid uncertainties. The associated Fokker-Planck-Kolmogorov equation is derived to describe the evolution of the voltage PDF, which enables the formulation of a PDF-shaping control strategy. To simplify the PDF control formulation, a B-spline neural network is introduced for real-time estimation and regulation of the voltage distribution. The proposed PDF control law updates voltage references for energy storage systems and synchronous generators using real-time PDF measurements and feedback signals. The proposed method is validated on a modified Kundur’s two-area system. Simulation results demonstrate that the controller can significantly improve the voltage stability under stochastic conditions, highlighting its effectiveness in modern inverter-rich grids.

Gui, Yonghao [ORNL] (ORCID:0000000250435534)↗
DOE OSTI2026

Nutrient limitation shapes functional traits of mycorrhizal fungi and phosphorus-cycling bacteria across an elevation gradient

In nutrient-limited high-elevation ecosystems, plants rely on arbuscular mycorrhizal (AM) fungi to provide mineral phosphorus (P) in the form of phosphate (PO43-). AM fungi gather these nutrients from phosphorus-cycling bacteria (PCBs) that can mineralize PO43- from organic matter and solubilize mineral-bound P. How climate, soil factors, and nutrient limitation influence AM fungi and PCB assembly remains unclear. We collected soil from montane meadows across a 1,000-m elevation gradient on three replicate mountainsides and analyzed AM fungal marker genes, P-cycling genes from shotgun metagenomes, and edaphic measurements. High-elevation soils had nearly 50-fold less soil PO₄³⁻ and 60% more AM fungal hyphae than low-elevation soils. AM fungal turnover was linked to changes in pH, organic carbon, and PO₄³-. The composition of 198 P-cycling genes was influenced by the AM fungal community structure. Drivers of individual PCB functional genes, including pH and organic carbon, varied with gene phylogeny. We found a trade-off in P-cycling strategies across elevation: P-rich, low-elevation soils supported root-colonizing AM fungi and organic P-mineralizing bacteria. P-poor, high-elevation soils were dominated by stress-tolerant AM fungi and mineral P-solubilizing bacteria. Our results suggest that AM fungi and PCB community turnover across elevation are both shaped by pH, organic carbon, and P availability. With continued climate warming, the structure and function of mountaintop ecosystems might shift to resemble lower elevations, disrupting long-established and specialized microbial assemblages, with consequences for P-cycling dynamics and the total P available to plant communities.IMPORTANCEPhosphorus (P) limits plant productivity in high-elevation ecosystems, yet the microbial networks that mobilize P, including arbuscular mycorrhizal (AM) fungi and phosphorus-cycling bacteria (PCBs), remain under-characterized in these nutrient-poor soils. We show that across a 10,00-m elevation gradient, AM fungi and P-cycling gene assemblages shift predictably with pH, organic carbon, and phosphate availability. Higher elevations, with less available P, select for stress-tolerant AM fungal taxa and PCB strategies geared toward mineral solubilization, while low-elevation sites favor root colonization by AM fungi and organic P mineralization. These results suggest that nutrient limitation can constrain microbial community assembly in consistent ways across landscapes. High mountain soils are low in P and rely on a network of underground AM fungi and PCB to deliver nutrients to plants. This study shows how those underground relationships reorganize with elevation and how climate change could collapse long-standing microbial strategies by pushing high-elevation ecosystems toward lowland conditions. As soils warm and dry, the microbial scaffolding that supports alpine plant life may become increasingly unstable.

arbuscular mycorrhizal fungi↗
DOE OSTI2026

Mapping the soil microbiome functions shaping wetland methane emissions

Accounting for only 8% of Earth’s land cover, freshwater wetlands remain the foremost contributors to global methane emissions. Yet the microorganisms and processes underlying methane emissions from wetland soils remain poorly understood. Over a five-year period, we surveyed the microbial membership and in situ methane measurements from over 700 samples in one of the most prolific methane-emitting wetlands in the United States. We constructed a catalog of 2,502 metagenome-assembled genomes (MAGs), with more than half of the 70 bacterial and archaeal phyla sampled containing novel lineages. Integration of these data with 133 soil metatranscriptomes provided a genome-resolved view of the biogeochemical specialization and versatility expressed over wetland soil spatial and temporal gradients. Centimeter-scale depth differences best explained patterns of microbial community structure and transcribed functionalities, even more than land cover or temporal information. Moreover, while extended flooding restructured soil redox, this perturbation failed to reconfigure the transcriptional profiles of methane-cycling microorganisms, contrasting with theoretically expected responses to hydrological perturbations. Co-expression analyses, coupled with depth-resolved methane measurements, revealed the metabolisms and trophic structures most predictive of methane hotspots. Mapping the spatiotemporal transcriptional patterns on this compendium of biogeochemically classified soil-derived genomes begins to untangle the microbial carbon, energy, and nutrient processing contributing to wetland methane production.

MAG↗
DOE OSTI2026

Topography and functional traits shape the distribution of key shrub plant functional types in low-Arctic tundra

The expansion of shrubs in the Arctic tundra fundamentally modifies land-atmosphere interactions. However, it remains unclear how shrub distribution and expansion differ across key species due to challenges with discriminating tundra plant species at regional scales. Here, we combined multi-scale, multi-platform remote sensing and in situ trait measurements to elucidate the distribution patterns and primary controls of two representative deciduous-tall-shrub (DTS) genera, Alnus and Salix, in low-Arctic tundra. We show that topographic features were a key control on DTSs, creating heterogeneous, but predictable distributions of Alnus and Salix fractional cover (fCover). Alnus was more tolerant of elevation and slope and was found on hilly uplands (slope >10°) within a specific elevational band (200–400 m above sea level [MSL]). In contrast, Salix occurred at lower elevations (50–300 m MSL) on gentler slopes (3-10°) and required adequate soil moisture associated with its profligate water use. We also show that niche differentiation between Alnus and Salix changed with patch size, where larger patches were more specialized in resource requirements than individual plants of Alnus and Salix. To understand what constrains the growth of DTSs at locations with low fCover, we developed environmental limiting factor models, which showed that topography limits the upper bound of Alnus and Salix fCover in 69.2% and 48.7% of the landscape, respectively. These findings highlight a critical need to better understand and represent topography-controlled processes and functional traits in regulating shrub distribution, as well as a need for more detailed species classification to predict shrubification in the Arctic.

alder↗
DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗
DOE OSTI2021

Independent normalization for γ -ray strength functions: The shape method

Here, the shape method, a novel approach to obtain the functional form of the γ-ray strength function (γSF), is introduced. In connection with the Oslo method the slope of the nuclear level density (NLD) and γSF can be obtained simultaneously even in the absence of neutron resonance spacing data. The foundation of the shape method lies in the primary γ-ray transitions which preserve information on the functional form of the γSF. The shape method has been applied to 56 Fe, 92 Zr, and 164 Dy, which are representative cases for the variety of situations encountered in typical NLD and γSF studies. The comparisons of results from the shape method to those from the Oslo method demonstrate that the functional form of the γSF is retained regardless of nuclear structure details or J π values of the states fed by the primary transitions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗
DOE OSTI2022

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗
DOE OSTI2022

A deformation-dependent coupled Lagrangian/semi-Lagrangian meshfree hydromechanical formulation for landslide modeling

Abstract The numerical modelling of natural disasters such as landslides presents several challenges for conventional mesh-based methods such as the finite element method (FEM) due to the presence of numerically challenging phenomena such as severe material deformation and fragmentation. In contrast, meshfree methods such as the reproducing kernel particle method (RKPM) possess unique features conducive to modelling extreme events such as the absence of a structured mesh and the ease of adaptive refinement, among others. While the semi-Lagrangian reproducing kernel (SL-RK) shape functions of RKPM defined in the current configuration have proven to be effective in extreme event modelling, the computational cost for the re-evaluation of the shape functions at every time step is costly. In this work, a deformation-dependent coupling of the Lagrangian reproducing kernel (L-RK) and SL-RK approximations is proposed for the solution of a hydro-mechanical formulation for effective simulations of landslides. The ramp function is constructed based on an equivalent plastic strain as a deformation-dependent transition from L-RK shape functions to SL-RK ones as the deformation progresses. The particular focus of the paper will be on modelling seepage-induced landslides with a mixed $$u$$ u – $$p$$ p formulation to couple the solid and fluid phases. Examples are presented to examine the effectiveness of this coupled Lagrangian/semi-Lagrangian reproducing kernel (L–SL RK) formulation and to highlight its performance in landslide modelling.

Baek, Jonghyuk↗
DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗