Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tree decomposition”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗
DOE OSTI2021

The Stellar decomposition: A compact representation for simplicial complexes and beyond

Here, we introduce the Stellar decomposition, a model for efficient topological data structures over a broad range of simplicial and cell complexes. A Stellar decomposition of a complex is a collection of regions indexing the complex’s vertices and cells such that each region has sufficient information to locally reconstruct the star of its vertices, i.e., the cells incident in the region’s vertices. Stellar decompositions are general in that they can compactly represent and efficiently traverse arbitrary complexes with a manifold or non-manifold domain. They are scalable to complexes in high dimension and of large size, and they enable users to easily construct tailored application-dependent data structures using a fraction of the memory required by a corresponding global topological data structure on the complex. As a concrete realization of this model for spatially embedded complexes, we introduce the Stellar tree, which combines a nested spatial tree with a simple tuning parameter to control the number of vertices in a region. Stellar trees exploit the complex’s spatial locality by reordering vertex and cell indices according to the spatial decomposition and by compressing sequential ranges of indices. Stellar trees are competitive with state-of-the-art topological data structures for manifold simplicial complexes and offer significant improvements for cell complexes and non-manifold simplicial complexes. We conclude with a high-level description of several mesh processing and analysis applications that utilize Stellar trees to process large datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Coarse woody debris decomposition assessment tool: Model development and sensitivity analysis

Coarse woody debris (CWD) is an important component in forests, hosting a variety of organisms that have critical roles in nutrient cycling and carbon (C) storage. We developed a process-based model using literature, field observations, and expert knowledge to assess woody debris decomposition in forests and the movement of wood C into the soil and atmosphere. The sensitivity analysis was conducted against the primary ecological drivers (wood properties and ambient conditions) used as model inputs. The analysis used eighty-nine climate datasets from North America, from tropical (14.2° N) to boreal (65.0° N) zones, with large ranges in annual mean temperature (26.5°C in tropical to -11.8°C in boreal), annual precipitation (6,143 to 181 mm), annual snowfall (0 to 612 kg m -2 ), and altitude (3 to 2,824 m above mean see level). The sensitivity analysis showed that CWD decomposition was strongly affected by climate, geographical location and altitude, which together regulate the activity of both microbial and invertebrate wood-decomposers. CWD decomposition rate increased with increments in temperature and precipitation, but decreased with increases in latitude and altitude. CWD decomposition was also sensitive to wood size, density, position (standing vs downed), and tree species. The sensitivity analysis showed that fungi are the most important decomposers of woody debris, accounting for over 50% mass loss in nearly all climatic zones in North America. The model includes invertebrate decomposers, focusing mostly on termites, which can have an important role in CWD decomposition in tropical and some subtropical regions. The role of termites in woody debris decomposition varied widely, between 0 and 40%, from temperate areas to tropical regions. Woody debris decomposition rates simulated for eighty-nine locations in North America were within the published range of woody debris decomposition rates for regions in northern hemisphere from 1.6° N to 68.3° N and in Australia.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Soil organic matter, tree communities, and fungal communities across mycorrhizal gradients in the Eastern United States

We collected these data to investigate how mycorrhizal associations are related to soil C and N across four sites in the eastern U.S. broadleaf forest biome, which capture broad variability in climate and tree species. Our sites represented each of the four ecoregions in the eastern United States temperate forest—warm continental (New Hampshire), hot continental (Wisconsin), Prairie (Illinios), and subtropical (Georgia). This breadth naturally provided tree species diversity and allowed us to also investigate how soil C and N are related to canopy tree and EcM fungal community composition. Specifically, we analyzed the effects of mycorrhizal association, canopy tree family, and EcM fungal taxa on the proportion of C and N associated with soil minerals (MAOM) and bulk soil C:N. These are plot level data from forests in New Hampshire, Georgia, Wisconsin, and Illinois. The plots are 10 m in radius and vary across gradients of ectomycorrhizal tree basal area. Within each plot, data include the basal area of all trees >2 cm in diameter identified to species, the abundance of different fungal taxa based on ITS sequences, and the %C, %N, C:N of bulk soil and of soil density fractions (free light, occluded light, and particulate) as well as the proportion of soil C and N within each fraction. The soil data (both organic matter and fungal communities) are from the top 10 cm of the mineral soil. This also includes a meta-analysis of leaf litter k and C:N values for the tree species present in the forests.The "Tree_BA_ALL.csv" files include the basal area of all trees >2 cm in diameter at breast height within our 10 cm radii plots. The "Soil_Data_Final.csv" file contains the soil carbon and nitrogen in bulk and density fractions (concentrations as well as proportion for free light fraction (fLF), occluded light fraction (oLF), and heavy or dense fraction (HF), oxalate extractable iron and aluminum concentrations, as well as the geographic coordinates of each plot's location. The "FungiNonRelative.csv" file contains the abundances of fungal taxa identified via ITS sequencing for each plot (with the exception of 4 plots from GA where there was not enough quality DNA to extract).The "R_ReadyK.csv" and "R_Ready_CN.csv" files contain leaf litter decomposition constant (k) and C:N values from the literature for the tree species found at our sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

FleCSPH notes

FleCSI is a compile-time con gurable framework designed to support multi-physics application development. As such, FleCSI provides a very general set of infrastructure design patterns that can be specialized and extended to suit the needs of a broad variety of solver and data requirements. FleCSI currently supports multi-dimensional mesh topology, geometry, and adjacency information, as well as n-dimensional hashed-tree data structures, graph partitioning interfaces, and dependency closures. FleCSI introduces a functional programming model with control, execution, and data abstractions that are consistent both with MPI and with state-of-the-art, task-based runtimes such as Legion and Charm++. The abstraction layer insulates developers from the underlying runtime, while allowing support for multiple runtime systems including conventional models like asynchronous MPI. The intent is to provide developers with a concrete set of user-friendly programming tools that can be used now, while allowing exibility in choosing runtime implementations and optimization that can be applied to future architectures and runtimes. FleCSI's control and execution models provide formal nomenclature for describing poorly understood concepts such as kernels and tasks. FleCSI's data model provides a low-buy-in approach that makes it an attractive option for many application projects, as developers are not locked into particular layouts or data structure representations. FleCSI currently provides a parallel but not distributed implementation of Binary, Quad and Oct-tree topology. This implementation is base on space lling curves domain decomposition, the Morton order. The current FleCSI version requires the implementation of a driver and a specialization driver. The role of the specialization driver is to provide the data distribution. This feature is not complete in FleCSI code and we provide it. The next step will be to incorporate it directly from FleCSPH to FleCSI as we reach a good level of performance. Then the driver represent the general execution of the resolution without worrying of the data locality and communications. As FleCSI is an On-Development code the structure may change in the future and we keep track of these changes in FleCSPH.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Monte Carlo Tree Search for Integrated Planning, Learning, and Execution in Nondeterministic Python

We present a novel use of Monte Carlo Tree Search (MCTS),adapted to explore a search space produced by the choice points embedded in Python code. The choice points are non-deterministic assignment statements and subroutine calls. We present MCTS extensions required for doing tree search in this context which includes control constructs like hierarchical decomposition (subroutine calls), iterative while loops and conditional statements. We demonstrate how the system works in a simulated rideshare scenario in an urban setting, and present preliminary experiments as a proof of concept.

Automatic planning↗
DOE OSTI2020

A comparison of lodgepole and spruce needle chemistry impacts on terrestrial biogeochemical processes during isolated decomposition

This study investigates the isolated decomposition of spruce and lodgepole conifer needles to enhance our understanding of how needle litter impacts near-surface terrestrial biogeochemical processes. Harvested needles were exported to a subalpine meadow to enable a discrete analysis of the decomposition processes over 2 years. Initial chemistry revealed the lodgepole needles to be less recalcitrant with a lower carbon to nitrogen (C:N) ratio. Total C and N fundamentally shifted within needle species over time with decreased C:N ratios for spruce and increased ratios for lodgepole. Differences in chemistry correlated with CO 2 production and soil microbial communities. The most pronounced trends were associated with lodgepole needles in comparison to the spruce and needle-free controls. Increased organic carbon and nitrogen concentrations associated with needle presence in soil extractions further corroborate the results with clear biogeochemical signatures in association with needle chemistry. Interestingly, no clear differentiation was observed as a function of bark beetle impacted spruce needles vs those derived from healthy spruce trees despite initial differences in needle chemistry. These results reveal that the inherent chemistry associated with tree species has a greater impact on soil biogeochemical signatures during isolated needle decomposition. By extension, biogeochemical shifts associated with bark beetle infestation are likely driven more by changes such as the cessation of rhizospheric processes than by needle litter decomposition.

Bark Beetle Disturbance↗
DOE OSTI2020

All-multiplicity one-loop amplitudes in Born-Infeld electrodynamics from generalized unitarity

We initiate a study of non-supersymmetric Born-Infeld electrodynamics in 4d at the quantum level. Explicit all-multiplicity expressions are calculated for the purely rational one-loop amplitudes in the self-dual (+ + … +) and next-to-self-dual (- + … +) helicity sectors. Using a supersymmetric decomposition, d-dimensional unitarity cuts of the integrand factorize into tree-amplitudes in a 4d model of Born-Infeld photons coupled to a massive complex scalar. The two-scalar tree-amplitudes needed to construct the Born-Infeld integrand are computed using two complimentary approaches: (1) as a double-copy of Yang-Mills coupled to a massive adjoint scalar with a dimensionally reduced form of Chiral Perturbation Theory, and (2) by imposing consistency with low-energy theorems under a reduction from 4d to 3d and T-duality. The Born-Infeld integrand is integrated in d = 4 - 2ϵ dimensions at order Ο(ϵ 0 ) using the dimension-shifting formalism. We comment on the implications for electromagnetic duality in quantum Born-Infeld theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2003

Enhancing and Archiving the APS Catalog of the POSS I

We have worked on two different projects: 1) Archiving the APS Catalog of the POSS I for distribution to NASA's NED at IPAC, SIMBAD in France, and individual astronomers and 2) The automated morphological classification of galaxies. We have completed archiving the Catalog into easily readable binary files. The database together with the software to read it has been distributed on DVD's to the national and international data centers and to individual astronomers. The archived Catalog contains more than 89 million objects in 632 fields in the first epoch Palomar Observatory Sky Survey. Additional image parameters not available in the original on-line version are also included in the archived version. The archived Catalog is also available and can be queried at the APS web site (URL: http://aps.umn.edu) which has been improved with a much faster and more efficient querying system. The Catalog can be downloaded as binary datafiles with the source code for reading it. It is also being integrated into the SkyQuery system which includes the Sloan Digital Sky Survey, 2MASS, and the FIRST radio sky survey. We experimented with different classification algorithms to automate the morphological classification of galaxies. This is an especially difficult problem because there are not only a large number of attributes or parameters and measurement uncertainties, but also the added complication of human disagreement about the adopted types. To solve this problem we used 837 galaxy images from nine POSS I fields at the North Galactic Pole classified by two independent astronomers for which they agree on the morphological types. The initial goal was to separate the galaxies into the three broad classes relevant to issues of large scale structure and galaxy formation and evolution: early (ellipticals and lenticulars), spirals, and late (irregulars) with an accuracy or success rate that rivals the best astronomer classifiers. We also needed to identify a set of parameters derived from the digitized images that separate the galaxies by type. The human eye can easily recognize complicated patterns in images such as spiral arms which can be spotty, blotchy affairs that are difficult for automated techniques. A galaxy image can potentially be described by hundreds of parameters, all of which may have some relation to the morphological type. In the set of initial experiments we used 624 such parameters, in two colors, blue and red. These parameters include the surface brightness and color measured at different radii, ratios of these parameters at different radii, concentration indices, Fourier transforms and wavelet decomposition coefficients. We experimented with three different classes of classification algorithms; decision trees, k-nearest neighbors, and support vector machines (SVM). A range of experiments were conducted and we eventually narrowed the parameters to 23 selected parameters. SVM consistently outperformed the other algorithms with both sets of features. By combining the results from the different algorithms in a weighted scheme we achieved an overall classification success of 86%.

Humphreys, Roberta M.↗
NASA NTRS2023

Search for UnderUtilized Airspace for Extensible Traffic Management Operations Based on Air Traffic Patterns

This paper presents a new method to facilitate integrating new vehicle traffic operations, operating with existing air traffic operations. The extensible traffic management (xTM) concept assumes that the new vehicles can operate in a dedicated Cooperative Area (CA) with minimal interaction with conventional air traffic and requiring minimal air traffic supervision. Our method assumes that a new xTM CA can be created when an underutilized airspace with little or no traffic can be identified. Our approach involves modeling airspace as a tree data structure and iteratively subdividing it into smaller cells, with underutilized airspace defined as any cells without flight tracks. The benefits of our approach include its applicability to all xTM scenarios, the ability to handle both 2D and 3D space using a unique tree data structure, and computational efficiency for key functions such as space decomposition, labeling of connected cells, and searching of cells containing a given point. By automatically searching for underutilized airspace based on operating air traffic patterns, we can optimize airspace utilization and improve air traffic management. Our proposed approach can quantitatively determine when and where to allow xTM operations in the National Airspace System.

extensible traffic management↗
DOE OSTI2018

Leaf litter decay rates differ between mycorrhizal groups in temperate, but not tropical, forests

Whereas the primary controls on litter decomposition are well established, we lack a framework for predicting interspecific differences in litter decay within and across ecosystems. Given previous research linking tree mycorrhizal association with carbon and nutrient dynamics, we hypothesized that the two dominant mycorrhizal groups in forests – arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi – differ in litter decomposition rates. We compiled leaf litter chemistry and decay data for AM- and ECM-associating angiosperms and gymnosperms (> 200 species) from temperate and tropical/subtropical, and investigated relationships among decay rates, mycorrhizal association, phylogeny and climate. In temperate forests, AM litters decayed faster than ECM litters, with litter nitrogen and phylogeny best explaining variation in litter decay. In sub/tropical forests, we found no significant difference in litter decay rate between mycorrhizal groups, and variation in decay rates was best explained by litter phosphorus. Finally, our results suggest that knowledge of tree mycorrhizal association may improve predictions of species effects on ecosystem processes, particularly in temperate forests where AM and ECM species commonly co-occur, providing a predictive framework for linking litter quality, organic matter dynamics and nutrient acquisition in forests.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Proof of a new colour decomposition for QCD amplitudes

Recently, Johansson and Ochirov conjectured the form of a new colour decom-position for QCD tree-level amplitudes. This note provides a proof of that conjecture. The proof is based on ‘Mario World’ Feynman diagrams, which exhibit the hierarchical Dyck structure previously found to be very useful when dealing with multi-quark amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Caravel: A C++ framework for the computation of multi-loop amplitudes with numerical unitarity

We present the first public version of Caravel, a C++17 framework for the computation of multi-loop scattering amplitudes in quantum field theory, based on the numerical unitarity method. Caravel is composed of modules for the D-dimensional decomposition of integrands of scattering amplitudes into master and surface terms, the computation of tree-level amplitudes in floating point or finite-field arithmetic, the numerical computation of one- and two-loop amplitudes in QCD and Einstein gravity, and functional reconstruction tools. Here, we provide programs that showcase Caravel's main functionalities and allow to compute selected one- and two-loop amplitudes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

An approach to solving large reliability models

This paper describes a unified approach to the problem of solving large realistic reliability models. The methodology integrates behavioral decomposition, state trunction, and efficient sparse matrix-based numerical methods. The use of fault trees, together with ancillary information regarding dependencies to automatically generate the underlying Markov model state space is proposed. The effectiveness of this approach is illustrated by modeling a state-of-the-art flight control system and a multiprocessor system. Nonexponential distributions for times to failure of components are assumed in the latter example. The modeling tool used for most of this analysis is HARP (the Hybrid Automated Reliability Predictor).

Boyd, Mark A.↗
DOE OSTI2025

Extremely Scalable Distributed Computation of Contour Trees via Pre-Simplification

Contour trees offer an abstract representation of the level set topology in scalar fields and are widely used in topological data analysis and visualization. However, applying contour trees to large-scale scientific datasets remains challenging due to scalability limitations. Recent developments in distributed hierarchical contour trees have addressed these challenges by enabling scalable computation across distributed systems. Building on these structures, advanced analytical tasks—such as volumetric branch decomposition and contour extraction—have been introduced to facilitate large-scale scientific analysis. Despite these advancements, such analytical tasks substantially increase memory usage, which hampers scalability. In this paper, we propose a pre-simplification strategy to significantly reduce the memory overhead associated with analytical tasks on distributed hierarchical contour trees. We demonstrate enhanced scalability through strong scaling experiments, constructing the largest known contour tree—comprising over half a trillion nodes with complex topology—in under 15 minutes on a dataset containing 550 billion elements.

Li, Mingzhe [University of Utah]↗
DOE OSTI2023

An Infinite Domain 3D Poisson Solver Based on the Barnes-Hut Algorithm

We present a domain decomposition method for the solution of the 3D Poisson equation with infinite domain boundary conditions. The method is based on an application of the Barnes-Hut tree particle scheme adapted to gridded data. Long range interactions are computed using the first two terms in the Cartesian multipole expansion of Green’s function convoluted with the charge while short range computations are performed using Hockney’s domain doubling algorithm. A standard domain decomposition strategy requires O(N 2 ) applications of Hockney’s algorithm, where N is the number of subdomains that intersect that charge support, while in the present approach only O(Nlog 2 N) such computations suffice. The discretization scheme employed is a sixth order Mehrstellen approximation of the 3D Laplace opera tor. The method exhibits satisfactory accuracy at a substantially reduced computational cost compared to the full domain decomposition Hockney’s algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

SPRUCE Root Tip and Ectomycorrhizal Fungi Colonization Measurements from Ingrowth Cores, 2017

This data set contains root tip and ectomycorrhizal fungi colonization measurements taken from ingrowth cores from the SPRUCE experiment (Hanson et al. 2017) that were deployed during the 2017 growing season (2017-06 to 2017-10-01). This study explored the relationship between warming treatments and fine-root growth. Increased fine-root growth may increase root exudates and accelerate turnover, representing an underlying mechanism for peat decomposition through priming, as exudates provide a labile carbon source to the microbial community. Roots of two tree species were studied: an evergreen conifer Picea mariana (black spruce) and a deciduous conifer Larix laricina (tamarack). Measurements include root tips counts and densities by tree species and the abundance of ectomycorrhizal colonization on root tips. This dataset contains one data file in comma separate (.csv) format. Additional metadata are provided: one data dictionary and a file-level metadata file in comma separate (.csv) format and a user guide in PDF (*.pdf) format.

black spruce [Picea mariana]↗
DOE OSTI2021

Phenolic acid-degrading Paraburkholderia prime decomposition in forest soil

Abstract Plant-derived phenolic acids are catabolized by soil microorganisms whose activity may enhance the decomposition of soil organic carbon (SOC). We characterized whether phenolic acid-degrading bacteria enhance SOC mineralization in forest soils when primed with 13C-labeled p-hydroxybenzoic acid (pHB). We further tested whether pHB-induced priming could explain differences in SOC content among mono-specific tree plantations in a 70-year-old common garden experiment. pHB addition primed significant losses of SOC (3–13 µmols C g−1 dry wt soil over 7 days) compared to glucose, which reduced mineralization (-3 to -8 µmols C g−1 dry wt soil over 7 days). The principal degraders of pHB were Paraburkholderia and Caballeronia in all plantations regardless of tree species or soil type, with one predominant phylotype (RP11ASV) enriched 23-fold following peak pHB respiration. We isolated and confirmed the phenolic degrading activity of a strain matching this phylotype (RP11T), which encoded numerous oxidative enzymes, including secretion signal-bearing laccase, Dyp-type peroxidase and aryl-alcohol oxidase. Increased relative abundance of RP11ASV corresponded with higher pHB respiration and expression of pHB monooxygenase (pobA), which was inversely proportional to SOC content among plantations. pobA expression proved a responsive measure of priming activity. We found that stimulating phenolic-acid degrading bacteria can prime decomposition and that this activity, corresponding with differences in tree species, is a potential mechanism in SOC cycling in forests. Overall, this study highlights the ecology and function of Paraburkholderia whose associations with plant roots and capacity to degrade phenolics suggest a role for specialized bacteria in the priming effect.

54 ENVIRONMENTAL SCIENCES↗