SEARCH · Search NASA

Results for “Separated flows”

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

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Hydrodynamics of a multicomponent vesicle under strong confinement

Herein we numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Effect of hollow fiber membrane packing on the performance of modules formed with fiber tows

Membrane processes are the preferred option for many gas separations. Hollow fiber membranes are commonly used in these processes. Commercial hollow fiber membrane modules are often fabricated from small groups of fibers (i.e., a tow) that are arranged into a larger cylindrical bundle. The existence of fiber tows can lead to packing non-uniformity. This nonuniformity can result in nonideal fluid distribution in the shell and be detrimental to module performance. In this work, full three-dimensional models of flow and mass transfer within fiber bundles formed from tows have been developed. These models account for the detailed fiber arrangement within modules and are solved to evaluate module separation performance in terms of stage cut and required membrane area as a function product gas purity for carbon capture applications. Results are presented illustrating the effects of varying fiber packing arrangements within a tow and varying tow packing arrangements. Further, the results suggest smaller tows are preferred over larger tows, at fixed overall fiber packing and lower intra- and inter-tow packings are preferred. To reduce the computational cost of simulating large fiber bundles, the use of an equivalent planar bundle is proposed (EPB). A simple procedure for determining the geometry of the EPB is described and validated through comparisons of calculated performance metrics with results obtained from full three-dimensional simulations.

42 ENGINEERING↗

DOE OSTI2023

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

A Novel Approach for Bounding the Stress Experienced by the Core of Utility-Scale Printed Circuit Heat Exchangers Under Thermohydraulic Loads

Printed circuit heat exchangers (PCHEs) have rapidly gained popularity since being introduced nearly three decades ago, and they are currently widely deployed in the petrochemical and aviation industry. Their compactness, thermohydraulic efficiency, inherent suitability for high temperature/pressure fluids containment, and demonstrated durability are some of the reasons the nuclear industry is seeking to adopt this technology as well. However, the relatively strict nuclear-related regulatory design code, especially when classified as critical to the safety of the reactors, is posing challenges to adopting the technology. From stress analysis point of view, one undesirable feature of PCHEs is their geometrical complexity, which is implied by their multilength-scale features. As a result, a full-scale model of a utility-scale exchanger cannot simply be solved on a computer because meshing such components results in a vast number of degrees-of-freedom. Here, this work seeks to address the challenge of stress analyses to PCHEs by presenting a method to simplify the geometry of PCHE designs. The models proposed by this work can be practically analyzed on a standard computer and provide a path for implementing ASME design rules. The analyses presented herein are divided into five separate investigations. Each is carried out to incrementally simplify the analyzed model by addressing features such as the shapes of the flow passages, the complex distribution of stress in large components, the three-dimensionality of the stress and strain, the thermal stresses caused by thermohydraulic operation observed experimentally and more.

42 ENGINEERING↗

DOE OSTI2022

Machine learning and atomic layer deposition: Predicting saturation times from reactor growth profiles using artificial neural networks

In this work, we explore the application of deep neural networks to the optimization of atomic layer deposition (ALD) processes. In particular, we focus on a one-shot optimization problem, where we try to predict the optimal dose time that leads to saturation everywhere in the reactor based on thickness values measured at different points of an ALD reactor after a single trial growth. In order to tackle this problem, we introduce a dataset designed to train neural networks to predict saturation times based on these inputs for a cross-flow ALD reactor. Here, we then explore the predictive ability of artificial neural networks of different depths and sizes using a separate testing dataset to evaluate their accuracies. The results obtained show that networks trained using stochastic gradient descent methods can accurately predict saturation times without requiring any additional information on the surface kinetics. This provides a viable approach to minimize the number of experiments required to optimize new ALD processes in a known reactor, and it highlights the way machine learning can be leveraged for thin film growth and manufacturing. While the datasets and training procedure depend on the reactor geometry, the trained neural networks provide a general surrogate model connecting thickness values and trial dose times with optimal saturation times that can be reused for different ALD processes within the same reactor.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Using Bayesian Methodology to Estimate Liquefied Natural Gas Leak Frequencies

This analysis provides estimates on the leak frequencies of nine components found in liquefied natural gas (LNG) facilities. Data was taken from a variety of sources, with 25 different data sets included in the analysis. A hierarchical Bayesian model was used that assumes that the log leak frequency follows a normal distribution and the logarithm of the mean of this normal distribution is a linear function of the logarithm of the fractional leak area. This type of model uses uninformed prior distributions that are updated with applicable data. Separate models are fit for each component listed. Five order-of-magnitude fractional leak areas are considered, based on the flow area of the component. Three types of supporting analyses were performed: sensitivity of the model to the data set used, sensitivity of the leak frequency estimates to differences in the model structure or prior distributions, and sufficiency of sample sized used for convergence. Recommended leak frequency distributions for all component types and leak sizes are given. These leak frequency predictions can be used for quantitative risk assessments in the future.

03 NATURAL GAS↗

DOE OSTI2026

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Hydrogen Enrichment for Solid Oxide Fuel Cell Gas Turbine (SOFC/GT) Hybrid Power Generation Systems Using Mechanical Gas Separation

This presentation, given at the 37 th annual Pittsburgh Coal Conference on September 8-11, 2020, includes information on the process, previous work, normal and oblique shocks model, diffuser design model, and calculation results (expansion, recompression with no losses, and recompression with mass depletion). Summary includes the concept of aerodynamic separation of water and carbon dioxide from a SYNGAS stream and that deposition was modeled by mass depletion from a compressible flow.

VanOsdol, John↗

DOE OSTI2022

High energy density redox flow device

Redox flow devices are described including a positive electrode current collector, a negative electrode current collector, and an ion-permeable membrane separating said positive and negative current collectors, positioned and arranged to define a positive electroactive zone and a negative electroactive zone; wherein at least one of said positive and negative electroactive zone comprises a flowable semi-solid composition comprising ion storage compound particles capable of taking up or releasing said ions during operation of the cell, and wherein the ion storage compound particles have a polydisperse size distribution in which the finest particles present in at least 5 vol % of the total volume, is at least a factor of 5 smaller than the largest particles present in at least 5 vol % of the total volume.

Chiang, Yet-Ming↗

DOE OSTI2024

High energy density redox flow device

Chiang, Yet-Ming↗

DOE OSTI2023

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Computationally efficient optimization models for preliminary distillation column design and separation energy targeting

We propose versatile shortcut distillation column and separation energy targeting models that are well-suited for superstructure-based process synthesis. The models are based on a novel reformulation of the Underwood equations to address systems where the components that are present in the feed can vary due to zero flow rates of some components. Also, we propose valid constraints, resulting in a significant enhancement of the computational performance of the models. Furthermore, the proposed distillation column model can automatically identify adequate key components and the energy requirement of a desired separation task, while considering a wide range of types of separations including non-sharp/sloppy splits. Also, the proposed separation energy targeting model can be used to estimate an energy requirement target for the separation of a mixture without finding detail network configurations. Due to their versatility and computational efficiency, the proposed models can be readily used for and expand the scope of superstructure-based process synthesis approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Solid–liquid separation of lignocellulosic sugars from biomass by rotating ceramic disc filtration

In many biomass conversion processes, the separation of cellulosic sugars from residual, lignin-rich solids is a critical step, and achieving high recovery yields of sugars by conventional tangential crossflow and vacuum filtration is challenged by the presence of biomass solids, which rapidly foul filters, resulting in decreased throughput. Considering the performance limitations of traditional filtration methods, dynamic filtration, which generates high shear at the membrane surface to decrease fouling, is emerging as a viable alternative for demanding solid–liquid separations. For high solids separations, there is little available information regarding the performance, limitations, and energy consumption of dynamic filtration. To that end, here we characterized the performance of a dynamic filtration module, specifically a rotating ceramic disc (RCD) filter, for the aseptic recovery of cellulosic sugars from biomass solids following pretreatment and enzymatic hydrolysis. We show how RCD rotational velocity and percent biomass solids impact the filter throughput. Additionally, we used computational fluid dynamics (CFD) simulations to estimate the shear rate at the membrane surface and to visualize hydrodynamic profiles within the module. With the combined CFD simulations and experimental results, we estimated the energy demand and operating expenses for a viable dynamic filtration system operating with a lignocellulosic feed slurry. Our results indicate that an RCD filter can achieve ≥95% recovery of sugars and produce a retentate slurry containing 12 wt% insoluble solids with low energy consumption (a 2.2-fold improvement over cross-flow filtration) and low operating costs ($\$$0.06 per kg sugars). These results show a viable path towards operationally reliable, energy efficient, and cost-effective separations of sterilized cellulosic sugars from biomass solids and highlight the potential of dynamic filtration systems for challenging solid–liquid separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2021

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗

DOE OSTI2022

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗

DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗