Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Partitioning”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impacts of divalent cations (Mg 2+ and Ca 2+ ) on PFAS bioaccumulation in freshwater macroinvertebrates representing different foraging modes

Per- and polyfluoroalkyl substances (PFAS) have extensively contaminated freshwater aquatic ecosystems where they can be transported in water and partition to sediment and biota. In this paper, three freshwater benthic macroinvertebrates with different foraging modes were exposed to environmentally relevant concentrations of eight perfluoroalkyl carboxylates (PFCA), three perfluoroalkyl sulfonates (PFSA), and three fluorotelomer sulfonates (FTS) at varying divalent cation concentrations of magnesium (Mg 2+ ) and calcium (Ca 2+ ). Divalent cations can impact PFAS partitioning to solids, especially to sediments, at higher concentrations. Sediment dwelling worms (Lumbriculus variegatus), epibenthic grazing snails (Physella acuta), and sediment-dwelling filter-feeding bivalves (Elliptio complanata) were selected due to their unique foraging modes. Microcosms were composed of synthetic sediment, culture water, macroinvertebrates, and PFAS and consisted of a 28-day exposure period. L. variegatus had significantly higher PFAS bioaccumulation than P. acuta and E. complanata, likely due to higher levels of interactions with and ingestion of the contaminated sediment. “High Mg 2+ ” (7.5 mM Mg 2+ ) and “High Ca 2+ ” (7.5 mM Ca 2+ ) conditions generally had statistically higher bioaccumulation factors (BAF) than the “Reference Condition” (0.2 mM Ca 2+ and 0.2 mM Mg 2+ ) for PFAS with perfluorinated chain lengths greater than eight carbons. Long-chain PFAS dominated the PFAS profiles of the macroinvertebrates for all groups of compounds studied (PFCA, PFSA, and FTS). Furthermore, these results indicate that the study organism has the greatest impact on bioaccumulation, although divalent cation concentration had observable impacts between organisms depending on the environmental conditions. Elevated cation concentrations in the microcosms led to significantly greater bioaccumulation in the test organisms compared to the experimental reference conditions for long-chain PFAS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Can activated sludge treatments and advanced oxidation processes remove organophosphorus flame retardants?

This study aims to determine the occurrence of 10 OPFRs (including chlorinated, nonchlorinated alkyl and aryl compounds) in influent, effluent wastewaters and partitioning into sludge of 5 wastewater treatment plants (WWTP) in Catalonia (Spain). All target OPFRs were detected in the WWTPs influents, and the total concentration ranged from 3.67 µg L{sup −1} to 150 µg L{sup −1}. During activated sludge treatment, most OPFRs were accumulated in the sludge at concentrations from 35.3 to 9980 ng g{sup −1} dw. Chlorinated compounds tris(2-chloroethyl) phosphate (TCEP), tris(2-chloroisopropyl) phosphate (TCIPP) and tris(2,3-dichloropropyl) phosphate (TDCPP) were not removed by the conventional activated sludge treatment and they were released by the effluents at approximately the same inlet concentration. On the contrary, aryl compounds tris(methylphenyl) phosphate (TMPP) and 2-ethylhexyl diphenyl phosphate (EHDP) together with alkyl tris(2-ethylhexyl) phosphate (TEHP) were not detected in any of the effluents. Advanced oxidation processes (UV/H{sub 2}O{sub 2} and O{sub 3}) were applied to investigate the degradability of recalcitrant OPFRs in WWTP effluents. Those detected in the effluent sample (TCEP, TCIPP, TDCPP, tributyl phosphate (TNBP), tri-iso-butyl phosphate (TIBP) and tris(2-butoxyethyl) phosphate (TBOEP)) had very low direct UV-C photolysis rates. TBOEP, TNBP and TIBP were degraded by UV/H{sub 2}O{sub 2} and O{sub 3}. Chlorinated compounds TCEP, TDCPP and TCIPP were the most recalcitrant OPFR to the advanced oxidation processes applied. The study provides information on the partitioning and degradability pathways of OPFR within conventional activated sludge WWTPs. - Highlights: • OPFRs were detected in wastewater and sludge of all studied WWTPs. • Alkyl and chloroalkyl phosphates were present in secondary treatment effluents. • TBOEP, TNBP and TIBP were degraded by UV/H{sub 2}O{sub 2} and O{sub 3} treatment. • TCEP, TCIPP and TDCPP were resistant to both secondary and tertiary treatment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Oral relative bioavailability of Dichlorodiphenyltrichloroethane (DDT) in contaminated soil and its prediction using in vitro strategies for exposure refinement

In this study, the bioavailability of DDTr (sum of DDT, DDD and DDE isomers) in pesticide-contaminated soil was assessed using an in vivo mouse model. DDTr relative bioavailability (RBA) ranged from 18.7±0.9 (As35) to 60.8±7.8% (As36) indicating that a significant portion of soil-bound DDTr was not available for absorption following ingestion. When DDTr bioaccessibility was assessed using the organic Physiologically Based Extraction Test (org-PBET), the inclusion of a sorption sink (silicone cord) enhanced DDTr desorption by up to 20-fold (1.6–3.8% versus 18.9–56.3%) compared to DDTr partitioning into gastrointestinal fluid alone. Enhanced desorption occurred as a result of the silicone cord acting as a reservoir for solubilized DDTr to partition into, thereby creating a flux for further desorption until equilibrium was achieved. When the relationship between in vivo and in vitro data was assessed, a strong correlation was observed between the mouse bioassay and the org-PBET+silicone cord (slope=0.94, y-intercept=3.5, r{sup 2}=0.72) suggesting that the in vitro approach may provide a robust surrogate measure for the prediction of DDTr RBA in contaminated soil. - Highlights: • An optimised mouse assay was used to quantify DDTr relative bioavailability in soil. • DDTr bioaccessibility was also determined using an in vitro sorption sink approach. • A strong correlation was observed between in vivo and in vitro data. • The sorption sink approach may be used to predict DDTr relative bioavailability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Quantification and mapping of urban fluxes under climate change: Application of WRF-SUEWS model to Greater Porto area (Portugal)

Climate change and the growth of urban populations are two of the main challenges facing Europe today. These issues are linked as climate change results in serious challenges for cities. Recent attention has focused on how urban surface-atmosphere exchanges of heat and water will be affected by climate change and the implications for urban planning and sustainability. In this study energy fluxes for Greater Porto area, Portugal, were estimated and the influence of the projected climate change evaluated. To accomplish this, the Weather Research and Forecasting Model (WRF) and the Surface Urban Energy and Water Balance Scheme (SUEWS) were applied for two climatological scenarios: a present (or reference, 1986–2005) scenario and a future scenario (2046–2065), in this case the Representative Concentration Pathway RCP8.5, which reflects the worst set of expectations (with the most onerous impacts). The results show that for the future climate conditions, the incoming shortwave radiation will increase by around 10%, the sensible heat flux around 40% and the net storage heat flux around 35%. In contrast, the latent heat flux will decrease about 20%. The changes in the magnitude of the different fluxes result in an increase of the net all-wave radiation by 15%. The implications of the changes of the energy balance on the meteorological variables are discussed, particularly in terms of temperature and precipitation. - Highlights: • Assessment of energy fluxes behaviour under past period and medium-term climate change projection. • Evaluation of climate change at urban scale. • Meteorological variables alters the partitioning of the energy fluxes. • Changes in the partition of the annual energy balance are found between the two analysed periods. • Increase in the magnitude of sensible and storage heat fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Earthquake-triggered submarine canyon flushing transfers young terrestrial and marine organic carbon into the deep sea

Submarine canyons transfer substantial amounts of sediment and organic carbon (OC) into the deep ocean, nourishing deep-sea ecosystems and contributing to the global carbon cycle through OC burial and sequestration. Tracking lateral OC transport through submarine canyon systems is challenged by the deep-ocean setting, difficulties with constraining episodic depositional events, and the need to assess the composition and age of marine and terrestrial organic matter. We apply innovative parallel ramped pyrolysis oxidation accelerator mass spectrometry and pyrolysis-gas chromatography-mass spectrometry with isotope analyses to track OC age and sources in the 2016 Kaikoura earthquake-triggered, canyon-flushing event that deposited along >1300 km of a submarine canyon-channel system, offshore Aotearoa New Zealand. Specifically, these techniques allow us to determine the ages, sources, and partitioning of OC within the Kaikoura turbidite deposit and test hypotheses of how submarine canyon systems contribute to lateral OC flux and burial. Our results show that, despite considerable canyon floor erosion, substantial amounts of young OC were flushed into the deep sea, with relatively little (~2%) pre-Holocene OC contributions. Even without a direct connection between rivers and submarine canyons, most (~55%) of the OC in the Kaikoura event bed is from terrestrial sources. However, the deposit also contains substantial amounts (~22%) of marine-derived OC and ~23% of the material is of unassignable origin. Particle sorting imparts variability on the age and composition of OC within turbidite deposits and along the turbidity current flow path. Terrestrial-derived OC is preferentially older than marine-derived OC and concentrated in coarser particle sizes found more commonly at the deposit base and in proximal settings. Young, marine-derived OC is concentrated at the surface of the deposits and tends to be enriched in finer particle sizes. Such OC partitioning in turbidites supports the relevance of depositional models for predicting and quantifying distribution of OC in deep-sea deposits. Earthquake-triggered, canyon flushing events and resulting turbidites enhance OC burial efficiency and can sequester OC effectively, contributing an important carbon sink to the sedimentary carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Distribution of rare earth elements in the pilot-scale processing of fly ashes derived from eastern Kentucky coals: Comparisons of the feed and processed ashes

Beneficiated fly ash from the combustion of Central Appalachian high volatile bituminous coals was extracted with HNO 3 in a pilot-scale processing plant. Several major oxides (notably CaO and SO 3 , but also including Fe 2 O 3 , MgO, K 2 O, and P 2 O 5 ) and minor elements (Mn, As, Sr, Ba, and Pb) are depleted in the post-HNO 3 -extraction spent ash. The total lanthanides, Y, and Sc concentration is reduced by about 20% in the spent ash, with Gd showing the greatest decrease. Along with Gd, Nd and Dy are also well differentiated between the feed and spent ashes, with La and Sm showing minimal partitioning. The Gd decrease is correlated with the depletion of Fe 2 O 3 . The heavy rare earth elements (REE heavier than Eu) and Y are disproportionately concentrated in the HNO 3 -leachate compared to the light REE. For the ashes studied, Sc did not partition between the feed and spent ashes. Pozzolanicity tests show that the compressive strength and strength activity indices of the spent ash + ordinary Portland cement (OPC) mixes are comparable to 100% OPC, indicating that the spent ashes produced in the pilot-scale runs have the potential to be sold as a Class F fly ash. Finally, the beneficiated ash chemistry influences the chemistry of the post-HNO 3 -extraction spent ash and the HNO 3 -leachate. A 500-ppm-REE fly ash will presumably be a more economically favorable feedstock than an ash with a significantly lesser concentration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2016

Theoretical estimates of equilibrium sulfur isotope effects in aqueous sulfur systems: Highlighting the role of isomers in the sulfite and sulfoxylate systems

We present theoretical calculations for all three isotope ratios of sulfur ( 33 S/ 32 S, 34 S/ 32 S, 36 S/ 32 S) at the B3LYP/6-31+G(d,p) level of theory for aqueous sulfur compounds modeled in 30–40H 2 O clusters spanning the range of sulfur oxidation state (S n , n = -2 to +6) for estimating equilibrium fractionation factors in aqueous systems. Computed 34 β values based on major isotope ( 34 S/ 32 S) reduced partition function ratios (RPFRs) scale to a first order with sulfur oxidation state and coordination, where 34 β generally increase with higher oxidation state and increasing coordination of the sulfur atom. Exponents defining mass dependent relationships based on β values ( x/34 κ = ln( x β)/ln( 34 β), x = 33 or 36) conform to tight ranges over a wide range of temperature for all aqueous compounds ( 33/34 κ ≈ 0.5148–0.5159, 36/34 κ ≈ 1.89–1.90 from T ≽ 0 °C). The exponents converge near a singular value for all compounds at the high temperature limit ( 33/34 κT→∞= 0.51587 ± 0.00003 and 36/34 κT→∞= 1.8905 ± 0.0002; 1 s.d. of all computed compounds), and typically follow trends based on oxidation state and coordination similar to those seen in 34 β values at lower temperatures. Theoretical equilibrium fractionation factors computed from these β-values are compared to experimental constraints for HSO 3 - T(aq)/SO 2(g, aq) , SO 2(aq) /SO 2(g) , H 2 S (aq) /H 2 S (g) , H 2 S (aq) /HS - (aq) , SO 4 2- (aq) /H 2 S T(aq) , S 2 O 3 2- (aq) (intramolecular), and S 2 O 3 2- (aq) /H 2 S T(aq) , and generally agree within a reasonable estimation of uncertainties. We make predictions of fractionation factors where other constraints are unavailable. Isotope partitioning of the isomers of protonated compounds in the sulfite and sulfoxylate systems depend strongly on whether protons are bound to either sulfur or oxygen atoms. The magnitude of the HSO 3 - T /SO 3 2- major isotope ( 34 S/ 32 S) fractionation factor is predicted to increase with temperature from 0 to 70 °C due to the combined effects of the large magnitude (HS)O 3 - /SO 3 2- fractionation factor (1000ln 34 α (HS)bisulfite-sulfite = 19.9‰, 25 °C) relative to the (HO)SO 2 - /SO 3 2- fractionation factor (1000ln 34 α (HO)bisulfite–sulfite = -2.2‰, 25 °C), and the increased stability of the (HS)O 3 - isomer with increasing temperature. Finally, we argue that isomerization phenomenon should be considered in models of the sulfur cycle, including models that describe the overall sulfur isotope fractionations associated with microbial metabolism (e.g., microbial sulfate reduction).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Refined views of ancient ocean chemistry: Tracking trace element incorporation in pyrite framboids using atom probe tomography

The trace element chemistry of pyrite can be used to determine the origin, timing, and conditions of formation of ore deposits; as a vectoring tool for mineral exploration; and to determine the evolution of the Earth’s oceanic and atmospheric chemistry. However, little is known about whether trace elements are held with the pyrite structure or within nano-inclusions of other phases. This distinction is important for two primary reasons. First, trace element incorporation into the pyrite structure can affect the partitioning of other trace elements. Second, if trace elements are held within nano-inclusions, the partition coefficient of the mineral phases that make up the nano-inclusion, rather than pyrite, are the critical consideration in related interpretations. Previous studies addressing this topic have used laser ablation inductively coupled plasma mass spectrometry; however, the resulting data do not provide sufficient spatial resolution to delineate trace element distributions unless the inclusions are large. Further, they use these flat element profiles in time resolved laser ablation output graphs to argue that pyrite trace element content provides a direct relationship between trace element content of seawater and pyrite trace element content. To improve resolution, we have used atom probe tomography to characterize trace element distributions in pyrite framboids from the Cariaco Basin and Demerara Rise in three-dimensions at sub-nanometer resolution. Manganese was found to be concentrated in the pyrite nanocrystal part of the original framboid structure. In contrast, Ni was mostly found along the grain boundary, though it still appeared to be contained within the pyrite structure. Copper was concentrated in later pyrite overgrowths, and As varied in its location. These observations suggest that some important trace elements are incorporated into pyrite during early diagenesis, even in euxinic settings dominated by water-column pyrite formation. Statistical analysis was used to determine whether trace elements were incorporated in the lattice or within nanoscale inclusions (referred to here are nano-inclusions). We found that As, Ni, Cu, and Mn were commonly held within the pyrite structure, but As, Ni, and Cu can also be held as nano-inclusions or within grain interfaces. Incorporation of As is known to enhance the incorporation of other trace elements and in this case appears to correlated to elevated Ni and Cu concentrations in the Cariaco Basin samples. Furthermore, understanding these relationships strongly impacts our ability to utilize pyrite trace metal concentrations to analyze and quantify early ocean chemistry and its evolution through time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Soil organic matter is principally root derived in an Ultisol under oak forest

To a large degree, the sources and stability of soil organic carbon remain poorly constrained. A clear understanding of links among the components of the soil C cycle is hampered by the complexity of the system as well as challenges associated with partitioning bulk soil C into meaningful fractions. A large accidental 14 CO 2 release at the Oak Ridge Reservation in Tennessee, USA provided a strong label pulse into adjacent, well-studied oak forests, resulting in highly elevated Δ 14 C values in leaf litter (~1000‰) and roots (~260–450‰). A four-year manipulative study was conducted to determine the relative contribution of litter versus roots to the bulk mineral soil C pool, as well as to free light, occluded light and heavy fractions. The heavy fraction was further split into fractions with densities of 1.7–2.4 g cm –3 and >2.4 g cm –3 to test the homogeneity of the mineral-associated fraction of C. Substantial concentrations of label were detected in all soil fractions within a year of the 14 CO 2 release, indicating rapid incorporation of newly fixed photosynthates in all fractions of soil organic C, regardless of differences in stability inferred by previous work. This rapid incorporation of label occurred only in treatments where roots were labeled, indicating that roots are the major source of inputs to mineral soil C stocks at these sites. Separation of the heavy fraction into subfractions of intermediate (1.7–2.4 g cm –3 ) and high (>2.4 g cm –3 ) density indicated that both subfractions incorporated label at similar rates, despite significant differences in degree of microbial processing. In general, the rate of label incorporation suggested a much faster turnover for all fractions than indicated by natural radiocarbon abundance values. This study suggests that within each soil fraction there are portions of slow-cycling and fast-cycling materials, and the determination of an average turnover time or mean age is dependent on experimental approach. The rapid incorporation of label into all fractions within a year regardless of inferred stability implies a high degree of heterogeneity in all fractions regardless of how finely the soils are partitioned. Further refinement of the nature and drivers of this heterogeneity could yield important insights into the soil C cycle.

58 GEOSCIENCES↗
DOE OSTI2016

Predictive Parameters of Symptomatic Hematochezia Following 5-Fraction Gantry-Based SABR in Prostate Cancer

Purpose: This study identified predictors of high-grade late hematochezia (HH) following 5-fraction gantry-based stereotactic ablative radiation therapy (SABR). Methods and Materials: Hematochezia data for 258 patients who received 35 to 40 Gy SABR in 5-fractions as part of sequential phase 2 prospective trials was retrieved. Grade 2 or higher late rectal bleeding was labeled HH. Hematochezia needing steroid suppositories, 4% formalin, or 1 to 2 sessions of argon plasma coagulation (APC) was labeled grade 2. More than 2 sessions of APC, blood transfusion, or a course of hyperbaric oxygen was grade 3 and development of visceral fistula, grade 4. Various dosimetric and clinical factors were analyzed using univariate and multivariate analyses. Receiver operating characteristic (ROC) curve analysis and recursive partitioning analysis were used to determine clinically valid cut-off points and identify risk groups, respectively. Results: HH was observed in 19.4%, grade ≥3 toxicity in 3.1%. Median follow-up was 29.7 months (interquartile range [IQR]: 20.6-61.7) Median time to develop HH was 11.7 months (IQR: 9.0-15.2) from the start of radiation. At 2 years, cumulative HH was 4.9%, 27.2%, and 42.1% in patients who received 35 Gy to prostate (4-mm planning target volume [PTV] margin), 40 Gy to prostate (5-mm PTV margin), and 40 Gy to prostate/seminal vesicles (5-mm PTV margin), respectively (P<.0001). In the ROC analysis, volume of rectum receiving radiation dose of 38 Gy (V38) was a strong predictor of HH with an area under the curve of 0.65. In multivariate analysis, rectal V38 (≥2.0 cm{sup 3}; odds ratio [OR]: 4.7); use of anticoagulants in the follow-up period (OR: 6.5) and presence of hemorrhoids (OR: 2.7) were the strongest predictors. Recursive partitioning analysis showed rectal V38 < 2.0 cm{sup 3}, and use of anticoagulants or rectal V38 ≥ 2.0 cm{sup 3} plus 1 other risk factor resulted in an HH risk of >30%. Conclusions: Rectal V38 and 2 clinical factors were strong predictors of HH following 5-fraction SABR. Planning constraints should keep rectal V38 below 2.0 cm{sup 3}.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

A new weak Galerkin finite element method for elliptic interface problems

We introduce and analyze a new weak Galerkin (WG) finite element method in this paper for solving second order elliptic equations with discontinuous coefficients and interfaces. Comparing with the existing WG algorithm for solving the same type problems, the present WG method has a simpler variational formulation and fewer unknowns. Moreover, the new WG algorithm allows the use of finite element partitions consisting of general polytopal meshes and can be easily generalized to high orders. Optimal order error estimates in both H1 and L2 norms are established for the present WG finite element solutions. We conducted extensive numerical experiments in order to examine the accuracy, flexibility, and robustness of the proposed WG interface approach. In solving regular elliptic interface problems, high order convergences are numerically confirmed by using piecewise polynomial basis functions of high degrees. Moreover, the WG method is shown to be able to accommodate very complicated interfaces, due to its flexibility in choosing finite element partitions. Finally, in dealing with challenging problems with low regularities, the piecewise linear WG method is capable of delivering a second order of accuracy in L∞ norm for both C1 and H2 continuous solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A robust hierarchical solver for ill-conditioned systems with applications to ice sheet modeling

A hierarchical solver is proposed for solving sparse ill-conditioned linear systems in parallel. The solver is based on a modification of the LoRaSp method, but employs a deferred-compression technique, which provably reduces the approximation error and significantly improves efficiency. Moreover, the deferred-compression technique introduces minimal overhead and does not affect parallelism. As a result, the new solver achieves linear computational complexity under mild assumptions and excellent parallel scalability. To demonstrate the performance of the new solver, we focus on applying it to solve sparse linear systems arising from ice sheet modeling. The strong anisotropic phenomena associated with the thin structure of ice sheets creates serious challenges for existing solvers. To address the anisotropy, we additionally developed a customized partitioning scheme for the solver, which captures the strong-coupling direction accurately. In general, the partitioning can be computed algebraically with existing software packages, and thus the new solver is generalizable for solving other sparse linear systems. Our results show that ice sheet problems of about 300 million degrees of freedom have been solved in just a few minutes using a thousand processors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A unified formulation of splitting-based implicit time integration schemes

Splitting-based time integration approaches such as fractional step, alternating direction implicit, operator splitting, and locally one dimensional methods partition the system of interest into components, and solve individual components implicitly in a cost-effective way. Here this work proposes a unified formulation of splitting time integration schemes in the framework of general-structure additive Runge–Kutta (GARK) methods. Specifically, we develop implicit-implicit (IMIM) GARK schemes, provide the order conditions for this class, and explain their application to partitioned systems of ordinary differential equations. We show that classical splitting methods belong to the IMIM GARK family, and therefore can be studied in this unified framework. New IMIM-GARK splitting methods are developed and tested using parabolic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗