SEARCH · Search NASA

Results for “computational geometry”

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.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗

DOE OSTI2024

Optimization of Fairing Geometry for ORPC Modular RivGen Power System: Cooperative Research and Development (Final Report)

The work will optimize the hydrodynamic performance and flow augmentation of Ocean Renewable Power Company's (ORPC) modular fairing for their cross-flow marine hydrokinetic (MHK) turbine, using a computational fluid dynamics study. The influence of the fairing cross-sectional shape and rotor-fairing spacing will be assessed, to maximize power production and minimize structural loads.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

DOE OSTI2026

Neural-Network Inverse Design of SRF Cavities and Transmons for Bosonic Quantum Computation

Three-dimensional superconducting radio-frequency (SRF) cavities provide exceptionally long-lived electromagnetic modes and, when coupled to nonlinear elements such as transmon qubits, become promising architectures for bosonic quantum information processing. The inverse design of such systems, i.e., recovering device geometries that produce specified electromagnetic and coupling targets, is generally a one-to-many problem. The qubit-cavity coupling strength depends sensitively on both the transmon geometry and its position within the cavity's electromagnetic field. As these systems scale up and their design parameter spaces grow, the cost of conventional iterative simulation becomes prohibitive. We present two deep neural network (DNN) approaches that address this inverse-design problem at complementary levels of the design stack. The first proposes SRF cavity geometries that produce target cavity observables. The second proposes transmon qubit designs that produce target qubit-cavity parameters - the coupling rate, qubit frequency, and anharmonicity $(g, ν_q, α)$. The recovered candidate designs match the targets to within ~5% (cavity) and ~2% (transmon), confirmed by end-to-end re-simulation. Both approaches map desired device behavior directly to candidate designs, a fast alternative to the iterative simulation studies usually required.

Yaker, Joseph [Fermilab; Northwestern U.]↗

DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗

DOE OSTI2025

Spin-on deposition of amorphous zeolitic imidazolate framework films for lithography applications

Amorphous zeolitic imidazolate framework (aZIF) films have been recently introduced as resists for electron beam and extreme ultraviolet lithography. aZIFs are also being considered for separation applications, including thin film membranes. However, the reported methods for aZIF deposition are currently based on highly empirical trial-and-error approaches that hinder control of film composition, thickness and uniformity as well as scale-up and transferability to different coating geometries. This work presents a method for depositing aZIF films with controllable thickness using dilute precursors mixed immediately before encountering the substrate. Importantly, the method is amenable to quantitative analysis by computational fluid dynamics to extract intrinsic deposition rates and limiting reactant transport diffusivities, enabling predictive physics-based modeling of the deposition process. This allows the deposition method to be adapted for spin coating on silicon wafers to prepare high-quality aZIF films with consistently controlled thickness. Using this approach, high-resolution resist performance and wafer-scale use for beyond extreme-ultraviolet lithography of aZIF films is demonstrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Transient Multiphysics Simulations with Pin Power Reconstruction in the Griffin Reactor Physics Code

This work introduces the pin power reconstruction capability available in the Griffin reactor physics code. This capability is implemented in an unstructured mesh framework, and the methods introduced are applied to the 2D SIMBA reactor core, which has assemblies and pins arranged in a hexagonal lattice. Since this reactor has a non-Cartesian geometry and also operates in the thermal spectrum, a general approach to pin power reconstruction is adopted, where SPH-based equivalence is leveraged to preserve assembly-wise reaction rates, while computing full-core form functions to preserve pin-wise fission production rates within the fuel pins of the reactor core. In a 2D microreactor benchmark problem, this pin power reconstruction approach was shown to reproduce pin powers compared to the Serpent2 Monte Carlo code for fixed temperature conditions and control drum rotation angles, yielding a core-wide RMS error level of 0.6\% and a maximum absolute pin error of 2.3\%. In addition, a tabulated library of multigroup cross sections, SPH factors, and form functions was generated to demonstrate the applicability of pin power reconstruction to a thermal feedback problem. Finally, a control drum transient was successfully simulated, showcasing the application of pin power reconstruction in a transient multiphysics feedback problem.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2025

LLRF System Analysis for the Fermilab PIP-II LINAC

Developing long-lived quantum processing units (QPUs) capable of supporting high-fidelity quantum operations is a crucial challenge on the path toward fault-tolerant quantum computing. TESLA-shaped superconducting RF (SRF) cavities, known for photon relaxation times on the order of seconds, provide an excellent foundation for 3D QPUs and quantum memory. This talk presents a novel design that leverages TESLA cavity modes coupled to ancillary transmon qubits, optimized to preserve coherence and control. By carefully engineering the package geometry, optimizing Hamiltonian parameters, and minimizing lossy participation ratios, we achieve photon relaxation times of over 16 ms and 20 ms for the two cavity modes, representing a significant improvement over previous multimode quantum memories. Despite the reduced coupling between the qubit and cavity modes, which is necessary to preserve long lifetimes, the platform supports robust and universal control schemes that are not limited by low coupling strength. We will also discuss how this architecture can lead to scalable, modular quantum computing systems.

Varghese, P. [Fermilab]↗

DOE OSTI2026

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗

DOE OSTI2025

Residual stress distribution in an additively manufactured complex structure by neutron diffraction measurement

Residual stress in an aerodynamically shaped Ni-based superalloy airfoil fabricated by laser powder bed fusion was measured by neutron diffraction. The experiment was conducted by considering the complex shape, implementing computer aided experiment planning, and automatic alignment at each rapid measurement. The 3-dimensional (3D) residual stress distribution in the airfoil is presented in this work, which lacks symmetry due to the complex geometry of the airfoil. In conclusion, the results provide theoretical thermal processing models a complete residual stress dataset of simulation validation on 3D shape complex structure.

Residual stress↗

DOE OSTI2024

The role of longitudinal decorrelations for measurements of anisotropic flow in small collision systems

Within a (3+1)D viscous hydrodynamic model we compute anisotropic flow in small system collisions as performed at the Relativistic Heavy Ion Collider and measured by the STAR and PHENIX Collaborations. We emphasize the importance of the rapidity dependence of the geometry for interpreting the differences encountered in measurements by the two collaborations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Design of Novel Hot Gas Component for Gas Turbine Engines Enabled by Materials and Additive Manufacturing Process Development

This CRADA project was the result of a project award under FOA-DOE-0001980. The overarching FOA project team consisted of researchers from Carpenter Technology Corporation (CTC), Solar Turbines Incorporated (Solar), Pennsylvania State University (PSU), University of California-Santa Barbara (UCSB), and Oak Ridge National Laboratory (ORNL). Evaluations were conducted on two high-γ’ superalloys that were designed by CTC and the UCSB. One alloy named GammaPrint-700 (GP-700) is a cobalt-base superalloy. The other alloy named GammaPrint-1100 (GP-1100) is a nickel-base (Ni-base) superalloy. PSU provided expertise and experimental testing of the thermal performance of AM micro-cooling architectures. ORNL provided expertise with the AM superalloy materials characterization and AM processing science. Solar provided turbine component design expertise. The focus of this CRADA report is to document the efforts between ORNL and CTC towards the development of superalloys designed for AM. The project goal was to use an AM processable high-temperature superalloy and design for Additive Manufacturing (DfAM) techniques to design an efficient turbine component (i.e. a turbine tip shoe) with enhanced cooling features that can only be fabricated through additive manufacturing (AM). The efficiencies of existing combined heat and power (CHP) engines are capped by both component design and materials limitations. However, AM of a tip shoe component from a γ’strengthened superalloy offers the design flexibility to increase the efficiency and power of an industrial gas turbine. This project brought about advancements in the DfAM tip shoe design space and in the area of high temperature superalloys processable through laser powder bed fusion (LPBF) AM. State of art computation design tools were utilized to optimize unique cooling features into a tip shoe component design. A two-prong materials development approach was taken to support development of the AM tip shoe geometry. The first approach centered on investigating the processability and the appropriate process science for the industry standard high-γ’ nickel-base (Ni-base) superalloy Mar-M247. This superalloy is typically cast and considered non-weldable by traditional welding standards. In the course of this work, the alloy was not deemed feasible for process scale-up due to significant cracking issues during printing. The second approach focused on the development and evaluation of a novel cobalt-base superalloy, GammaPrint™-700 (GP-700 and a Ni-base superalloy, GammaPrint™-1100 (GP-1100) designed to mitigate the significant AM processing issues with Mar-M247. The processability of these two alloys were investigated through electron beam melting (EBM) binder-jet AM (BJAM), and LPBF as a risk mitigation for manufacturability. To be considered a candidate material for down-selection to proceed to full-scale AM tip shoe engine testing trials, the high temperature creep rupture strength was required to achieve at a minimum, a Larsen Miller Parameter (LMP) increase of 10.9% over the baseline material LPBF AM Hastelloy X.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Design of Novel Hot Gas Component for Gas Turbine Engines Enabled by Materials and Additive Manufacturing Process Development

This CRADA project was the result of a project award under FOA-DOE-0001980. The overarching FOA project team consisted of researchers from Carpenter Technology Corporation (CTC), Solar Turbines Incorporated (Solar), Pennsylvania State University (PSU), University of California-Santa Barbara (UCSB), and Oak Ridge National Laboratory (ORNL). Evaluations were conducted on two high-γ’ superalloys that were designed by CTC and the UCSB. One alloy named GammaPrint-700 (GP-700) is a cobalt-base superalloy. The other alloy named GammaPrint-1100 (GP-1100) is a nickel-base (Ni-base) superalloy. PSU provided expertise and experimental testing of the thermal performance of AM micro-cooling architectures. ORNL provided expertise with the AM superalloy materials characterization and AM processing science. Solar provided turbine component design expertise. The focus of this CRADA report is to document the efforts between ORNL and CTC towards the development of superalloys designed for AM. The project goal was to use an AM processable high-temperature superalloy and design for Additive Manufacturing (DfAM) techniques to design an efficient turbine component (i.e. a turbine tip shoe) with enhanced cooling features that can only be fabricated through additive manufacturing (AM). The efficiencies of existing combined heat and power (CHP) engines are capped by both component design and materials limitations. However, AM of a tip shoe component from a γ’strengthened superalloy offers the design flexibility to increase the efficiency and power of an industrial gas turbine. This project brought about advancements in the DfAM tip shoe design space and in the area of high temperature superalloys processable through laser powder bed fusion (LPBF) AM. State of art computation design tools were utilized to optimize unique cooling features into a tip shoe component design. A two-prong materials development approach was taken to support development of the AM tip shoe geometry. The first approach centered on investigating the processability and the appropriate process science for the industry standard high-γ’ nickel-base (Ni-base) superalloy Mar-M247. This superalloy is typically cast and considered non-weldable by traditional welding standards. In the course of this work, the alloy was not deemed feasible for process scale-up due to significant cracking issues during printing. The second approach focused on the development and evaluation of a novel cobalt-base superalloy, GammaPrint™-700 (GP-700 and a Ni-base superalloy, GammaPrint™-1100 (GP-1100) designed to mitigate the significant AM processing issues with Mar-M247. The processability of these two alloys were investigated through electron beam melting (EBM) binder-jet AM (BJAM), and LPBF as a risk mitigation for manufacturability. To be considered a candidate material for down-selection to proceed to fullscale AM tip shoe engine testing trials, the high temperature creep rupture strength was required to achieve at a minimum, a Larsen Miller Parameter (LMP) increase of 10.9% over the baseline material LPBF AM Hastelloy X.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Simulated Evaluation of Powder Flowability Through a Partially Obstructed Consumable in Blown Powder Directed Energy Deposition Systems

Abstract In the interest of continued industrialization of metal additive manufacturing in modern production environments, cost is often referenced as a primary deterrent to new adopters. Conventional economic models for additive systems, processes, and supply chains often focus on specific process applications with little generalizability, or they neglect significant costs associated with production such as machine maintenance and consumable part replacement. Compounding the latter issue are substantial knowledge gaps in consumable part wear characterization for additive and other convergent manufacturing systems. In coaxial blown powder directed energy deposition systems, gas atomized metal powder is wasted during material deposition at a rate that is partly dependent on present wear phenomena in a consumable nozzle housed in the cladding head assembly. The price and lead time required to replace the nozzle incentivizes its reuse even when visibly worn. Often this initiates a process quality decline in the form of underbuilt geometry and internal defects due to losses in powder catchment efficiency. While depositing H13 steel using a hybrid manufacturing machine tool equipped with such a deposition system, a unique partial clog with a bridge-like structure formed at the consumable nozzle exit when supporting argon gas flows failed mid-process. To further understand coaxial multi-phase powder flow in the event of support gas failure, a computational fluid dynamics simulation is tailored to relevant process parameters, H13 powder material profile, and machine operator observations collected after the incident. The resulting differences in powder flow compared to control gas flow parameters is presented and discussed. The powder flowability and performance of the clogged nozzle is then assessed by using an optical profilometer to extract the profile of the clog and recreate the clog geometry within the simulation environment. In past work this simulation has been experimentally validated for a 316L steel powder material profile and used specifically for analyzing powder stream geometry and catchment efficiency. After the initial powder flow characterization, the clog is removed, and the nozzle is reprofiled. After removing the obstructing clog, the newly unobstructed nozzle geometry, the original off the shelf nozzle geometry, and additional nozzle profiles exploring different consumable refurbishment strategies are reevaluated in the simulation. Powder catchment efficiency for all variant nozzle geometries and relevant flow variables are compared and discussed, along with potential mitigation strategies for optimizing powder flowability with worn consumables. This work expands on the known morphology of blown powder obstructions and wear defects present in consumable coaxial nozzles while discussing pragmatic simulation driven responses to unanticipated subsystem failure in hybrid manufacturing machining platforms.

DeWitte, Lisa↗

DOE OSTI2025

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗