SEARCH · Search NASA

Results for “Volume of fluid”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Quantitative Assessment of Parent Well Effect on Hydraulic Fracture Propagation at HFTS2: Insights from Cross-Well Strain Measurements and Microseismic Data

Understanding fracture propagation behavior is essential for optimizing hydraulic fracturing in unconventional reservoirs. This study demonstrates the value of integrating Low-Frequency Distributed Acoustic Sensing (LF-DAS) and microseismic data, which together provide a more complete picture of fracture growth. Using data from Hydraulic Fracturing Test Site 2 (HFTS2), we identify stress changes in depletion zones induced by parent wells as a key factor influencing fracture propagation. This result is shown by new measurements of in-situ fracture propagation velocity and fracture-hit volume (fluid volume at fracture hit?) from LF-DAS and event density from microseismic. These findings highlight the importance of considering parent well effects, well spacing, and stimulation sequencing in completion design to improve reservoir development and production efficiency.

depletion zones↗

DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗

DOE OSTI2022

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Constraining maximum event magnitude during injection-triggered seismicity

Abstract Understanding mechanisms controlling fluid injection-triggered seismicity is key in defining strategies to ameliorate it. Recent triggered events (e.g. Pohang, Mw 5.5) have exceeded predictions of average energy release by a factor of >1000x, necessitating robust methodologies to both define critical antecedent conditions and to thereby constrain anticipated event size. We define maximum event magnitudes resulting from triggering as a function of pre-existing critical stresses and fluid injection volume. Fluid injection experiments on prestressed laboratory faults confirm these estimates of triggered moment magnitudes for varied boundary conditions and injection rates. In addition, observed ratios of shear slip to dilation rates on individual faults signal triggering and may serve as a measurable proxy for impending rupture. This new framework provides a robust method of constraining maximum event size for preloaded faults and unifies prior laboratory and field observations that span sixteen decades in injection volume and four decades in length scale.

58 GEOSCIENCES↗

DOE OSTI2020

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗

DOE OSTI2022

Vorticity-induced flow-focusing leads to bubble entrainment in an inkjet printhead: Synchrotron x-ray and volume-of-fluid visualizations

The oscillatory flows present in an inkjet printhead can lead to strong deformations of the air-liquid interface at the nozzle exit. Such deformations may lead to an inward directed air jet with bubble pinch-off and the subsequent entrainment of an air bubble, which is highly detrimental to the stability of inkjet printing. Understanding the mechanisms of bubble entrainment is therefore crucial to improving print stability. Here, in the present work, we use ultrafast x-ray phase-contrast imaging and direct numerical simulations based on the volume-of-fluid method to study the mechanisms underlying the bubble entrainment in a piezoacoustic printhead. We first demonstrate good agreement between experiments and numerics. We then show the different classes of bubble pinch-off obtained in experiments, and that those were also captured numerically. The numerical results are then used to show that the baroclinic torque, which is generated at the gas-liquid interface due to the misalignment of density and pressure gradients, results in a flow-focusing effect that drives the formation of the air jet from which a bubble can pinch off.

42 ENGINEERING↗

DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗

DOE OSTI2023

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗

DOE OSTI2021

Analysis of Completion Design Impact on Cluster Efficiency and Pressure-Based Well Communication in HFTS-2 Delaware Basin

The Hydraulic Fracturing Test Site 2 (HFTS-2) is a joint industry project in the Delaware basin to advance hydraulic fracturing understanding and improve productivity in shale reservoirs. The project integrates multi-disciplinary approaches to evaluate different completion designs, well spacing, inter-well communication, and stimulated rock volume, among other factors. This paper focused on two major areas related to hydraulic fracture performance. First, an analysis of different completion designs on cluster efficiency based on near well Distributed Acoustic Sensing (DAS). Second, an evaluation of well and completion designs on inter-well fracture driven interactions (FDIs) based on downhole pressure monitoring. Fluid/sand distribution and cluster efficiency analyses were based on near wellbore DAS data collected from two adjacent horizontal wells completed in two different landing zones in the Wolfcamp formation. These wells had different completion designs aiming to evaluate the effect of normal vs. extended stage lengths, perforation hole tapering and limited entry. Standard deviation from ideal fluid/sand distribution and waterfall plots were used to evaluate cluster efficiency for each design and stage. Inter-well FDIs analysis was conducted among the horizontal wells and a vertical monitor well. One horizontal well served as the monitor well while the other horizontal well was being treated. The vertical well was instrumented with downhole pressure and temperature gauges to aid monitoring fracture height growth. The pressure response during and after fracturing was characterized based on maximum pressure increase value and slope. Pressure response vs. FDIs trigger factors such distance, cluster efficiency and stage fluid volume were also analyzed. Based on the different completion and perforation designs tested, DAS analysis suggests that limited entry design worked best. Extreme limited entry showed the potential of high perforation erosion and reduced cluster efficiency. Here, the limited entry and tapered perforation design demonstrated potential to improve the cluster efficiency for extended stage lengths. Pressure monitoring across formation units proved to be critical to understand fracture interactions and fracture vertical growth. Pressure communication across different formation units during hydraulic fracturing operation indicate fractures grew upwards during Wolfcamp wells fracturing. However, this pressure communication dissipated over time. High intensity FDIs were recorded when the frac stages Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R001/2477607/urtec-2021-5289-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5289 were closer to the pressure gauge location in the monitor wells. Some of these stages that produced high intensity FDIs also had high fluid volume per cluster and low cluster efficiency. The multi-disciplinary and high-quality data collected from HFTS-2 helped to further understand why completion approaches such as limited entry and tapered perforation design are successful in improving cluster efficiency. The DAS data combined with downhole high-resolution pressure measurements also helped to quantify the effect of lower cluster efficiency data on the incidence and intensity of FDIs.

58 GEOSCIENCES↗

DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Finding a contaminant source in a volume of flowing fluid

A method includes: mapping a boundary of a volume of flowing fluid; partitioning the volume by a computational mesh; finding a contaminant location at a first sensor that is disposed within the volume; obtaining a measured velocity of the flowing fluid within the volume; generating a reversed velocity vector field within the mesh, in response to the measured velocity; time stepping the contaminant location from the first sensor along the reversed velocity vector field until the contaminant location intersects the boundary of the volume; and finding a contaminant source at the intersection of the time stepped contaminant location with the boundary of the volume.

Lopez-Marrero, Vanessa↗

DOE OSTI2022

Comprehensive single-shot biophysical cytometry using simultaneous quantitative phase imaging and Brillouin spectroscopy

Abstract Single-cell analysis, or cytometry, is a ubiquitous tool in the biomedical sciences. Whereas most cytometers use fluorescent probes to ascertain the presence or absence of targeted molecules, biophysical parameters such as the cell density, refractive index, and viscosity are difficult to obtain. In this work, we combine two complementary techniques—quantitative phase imaging and Brillouin spectroscopy—into a label-free image cytometry platform capable of measuring more than a dozen biophysical properties of individual cells simultaneously. Using a geometric simplification linked to freshly plated cells, we can acquire the cellular diameter, volume, refractive index, mass density, non-aqueous mass, fluid volume, dry volume, the fractional water content of cells, both by mass and by volume, the Brillouin shift, Brillouin linewidth, longitudinal modulus, longitudinal viscosity, the loss modulus, and the loss tangent, all from a single acquisition, and with no assumptions of underlying parameters. Our methods are validated across three cell populations, including a control population of CHO-K1 cells, cells exposed to tubulin-disrupting nocodazole, and cells under hypoosmotic shock. Our system will unlock new avenues of research in biophysics, cell biology, and medicine.

Steelman, Zachary A.↗

DOE OSTI2024

Bubbler Design Updates for Nuclear Safeguards Applications

The accurate monitoring of molten salt is crucial for nuclear safeguards, particularly in the context of mixtures used in molten salt reactors and pyropocessing. To attain material accountancy in actinide bearing molten salts, fluid volume and density are needed. Previous work utilized a triple bubbler sensor to measure the fluid level and density in processing tanks. This triple bubbler system utilized three gas dip-tubes strategically positioned at different heights within the molten salt such that the molten salt level, density, and surface tension were determined simultaneously. This investigation explores the effectiveness of an updated bubbler design with larger internal diameters for the dip-tubes to reduce plugging. The findings of this study will provide valuable insights into a bubbler design and configuration for applications in nuclear safeguards. This will contribute to enhanced accountability and transparency in the monitoring of actinide bearing molten salts, promoting objectives in nuclear safeguards.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Bubbler Design Updates for Nuclear Safeguards Application

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗

DOE OSTI2023

Integration of Biorelevant Pediatric Dissolution Methodology into PBPK Modeling to Predict In Vivo Performance and Bioequivalence of Generic Drugs in Pediatric Populations: a Carbamazepine Case Study

This study investigated the impact of gastro-intestinal fluid volume and bile salt (BS) concentration on the dissolution of carbamazepine (CBZ) immediate release (IR) 100 mg tablets and to integrate these in vitro biorelevant dissolution profiles into physiologically based pharmacokinetic modelling (PBPK) in pediatric and adult populations to determine the biopredictive dissolution profile. Dissolution profiles of CBZ IR tablets (100 mg) were generated in 50–900 mL biorelevant adult fasted state simulated gastric and intestinal fluid (Ad-FaSSGF and Ad-FaSSIF), also in three alternative compositions of biorelevant pediatric FaSSGF and FaSSIF medias at 200 mL. This study found that CBZ dissolution was poorly sensitive to changes in the composition of the biorelevant media, where dissimilar dissolution (F2 = 46.2) was only observed when the BS concentration was changed from 3000 to 89 μM (Ad-FaSSIF vs Ped-FaSSIF 50% 14 BS). PBPK modeling demonstrated the most predictive dissolution volume and media composition to forecast the PK was 500 mL of Ad-FaSSGF/Ad-FaSSIF media for adults and 200 mL Ped-FaSSGF/FaSSIF media for pediatrics. A virtual bioequivalence simulation was conducted by using Ad-FaSSGF and/or Ad-FaSSIF 500 mL or Ped-FaSSGF and/or Ped-FaSSIF 200 mL dissolution data for CBZ 100 mg (reference and generic test) IR product. The CBZ PBPK models showed bioequivalence of the product. This study demonstrates that the integration of biorelevant dissolution data can predict the PK profile of a poorly soluble drug in both populations. Further work using more pediatric drug products is needed to verify biorelevant dissolution data to predict the in vivo performance in pediatrics.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

High Performance Adaptive Physics Refinement to Enable Large-Scale Tracking of Cancer Cell Trajectory

The ability to track simulated cancer cells through the circulatory system, important for developing a mechanistic understanding of metastatic spread, pushes the limits of today's supercomputers by requiring the simulation of large fluid volumes at cellular-scale resolution. To overcome this challenge, we introduce a new adaptive physics refinement (APR) method that captures cellular-scale interaction across large domains and leverages a hybrid CPU-GPU approach to maximize performance. Through algorithmic advances that integrate multi-physics and multi-resolution models, we establish a finely resolved window with explicitly modeled cells coupled to a coarsely resolved bulk fluid domain. In this work we present multiple validations of the APR framework by comparing against fully resolved fluid-structure interaction methods and employ techniques, such as latency hiding and maximizing memory bandwidth, to effectively utilize heterogeneous node architectures. Collectively, these computational developments and performance optimizations provide a robust and scalable framework to enable system-level simulations of cancer cell transport.

59 BASIC BIOLOGICAL SCIENCES↗