Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical simulation”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2025

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2024

The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equations

We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.

Bailo, Rafael (ORCID:0000000180183799)↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2025

The CHIMERAS project: design framework for the Collisionless HIgh-beta Magnetized Experiment Researching Astrophysical Systems

From the near-Earth solar wind to the intracluster medium of galaxy clusters, collisionless, high-beta, magnetized plasmas pervade our universe. Energy and momentum transport from large-scale fields and flows to small-scale motions of plasma particles is ubiquitous in these systems, but a full picture of the underlying physical mechanisms remains elusive. The transfer is often mediated by a turbulent cascade of Alfvénic fluctuations as well as a variety of kinetic instabilities; these processes tend to be multi-scale and/or multi-dimensional, which makes them difficult to study using spacecraft missions and numerical simulations alone. Meanwhile, existing laboratory devices struggle to produce the collisionless, high ion beta (β i ≳ 1), magnetized plasmas across the range of scales necessary to address these problems. As envisioned in recent community planning documents, it is therefore important to build a next generation laboratory facility to create a β i ≳ 1, collisionless, magnetized plasma in the laboratory for the first time. A working group has been formed and is actively defining the necessary technical requirements to move the facility towards a construction-ready state. Recent progress includes the development of target parameters and diagnostic requirements as well as the identification of a need for source-target device geometry. As the working group is already leading to new synergies across the community, we anticipate a broad community of users funded by a variety of federal agencies (including National Aeronautics and Space Administration, Department of Energy and National Science Foundation) to make copious use of the future facility.

astrophysical plasmas↗
DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

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

Revealing the Influence of Binding Motifs on Electron Transfer and Recombination Kinetics for CdSe Quantum Dots Functionalized with a Modified Viologen

Anchoring of molecules to the surfaces of semiconductor nanocrystals (NCs) presents an opportunity to leverage the precise synthetic tunability of molecular function and the remarkable light harvesting properties of NCs to drive photochemical reactions. However, charge transfer between the two species depends not only on the energy level alignments but also on the details of their binding interactions, which are difficult to probe. Here, in this work, we characterize the binding between CdSe quantum dots (QDs) and a new phosphonated derivative of the electron acceptor methyl viologen, designed to attach to the QD surface via the phosphonate group. We use isothermal titration calorimetry to probe the thermodynamics of the QD–molecule interaction and use the parameters determined therein to analyze transient absorption spectroscopy measurements of forward and back electron transfer from QDs to the viologen. We find that the ligand-like phosphonate binding leads to an electron-transfer rate constant that is 3 orders of magnitude smaller than that for the face-on binding of the bipyridine ring of methyl viologen. Back electron transfer is also significantly slower in the derivative. Interestingly, a minor fraction of the phosphonated derivative also binds in the face-on configuration, with similar forward and back electron transfer kinetics as methyl viologen. Numerical simulations show that the ligand-like binding will lead to significantly improved quantum yields of photocatalysis over a wide range of reaction rates. By independently characterizing binding thermodynamics and charge transfer kinetics, this work reveals how the complexities underlying electron transfer at the NC–molecule interface determine photocatalytic outcomes. This work also represents a step toward controlling forward and back electron transfer kinetics via rational molecular design.

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

Controlling Vibronic Coupling in Chlorophyll Proteins: The Effects of Excitonic Delocalization and Vibrational Localization

Vibrational-electronic (vibronic) coupling plays a critical role in excitation energy transfer in molecular aggregates and pigment-protein complexes (PPCs). But the interplay between excitonic delocalization and vibronic interactions is complex, often leaving even qualitative questions as to what conceptual framework (e.g., Redfield versus Forster theory) should be used to interpret experimental results. To shed light on this issue, we report here on the interplay between excitonic delocalization and vibronic coupling in site-directed mutants of the water soluble chlorophyll protein (WSCP), as reflected in 77 K fluorescence spectra. Experimentally, we find that in PPCs where excitonic delocalization is disrupted (either by mutagenesis or heterodimer formation), the relative intensity of the vibrational sideband (VSB) in fluorescence spectra is suppressed by up to 37% compared to the native protein. Furthermore, numerical simulations reveal that this effect results from the localization of high-frequency vibrations in the coupled system; while excitonic delocalization suppresses the purely electronic transition due to H-aggregate-like dipole-dipole interference, high-frequency vibrations are unaffected, leading to a relative enhancement of the VSB. By comparing VSB intensities of PPCs both in the presence and absence of excitonic delocalization, we extract a set of “local” Huang-Rhys (HR) factors for Chl a in WSCP. More generally, our results suggest a significant role for geometric effects in controlling energy-transfer rates (which depend sensitively on absorption/fluorescence line shapes) in molecular aggregates and PPCs.

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

Strong Surface-Enhanced Coherent Phonon Generation in van der Waals Materials

Terahertz (THz) coherent phonons have emerged as promising candidates for the next generation of high-speed, low-energy information carriers in atomically thin phononic or phonon-integrated on-chip devices. However, effectively manipulating THz coherent phonons remains a significant challenge. Here, in this study, we investigated THz coherent phonon generation in exfoliated van der Waals (vdW) flakes of Fe 3 GeTe 2 , Fe 5 GeTe 2 , and FePS 3 . We successfully generated the THz A 1g coherent phonon mode in these vdW flakes. An innovative approach involved partially exfoliating vdW flakes on a gold substrate and partially on a silicon (Si) substrate to compare the THz coherent phonon generation between both sides. Interestingly, we observed a significantly enhanced THz coherent phonon in the vdW/gold area compared with that in the vdW/Si area. Frequency-domain Raman mapping across the vdW flakes corroborated these findings. Numerical simulations further indicated a stronger enhanced surface field in vdW/gold structures than in vdW/Si structures. Consequently, we attribute the observed enhancement in THz coherent phonon generation to the increased surface field on the gold substrate. This enhancement was consistent across the three different vdW materials studied, suggesting the universality of this strategy. Our results hold promise for advancing the design of THz phononic and phonon-integrated devices.

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

Cavity-Assisted Coherent Phonon Generation and Control in a WSe 2 /Au Structure

Coherent phonons in the Terahertz (THz) regime have gained attention as potential candidates for next-generation high-speed, low-energy information carriers in atomically thin phononic or phonon-integrated on-chip devices. Nevertheless, achieving efficient control over THz coherent phonons continues to pose a considerable challenge. In this work, we explore THz coherent phonon generation in exfoliated van der Waals (vdW) flakes of WSe 2 on Au (WSe 2 /Au) and Si (WSe 2 /Si) using time-resolved pump-probe spectroscopy. The generation of THz coherent phonons was studied as a function of WSe 2 layer thickness and laser wavelength. Notably, a significant enhancement in THz coherent phonon generation was observed in the WSe 2 /Au structure, but only within specific ranges of WSe 2 thickness and laser wavelength. Further, the results from numerical simulations, which consider a self-hybridized optical cavity depending on WSe 2 thickness, along with optical reflectance and Raman spectroscopy measurements, aligned well with the time-domain observations of THz coherent phonon generation. We propose that the observed enhancement in THz coherent phonon generation is strongly influenced by light-matter interaction in the WSe 2 cavity, a mechanism that may be applicable to a broader range of vdW materials. These findings offer promising insights for the development of THz phononic or phonon-integrated devices.

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

Spin Polarization from Circularly Polarized Light Induced Charge Transfer

Here, we show how a spin polarization can be generated through the photoinduced electron transfer of an achiral donor–acceptor complex following chiral light excitation. In particular, we illustrate the basic energetic and symmetry requirements for chirality induced spin selectivity where the chirality emerges from the electronic degrees of freedom following excitation with circularly polarized light. We study this effect in a simple model of a metalloporphyrin complex with an axial acceptor ligand using quantum mechanical rate theories and numerical simulations. We find that the spin polarization emerges due to the selective excitation of a ring current within the porphryin, breaking the degeneracy of the two degenerate spin states. The resultant spin polarization increases with the spin orbit coupling between the metal in the porphyrin and the axial ligand, and is transient, with a lifetime dependent on the rate of dephasing from the Jahn–Teller distortion mode. This proposed effect should be observable in spin-resolved photoemission spectroscopy.

charge transfer↗
DOE OSTI2024

Automated Gold Nanorod Spectral Morphology Analysis Pipeline

The development of a colloidal synthesis procedure to produce nanomaterials with high shape and size purity is often a time-consuming, iterative process. This is often due to quantitative uncertainties in the required reaction conditions and the time, resources, and expertise intensive characterization methods required for quantitative determination of nanomaterial size and shape. Absorption spectroscopy is often the easiest method for colloidal nanomaterial characterization. However, due to the lack of a reliable method to extract nanoparticle shapes from absorption spectroscopy, it is generally treated as a more qualitative measure for metal nanoparticles. This work demonstrates a gold nanorod (AuNR) spectral morphology analysis tool, called AuNR-SMA, which is a fast and accurate method to extract quantitative structural information from colloidal AuNR absorption spectra. To demonstrate the practical utility of this model, we apply it to three distinct applications. First, we demonstrate this model's utility as an automated analysis tool in a high-throughput AuNR synthesis procedure by generating quantitative size information from optical spectra. Second, we use the predictions generated by this model to train a machine learning model to predict the resulting AuNR size distributions under specified reaction conditions. Third, we apply this model to spectra extracted from the literature where no size distributions are reported and impute unreported quantitative information on AuNR synthesis. This approach can potentially be extended to any other nanocrystal system where absorption spectra are size dependent, and accurate numerical simulation of absorption spectra is possible. In addition, this pipeline could be integrated into automated synthesis apparatuses to provide interpretable data from simple measurements, help explore the synthesis science of nanoparticles in a rational manner, or facilitate closed-loop workflows.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Single-Molecule Fluorescence Microscopy Reveals Energy Transfer Active versus Inactive Nanocrystal/Dye Conjugate Pairs

Defect-mediated energy transfer (EnT) is a radiative process that occurs between donor defect states in the forbidden bandgap of semiconductor nanocrystals (NCs) and dye molecules bound to their surfaces. The EnT efficiency depends on the number of dye molecules attached to each NC, the donor–acceptor distance, and the dipole orientation factor between the donor and acceptor, all of which vary across all individual NCs in a sample. While ensemble-level fluorescence spectroscopy measurements have provided average values for donor–acceptor distances, dye-to-NC ratios, and EnT rate constants, questions remain about the impact of donor/acceptor heterogeneity on observed EnT efficiencies. Notably, ensemble-level measurements cannot distinguish between bare NCs and EnT-active versus inactive NC/dye pairs in the same sample batch, limiting the ability to design systems with 100% EnT efficiency. To address this, we studied defect-mediated EnT between AlexaFluor 555 dye acceptors chemically bound to ZnO NC donors at the level of single molecules and single NCs. Interestingly, 20% of bound NC/dye pairs are EnT-inactive, likely contributing to residual defect photoluminescence (PL) observed in ensemble-level measurements and reducing overall EnT efficiency. Single particle-level ZnO defect PL and acceptor fluorescence trajectories exhibited distinct microfluctuations, which are absent in bare ZnO NCs. We hypothesized that our observations can be explained with a competitive dye fluorescence quenching pathway, possibly due to charge transfer between the excited state dye and the ZnO NC. Numerical simulations of single-molecule PL traces for this scenario produced microfluctuations consistent with the experimental results. These findings highlight the impact of sample heterogeneity on EnT processes and provide insights for designing light-harvesting systems with optimized EnT efficiency.

25 ENERGY STORAGE↗
DOE OSTI2026

Atomic Alignment in PbS Nanocrystal Superlattices with Compact Inorganic Ligands via Reversible Oriented Attachment of Nanocrystals

Nanocrystals (NCs) serve as versatile building blocks for the creation of functional materials with NC self-assembly offering opportunities to enable novel material properties. Here, in this work, we demonstrate that PbS NCs functionalized with strongly negatively charged metal chalcogenide complex (MCC) ligands, such as Sn 2 S 6 4– and AsS 4 3– , can self-assemble into all-inorganic superlattices with both long-range superlattice translational and atomic-lattice orientational order. Structural characterizations reveal that the NCs adopt an unexpected edge-to-edge alignment, and numerical simulation clarifies that orientational order is thermodynamically stabilized by many-body ion correlations originating from the dense electrolyte. Furthermore, we show that the superlattices of Sn 2 S 6 4– -functionalized PbS NCs can be fully disassembled back into the colloidal state, which is highly unusual for orientationally attached superlattices with atomic-lattice alignment. The reversible oriented attachment of NCs, enabling their dynamic assembly and disassembly into effectively single-crystalline superstructures, offers a pathway toward designing reconfigurable materials with adaptive and controllable electronic and optoelectronic properties.

Landau-Ginzburg theory↗
DOE OSTI2025

DeFault: DEep‐Learning‐Based FAULT Delineation Using the IBDP Passive Seismic Data at the Decatur CO2 Storage Site

Abstract The carbon capture, utilization, and storage (CCUS) framework is an essential component in reducing greenhouse gas emissions, with its success hinging on the comprehensive knowledge of subsurface geology and geomechanics. Passive seismic event relocation and fault detection offer vital insights into subsurface structures and the ability to monitor fluid migration pathways. Accurate identification and localization of seismic events, however, face significant challenges, including the necessity for high‐quality seismic data and advanced computational methods. To address these challenges, we introduce a novel deep learning method, , specifically designed for passive seismic source relocation and fault delineating for passive seismic monitoring projects. By leveraging data domain‐adaptation, allows us to train a neural network with labeled synthetic data and apply it directly to field data. Using , the passive seismic sources are automatically clustered based on their recording time and spatial locations, and subsequently, faults and fractures are delineated accordingly. We demonstrate the efficacy of on a field case study involving injection related microseismic data from Decatur, Illinois area. Our approach accurately and efficiently relocated passive seismic events, identified faults and could aid in potential damage induced by seismicity. Our results highlight the potential of as a valuable tool for passive seismic monitoring, emphasizing its role in ensuring CCUS project safety. This research bolsters the understanding of subsurface characterization in CCUS, illustrating machine learning’s capacity to refine these methods. Ultimately, our work has significant implications for CCUS technology deployment, an essential strategy in combating climate change. Plain Language Summary In our quest to tackle climate change, we use a strategy known as carbon capture, utilization, and storage (CCUS) to keep greenhouse gases out of the atmosphere. This strategy relies heavily on our ability to understand what's happening deep under the earth's surface. To make sure we store super critical safely, we need to accurately map out the geological structure, especially faults, but this is tough without high‐quality data and complex computer programs. We've developed a new tool called “DeFault,” which uses advanced machine learning to improve how we find and map these underground features. “DeFault” is smart enough to learn from numerically simulated data and then apply what it’s learned to real‐world situations. It groups together seismic activity—tiny tremors and shifts in the earth—based on when and where they happen, which helps us spot where there might be cracks or faults. We tested “DeFault” in Illinois, where CO 2 is injected underground, and it successfully pinpointed where these tremors occurred and mapped out the faults, helping to prevent accidents accurately in the future. Our study shows that “DeFault” will be a powerful ally in making CCUS safer and more effective, especially for the Illinois Basin Decatur Project. Key Points Faults and fractures introduced by carbon storage can be monitored by passive seismicity DeFault algorithm enables an automatic process for accurate and efficient passive seismic event locating and clustering

58 GEOSCIENCES↗