SEARCH · Engineering Papers

Results for “adaptive methods”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Understanding the evolution of double perovskite band structure upon dimensional reduction

Recent investigations into the effects of dimensional reduction on halide double perovskites have revealed an intriguing change in band structure when the three-dimensional (3D) perovskite is reduced to a two-dimensional (2D) perovskite with inorganic sheets of monolayer thickness (n = 1). The indirect bandgap of 3D Cs 2 AgBiBr 6 becomes direct in the n = 1 perovskite whereas the direct bandgap of 3D Cs 2 AgTlBr 6 becomes indirect at the n = 1 limit. Here, we apply a linear combination of atomic orbitals approach to uncover the orbital basis for this bandgap symmetry transition with dimensional reduction. We adapt our previously established method for predicting band structures of 3D double perovskites for application to their 2D congeners, emphasizing new considerations required for the 2D lattice. In particular, we consider the inequivalence of the terminal and bridging halides and the consequences of applying translational symmetry only along two dimensions. The valence and conduction bands of the layered perovskites can be derived from symmetry adapted linear combinations of halide p orbitals propagated across the 2D lattice. The dispersion of each band is then determined by the bonding and antibonding interactions of the metal and halide orbitals, thus affording predictions of the essential features of the band structure. We demonstrate this analysis for 2D Ag–Bi and Ag–Tl perovskites with sheets of mono- and bilayer thickness, establishing a detailed understanding of their band structures, which enables us to identify the key factors that drive the bandgap symmetry transitions observed at the n = 1 limit. Importantly, these insights also allow us to make the general prediction that direct → indirect or indirect → direct bandgap transitions in the monolayer limit are most likely in double perovskite compositions that involve participation of metal d orbitals at the band edges or that have no metal-orbital contributions to the valence band, laying the groundwork for the targeted realization of this phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Machine learning for international freight transportation management: A comprehensive review

Machine learning (ML) offers a promising avenue for international freight transportation management (IFTM) given its capability to harness the power of data that have become increasingly available to freight transportation researchers and practitioners. This paper conducts a comprehensive investigation of the state-of-the-art in developing ML models for applications to different aspects of IFTM. We start by giving an overview of various fundamental ML methods. Then, how different ML methods have been employed, adapted, and applied to a multitude of subject areas in IFTM are discussed, including demand forecast, operation and asset maintenance, and vehicle trajectory and on-time performance prediction. The potential data sources that may be used to develop ML models are further examined. Subsequently, a synthesis of the exiting work is performed to identify the specific topics addressed in the existing research, ML methods used, the trends of research, and opportunities for further explorations. Four directions for future research are proposed in the end.

Data sources↗

DOE OSTI2023

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗

DOE OSTI2024

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Hydrogen Detection Strategies to Support H2@SCALE - The NREL Sensor Laboratory

Hydrogen represents a major pathway to decarbonize and stabilize the national and international energy industry and select manufacturing markets. To facilitate the development of hydrogen markets, the US Department of Energy initiated H2@Scale to bring together stakeholders to advance affordable hydrogen production, transport, storage, and utilization to increase revenue opportunities in multiple energy sectors. One major impediment to hydrogen implementation is cost. To expedite the use of hydrogen in energy and other markets, the United States announced in 2021 the Hydrogen Shot, which seeks to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). As the cost of hydrogen drops, new applications will emerge that will require unique configurations of existing equipment and infrastructure, and eventually lead to advances in the generation and utilization of hydrogen. As the hydrogen economy expands, sensors and detection methods will need to adapt to changing infrastructure demands to address the primary targets of health & safety, emissions monitoring, and process control. The NREL Sensor Laboratory is playing a pivotal role in advancing the use of hydrogen sensors and detection methodologies in each of these categories to support DOE's mission for safe and efficient utilization in emerging markets. Health & safety monitors are required to ensure that operators and facilities can react to unintended hydrogen releases, either as GH2, LH2, or as a constituent of blends (e.g., natural gas or ammonia). Current detection methodologies focus on safety applications to detect near its lower flammable limit (4 vol %), and typically include point sensors in applications such as fixed or mobile detectors (e.g., personal gas monitors). Methodologies amenable for area detection include acoustic, emerging optical imaging methods, and flame detectors. Comparable detection strategies can be utilized for emissions monitoring and quantization, however few methods can simultaneously cover both low (emissions) and high (health & safety) levels. Deployment of emission level detectors will be required to 1) reduce product loss through small but potentially significant leaks from an environmental or cost perspective, 2) reduce downtime of high demand systems by early identification of eminent system failures (leaks through pump or compressor seals indicative of impending failure), and 3) address potential emission monitoring requirements that may be set by regulating bodies. The first two points should be adopted by industry to reduce the cost-of-goods-sold. The third main category for hydrogen detection relates to process control and may be advantageous for many existing applications. Two main applications are emerging. For example, the purity requirements for hydrogen that is dispensed from refueling systems for hydrogen fuel cell electric vehicles (FCEV) is rigorously regulated by the Standard SAE J2719, which prescribes maximum allowable levels of multiple impurities in the hydrogen fuel and must be verified by a regulatory body. Hydrogen contaminant detectors (HCD) integrated to the fueling station can assure this compliance. HCDs must be able operate in 100% H2 backgrounds and be able to distinguish between multiple contaminants at low ppm to low ppb levels. Secondly, as a strategy to decarbonize the natural gas grid, there are proposals to blend hydrogen with natural gas. This blending will affect transport applications (pipeline infrastructure), stationary combustion systems (turbines), and consumer and commercial appliances. In the short-term, hydrogen levels up to 20% are proposed. Variations in the hydrogen level can have dramatic impact on the combustion process and on the potential response of safety sensors. These mixtures may be regulated so that the concentration at a delivery point must be monitored with high precision. However, routine maintenance may introduce background gases such as ambient air (with water) or maintenance gases (introduced with welding processes or adhesive outgassing.) Therefore, the detection methodology must be robust enough to recover or respond to various contaminants. Several reviews can be found in literature addressing sensing and detection technologies, including their limitations and applications. However, for most applications, limitations can be alleviated by combining various detection techniques either through system integration or implementation of machine learning methods (artificial intelligence). In this presentation, we will discuss several applications, highlight their current approach for hydrogen detection, and suggest detection strategies to supplement their limitations.

ENERGY STORAGE,HYDROGEN↗

DOE OSTI2022

A probabilistic inverse prediction method for predicting plutonium processing conditions

In the past decade, nuclear chemists and physicists have been conducting studies to investigate the signatures associated with the production of special nuclear material (SNM). In particular, these studies aim to determine how various processing parameters impact the physical, chemical, and morphological properties of the resulting special nuclear material. By better understanding how these properties relate to the processing parameters, scientists can better contribute to nuclear forensics investigations by quantifying their results and ultimately shortening the forensic timeline. This paper aims to statistically analyze and quantify the relationships that exist between the processing conditions used in these experiments and the various properties of the nuclear end-product by invoking inverse methods. In particular, these methods make use of Bayesian Adaptive Spline Surface models in conjunction with Bayesian model calibration techniques to probabilistically determine processing conditions as an inverse function of morphological characteristics. Not only does the model presented in this paper allow for providing point estimates of a sample of special nuclear material, but it also incorporates uncertainty into these predictions. This model proves sufficient for predicting processing conditions within a standard deviation of the observed processing conditions, on average, provides a solid foundation for future work in predicting processing conditions of particles of special nuclear material using only their observed morphological characteristics, and is generalizable to the field of chemometrics for applicability across different materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

High-throughput, dynamic, multi-dimensional: an expanding repertoire of plant respiration measurements

A recent burst of technological innovation and adaptation has greatly improved our ability to capture respiration rate data from plant sources. At the tissue level, several independent respiration measurement options are now available, each with distinct advantages and suitability, including high-throughput sampling capacity. These advancements facilitate the inclusion of respiration rate data into large-scale biological studies such as genetic screens, ecological surveys, crop breeding trials, and multi-omics molecular studies. As a result, our understanding of the correlations of respiration with other biological and biochemical measurements is rapidly increasing. Difficult questions persist concerning the interpretation and utilization of respiration data; concepts such as allocation of respiration to growth versus maintenance, the unnecessary or inefficient use of carbon and energy by respiration, and predictions of future respiration rates in response to environmental change are all insufficiently grounded in empirical data. However, we emphasize that new experimental designs involving novel combinations of respiration rate data with other measurements will flesh-out our current theories of respiration. Furthermore, dynamic recordings of respiration rate, which have long been used at the scale of mitochondria, are increasingly being used at larger scales of size and time to reflect processes of cellular signal transduction and physiological response to the environment. Herein, we also highlight how respiratory methods are being better adapted to different plant tissues including roots and seeds, which have been somewhat neglected historically.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Probabilistic Assessment of High-Renewables Power Systems: Current Work and Future Directions

In a future with higher penetrations of variable and energy-limited resources, the dominant sources of resource adequacy risk begin to shift from (assumed) independent thermal unit outages to spatially- and temporally-correlated weather-driven phenomena. Traditional resource adequacy assessment and planning methods can be ill-suited to such scenarios. This presentation will highlight aspects of NREL's work in adapting and developing new probabilistic methods and tools in anticipation of the needs of the power system planning exercises of tomorrow.

17 WIND ENERGY↗

DOE OSTI2021

Quantifying Properties for a Mechanistic, Predictive Understanding of Aqueous Impact on Aging of Medium and Low Voltage AC and DC Cabling in Nuclear Power Plants

To address the gap in knowledge in understanding degradation in relevant service conditions, this project aimed to develop a mechanistic, predictive model of medium and low voltage cable failure based on the primary environmental degradation parameters of aqueous immersion time, temperature, and the oxidation extent. To do so, we took a two-fold approach toward evaluating degradation as related to the aqueous condition. First, we evaluated the chemical, mechanical, and electrical properties of polymers that comprise the cabling insulation under varied aqueous conditions at different temperatures. Secondly, we evaluated the performance of insulation materials under similar accelerated aging conditions in addition to the dielectric breakdown of the cabling under submersion conditions. Thermal oxidation and immersion of LDPE thin films were done using a Parr vessel and the extent of oxidation was monitored using the carbonyl index from ATR-FTIR spectrum, a ratio quantifying the carbonyl functionality produced on the LDPE film surface. Increasing the temperature, oxidation time, and oxygen pressure increased the water-vapor permeability and the carbonyl content. Thermogravimetric analysis of the films confirmed that there is indeed an initial weight loss of the smaller molecular weight molecules and volatiles. At higher temperatures, there is a secondary mass loss. Up to 70°C, the mass loss curves were similar until 80°C, when the aged films decreased at a much greater rate with respect to increasing temperature. This was seen in the 10, 50, 90% weight loss temperatures. The mechanical properties of sheet PP and HDPE as well as injection molded LDPE, HDPE, HDPP dog bones were plotted with respect to the time spent in the Parr vessel during thermal oxidation. For the sheet PP and HDPE, the UTS and the modulus of elasticity decreased while the elongation increased. The injection-molded LDPE, HDPE, and HDPP dog bones experienced a similar trend where the changes in the properties were within experimental uncertainty. This analysis of LDPE in dry and immersive oxygen-rich environments is an important baseline for future experiments looking at other degradation mechanisms. A predictive aging model for low-density polyethylene insulative cable housings was developed. This model revealed a combination of physical and chemical processes which lead to an increase in permeability and a decrease in strength of the insulator. In creating this predictive model, a gap of knowledge in the literature was revealed in two areas: understanding how the crystallinity of a polymer changes with time and visualizing the predicted pores formed through the insulator which leads to failure. Developing and using an adapted ATR-FTIR method, crystallinity was monitored and compared to bulk crystallinity found via DSC. Inhomogeneity in crystallinity changes were observed but the FTIR method but was found to be less reliable. Pore visualization was found utilizing electrical impedance spectroscopy saturated aged polyethylene films with synthesized citrate-capped gold nanoparticles, and chloroauric acid precursor. results indicate a decreasing impedance of the polyethylene films caused by the transport of ions through the film. SEM imaging was then performed for elemental analysis of the films and counter electrodes for the presence of ions. Chloride and gold ions were detected; however, no nanoparticles were found. This gives an estimated pore size of at least 0.3 nm through the aged film. Cyclic submergence of polyethylene and polypropylene in aqueous solutions of copper sulfate and Harrison’s solution were examined. It was discovered that aging in these mixed conditions and at 90°C did cause a small, yet significant increase to tensile strength for both unaged polyethylene and polypropylene. The PE and PP in cycled and submerged conditions did not have tensile strengths significantly different from dry-aged after 16 weeks. For both solutions, it was determined that capacitance increases with both water tree depth and cable temperature and is relatively independent of changing water tree AR. As for resistance, there is no apparent change between depth percentages of 10 and 80 for both solutions. This is because the water tree has not yet entered the conductor and thus there is no shorting current flowing through the tree yet. Between 80 and 100% water tree depth, it is observed that for both solutions there is a rapid decrease in resistance because the tree is now impacting the conductor shield and allowing shorting current to flow through it. The relationship between resistance and temperature was found to be inversed. As temperature increases, resistance decreases for both water tree solutions with distilled water having the most apparent difference across the range of temperatures simulated. With regards to water tree AR, resistance was found to be relatively independent at depth percentages of 10, 50, and 70, but there was some relationship at 100% depth for both solutions. At this depth, as water tree AR increased, resistance was found to also increase. Through the examination of voltage and electric field distribution plots, it was confirmed for EPR cables that increasing water tree depth leads to increasing distortion. The increase in shorting current at 100% depth for aqueous copper sulfate compared to distilled water was also visualized. Lastly, the relationship between shorting current and temperature for narrow water trees was visualized and validated using the distribution plots. This study has led to a better understanding of the effect of cable temperature, water tree solution, and geometry of the tree region on cable degradation. Furthermore, this study has documented simulation results of water tree degradation in EPR MV cables, which had been previously lacking in information. Future work in this area will add frequency into the mix and examine what effect it has on the rate of cable degradation as a result of water treeing. To compliment the work being carried by UMD on different polymer chemistries, harvested medium voltage (MV) cables were selected for accelerated aging study at ORNL. These medium voltage cables are representative of the vintages and materials that are currently in use at existing nuclear power plants (NPPs). These cables were immersed in water at 90°C and energized for a period of two years at elevated voltage. Partial discharge was measured during this period of time to track potential degradation. Unfortunately, due to the absence of failure, the extent of integration between the UMD techniques and the harvested MV insulation was limited. However, based on the findings and from the UMD research in the previous sections on polyethylene and polypropylene, follow-on characterization for harvested MV cable insulation should focus on sample preparation to take advantage of the UMD techniques to better understand mechanistic degradation in MV cable insulations. This would include: 1.) Utilization of solutions and impedance measurement insulation to track permeability in harvested and accelerated aging insulation samples, 2.) Adaptation of gold nanoparticle synthesis for electrical impedance spectroscopy with supporting SEM characterization to study pore formation in harvested and accelerated aged insulation samples, & 3.) ATR-FTIR crystallinity with multiple temperature dependent harvested MV cable insulation under different air and submerged accelerated aging. In addition, low voltage dielectric spectroscopy and high voltage dissipation factor, tan δ, could be implemented as a condition monitoring tool for harvested MV cable insulation in submerged environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Analysis of Hydrocyclone Geometry via Rapid Optimization Based on Computational Fluid Dynamics

Abstract Hydrocyclones exploit density gradients for the centrifugal separation of dispersions in a continuous liquid. Selection of the geometrics for optimal separation is case specific, like the media characteristics. The existing optimization method based on computational fluid dynamics (CFD) provides a powerful analytical tool but requires long computational times. The most common praxis for CFD optimization is via the single‐factor optimization method (SFOM). In this study, a novel approach is presented as an improved rapid optimization method that implements a dynamic‐mesh and user‐defined function optimization method (DUOM). The DUOM adapts the dynamic‐mesh approach from other applications to the optimization analysis of hydrocyclones. The DUOM reduced the computational time by 31.1 %, compared to the SFOM.

Liu, Lin↗

DOE OSTI2022

Adaptive laboratory evolution for increased temperature tolerance of the diatom Nitzschia inconspicua

Abstract Outdoor microalgal cultivation for the production of valuable biofuels and bioproducts typically requires high insolation and strains with high thermal (>37°C) tolerance. While some strains are naturally thermotolerant, other strains of interest require improved performance at elevated temperatures to enhance industrial viability. In this study, adaptive laboratory evolution (ALE) was performed for over 300 days using consecutive 0.5°C temperature increases in a constant temperature incubator to attain greater thermal tolerance in the industrially relevant diatom Nitzschia inconspicua str. Hildebrandi. The adapted strain was able to grow at a constant temperature of 37.5°C; whereas this constant temperature was lethal to the parental control, which had an upper‐temperature boundary of 35.5°C before adaptive evolution. Several high‐temperature clonal isolates were obtained from the evolved population following ALE, and increased temperature tolerance was observed in the clonal, parent, and non‐clonal adapted cultures. This ALE method demonstrates the development of enhanced industrial algal strains without the production of genetically modified organisms (GMOs).

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

High-Resolution Imaging of Cellulose Organization in Cell Walls by Field Emission Scanning Electron Microscopy

Field emission scanning electron microscopy (FESEM) is a powerful tool for analyzing surface structures of biological and nonbiological samples. However, when it is used to study fine structures of nanometer-sized microfibrils of epidermal cell walls, one often encounters tremendous challenges to acquire clear and undistorted images because of two major issues: (1) Preparation of samples suitable for high resolution imaging; due to the delicateness of some plant materials, such as onion epidermal cell walls, many things can happen during sample processing, which subsequently result in damaged samples or introduce artifacts. (2) Difficulties to acquire clear images of samples which are electron-beam sensitive and prone to charging artifacts at magnifications over 100,000×. In this chapter we described detailed procedures for sample preparation and conditions for high-resolution FESEM imaging of onion epidermal cell walls. Lastly, the methods can be readily adapted for other wall materials.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Amplitude/operator basis in chiral perturbation theory

We establish a systematic construction of the on-shell amplitude/operator basis for Chiral Perturbation Theory (ChPT) in D = 4 spacetime dimensions and with an arbitrary number of flavors N f . For kinematic factors, we employ spinor-helicity variables to construct the soft blocks, which are local amplitudes satisfying the Adler’s zero condition, as well as to take into account the reduction in the kinematic basis due to the Gram determinant, which arises at O(p 10 ) when the number of multiplicity N in an amplitude becomes large: N > D. For flavor factors, we include group-theoretic relations at small N f , N f ≤ N, which decreases the flavor basis. The result is obtained by adapting the Young tensor method of constructing the operator basis for generic effective field theories to the case of non-linearly realized symmetries. Working in the massless quark limit, we present purely mesonic operators for both even- and odd-parity at O(p 6 ) and O(p 8 ) for N = 6 and arbitrary N f , and establish a direct correspondence between the amplitude basis and the operator basis. Furthermore, the redundancy due to the Gram determinant is studied at O(p 10 ) for N = 6, 8 and 10.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

FIRM: federated image reconstruction using multimodal tomographic data

Here, we propose a federated algorithm for reconstructing images using multimodal tomographic data sourced from dispersed locations, addressing the challenges of traditional unimodal approaches that are prone to noise and reduced image quality, as well as the limitations of centralized multimodal approaches that require extensive data transfer, leading to significant communication overhead, storage demands, and potential data privacy concerns. Our approach formulates a joint inverse optimization problem incorporating multimodality constraints and solves it in a federated framework through local gradient computations complemented by lightweight central operations, thereby ensuring data decentralization. Leveraging the connection between our federated algorithm and the quadratic penalty method, we introduce an adaptive step-size rule with guaranteed sublinear convergence. Numerical results demonstrate superior computational efficiency and improved image reconstruction quality compared to existing approaches.

federated algorithm↗

DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗