Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Heat accelerates degradation of β -lactoglobulin fibrils at neutral pH

Fibrous aggregates of β-lactoglobulin display superior mechanical and interfacial properties compared to the native protein. These properties directly link to the protein morphology and structure. When incorporated into food matrices, during processing protein fibrils are exposed to pH shifts and high temperature conditions, which accelerate their degradation. In the present study, neutralized β-lactoglobulin fibrils were heated at 100 °C and 121 °C for various times to assess their degradation. Fibril morphology, structure, and viscosity in solution were examined by microscopy, scattering, spectroscopy, and rheology. Atomic force microscopy showed the contour length of the protein fibrils decreased gradually with heating at 100 °C and 121 °C, with greater decreases at 121 °C. Increased fibril diameters (~15–25 nm) were observed at 121 °C for 5–15 min heating and were disrupted upon further heating. Small-angle x-ray scattering indicated an increase in fibril radius with heating at pH 7 followed by a decrease at prolonged heating, whereas fibril length decreased continuously with heating. Thioflavin T fluorescence, circular dichroism and Fourier transform infrared spectroscopy confirmed the conversion of β-sheet to random coils as fibrils were degraded during thermal treatment at pH 7. Surface hydrophobicity of fibrils decreased with increase in heating temperature and time, coinciding with an increase in the content of non-aggregated proteins. Viscosity of fibril solutions increased when fibrils were heated at 100 °C, whereas at 121 °C their viscosity first increased and then decreased. Furthermore, these findings imply heating at 100 °C and 121 °C facilitates degradation and depolymerisation of β-lactoglobulin fibrils with aggregation as an intermediate step.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗
DOE OSTI2026

Comparison of automated chemical-guided segmentation and human annotation of soil organic matter in X-ray microcomputed tomography imaging in contrasted soil types

Soil organic matter (OM) formation and persistence is strongly influenced by the spatial distribution of organic substrates and microscale soil heterogeneity by dictating OM accessibility to microorganisms. However, traditional size and/or density fractionation techniques disrupt aggregate architecture, eliminating spatial information needed to fully understand intra-aggregate OM distribution. To quantify three-dimensional OM spatial distribution and automate segmentation in X-ray microcomputed tomography (µCT) imaging without human annotation bias, we developed an iodine gas vapor (I2) based staining workflow that eliminates labor-intensive manual annotation while maintaining segmentation accuracy, using aggregates from four taxonomically diverse soils (Xerofluvent, Haploxeroll Sphagnofibrist, Palehumult) with an 8-fold range of soil organic carbon. Human annotation of 10 µCT slices by the experienced and inexperienced annotators resulted in variations up to 3% in the Dice similarity coefficient (DSC), reflecting a degree of inherent subjectivity of manual labeling. Such inconsistencies are expected to compound as the number of manually annotated slices increases. Dual-energy µCT imaging at 33.1 keV (below the iodine (I) K-edge) and 33.2 keV (above the I K-edge) was used to resolve aggregate microstructure following I2 staining. The automated image subtraction pipeline identified OM regions by the I Kedge induced brightness increases, achieving DSC values of 0.58–0.83 relative to an experienced annotator. Sensitivity analyses revealed that the reconstruction alpha value—optimized via the open-source tool TomocuPy—and the 3D registration slice count were the primary determinants of accuracy, providing a novel benchmark for dual-energy soil imaging. The pipeline without GPU acceleration achieved 9.6 to 43.2 times faster than manual annotation. Using GPU-accelerated image post-processing and affine transformation matrices, the pipeline successfully segmented OM elements for large-scale datasets (3232×3232 pixel, 2048 slices) within ~5200 s from raw file acquisition to segmented output. The high-throughput approach enables the quantification of OM spatial distribution across diverse and heterogeneous soil.

Soil microbial biomass↗
DOE OSTI2022

Branches, quivers, and ideals for knot complements

We generalize the F K invariant, i.e. Z^ for the complement of a knot K in the 3-sphere, the knots-quivers correspondence, and A-polynomials of knots, and find several interconnections between them. We associate an F K invariant to any branch of the A-polynomial of K and we work out explicit expressions for several simple knots. We show that these F K invariants can be written in the form of a quiver generating series, in analogy with the knots-quivers correspondence. We discuss various methods to obtain such quiver representations, among others using R-matrices. We generalize the quantum a-deformed A-polynomial to an ideal that contains the recursion relation in the group rank, i.e. in the parameter a, and describe its classical limit in terms of the Coulomb branch of a 3d-5d theory. We also provide t-deformed versions. Furthermore, we study how the quiver formulation for closed 3-manifolds obtained by surgery leads to the superpotential of 3d N = 2 theory T[M 3 ] and to the data of the associated modular tensor category MTC[M 3 ].

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Exploring the binding properties and activities of ancestral expansins

Bacterial expansins are non-lytic proteins capable of loosening cellulose networks, offering promising applications in agriculture, biotechnology, and material science. Their ability to disrupt noncovalent interactions in biopolymer matrices such as cellulose and chitin positions them as valuable tools for upgrading abundant natural materials. However, their industrial use remains limited due to their relatively low wall-loosening activity compared to plant expansins. To address this limitation, we applied Ancestral Sequence Resurrection (ASR) to reconstruct and characterize ancient variants of the Bacillus subtilis expansin BsEXLX1. ASR is a powerful evolutionary tool that enables the inference and synthesis of ancestral proteins, allowing researchers to explore functional traits that may have been lost over time. This approach not only provides insights into protein evolution but also facilitates the design of proteins with enhanced properties, such as improved substrate affinity or structural stability. In this study, we combined biochemical and biophysical assays to evaluate the activity and binding behavior of ancestral expansins. Our results reveal that ancestral variants exhibit increased cellulose affinity, reduced binding to acidic polysaccharides, and greater salt resistance. Furthermore, these traits enhance their wall-loosening activity and demonstrate the utility of ASR in engineering surface-active proteins for industrial applications, particularly in biomass processing and cellulose modification.

09 BIOMASS FUELS↗
DOE OSTI2020

Analytical Determination of Radioactive Strontium and Cesium by Thermal Ionization Mass Spectrometry

Thermal ionization mass spectrometry (TIMS) has been evaluated for the detection of the radioactive isotopes of strontium and cesium. The instrument was investigated for the analysis of isotopic ratios of 89Sr to 90Sr and 135Cs to 137Cs in the presence of atomic isobars (Y and Zr for Sr analysis and Ba for Cs analysis). Signicant discrimination using TIMS was found for both Sr and Cs through evaluation using stable isotopes. Fusion procedures with different flux were evaluated for preparation of sample matrices. Chemical separations were developed and used after fusion to isolate the Sr and Cs from a prepared background matrix (i.e., Montana Soil). A graded approach, increasing in stable background isotopes, was performed to monitor the chemistry and instrument response. A thermal neutron irradiation was performed to produce a sample that contained the radioisotopes of Sr and Cs. Even though the irradiated sample was not a certified standard it does provide accurate expectation values via published, well understood cumulative fission yield data. The results presented in this work correlate extremely well with the published data. Furthermore, a significant decrease in uncertainty (by a factor of 5 { 10) was found for the isotopic ratios measured with TIMS compared to the published values.

TIMS, environmental contamination, radiostrontium,↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Studying the effect of drug-to-excipient ratio on drug release profile for drug coated balloons

Drug-coated balloons (DCB) have emerged as the alternative procedure for restenosis because of their ability to treat a variety of occlusion types with a uniform dose of anti-proliferative drugs. DCB are balloons coated with antiproliferative drugs encapsulated in a polymer matrix. There are several types of coating matrices used to produce DCB. In this study, the relationship between coating composition and drug release under physiologically relevant conditions was examined to understand how differences in coating composition impacts the drug transfer from the balloon surface to the simulated body fluids. To conduct the experiments, the balloons were coated with different paclitaxel (drug)-to-iopromide (excipient) ratios (3:1, 3:2 and 1:2) using an in-house developed micro-pipetting method. Scanning electron microscopy (SEM) images showed that the 3:1 PTX:IOP ratio produced a more uniform, crystalline microstructure with a thinner coating throughout the balloon surface compared to the other drug-to-excipient ratios. The 1:2 PTX:IOP ratio showed the least crystalline microstructure among the three ratios evaluated in this study. Three different drug elution conditions were tested. The amount of drug released to the medium was quantified by high performance liquid chromatography (HPLC). Our soaking study and submerge & deploy study showed that ~20% of the drug transferred to the target site under physiological conditions. A track and deploy method was performed using a “mock” artery, to simulate an in vitro environment. Coated balloons were passed through the mock artery to mimic tracking turns the balloon within the arteries during the angioplasty procedures. Seven elution samples were collected at different stages of the procedure. Finally, drug release results suggest that the higher excipient ratio helps to deliver the lipophilic drug to the target site under simulated conditions but causes higher drug loss during the balloon transfer process.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

A characterization of recombinant Arabidopsis FRIABLE1 (FRB1) reveals robust rhamnogalacturonan-I rhamnosyltransferase activity and critical catalytic residues

Plant cell walls are glycan-rich extracellular matrices that fundamentally impact essential cellular processes, such as growth, adhesion, and cell shape acquisition. Understanding plant cell wall glycans requires the identification and characterization of the biosynthetic enzymes that produce these polymers. Most successful in vitro protein expression studies of plant cell wall glycosyltransferases have relied on insect, fungal/yeast, or human cell expression systems, whereas prokaryotic expression systems have been generally unsuccessful. Here, we show that Arabidopsis FRIABLE1 (FRB1)/rhamnogalacturonan-I rhamnosyltransferase 8 (RRT8) can be produced in Escherichia coli RosettaGami2 cells as N-terminal maltose-binding protein fusion proteins containing C-terminal 6X-His-tags. We also report the catalytic constants of FRB1/RRT8 with apparent K M and K cat values of 226 μM and 33 min -1 for UDP-Rhamnose and 117 μM and 28.7 min -1 for rhamnogalacturonan-I (RG-I), respectively. We examine the catalytic activities of mutated FRB1/RRT8 proteins based on an AlphaFold 3-generated FRB1/RRT8 protein structural model with a virtually docked UDP-Rha donor. Enzymatic characterization of the mutated and wildtype FRB1/RRT8 protein confirmed that mutation of predicted catalytic site amino acid residues resulted in a 20-fold reduction in RRT activity. FRB1 also robustly polymerizes RG-I in combination with RG-I galacturonosyltransferase 1. These results show how a robust E. coli expression system combined with artificial intelligence tools can be used to increase understanding of plant cell wall glycosyltransferase structure and function.

glycosyltransferase↗
DOE OSTI2020

Multigrid deflation for Lattice QCD

Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enhanced multiscale restriction-smoothed basis (MsRSB) preconditioning with applications to porous media flow and geomechanics

A novel method to enable application of the Multiscale Restricted Smoothed Basis (MsRSB) method to non M-matrices is presented. The original MsRSB method is enhanced with a filtering strategy enforcing M-matrix properties to enable the robust application of MsRSB as a preconditioner. Through applications to porous media flow and linear elastic geomechanics, the method is proven to be effective for scalar and vector problems with multipoint finite volume (FV) and finite element (FE) discretization schemes, respectively. Realistic complex (un)structured two- and three-dimensional test cases are considered to illustrate the method's performance.

58 GEOSCIENCES↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

GenMod: A generative modeling approach for spectral representation of PDEs with random inputs

Here, we propose a method for quantifying uncertainty in high-dimensional PDE systems with random parameters, where the number of solution evaluations is small. Parametric PDE solutions are often approximated using a spectral decomposition based on polynomial chaos expansions. For the class of systems we consider (i.e., high dimensional with limited solution evaluations) the coefficients are given by an underdetermined linear system in a regression formulation. This implies additional assumptions, such as sparsity of the coefficient vector, are needed to approximate the solution. Here, we present an approach where we assume the coefficients are close to the range of a generative model that maps from a low to a high dimensional space of coefficients. Our approach is inspired be recent work examining how generative models can be used for compressed sensing in systems with random Gaussian measurement matrices. Using results from PDE theory on coefficient decay rates, we construct an explicit generative model that predicts the polynomial chaos coefficient magnitudes. The algorithm we developed to find the coefficients, which we call GenMod, is composed of two main steps. First, we predict the coefficient signs using Orthogonal Matching Pursuit. Then, we assume the coefficients are within a sparse deviation from the range of a sign-adjusted generative model. This allows us to find the coefficients by solving a nonconvex optimization problem, over the input space of the generative model and the space of sparse vectors. We obtain theoretical recovery results for a Lipschitz continuous generative model and for a more specific generative model, based on coefficient decay rate bounds. We examine three high-dimensional problems and show that, for all three examples, the generative model approach outperforms sparsity promoting methods at small sample sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries

Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.

97 MATHEMATICS AND COMPUTING↗
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↗