Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimensional analysis”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SpectraCodec: A Hilbert curve-based method for encoding metadata in mass spectra for machine learning applications (SpectraCodec) v1

Machine learning approaches to mass spectrometry (MS) data analysis require structured metadata for optimal performance. However, current MS file formats necessitate external metadata sources, creating integration challenges that impede analytical workflows. Here, we present a novel approach for encoding metadata directly within mzML files using one-hot encoding of ASCII characters mapped via Hilbert space-filling curves. This strategy embeds metadata in the first spectrum's m/z-intensity space, ensuring persistence with the primary data, eliminating the need for external metadata files, and maintaining compatibility with existing MS software. We demonstrate that the Hilbert curve mapping efficiently utilizes the two-dimensional spectral space while maintaining robust data recovery. This method offers a practical solution for machine learning applications in mass spectrometry by ensuring metadata and spectral data remain unified through all stages of analysis.

Bowen, Benjamin [Lawrence Berkeley National Labora↗
DOE OSTI2025

Transverse orbital angular momentum and polarization entangled spatiotemporal structured light

Intra-system entanglement occurs between non-separable modes within the same system. For optical systems, the various degrees of freedom of light represent different modes, and the potential use of light to create higher dimensional classical entangle states offers a promising potential to drive new technological developments. In this work, we present experimental results demonstrating the orthogonality between transverse orbital angular momentum (t-OAM) of different spatiotemporal topological charges, a previously unverified property of t-OAM. Based on those results, we developed methods to create and characterize a novel family of t-OAM and polarization entangled spatiotemporal structured light. We further provide theoretical analysis to support our study of the entanglement between those modes. By demonstrating the feasibility of leveraging t-OAM as a new family of modes for classical entanglement, our work represents a new advancement towards higher dimensional classical entanglement strategies.

42 ENGINEERING↗
DOE OSTI2022

Comparative Effectiveness Analysis of 3D-Conformal Radiation Therapy Versus Intensity Modulated Radiation Therapy (IMRT) in a Prospective Multicenter Cohort of Patients With Breast Cancer

Simple intensity modulation of radiation therapy reduces acute toxicity compared with 2-dimensional techniques in adjuvant breast cancer treatment, but it remains unknown whether more complex or inverse-planned intensity modulated radiation therapy (IMRT) offers an advantage over forward-planned, 3-dimensional conformal radiation therapy (3DCRT).

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Three-Dimensional Grid Visualization for Planning Activities: A Dubai Case Study

National Laboratory of the Rockies (NLR), in collaboration with the Dubai Electricity and Water Authority (DEWA) and Infra-X, has undertaken the Energy Visualization Analysis Project. The aim of this project is to enhance analytical and 3D visualization capabilities for distribution network planning and renewable energy integration. As modern grid continues to evolve with large-scale solar PV deployment and emerging distributed energy resources (DERs), the ability to effectively analyze, visualize, and communicate complex grid behaviors has become increasingly critical. The project focuses on developing empirical use cases based on real distribution feeder data and engineering workflows, ensuring the outcomes are directly aligned with operational environment. Through time-series power flow simulations and nodal hosting capacity analysis, the study quantifies the impacts of high PV penetration on voltage and thermal limits within representative 11 kV feeders. These analyses identify specific nodes and conditions where DER integration challenges arise. Furthermore, a Battery Energy Storage System (BESS) optimization algorithm was applied to determine the optimal size and placement of storage systems that can mitigate network constraints and enhance hosting capacity. The comparative results between base-case and BESS-augmented scenarios clearly demonstrate improvements in network stability and load management efficiency. In parallel, the NLR team developed an immersive 3D visualization framework, enabling interactive exploration of grid simulations using commodity head-mounted display (HMD) systems. This framework transforms conventional 2D simulation data into spatially intuitive visual environments - allowing engineers to analyze feeder conditions, PV hosting potential, and BESS effects in real time. This report represents the first foundational phase in establishing a visualization-driven analytical ecosystem. It provides a methodological foundation for data integration, visualization architecture, and simulation-based decision support, paving the way for large-scale adoption of immersive visualization across DEWA's Smart Grid Initiative, R&D activities, and future network resilience studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Beyond microbial abundance: metadata integration enhances disease prediction in human microbiome studies

Multiple studies have highlighted the interaction of the human microbiome with physiological systems such as the gut, immune, liver, and skin, via key axes. Advances in sequencing technologies and high-performance computing have enabled the analysis of large-scale metagenomic data, facilitating the use of machine learning to predict disease likelihood from microbiome profiles. However, challenges such as compositionality, high dimensionality, sparsity, and limited sample sizes have hindered the development of actionable models. One strategy to improve these models is by incorporating key metadata from both the human host and sample collection/processing protocols. This remains challenging due to sparsity and inconsistency in metadata annotation and availability. In this paper, we introduce a machine learning-based pipeline for predicting human disease states by integrating host and protocol metadata with microbiome abundance profiles from 68 different studies, processed through a consistent pipeline. Our findings indicate that metadata can enhance machine learning predictions, particularly at higher taxonomic ranks like Kingdom and Phylum, though this effect diminishes at lower ranks. Our study leverages a large collection of microbiome datasets comprising 11,208 samples, therefore enhancing the robustness and statistical confidence of our findings. This work is a critical step toward utilizing microbiome and metadata for predicting diseases such as gastrointestinal infections, diabetes, cancer, and neurological disorders.

Mathematics and Computing↗
DOE OSTI2025

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗
DOE OSTI2021

IMS Rapid Response FY21 Summary Report for: Integrating Patterned Probes with Four-Dimensional Scanning Transmission Electron Microscopy for Unrivaled Crystallographic Structure Determination in Nanomaterials

The initial goal of our 4-dimensional scanning transmission electron microscopy (4D-STEM)-based project was to develop strain resolution two orders of magnitude better than what is now currently possible with electron-based scattering techniques, all while collecting scattering information from 7 different tilt axes at one time [multi-beam electron diffraction (MBED)1 ] through the development of a new electron probe-forming aperture with non-circular features (patterned probes 2 ). We set out to accomplish this through a collaboration with Drs. Colin Ophus and Ben Savitsky at Lawrence Berkeley Laboratory (they are the world-leading experts in developing the complex computational codes required to perform orientation analysis and quantitative strain mapping on our 4D-STEM data sets. We are motivated to invest in this area as it will be the only technique sensitive enough to perform three- dimensional automated crystallographic orientation mapping (ACOM) and strain mapping for materials exposed to external stimulus (a focus of our larger efforts).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chimera Spectrum Diagnostics for Peptides Using Two-Dimensional Partial Covariance Mass Spectrometry

The rate of successful identification of peptide sequences by tandem mass spectrometry (MS/MS) is adversely affected by the common occurrence of co-isolation and co-fragmentation of two or more isobaric or isomeric parent ions. This results in so-called `chimera spectra’, which feature peaks of the fragment ions from more than a single precursor ion. The totality of the fragment ion peaks in chimera spectra cannot be assigned to a single peptide sequence, which contradicts a fundamental assumption of the standard automated MS/MS spectra analysis tools, such as protein database search engines. This calls for a diagnostic method able to identify chimera spectra to single out the cases where this assumption is not valid. Here, we demonstrate that, within the recently developed two-dimensional partial covariance mass spectrometry (2D-PC-MS), it is possible to reliably identify chimera spectra directly from the two-dimensional fragment ion spectrum, irrespective of whether the co-isolated peptide ions are isobaric up to a finite mass accuracy or isomeric. We introduce ‘3-57 chimera tag’ technique for chimera spectrum diagnostics based on 2D-PC-MS and perform numerical simulations to examine its efficiency. We experimentally demonstrate the detection of a mixture of two isomeric parent ions, even under conditions when one isomeric peptide is at one five-hundredth of the molar concentration of the second isomer.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

Magnetization dynamics of weakly interacting sub-100 nm square artificial spin ices

Artificial Spin Ice (ASI), consisting of a two dimensional array of nanoscale magnetic elements, provides a fascinating opportunity to observe the physics of out-of-equilibrium systems. Initial studies concentrated on the static, frozen state, whilst more recent studies have accessed the out-of-equilibrium dynamic, fluctuating state. This opens up exciting possibilities such as the observation of systems exploring their energy landscape through monopole quasiparticle creation, potentially leading to ASI magnetricity, and to directly observe unconventional phase transitions. In this work we have measured and analysed the magnetic relaxation of thermally active ASI systems by means of SQUID magnetometry. We have investigated the effect of the interaction strength on the magnetization dynamics at different temperatures in the range where the nanomagnets are thermally active. We have observed that they follow an Arrhenius-type Néel-Brown behaviour. An unexpected negative correlation of the average blocking temperature with the interaction strength is also observed, which is supported by Monte Carlo simulations. The magnetization relaxation measurements show faster relaxation for more strongly coupled nanoelements with similar dimensions. The analysis of the stretching exponents obtained from the measurements suggest 1-D chain-like magnetization dynamics. This indicates that the nature of the interactions between nanoelements lowers the dimensionality of the ASI from 2-D to 1-D. Finally, we present a way to quantify the effective interaction energy of a square ASI system, and compare it to the interaction energy computed with micromagnetic simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Effects of high-dose neutron irradiation at light-water reactor relevant temperature on the mechanical properties of SiC/SiC composites

For this study, the neutron dose-dependent evolutions and the underlying mechanisms of properties of SiC fiber-reinforced SiC matrix (SiC/SiC) composites at a temperature relevant to light-water reactors (∼600 K) were investigated and analyzed. Chemical vapor infiltrated (CVI) SiC/SiC composites reinforced with Hi-Nicalon Type S or Tyranno SA3 fiber were neutron-irradiated to doses up to 30 dpa. The irradiated composites retained their flexural strengths. The thermal diffusivity and dimensional changes were mostly retained from 2.0 to 30.2 dpa. Discrepancies in irradiation responses among CVI SiC/SiC composites from different sources were found. Additional microstructural analysis using Raman spectroscopy and numerical analysis on irradiation effect on residual stress were used to explain how the microstructural variables, especially of carbon interphases, affect the mechanical properties in the 30–40 dpa dose range.

Continuous fiber-reinforced ceramic matrix composi↗
DOE OSTI2024

Antiferromagnetic Ordering in Quasi-One-Dimensional FeBi 4 S 7

Here, we report a detailed study of the synthesis, composition, magnetic structure, and transport properties of a quasi-one-dimensional antiferromagnet FeBi 4 S 7 that contains chains of edge-sharing FeS 6 octahedra. High-resolution powder X-ray diffraction (PXRD) analysis, aided by variation of synthetic conditions, suggests that the true formula of the material is Fe 1.2 Bi 3.8 S 7 , due to the minor substitution of Fe into Bi sites. This finding is in agreement with crystal structure refinement from neutron powder diffraction data as well as with the small band gap of 0.23 eV determined from electrical transport measurements. Analysis of the neutron diffraction pattern collected below the antiferromagnetic ordering temperature of 64 K revealed ferromagnetic coupling between the Fe moments in the chains of FeS 6 octahedra. The overall ordering, however, is antiferromagnetic due to the antiparallel arrangement of moments on neighboring chains. The collinear spin arrangement is described by a k-vector (1, 0, 1/2), which indicates doubling of the unit cell in the c direction and the loss of the C-centering translation as compared to the nuclear cell. The ferromagnetic nature of the sulfide-bridged chains of Fe 2+ ions in FeBi 4 S 7 , in contrast to the antiferromagnetic coupling between Fe moments in compounds with similar structural fragments, can be justified by the analysis of metric parameters that characterize the Fe–S bonding in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Revealing Variable Dependences in Hexagonal Boron Nitride Synthesis via Machine Learning

Wafer-scale monolayer two-dimensional (2D) materials have been realized by epitaxial chemical vapor deposition (CVD) in recent years. To scale up the synthesis of 2D materials, a systematic analysis of how the growth dynamics depend on the growth parameters is essential to unravel its mechanisms. However, the studies of CVD-grown 2D materials mostly adopted the control variate method and considered each parameter as an independent variable, which is not comprehensive for 2D materials growth optimization. Herein, we synthesized a representative 2D material, monolayer hexagonal boron nitride (hBN), on single-crystalline Cu (111) by epitaxial chemical vapor deposition and varied the growth parameters to regulate the hBN domain sizes. Furthermore, we explored the correlation between two growth parameters and provided the growth windows for large flake sizes by the Gaussian process. Here, this new analysis approach based on machine learning provides a more comprehensive understanding of the growth mechanism for 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Amino acid interacting network in the receptor-binding domain of SARS-CoV-2 spike protein

The relation between amino acid (AA) sequence and biologically active conformation controls the process of polypeptide chains folding into three-dimensional (3d) protein structures. The recent achievements in the resolution achieved in cryo-electron microscopy coupled with improvements in computational methodologies have accelerated the analysis of structures and properties of proteins. However, the detailed interaction between AAs has not been fully elucidated. Herein, we present a de novo method to evaluate inter-amino acid interactions based on the concept of accurately evaluating the amino acid bond pairs (AABP). The results obtained enabled the identification of complex 3d long-range interconnected AA interacting network in proteins. The method is applied to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. We show that although nearest-neighbor AAs in the primary sequence have large AABP, other nonlocal AAs make substantial contribution to AABP with significant participation of both covalent and hydrogen bonding. Detailed analysis of AABP in RBD reveals the pivotal role they play in sequence conservation with profound implications on residue mutations and for therapeutic drug design. This approach could be easily applied to many other proteins of biomedical interest in life sciences.

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

Optimization techniques in self-similar compressible flow

We investigate the one-dimensional (1D) inviscid compressible flow equations for an ideal gas through the lens of optimization techniques. It is the case that, to our knowledge, optimization analysis applied to the so-called “linear velocity” solutions of the Euler compressible flow equations has not been previously conducted. Through both gradient-based and variational techniques, new variants of well-studied flow scenarios, i.e., self-similar, 1D, linear velocity solution class to idealized inviscid compressible flow equations, are determined, as encoded in both the kinematic and thermodynamic properties of this self-similar solution class. With the kinematics of the said solutions being driven by a self-similar “scale radius” and the thermodynamics being driven separately through the appearance of an arbitrary function, a myriad of new solution classes is possible. Acting as a guide to more realistic physical circumstances as well as discovery, it is the hope that the presented cases serve as the framework for future investigations into the intersection of self-similarity and optimization techniques. Fields of study that may find this work to be of interest include aerodynamic design, flow control, inertial confinement fusion, physics-informed neural networks, and other related areas of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗
DOE OSTI2020

Experimental Investigation of Polymer Induced Fouling of Heater Tubes in the First-ever Polymer Flood Pilot on Alaska North Slope

Polymer flooding has been validated in the lab and by successful applications in other countries, such as Canada and China, to be an effective EOR technique for heavy oil reservoirs. Currently, polymer flooding is being pilot tested for the first time in the Schrader Bluff viscous oil reservoir at Milne Point field on Alaska North Slope (ANS). One of the major concerns of the operator is the impact of polymer on the oil production system after polymer breakthrough, especially the polymer induced fouling issues in the heat exchanger. This study investigates the propensity of polymer fouling on the heater tubes as a function of different variables, with the ultimate goal of determining safe and efficient operating conditions. A unique experimental set-up was indigenously designed and developed to simulate the fouling process on the heating tube. The influence of heating tube skin temperature, tube material, and polymer concentration on fouling tendency was investigated. Under each test condition, the test was run five times with the same tube, and in each run, the freshly prepared synthetic brine and polymer solution were heated from 77°F to 122°F to mimic field operating conditions. The heating time and fouling amount were recorded for each run. The morphology and composition of the deposit samples were analyzed by environmental scanning electron microscopy (ESEM) and X-ray diffraction (XRD), respectively. It has been found that, in general, the presence of polymer in produced fluids would aggravate the fouling issues on both carbon steel and stainless steel surfaces at all tested skin temperatures (165°F, 250 °F and 350°F), but only higher skin temperatures of 250°F and 350°F could cause polymer induced fouling issues on the copper tube surface and the fouling tendency increased with polymer concentration. At the lower skin temperatures of 165°F, no polymer-induced fouling was identified on the copper tube. A critical temperature that is related to the cloud point of the polymer solution was believed to exist, below which polymer induced fouling would not occur and only mineral scale deposited, but above which the polymer would obviously aggravate the fouling issue. The heating efficiency of the tube would be reduced gradually as more fouling material accumulates on its surface. The ESEM results indicate that if polymer precipitated and deposited on the surface, it would bond to the mineral crystals to form a stronger three- dimensional network structure, and that is why the polymer induced fouling was tougher and more difficult to be removed from the tube surface. The XRD analysis results confirm that the presence of polymer in the fluids can enhance the mineral scale propensity. Overall, the study results have provided practical guidance to the field operator for the ongoing polymer flooding pilot test on ANS. This study may also prove to be valuable for other chemical EOR projects around the world.

Alaska North Slope↗
DOE OSTI2020

Novel Geometric Operations for Linear Programming

This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗