Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initialization method”

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

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗
DOE OSTI2025

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Roadmap for a Lightning Modeling Grand Challenge

This document is a roadmap for building an interconnected model of the physical processes that produce a lightning discharge, and its observable optical and radio signals. We call this a Lightning Modeling Grand Challenge, recognizing that significant effort and coordination of human and financial resources is required to realize the capability. The roadmap serves to outline the coordination of resources necessary to enable stitching together existing knowledge and model components to make a lightning prediction, and to test these predictions with observations. Such a capability does not currently exist. The roadmap is motivated not only by a spirit of scientific inquiry, but by practical challenges faced by US Federal and societal stakeholders. Advancements in lightning observations have outpaced our tests of integrated understanding, leaving many stakeholders unsure how to design their missions to properly detect and discriminate lightning, and unsure how to apply the sometimes-disagreeing lightning signals from diverse instruments. The time is right to connect existing theories and models to support stakeholders in understanding the signals they observe, for needs as diverse as climate monitoring, national security, weather forecasting, public safety, and protection of natural and built environments. The roadmap’s two main technical sections describe the components of a linked physical model, followed by a description of models of lightning signals and sensors that are driven by outputs from the physical model. The goal is to predict the time-varying physical properties of lightning that are self-consistent with the thunderstorm’s structure and dynamics. These lightning signals then propagate through the storm, with realistic dispersion and attenuation, to receivers on the ground or in space. At a high level, the model begins with weather (cloud) model output, including explicit prediction of the electrification of cloud particles. The cloud’s electrical structure drives a model of lightning physics, from initiation, through channel development, and discharges along those channels. Key lightning parameters, such as the temperature and currents in the channel, and their space and time distribution, are then used to produce optical and electromagnetic signal sources that propagate to modeled receivers. This architecture therefore generates a dataset suitable for comparison to existing and envisioned observing systems. The need for additional measurements and field campaigns to support model development is described. In each model sub-component, inputs, outputs, uncertainties, evaluation methods, and next steps are summarized, interleaved with references to the scientific literature. Identifying boundaries between the model sub-components aids in segmenting an integrated, complex model into practical work packages and system sub-components, allowing a diverse team to contribute and maintain the system. We estimate that at least five years of effort and a $\$$10M initial investment is necessary to make a significant step forward. Mechanisms to facilitate community coordination, including annual workshops and open-source code repositories, are described.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Ductility of 49Fe-49Co-2V Soft Magnetic Alloy Bar: Surface Effects and Test Methods

The tensile ductility of 49Fe-49Co-2V (Hiperco® 50A) bar was investigated in both as-received and heat-treated conditions. The as-received/machined specimens exhibit very low ductility compared to samples where heat treatment was the final step prior to testing. Microstructural characterization showed that internal residual strain from bar processing and, most importantly, surface machining damage, cause lower elongation in the as-received material. Because fracture of this intermetallic alloy initiates at the surface, it is particularly susceptible to surface machining damage, i.e., the near-surface region has already exhausted most of its ability to accumulate tensile strain. During heat treatment, the internal residual strain and near-surface machining damage are eliminated and ductility is improved, despite a higher degree of crystallographic ordering in the heat-treated condition (which typically lowers ductility). Furthermore, if machining is again performed after heat treatment, the material again exhibits brittle behavior, even with only light touch-up machining passes. Here, in this work, methods of tensile strain measurement were investigated, namely conventional knife-edge extensometry and noncontact digital image correlation (DIC) on heat-treated material. For clip-on knife-edge extensometry, the range of failure strain was 2.5-5.5% for heat-treated Hiperco. For noncontact methods, ductility up to 7% was observed. The results highlight the tendency for the alloy to fail at surface imperfections, even those produced by application of the extensometer itself. Noncontact laser extensometry is recommended for determining the intrinsic ductility of the alloy. A method of laser surface modification was developed which increased ductility by ~ 100% compared to unmodified samples. The high cooling rates achieved during laser surface processing can bypass the ordering reaction and produce a ductile disordered structure at the surface that exhibits ductile fracture characteristics.

EBSD↗
DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Pore-Selective Fullerene Loading in a Porphyrin-Based Metal–Organic Framework Controls Photoinduced Charge-Transfer Dynamics

Building porous donor−acceptor networks based on host−guest interactions in metal organic frameworks (MOFs) provides unique opportunities for tuning charge separation in highly tailorable materials. Here we focus on installing electron-rich porphyrins and electron-deficient fullerene derivatives in the PCN-222 MOF using a solvent assisted ligand insertion (SALI) method. The fullerene is primarily bound in the large pore, where it is subject to distinct dielectric environments through dimethylformamide (DMF) and 1,4-dioxane solvent exchange. Following photoexcitation, sub-picosecond charge transfer involving initial exciplex population is observed, with different charge recombination pathways and lifetimes depending on solvent polarity through modulation of charge-transfer state energies. While the 1,4-dioxane environment yields charge recombination within 1 ns via local fullerene and porphyrin triplet state population, DMF results in charge recombination directly to the ground state on much longer time scales, including some lifetime components in the microsecond range. Fullerene loading influences these kinetics, and the potential for charge delocalization due to fullerene aggregation within the pores is evaluated by using molecular dynamics simulations.

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

Data-driven design of electrolyte additives supporting high-performance 5 V LiNi 0.5 Mn 1.5 O 4 positive electrodes

LiNi 0.5 Mn 1.5 O 4 (LNMO) is a high-capacity spinel-structured material with an average lithiation/de-lithiation potential at ca. 4.6–4.7 V vs Li + /Li, far exceeding the stability limits of electrolytes. An efficient way to enable LNMO in lithium-ion batteries is to reformulate an electrolyte composition that stabilizes both graphitic (Gr) negative electrode with solid-electrolyte-interphase and LNMO with cathode-electrolyte-interphase. In this study, we select and test a diverse collection of 28 single and dual additives for the Gr||LNMO battery system. Subsequently, we train machine learning models on this dataset and employ the trained models to suggest 6 binary compositions out of 125, based on predicted final area-specific-impedance, impedance rise, and final specific-capacity. Such machine learning-generated new additives outperform the initial dataset. This finding not only underscores the efficacy of machine learning in identifying materials in a highly complicated application space but also showcases an accelerated material discovery workflow that directly integrates data-driven methods with battery testing experiments.

batteries↗
DOE OSTI2025

Criticality analysis of nuclear binding energy neural networks

Machine learning methods, in particular deep learning methods such as artificial neural networks (ANNs) with many layers, have become widespread and useful tools in nuclear physics. However, these ANNs are typically treated as ‘black boxes’, with their architecture (width, depth, and weight/bias initialization) and the training algorithm and parameters chosen empirically by optimizing learning based on limited exploration. We test a non-empirical approach to understanding and optimizing nuclear physics ANNs by adapting a criticality analysis based on renormalization group flows in terms of the hyperparameters for weight/bias initialization, training rates, and the ratio of depth to width. This treatment utilizes the statistical properties of neural network initialization to find a generating functional for network outputs at any layer, allowing for a path integral formulation of the ANN outputs as a Euclidean statistical field theory. We use a prototypical example to test the applicability of this approach: a simple ANN for nuclear binding energies. We find that with training using a stochastic gradient descent optimizer, the predicted criticality behavior is realized, and optimal performance is found with critical tuning. However, the use of an adaptive learning algorithm leads to somewhat superior results without concern for tuning and thus obscures the analysis. Nevertheless, the criticality analysis offers a way to look within the black box of ANNs, which is a first step towards potential improvements in network performance beyond using adaptive optimizers.

artificial neural network↗
DOE OSTI2024

Nontrivial fusion of Majorana zero modes in interacting quantum-dot arrays

Motivated by recent experimental reports of Majorana zero modes (MZMs) in quantum-dot systems at the “sweet spot,” where the electronic hopping t ℎ is equal to the superconducting coupling Δ, we study the time-dependent spectroscopy corresponding to the nontrivial fusion of MZMs. The term “nontrivial” refers to the fusion of Majoranas from different original pairs of MZMs, each with well-defined parities. We employ an experimentally accessible time-dependent real-space local density-of-states (LDOS) method to investigate the nontrivial MZM fusion outcomes in canonical chains and in a Y-shaped array of interacting electrons. In the case of quantum-dot chains where two pairs of MZMs are initially disconnected, after fusion we find equal-height peaks in the electron and hole components of the LDOS, signaling nontrivial fusion into both the vacuum I and fermion Ψ channels with equal weight. For π-junction quantum-dot chains, where the superconducting phase has opposite signs on the left and right portions of the chain, after the nontrivial fusion we observed the formation of an exotic two-site MZM near the center of the chain, coexisting with another single-site MZM. Furthermore, we also studied the fusion of three MZMs in the Y-shaped geometry. In this case, after the fusion we observed the novel formation of another exotic multisite MZM, with properties depending on the connection and geometry of the central region of the Y-shaped quantum-dot array.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

An Open Benchmark of One Million High-Fidelity Cislunar Trajectories

Cislunar space spans from geosynchronous altitudes to beyond the Moon and will underpin future exploration, science, and security operations. We describe and release an open dataset of one million numerically propagated cislunar trajectories generated with the open-source Space Situational Awareness Python package (SSAPy). The model includes high-degree Earth/Moon gravity, solar gravity, and Earth/Sun radiation pressure; other planetary gravities are omitted by design for computational efficiency. Initial conditions uniformly sample commonly used osculating-element ranges, and each trajectory is propagated for up to six years under a single, fixed start epoch. The dataset is intended as a reusable benchmark for method development (e.g., space domain awareness, navigation, and machine-learning pipelines), a reference library for statistical studies of orbit families, and a starting point for community-driven extensions (e.g., alternative epochs). We report empirically observed stability trends (e.g., a band near ~5 GEO and persistence of some co-orbital classes including L4/L5 librators) as dataset descriptors rather than new dynamical results. The chief contribution is the scale, fidelity, organization (CSV/HDF5 with full state time series and metadata), and open availability, which together lower the barrier to comparative and data-driven studies in the cislunar regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dashboard for Marine Energy Site Assessment and Monitoring

The marine energy (ME) industry presently relies upon fragmented site assessment solutions that require high resource expenditure for deployment at each site and do not leverage the wealth of readily available tools and information. A wave energy resource assessment dashboard, currently in development, will substantially improve siting, permitting, operations, and maintenance of ME projects by providing an integrated solution that is a one-stop-shop for a developer’s needs. The Site Energy Assessment and MOnitoring Dashboard (SEAMOD) will be of commercial interest to anyone seeking to deploy an ME project and is easily expandable to include tidal and wind energy site assessments. The integrated dashboard is being developed using state-of-the-art database and cloud computing methods and data-assimilative modeling tools that can be coupled with low-cost, rapidly deployable wave buoys and environmental sensing hardware. The combined software and hardware dashboard will reduce wave energy site characterization and wave climate monitoring costs by more than 60 percent and provide assessments that meet international industry standards. To realize a thriving global ME industry, the physical environment at a potential deployment site must be understood, not only for resource characterization, but also for optimization of device and power conversion performance. SEAMOD directly addresses these needs with a commercially marketable product. SEAMOD is a low-cost solution that provides comprehensive ME resource assessments, baseline environmental monitoring, and offshore characterizations required for successful ME development. The key technical objectives for Phase I were a series of software development goals, which when implemented with monitoring solutions, produced an initial proof-of-concept low-cost wave energy resources dashboard. In Phase II, the development of the prototype SEAMOD continued. The basic framework employed was the development of a revised dashboard and monitoring tool customized for ME applications by focusing on IEC site assessment and method requirements. Development was focused on the integration of full hindcast metocean products to provide hindcast resource characterization and environmental information. The final integrated dashboard provides a low-cost solution that delivers comprehensive, scalable, industry-standard energy resource assessments and offshore characterizations required for successful ME development. The integrated dashboard offers visibility of the most recent site modeling, measurements, and historical data. The application and integration of consensus-based standards for wave energy resource assessment, as determined by the International Electrotechnical Commission (IEC), are crucial for the impact and value of SEAMOD. SEAMOD includes monthly, seasonal, and yearly statistics, as well as the total 30-year record, offering temporal resolution of the IEC parameters to aid potential developers in determining the available wave energy resources in their area of interest.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Synthesizing Highly Crystalline Graphite Powder from Bulk Polyethylene Waste for Lithium-Ion Battery Anodes

Upcycling plastic waste into graphite can potentially be used, in conjunction with other methods, to manage existing waste materials and diversify graphite supply chains. However, synthesizing large quantities of crystalline graphite powder from plastic waste, particularly polyethylene (PE), remains a challenge because PE decomposes into light gases during thermal processing and simple methods do not exist at any appreciable size scale to address this challenge. In this work, a method is developed for air processing bulk forms of PE waste to create stable carbon char that does not readily decompose during high-temperature processing. This method employs solid additives in the form of salts, which are combined with the PE melt during air processing to increase the effective surface area of the melt and improve the oxygen-driven chemistry that stabilizes PE for high-temperature processing. After removal of the solid salt additives from the PE-derived char, it is converted into a highly crystalline bulk graphite powder using an Fe-based catalytic process. The PE-derived graphite anode in a lithium-ion coin cell showed a specific capacity of 345 mAh/g at 0.05C with an initial Coulombic efficiency of 87% and reversible capacity retention of ~100% at different current rates. It also showed a specific capacity of up to 313 mAh/g at 0.5 discharge/charge cycles per hour (0.5C) and Coulombic efficiency of 99.9% after 250 cycles, indicating excellent electrochemical performance as an anode material for lithium-ion batteries. This method illustrates that there are opportunities for upcycling large quantities of PE waste to produce graphite powders suitable for use in LIBs.

25 ENERGY STORAGE↗
DOE OSTI2025

Performance Evaluation of LPBF Manufactured 316H Components

This work represents the continuation of a benchmark study that includes modeling, fabrication and characterization as demonstration to support industry’s adoption of advanced manufacturing processes in a variety of structures. This comprehensive study investigated the feasibility of using additive manufacturing (AM) technologies, specifically Laser Powder Direct Energy Deposition (LP-DED) and Laser Powder Bed Fusion (LPBF), to produce complex nuclear microreactor components using 316H stainless steel. The research focused on manufacturing an expanded elbow pipe component with transitioning sections, which are traditionally difficult and costly to produce through conventional manufacturing methods. The overall study’s primary objectives are therefore demonstrating AM viability for nuclear applications, optimizing process parameters, developing comprehensive material characterization protocols, validating computational modeling approaches, and establishing manufacturing guidelines for complex geometries. Although the initial work included the phased approach of cubical, upscaled cylindrical components, it is to enable to obtain more knowledge for the printing of the expanded elbow structure. The project achieved significant progress in process development by successfully optimizing LP-DED parameters to achieve 99.16-99.97% relative density in 316H stainless steel components. Through systematic evaluation of sixteen cube samples with varied laser powers (400-700W) and scan speeds (600-900 mm/min), optimal processing windows were identified at 500-550W with 600-700 mm/min or 650-700W with 650-900 mm/min scan speeds. The DED manufactured 316H demonstrated mechanical properties comparable or superior to wrought materials, with Young's modulus ranging from 153-208 GPa and controlled microstructural characteristics including greater than 95% face-centered cubic (FCC) phases and engineered cellular structures with sizes between 3.23-6.17 µm.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Diurnal evolution of non-precipitating marine stratocumuli in a large-eddy simulation ensemble

Abstract. We explore the cloud system evolution of non-precipitating marine stratocumuli with a focus on the impacts of the diurnal cycle and free-tropospheric (FT) humidity based on an ensemble of 244 large-eddy simulations generated by perturbing initial thermodynamic profiles and aerosol conditions. Cases are categorized based on their degree of decoupling and the cloud liquid water path (LWPc, based on model columns with cloud optical depths greater than one). A budget analysis method is proposed to analyze the evolution of cloud water in both coupled and decoupled boundary layers. More coupled clouds start with a relatively low LWPc and cloud fraction (fc) but experience the least decrease in LWPc and fc during the daytime. More decoupled clouds undergo greater daytime reduction in LWPc and fc, especially those with higher LWPc at sunrise because they suffer from faster weakening of net radiative cooling. During the nighttime, a positive correlation between FT humidity and the LWPc emerges, consistent with higher FT humidity reducing both radiative cooling and the humidity jump, both of which reduce entrainment and increase LWPc. The LWPc is more likely to decrease during the nighttime for a larger LWPc and greater inversion base height (zi), conditions under which entrainment dominates as turbulence develops. In the morning, the rate of the LWPc reduction depends on the LWPc at sunrise, zi, and the degree of decoupling, with distinct contributions from subsidence and radiation.

Chen, Yao-Sheng (ORCID:0000000208354132)↗
DOE OSTI2025

Initial Testing of a Small Gamma Detector (GR1+) for Detection of Holdup in HEPA Filters

Measuring and assessing holdup is a key component of proper inventorying of a Material Balance Area (MBA). The Nuclear Material Control and Accountability (NMC&A) group at the Savannah River Site is responsible for ensuring that the methods used to determine holdup in the MBAs are valid for the conditions of that MBA. Holdup can be contained inside processing equipment or inside ventilation systems that that pulls air out of the glovebox through HEPA filters and stainless steel piping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Durable Thin‐Film Porous Transport Electrodes for High Current Density PEM Water Electrolysis

Proton exchange membrane water electrolyzers rely on relatively expensive Ir-based catalysts for efficient and durable hydrogen production. To reduce system costs, Ir loadings can be reduced if performance and durability are maintained. Sputter deposition is a readily scalable method to synthesize uniform, low-loading catalyst layers with controlled composition. A catalyst applied directly to the porous transport layer can have advantages for performance, manufacturing simplicity, and catalyst recovery. Suitable porous transport layer porosity can minimize activity losses when reducing loadings. Here, methods are presented to deposit metallic Ir as well as amorphous and rutile Ir oxides. The activity and durability of these materials in the porous transport electrode architecture is evaluated. The metallic and amorphous forms have better initial activity, however, operation at 3 A cm −2 with 0.1 mg Ir cm −2 shows that only rutile IrO 2 maintains performance beyond 100 h with a 50 mV improvement after 700 h. A >10x reduced dissolution rate is shown for rutile IrO 2 . With a low-porosity transport layer and 0.4 mg Ir cm −2 , a steady-state voltage decay rate of 6 µV h −1 is achieved. The results demonstrate that sputter-deposited rutile IrO 2 porous transport electrodes with low Ir loading can be operated at high current density to reduce hydrogen production costs.

08 HYDROGEN↗
DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗
DOE OSTI2024

Temporal and Spatial Evolution of Non-Elastic Strain Accumulation in Stanstead Granite During Brittle Creep

Understanding the long-term behavior of brittle rocks requires fundamental consideration of time-dependent strain evolution and brittle creep processes. Previous studies have evaluated sub-critical crack growth during time-dependent deformation and damage evolution in brittle rocks; however, there is an incomplete knowledge of how damage evolves spatially and temporally within the body of intact rocks, where distributed regions of damage interact and coalesce during creep. This paper presents laboratory research focusing on evaluating brittle creep damage processes in Stanstead granite (SG) using 2-dimensional digital image correlation (2D-DIC). In the laboratory, the prismatic SG specimens were loaded beyond an estimated Crack Damage stress threshold (CD) level and then maintained a constant stress to initiate the creep process. DIC was used to characterize full-field spatiotemporal strain evolution, which was then interpreted in the context of local regions of “damage”, determined according to a strain-based criterion. Here, a method was proposed for identifying “existing” and “new” damage regions over specified intervals during the test, followed by spatial clustering of these regions to assess their spatiotemporal evolution. The clustering analysis results demonstrated the extension of existing damage regions was the main damage process during brittle creep, which is consistent with existing models of sub-critical crack growth. In addition, temporal analysis of tensile and shear strains on a point-by-point basis revealed both new damage formation and the strain concentration within existing damaged regions significantly contribute to overall specimen strain during primary creep. In contrast, during secondary creep, increases in specimen deformation are influenced by the accumulation of strains within already damaged regions.

58 GEOSCIENCES↗