Search NASA⌕ Search

SEARCH · Search NASA

Results for “force matrix”

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

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Long-Term Aging in Miscible Polymer Nanocomposites

Here we find that the initial, solvent-cast state of nanoparticles (NPs) in a polymer matrix temporally evolves during thermal annealing such that, at steady state, NPs maximize their distance from each other subject to mass balance constraints. The observed timescales for this unexpected structural reorganization, as probed by small-angle X-ray scattering, are temperature-dependent and can be prohibitively large, especially at temperatures around and below 1.2Tg. X-ray photon correlation spectroscopy measurements during reorganization reveal that the collective NP dynamics slow down with annealing at constant temperature; this is accompanied by changes in the low-frequency regime in macroscopic viscoelastic measurements in equilibrated materials. By ruling out other potential sources for these effects (i.e., electrostatic interactions, adsorbed layers), we attribute these results to a long-ranged repulsive force between the NPs caused by fluctuations in the polymer phase, i.e., the "anti-Casimir" effect proposed by Obhukhov and Semenov [Long-range interactions in polymer melts: The anti-Casimir effect. Phys Rev Lett 2005, 95 (3), 038305]. Thus, our results highlight the important role of long-term, slow NP reorganization on the structure and, subsequently, the properties of polymer nanocomposites (PNCs), even in the case of nominally miscible polymer nanoparticle hybrids.

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

Towards underwater additive manufacturing via additive friction stir deposition

Given the challenges in feed material supply and quality control, metal additive manufacturing has rarely been implemented in austere environments, especially underwater. This paper explores the underwater operation potential of an emerging solid-state additive technology: additive friction stir deposition, wherein material feeding and bonding are enabled by mechanical forces with minimal influences from water. It is demonstrated that additive friction stir deposition of 304 stainless steel can be successfully performed with the print head and substrate immersed in water. High temperature is reached in the deposition zone (>60% melting temperature); the material deposition behavior is similar to that in typical open-air operation. The as-deposited material is fully-dense, having fewer annealing twins and a substantially smaller grain size than the feed material (4.98 μm vs. 31.44 μm ). Such microstructural changes stem from dynamic recrystallization caused by the large strain and high temperature introduced during deposition. In addition to grain refinement, small equiaxed dispersoids (~2–3 μm or less) are formed and evenly distributed in the austenite steel matrix. Rich in Cr, Mn, and O, these particles likely result from the reaction between the elements in stainless steel and water at elevated temperatures.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Relativistic spin dynamics for vector mesons

We propose a relativistic theory for spin density matrices of vector mesons based on Kadanoff-Baym equations in the closed-time-path formalism. The theory puts the calculation of spin observables such as the spin density matrix element ρ 00 for vector mesons on a solid ground. Within the theory we formulate ρ 00 for ϕ mesons into a factorization form in separation of momentum and spacetime variables. We argue that the main contribution to ρ 00 at lower energies should be from the ϕ fields that can polarize the strange quark and antiquark in the same way as electromagnetic fields. The key observation is that there is correlation inside the ϕ meson wave function between the ϕ field that polarizes the strange quark and that polarizes the strange antiquark. This is reflected by the fact that the contributions to ρ 00 are all in squares of fields that are nonvanishing even if the fields may strongly fluctuate in spacetime. The fluctuation of strong force fields can be extracted from ρ 00 of unflavored vector mesons as links to fundamental properties of quantum chromodynamics. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2025

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗
DOE OSTI2020

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Surviving nanoscale interfacial stability in extreme thermal expansion contrast Zn(CN) 2 -epoxy resin matrix composites

Here, we report the remarkable stability at the nanoscale matrix-filler interface in a series of overall low coefficient of thermal expansion (CTE) zinc cyanide (ZC)/epoxy resin composites. These interfaces demonstrate the stability of highly contrasting materials. These composites are designed to optimize the properties of both materials: epoxy resins are valued for their high strength, superb chemical resistance, low dielectric properties, and adhesive abilities. ZC, a model flexible framework, exhibits negative thermal expansion (NTE) behavior and compensates for the epoxy resin’s detrimentally large thermal expansion values. Filled resins minimize the native 65 ppm/°C CTE values to 19 ppm/° C over a wide temperature range while maintaining linearity, yielding composites compatible with typical metal and ceramic substrates’—commonly bonded to the resin— 0–20 ppm/° C CTE values. Additionally, with significant intrinsic differences in the materials, we use wide-angle X-ray diffraction and atomic force microscopy to elucidate thermal movement. The ZC/epoxy interface is surprisingly resilient over 1,000 cycles between -55° C and 75° C, extremes typical to aerospace applications. Further, we show that the resin viscosity remains less than 21 Pa*s at up to 30 vol% ZC, an important consideration for workability. We also discuss settling distribution, glass transition temperature, flexural strength, and longevity. All results point favorably towards practical implementation in low thermal expansion application needs: filling the gap in component design.

42 ENGINEERING↗
DOE OSTI2021

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Block Length on the Properties of Multiblock Polysulfone-Poly(diallylpiperidinium hydroxide) Anion Exchange Membranes

Research in anion exchange membranes (AEM)s has continued with the development of materials bearing base stable cations. Designing AEMs that microphase separate into hydroxide conductive, hydrophilic domains within a hydrophobic and mechanically robust matrix has been shown to be successful for improving AEM performance. A series of multiblock polysulfone-poly(diallylpiperidinium hydroxide) copolymers (PSf-PDApipOH) of similar hydrophobic/hydrophilic composition was prepared in which the molecular weight of the hydroxide conducting PDApipOH segments was varied. The variable hydrophilic segment molecular weight was designed to assess the impact on microphase separation, hydroxide conductivity, and water management. The multiblock copolymers investigated were prepared by condensation polymerization of preformed4-fluorophenyl sulfone terminated poly(diallylpiperidinium hexafluorophosphate) (PDApipPF6) oligomers with polysulfone monomers. The structure-property relationship between the molecular weight of the conductive PDApipOH segments and AEM performance was demonstrated by evaluation of the microphase separation, water uptake, and hydroxide conductivity. Further, membranes fabricated from the polysulfone-poly(diallylpiperidinium hexafluorophosphate) (PSf-PDApipPF6) multiblock copolymers were shown to form well-connected conductive domains by SAXS and atomic force microscopy experiments. The PSf-PDApipOH membranes were highly conductive with the maximum hydroxide conductivity reaching 62.9 mS∙cm -1 at 60°C and 95% relative humidity. Furthermore, it was demonstrated that the conductivity increased with increasing PDApipOH segment molecular weight.

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

Microstructural evolution in Cu-Nb processed via friction consolidation

Immiscible alloys, whether in well-mixed or layered forms, are of increasing interest based on their novel structural and functional properties, such as enhanced thermal stability against grain growth or radiation-induced defect trapping at the interfaces. To address the need for new approaches to tailor microstructures, the microstructural development of an immiscible Cu-Nb alloy processed via friction consolidation of elemental powders is investigated. Friction consolidation is a solid-phase processing technique that imparts severe plastic strain into a deforming volume resulting in elevated temperatures below the melting temperature of the alloy. Two distinct processing pathways were chosen to understand the effect of thermomechanical conditions on the final microstructure. The microstructure was characterized using scanning electron microscopy, scanning transmission electron microscopy, and X-ray diffraction techniques. Path 1 exhibited larger strain, strain rate, and temperature as compared with path 2. In path 1, agglomerated Nb particles were present in the recrystallized ultrafine-grained Cu matrix, while in path 2 extremely fine and dispersed Nb particles were present in a highly deformed Cu matrix. In both pathways, supersaturation of Cu in Nb lattices was noted, but not vice-versa. The asymmetry in mixing is explained based on deformation-based, thermodynamic and kinetic factors. These findings provide a pathway for creation of novel tailored microstructures and improved properties in any number of binary immiscible alloy systems.

Cu-Nb alloys, Microstructural refinement, Friction↗
DOE OSTI2023

Plasma production and ion heating systems for the Material Plasma Exposure eXperiment

Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.

Goulding, Richard↗
DOE OSTI2021

Effects of varying gelatin coating concentrations on RANKL induced osteoclastogenesis

Here, we assessed the effects of varying concentrations of gelatin coating on Receptor Activator of Nuclear Factor κ-B Ligand (RANKL)-induced RAW264.7 murine macrophage differentiation into osteoclast (OC) via osteoclastogenesis. The microstructures of coating surfaces with different concentrations of gelatin were examined by scanning electron microscopy and atomic force microscopy. Increased gelatin coating concentrations led to decreased gel rigidity but increased surface adhesion force attenuated OC differentiation and the decreased actin ring formation in RANKL-induced osteoclastogenesis. The decreased actin ring formation is associated with decreased lysosomal-associated membrane protein 1 (LAMP1) activity and bone resorption in the differentiated OCs with different gelatin coating concentrations as compared to the cells differentiated without gelatin coatings. In addition, increasing concentrations of gelatin coating attenuated the medium TGF-β1 protein levels and the expression levels of TGF-β and type-I (R1) and type-II (R2) TGF-β receptors in OCs, suggesting the gelatin-induced suppression of TGF-β signaling for the regulation of RNAKL-induced OC differentiation. Taken together, these findings showed that changes in gelatin coating concentrations, which were associated with altered gel thickness and substrate rigidity, might attenuate TGF-β signaling events to modulate OC differentiation and concomitant actin ring formation and bone matrix resorption in RANKL-induced osteoclastogenesis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Programmable soft valves for digital and analog control

In soft devices, complex actuation sequences and precise force control typically require hard electronic valves and microcontrollers. Existing designs for entirely soft pneumatic control systems are capable of either digital or analog operation, but not both, and are limited by speed of actuation, range of pressure, time required for fabrication, or loss of power through pull-down resistors. Using the nonlinear mechanics intrinsic to structures composed of soft materials—in this case, by leveraging membrane inversion and tube kinking—two modular soft components are developed: a piston actuator and a bistable pneumatic switch. These two components combine to create valves capable of analog pressure regulation, simplified digital logic, controlled oscillation, nonvolatile memory storage, linear actuation, and interfacing with human users in both digital and analog formats. Three demonstrations showcase the capabilities of systems constructed from these valves: 1) a wearable glove capable of analog control of a soft artificial robotic hand based on input from a human user’s fingers, 2) a human-controlled cushion matrix designed for use in medical care, and 3) an untethered robot which travels a distance dynamically programmed at the time of operation to retrieve an object. This work illustrates pathways for complementary digital and analog control of soft robots using a unified valve design.

42 ENGINEERING↗
DOE OSTI2023

Mixed imbibition controls the advance of wetting fluid in multiscale geological media

Imbibition plays a central role during multiphase flows in geological media; yet, our understanding of imbibition dynamics across the multiscale features of rocks, including fractures, microcracks, and heterogeneous matrices is limited. Here, we deployed a unique combination of micro- and clinical X-ray CT imaging to probe directly the spatiotemporal evolution of imbibition in a natural, multiscale shale sample. The method enabled unraveling of mixed imbibition behavior, with spontaneous imbibition rates in the matrix pores and microcracks driven by capillarity proceeding as $\propto t^{-1/2}$, whereas pressure-driven flows in the macroscopic fractures yield apparent imbibition rates with $\propto t^{-1}$ scaling. Regardless of the aperture, we showed that spontaneous imbibition proceeds uniformly away from the fractures. The effects of microcracks are further investigated by comparing the dynamics of the imbibition fronts with and without microcracks. Finally, these results highlight the paramount role of microcracks on anomalous front roughening and the need to account for microcracks appropriately in the development of predictive theoretical models.

58 GEOSCIENCES↗