Search NASA⌕ Search

SEARCH · Search NASA

Results for “mathematical methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Universal Nuclear Accident Dosimeter

The Lawrence Livermore National Laboratory (LLNL) Universal Nuclear Accident Dosimetry (UNAD) project is a four-year initiative aimed at advancing nuclear accident dosimetry methods. This article presents an overview of the research, key findings, and the progress made throughout the project. The primary goals included a background into the history of nuclear accident dosimetry, consolidating current dosimetry techniques within the NNSA/DOE complex, fostering collaboration among subject matter experts, and exploring novel technologies for potential implementation. The technical focus centered on investigating new and novel technologies, instrumentation methods, and analysis methods to develop recommendations for a potential nuclear accident dosimeter (NAD) to be universally deployed through the DOE complex. A multilaboratory and multinational Usergroup was established, conducting periodic meetings to facilitate knowledge exchange. The UNAD team has participated in two international nuclear accident dosimetry intercomparison exercises and one characterization exercise, where the existing LLNL NAD and a prototype alanine electron paramagnetic dosimeter NAD were deployed. Ongoing improvements are being made to the prototype NAD based on results from the exercises, laboratory studies, and collaboration with other laboratories. A machine learning algorithm to optimize the geometry and conversion factors of the current LLNL NAD is being implemented, and the resulting design will be tested in the next exercise. In conclusion, key lessons learned and future directions for the project are discussed.

Electron paramagnetic resonance spectroscopy↗
DOE OSTI2025

Thermal Weight Determination and Interstate Coupling in State-Averaged ADAPT-VQE

Characterizing electronic thermal states at low temperatures is an important but challenging task in quantum chemistry and condensed matter physics, making it a prime candidate for a useful application in quantum computing. One of the most successful methods for state preparation on quantum computers is the Adaptive, Problem-Tailored (ADAPT) Variational Quantum Eigensolver (VQE), which has recently been generalized to treat excited states within a state-averaged framework as well as Gibbs states. In this work, we introduce Helmholtz-Optimized Thermal (HOT) ADAPT-VQE, an ancilla-free strategy for preparing Gibbs states that directly minimizes the Helmholtz free energy by targeting the dominant eigenstates of the thermal ensemble. We demonstrate the usefulness of HOT-ADAPT-VQE by predicting the free energy of two model systems with strongly correlated ground states: (1) the Fe 2+ cation in a magnetic field and (2) a [Cu 2 O 7 ] 10– fragment of the Mott insulator La 2 CuO 4 . Our results demonstrate that HOT-ADAPT-VQE significantly improves upon Gibbs-state estimates from multistate variants of ADAPT-VQE, often with substantially shallower quantum circuits, making it a promising candidate for thermal-state calculations.

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

Numerical Investigation of Susceptor-Catalyst Design for Ethylene Generation in Radio Frequency Based Reactors

Ethylene is the most widely produced petrochemical component in the world. Whether reactors are heated directly or indirectly via steam, manufacturers use economies of scale to overcome inherent thermodynamic inefficiencies when burning fossil fuels. While allowing large-scale operations to use alternative sources of energy and raw materials, new methods of supplying energy to chemical reactor systems can reduce the energy waste produced by conventional processes. One viable method for effectively supplying energy to reactor systems is electromagnetic (EM) induction heating. Using the properties of radio frequency (RF) waves, heterogeneous catalysts can be precisely targeted for heating inside reactors. Site-selective heating can greatly lower the energy requirements of the process by supplying heat to reaction sites while reducing needless heat transfer elsewhere. In this study, a microscale model was used to help create guidelines for susceptors and catalysts to improve ethylene production. The oxidative dehydrogenation of ethane is investigated by utilizing several catalysts and potential catalyst/susceptor combinations, with heat provided by an EM susceptor. Having the susceptors and catalyst function separately results in higher gradients in both heat and mass transfer, which drives transport through the catalyst region, while still providing adequate heating for endothermic reactions. However, using a susceptive core with a catalyst covering produces the maximum ethylene concentration.

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

Discovery of an Interlocked and Interwoven Molecular Topology in Nanocarbons via Dynamic C–C Bond Formation

Topologically complex carbon nanostructures are an exciting but largely unexplored class of materials due to their challenging synthesis. Previous methods are low yielding because they rely on irreversible C sp 2 −C sp 2 bond formation, which necessitates complex templating strategies to enforce entanglement. Here, reversible zirconocene coupling of alkynes is developed as a new method to access complex molecular topologies, where dynamic C−C bond formation facilitates entanglement under thermodynamic control, allowing the use of simple precursors without the need for preassembly. This strategy enables the scalable, high-yield synthesis of three topologically distinct nanocarbons, including the serendipitous discovery of a structure containing a new topological motif that was not previously identified or realized synthetically. This motif, consisting of an unusual combination of interlocking and interweaving, was recognized to be generalizable to a new topological class of molecules, introduced here as perplexanes.

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

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Lastly, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

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

Classical Preoptimization Approach for ADAPT-VQE: Maximizing the Potential of High-Performance Computing Resources to Improve Quantum Simulation of Chemical Applications

The ADAPT-VQE algorithm is a promising method for generating a compact ansatz based on derivatives of the underlying cost function, and it yields accurate predictions of electronic energies for molecules. In this work, we report the implementation and performance of ADAPT-VQE with our recently developed sparse wave function circuit solver (SWCS) in terms of accuracy and efficiency for molecular systems with up to 52 spin orbitals. The SWCS can be tuned to balance computational cost and accuracy, which extends the application of ADAPT-VQE for molecular electronic structure calculations to larger basis sets and a larger number of qubits. Using this tunable feature of the SWCS, we propose an alternative optimization procedure for ADAPT-VQE to reduce the computational cost of the optimization. Furthermore, by preoptimizing a quantum simulation with a parametrized ansatz generated with ADAPT-VQE/SWCS, we aim to utilize the power of classical high-performance computing in order to minimize the work required on noisy intermediate-scale quantum hardware, which offers a promising path toward demonstrating quantum advantage for chemical applications.

ADAPT-VQE↗
DOE OSTI2026

Dara: Automated Multiple-Hypothesis Phase Identification and Refinement from Powder X-ray Diffraction

Powder X-ray diffraction (XRD) is a foundational technique for characterizing crystalline materials. However, the reliable interpretation of XRD patterns, particularly in multiphase systems, remains a manual and expertise-demanding task. As a characterization method that only provides structural information, multiple reference phases can often be fit to a single pattern, leading to potential misinterpretation when alternative solutions are overlooked. To ease humans’ efforts and address the challenge, we introduce Dara (data-driven automated Rietveld analysis), a framework designed to automate the robust identification and refinement of multiple phases from powder XRD data. Dara performs an exhaustive tree search over all plausible phase combinations within a given chemical space and validates each hypothesis using the BGMN Rietveld refinement routine. Key features include structural database filtering, automatic clustering of isostructural phases during tree expansion, and peak-matching-based scoring to identify promising phases for refinement. When ambiguity exists, Dara generates multiple hypothesis which can then be decided between by human experts or with further characterization tools. By enhancing the reliability and accuracy of phase identification, Dara enables scalable analysis of realistic complex XRD patterns and provides a foundation for integration into multimodal characterization workflows, moving toward fully self-driving materials discovery.

Biological databases↗
DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗
DOE OSTI2024

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

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

Capturing Ring Opening in Photoexcited Enolic Acetylacetone upon Hydrogen Bond Dissociation by Ultrafast Electron Diffraction

Photoinduced biological and chemical reactions are often based on key structural transformations of a molecule driven across multiple electronic states. Acetylacetone (AcAc) is a prototypical system for complex chemical pathways involving several conical intersections (CI) and singlet–triplet intersystem crossings (ISC) characterized by distinct geometries. In the gas phase, AcAc is predominantly in a planar ring-like enolic form stabilized by a strong intramolecular O–H···O hydrogen bond. Following excitation into the S 2 (ππ*) state at 266 nm, acetylacetone undergoes rapid internal conversion followed by intersystem crossing. Such relaxation pathways are associated with structural changes including ring opening, deplanarization, and bond elongation. In this work, ultrafast electron diffraction (UED) at the SLAC MeV-UED setup is employed as a direct structural probe with a time resolution of 160 fs. Together with trajectory surface hopping simulations, analysis of the UED data provides a new perspective on the early time nuclear dynamics in acetylacetone. Specifically, AcAc is observed to undergo ring opening, deplanarization, and bond elongation all within the first 700 fs after photoexcitation. The monitored dynamics is associated mainly with the nuclear motion on the S 1 potential energy surface, formed after very rapid transfer from S 2 to S 1 , allowing AcAc to reach the conical intersection to intersystem crossing. Such time scales of nuclear motion are contrasted with the time scales of electronic transitions in AcAc that were previously characterized with spectroscopic methods, specifically internal conversion (<100 fs) and intersystem crossing (∼1.5 ps).

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

Spatially Patterned Architectures to Modulate CO 2 Reduction Cascade Catalysis Kinetics

Electrochemical CO 2 reduction using renewable sources of electrical energy holds promise for converting CO 2 into fuels and chemicals. The complex interactions among chemical/electrochemical reactions and mass transport make it difficult to analyze the effect of an individual process on electrode performance based only on experimental methods. Here, we developed a generalized steady-state simulation to describe an electrode surface in which sequential cascade catalysts are patterned in a periodic trench design. If appropriately constructed, this trench geometry is hypothesized to be able to yield a higher net current density for a CO 2 reduction (CO 2 R) cascade reaction. We have used realistic experimental reaction kinetics to investigate the role of trench geometry in mass transport, local microenvironments, and selectivity for a model CO 2 R cascade reaction. The model considers local concentration gradients of bicarbonate species at quasi-equilibrium and catalytic surface reactions based on concentration-dependent Butler–Volmer kinetics. Our results suggest that varying the spatial distribution of active sites plays a significant role in facilitating effective mass transport between active sites, modulating selectivity for the cascade reaction, and enhancing the yield of desirable cascade products. Moreover, we observe that this trench geometry significantly alters the cascade reaction rate by affecting the local pH, which can cause inadvertent depletion of available aqueous CO 2 to limit the CO 2 R cascade kinetics and modest suppression of the hydrogen evolution reaction (HER). The results highlight the trade-offs between mass transport, pH, and reaction kinetics that become apparent only when considering the coupled physics of all processes at the electrode surface. Here, this model can thus serve as a primary tool to build more selective and efficient patterned architectures for the CO 2 R cascade catalysis.

CO2 reduction↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2025

Computational Investigation of the Chemical Bond between An(III) Ions and Soft-Donor Ligands

The chemical bonding of actinide ions with arene and borohydride ligands is explored via quantum chemical methods to understand how the transuranium elements interact with softdonor ligands. Specifically, the [An(C 6 Me 6 )(BH 4 ) 3 ] complexes (An = U, Np, and Pu) and their reduced congeners are studied. Density functional theory (DFT) shows that the metal–ligand interactions in the neutral complexes are governed by electrostatic interactions. Both DFT and complete active space (CASSCF) results show that as one moves from U to Pu, the 5f-orbitals are stabilized leading to a poorer energy match with the ligand orbitals. This contributes to progressively weaker metal-arene and metal-borohydride interactions across the series due to a decrease in energy-driven covalency. A reduction in orbital contributions to bonding is obtained for the transuranium-arene interactions as well. Upon reduction, the arene is reduced, forming a δ-bond. This causes the An–arene distances to contract by 0.1–0.2 Å compared to the neutral complexes. The ground state is assigned as the intermediate-spin state where the arene radical is antiferromagnetically coupled to the metal-centered f-electrons in Np and Pu. On the other hand, the ferromagnetically and antiferromagnetically coupled states are close in energy in the uranium complex, but do not mix when spin– orbit coupling is included using a state-interaction approach (SO-CASPT2). The population of the CASSCF δ*-antibonding natural orbital increases from U to Pu consistent with the increased An–arene distances, weaker interactions, and decreasing covalency across the series. Although the An–B distance increases by ca. 0.06 Å upon reduction, both the neutral and reduced species involve an An(III)–borohydride bond and as such are qualitatively similar. The Np complexes can be assigned to have slightly weaker bonding than the uranium analogs but are overall “uranium-like”. The Pu complexes are predicted to have less covalent contributions to bonding in both the Pu–arene and Pu–borohydride interactions; however, the Pu–arene interaction is predicted to be particularly weak.

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

Emerging Per- and Polyfluoroalkyl Substances in Tap Water from the American Healthy Homes Survey II

Humans experience widespread exposure to anthropogenic per- and polyfluoroalkyl substances (PFAS) through various media, which can lead to a wide range of negative health impacts. Tap water is an important source of exposure in communities with any degree of contamination but routine or large-scale PFAS monitoring often depends on targeted analytical methods limited to measuring specific PFAS. We analyzed 680 tap water samples from the American Healthy Homes Survey II for PFAS using non-targeted analysis (NTA) to expand the range of detectable PFAS. Based on detection frequency and relative abundance, about half of the identified PFAS were found only by NTA. We identified (with varying degrees of confidence) 75 distinct PFAS, including 57 exclusively detected by NTA. The identified PFAS are members of seven structural subclasses differentiated by their head groups and degree of fluorination. Clustering analysis categorized the PFAS into four coabundance groups dominated by specific PFAS subclasses. One group uniquely identified by NTA contains zwitterionic PFAS and other PFAS transformation products which are likely associated with aqueous firefighting foam contaminants in a small number of spatially correlated samples. These results help further characterize the scope of exposure to emerging PFAS experienced by the U.S. population via tap water and augment nationwide targeted-PFAS monitoring programs.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

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

Ladder symmetries and Love numbers of Reissner-Nordström black holes

It is well known that asymptotically flat black holes in general relativity have vanishing tidal Love numbers. In the case of Schwarzschild and Kerr black holes, this property has been shown to be a consequence of a hidden structure of ladder symmetries for the perturbations. In this work, we extend the ladder symmetries to non-rotating charged black holes in general relativity. As opposed to previous works in this context, we adopt a more general definition of Love numbers, including quadratic operators that mix gravitational and electromagnetic perturbations in the point-particle effective field theory. We show that the calculation of a subset of those couplings in full general relativity is affected by an ambiguity in the split between source and response, which we resolve through an analytic continuation. As a result, we derive a novel master equation that unifies scalar, electromagnetic and gravitational perturbations around Reissner-Nordström black holes. The equation is hypergeometric and can be obtained from previous formulations via nontrivial field redefinitions, which allow to systematically remove some of the singularities and make the presence of the ladder symmetries more manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2025

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗