Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “channel calibration”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Calibration strategy of the PROSPECT-II detector with external and intrinsic sources

This paper presents an energy calibration scheme for an upgraded reactor antineutrino detector for the Precision Reactor Oscillation and Spectrum Experiment (PROSPECT). The PROSPECT collaboration is preparing an upgraded detector, PROSPECT-II (P-II), to advance capabilities for the investigation of fundamental neutrino physics, fission processes and associated reactor neutrino flux, and nuclear security applications. P-II will expand the statistical power of the original PROSPECT (P-I) dataset by at least an order of magnitude. The new design builds upon previous P-I design and focuses on improving the detector robustness and long-term stability to enable multi-year operation at one or more sites. The new design optimizes the fiducial volume by elimination of dead space previously occupied by internal calibration channels, which in turn necessitates the external deployment. In this paper, we describe a calibration strategy for P-II. Here, the expected performance of externally deployed calibration sources is evaluated using P-I data and a well-benchmarked simulation package by varying detector segmentation configurations in the analysis. The proposed external calibration scheme delivers a compatible energy scale model and achieves comparable performance with the inclusion of an additional AmBe neutron source, in comparison to the previous internal arrangement. Most importantly, the estimated uncertainty contribution from the external energy scale calibration model meets the precision requirements of the P-II experiment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Measurement of CP asymmetries in ${D}_{(s)}^{+}\to \eta {\pi}^{+}$ and ${D}_{(s)}^{+}\to {\eta}^{\prime }{\pi}^{+}$ decays

Searches for CP violation in the decays ${D}_{(s)}^{+}\to \eta {\pi}^{+}$ and ${D}_{(s)}^{+}\to {\eta}^{\prime }{\pi}^{+}$ are performed using pp collision data corresponding to 6 fb –1 of integrated luminosity collected by the LHCb experiment. The calibration channels ${D}_{(s)}^{+}\to \phi {\pi}^{+}$ are used to remove production and detection asymmetries. The resulting CP-violating asymmetries are $\mathcal{A}^{CP} = (D^+ \to ηπ^+) = (0.34 ± 0.66 ± 0.16 ± 0.05)$ %, $\mathcal{A}^{CP} = (D^+_s \to ηπ^+) = (0.32 ± 0.51 ± 0.12)$ %, $\mathcal{A}^{CP} = (D^+ \to η'π^+) = (0.49 ± 0.18 ± 0.06 ± 0.05)$ %, $\mathcal{A}^{CP} = (D^+_s \to η'π^+) = (0.01 ± 0.12 ± 0.08)$ %, where the first uncertainty is statistical, the second is systematic and the third, relevant for the D + channels, is due to the uncertainty on ${\mathcal{A}}^{CP}=\left({D}^{+}\to \phi {\pi}^{+}\right)$. These measurements, currently the most precise for three of the four channels considered, are consistent with the absence of CP violation. A combination of these results with previous LHCb measurements is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Optically multiplexed neutron time-of-flight technique for inertial confinement fusion

Neutron time-of-flight (nTOF) detectors are crucial in diagnosing the performance of inertial confinement fusion (ICF) experiments, which implode targets of deuterium–tritium fuel to achieve thermonuclear conditions. These detectors utilize the fusion neutron energy spectrum to extract key measurements, including the hotspot ion temperature and fuel areal density. Previous work [Danly et al., Rev. Sci. Instrum. 94, 043502 (2023)] has demonstrated adding 1D spatial resolution to an nTOF-like detector using a neutron aperture and streak camera to measure the ion temperature profile of an ICF implosion. By contrast, the study presented herein explores modifying the 1D detector to use a fast photomultiplier tube (PMT) to validate the design of a 2D spatially resolved instrument based on reconstruction from 1D profiles. The modification would collect time-of-flight traces from separate scintillators in an imaging array with one PMT using optical fibers of varying lengths to time-multiplex the signals. This technique has been demonstrated in ride-along experiments on the OMEGA laser with 20 fiber-coupled scintillator channels connected to a Photek PMT210. Results provide constraints on the fiber lengths and PMT gating requirements to promote pulse fidelity throughout all channels. Calibration of the detector to fixed nTOFs can provide a preliminary estimate of the instrument response function (IRF), although measurement of the IRF is currently under way. These results suggest that nTOF signals can potentially be time-multiplexed with fibers so long as the design is strategic to mitigate signal-to-noise reduction, modal dispersion, and charge build-up in the PMT, which has implications beyond ion temperature imaging.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

First-ever deployment of a SiPM-on-tile calorimeter in a collider: a parasitic test with 200 GeV p p collisions at RHIC.

We describe the testing of a prototype SiPM-on-tile iron-scintillator calorimeter at the Relativistic Heavy Ion Collider (RHIC) during its 200 GeV p p run in 2024. The prototype, measuring 20 × 20 cm 2 and 24 radiation lengths in depth, was positioned in the STAR experimental hall, approximately 8 m from the interaction point and 65 cm from the beam line, covering a pseudorapidity range of about 3.1 < η < 3.4. By using the dark current of a reference SiPM as a radiation monitor, we estimate that the prototype was exposed to a fluence of about 10 10 1-MeV n eq /cm 2 . Channel-by-channel calibration was performed in a data-driven way with the signature from minimum-ionizing particles during beam-on conditions. AGeant4 detector simulation, with inputs from thePythia8 event generator, describes measurements of energy spectra and hit multiplicities reasonably well. These results mark the first deployment, commissioning, calibration, and long-term operation of a SiPM-on-tile calorimeter in a collider environment. This experimental campaign will guide detector designs and operational strategies for the ePIC detector at the future EIC, as well as other applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Calibration of the underground muon detector of the Pierre Auger Observatory

To obtain direct measurements of the muon content of extensive air showers with energy above 10 16.5 eV, the Pierre Auger Observatory is currently being equipped with an underground muon detector (UMD), consisting of 219 10 m 2 -modules, each segmented into 64 scintillators coupled to silicon photomultipliers (SiPMs). Direct access to the shower muon content allows for the study of both of the composition of primary cosmic rays and of high-energy hadronic interactions in the forward direction. As the muon density can vary between tens of muons per m 2 close to the intersection of the shower axis with the ground to much less than one per m 2 when far away, the necessary broad dynamic range is achieved by the simultaneous implementation of two acquisition modes in the read-out electronics: the binary mode, tuned to count single muons, and the ADC mode, suited to measure a high number of them. In this work, we present the end-to-end calibration of the muon detector modules: first, the SiPMs are calibrated by means of the binary channel, and then, the ADC channel is calibrated using atmospheric muons, detected in parallel to the shower data acquisition. The laboratory and field measurements performed to develop the implementation of the full calibration chain of both binary and ADC channels are presented and discussed. The calibration procedure is reliable to work with the high amount of channels in the UMD, which will be operated continuously, in changing environmental conditions, for several years.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Introducing a Markov chain-based time calibration procedure for multi-channel particle detectors: application to the SuperFGD and ToF detectors of the T2K experiment

Inter-channel mis-synchronisation can be a limiting factor to the time resolution of high performance timing detectors with multiple readout channels and independent electronics units. In these systems, time calibration methods employed must be able to efficiently correct for minimal mis-synchronisation between channels and achieve the best detector performance. We present an iterative time calibration method based on Markov Chains, suitable for detector systems with multiple readout channels. Starting from correlated hit pairs alone, and without requiring an external reference time measurement, the method solves for fixed per-channel offsets, with precision limited only by the intrinsic single-channel resolution. A mathematical proof that the method is able to find the correct time offsets to be assigned to each detector channel in order to achieve inter-channel synchronisation is given, and it is shown that the number of iterations to reach convergence within the desired precision is controllable with a single parameter. Numerical studies are used to confirm unbiased recovery of true offsets. Finally, the application of the calibration method to the Super Fine-Grained Detector (SuperFGD) and the Time of Flight (TOF) detector at the upgraded T2K near detector (ND280) shows good improvement in overall timing resolution, demonstrating the effectiveness in a real-world scenario and scalability.

calibration and fitting methods↗
DOE OSTI2026

Fast Muon Capture Monitoring in Mu2e with the CAPHRI Detector

The Mu2e experiment at Fermilab will search for the charged lepton flavor violating (CLFV) process of a neutrinoless muon-to-electron conversion in the field of an aluminum nucleus. Reaching the experiment’s target sensitivity requires precise normalization of the physics signal through accurate monitoring of the muon capture rate on the stopping target. For this purpose, the Calorimeter Precise High-Resolution Intensity detector (CAPHRI) has been developed. The detector is composed of four LYSO crystals installed in the upstream disk of the Mu2e calorimeter and read out with the standard calorimeter readout. CAPHRI measures the muon capture rate by detecting the characteristic 1.8~MeV gamma emission line of the $^{27}Al(\mu^−, \nu n \gamma) ^{26}Mg$ nuclear reaction. The fast, precise response enables injection-by-injection monitoring of proton beam intensity fluctuations. We report on the commissioning and performance characterization of CAPHRI. The response of each channel is calibrated at two SiPM overvoltages using both the intrinsic self-emission of the LYSO crystals and cosmic ray signals. In parallel, Monte Carlo simulations are used to evaluate the detector acceptance and the expected signal-to-background ratio under realistic running conditions. Preliminary results show a crystal light yield consistent with expectations and a channel inter-calibration at the 2--4% level. Simulation studies indicate that the detector acceptance and background rejection satisfy the requirements for physics operations, with about 1000 detected events per beam injection at a beam power of 1.5~kW. These results demonstrate that CAPHRI is an effective tool for beam monitoring and signal normalization in Mu2e.

Ciccarella, V. [Frascati; U. Rome La Sapienza (mai↗
DOE OSTI2026

Calorimeter calibration and performance for the Mu2e experiment

The Mu2e experiment at Fermilab will search for the charged lepton flavour-violating conversion of a muon into an electron, aiming to reach a sensitivity of $R_{\mu e} \sim 10^{-17}$, an improvement of four orders of magnitude over previous limits. To reach this goal, Mu2e will use an intense pulsed muon beam and a detector system composed of a high-precision straw tube tracker and a pure CsI crystal calorimeter. The calorimeter plays a crucial role in the experiment, as it provides particle identification capabilities that are necessary for background suppression. To perform its tasks, the detector must achieve an energy resolution better than 10% and a timing resolution below 500 ps for 100 MeV electrons. Cosmic-ray data and laser pulses are used to equalize the response of each channel, to calibrate the energy scale and to monitor the system's stability over time. This poster reports on the calibration and analysis techniques developed to ensure that the calorimeter requirements for precise energy and time measurements are met. Results for the calorimeter performance obtained during the commissioning phase will be discussed, and an overview of the current status in the Mu2e experimental hall will be presented.

Salamino, Sabrina [Frascati]↗
DOE OSTI2021

Sensitivity-Informed Bayesian Inference for Home PLC Network Models with Unknown Parameters

Bayesian inference is used to calibrate a bottom-up home PLC network model with unknown loads and wires at frequencies up to 30 MHz. A network topology with over 50 parameters is calibrated using global sensitivity analysis and transitional Markov Chain Monte Carlo (TMCMC). The sensitivity-informed Bayesian inference computes Sobol indices for each network parameter and applies TMCMC to calibrate the most sensitive parameters for a given network topology. A greedy random search with TMCMC is used to refine the discrete random variables of the network. This results in a model that can accurately compute the transfer function despite noisy training data and a high dimensional parameter space. The model is able to infer some parameters of the network used to produce the training data, and accurately computes the transfer function under extrapolative scenarios.

42 ENGINEERING↗
DOE OSTI2022

Calibration and Localization of Optically Pumped Magnetometers Using Electromagnetic Coils

In this paper, we propose a method to estimate the position, orientation, and gain of a magnetic field sensor using a set of (large) electromagnetic coils. We apply the method for calibrating an array of optically pumped magnetometers (OPMs) for magnetoencephalography (MEG). We first measure the magnetic fields of the coils at multiple known positions using a well-calibrated triaxial magnetometer, and model these discreetly sampled fields using vector spherical harmonics (VSH) functions. We then localize and calibrate an OPM by minimizing the sum of squared errors between the model signals and the OPM responses to the coil fields. We show that by using homogeneous and first-order gradient fields, the OPM sensor parameters (gain, position, and orientation) can be obtained from a set of linear equations with pseudo-inverses of two matrices. The currents that should be applied to the coils for approximating these low-order field components can be determined based on the VSH models. Computationally simple initial estimates of the OPM sensor parameters follow. As a first test of the method, we placed a fluxgate magnetometer at multiple positions and estimated the RMS position, orientation, and gain errors of the method to be 1.0 mm, 0.2°, and 0.8%, respectively. Lastly, we calibrated a 48-channel OPM array. The accuracy of the OPM calibration was tested by using the OPM array to localize magnetic dipoles in a phantom, which resulted in an average dipole position error of 3.3 mm. The results demonstrate the feasibility of using electromagnetic coils to calibrate and localize OPMs for MEG.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Development of Land‐River Two‐Way Hydrologic Coupling for Floodplain Inundation in the Energy Exascale Earth System Model

Abstract Floodplain inundation links river and land systems through significant water, sediment, and nutrient exchanges. However, these two‐way interactions between land and river are currently missing in most Earth System Models. In this study, we introduced the two‐way hydrological coupling between the land component, E3SM Land Model, and the river component, Model for Scale Adaptive River Transport, in Energy Exascale Earth System Model (E3SM) to study the impacts of floodplain inundation on land and river processes. We calibrated the river channel geometry and developed a new data‐driven inundation scheme to improve the simulation of inundation dynamics in E3SM. The new inundation scheme captures 96% of the spatial variation of inundation area in a satellite inundation product at global scale, in contrast with 7% when the default inundation scheme of E3SM was used. Global simulations including the new inundation scheme performed at resolution with and without land‐river two‐way coupling was used to quantify the impact of coupling. Comparisons show that two‐way coupling modifies the water and energy cycle in 20% of the global land cells. Specifically, riverine inundation is reduced by two‐way coupling, but inland inundation is intensified. Wetter periods are more impacted by the two‐way coupling at the global scale, while regions with different climates exhibit different sensitivities. The two‐way exchange of water between the land and river components of E3SM provides the foundation for enabling two‐way coupling of land‐river sediment and biogeochemical fluxes. These capabilities will be used to improve understanding of the interactions between water and biogeochemical cycles and their response to human perturbations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

In-flight polarization angle calibration for LiteBIRD: blind challenge and cosmological implications

We present a demonstration of the in-flight polarization angle calibration for the JAXA/ISAS second strategic large class mission, LiteBIRD, and estimate its impact on the measurement of the tensor-to-scalar ratio parameter, r, using simulated data. We generate a set of simulated sky maps with CMB and polarized foreground emission, and inject instrumental noise and polarization angle offsets to the 22 (partially overlapping) LiteBIRD frequency channels. Our in-flight angle calibration relies on nulling the EB cross correlation of the polarized signal in each channel. This calibration step has been carried out by two independent groups with a blind analysis, allowing an accuracy of the order of a few arc-minutes to be reached on the estimate of the angle offsets. Both the corrected and uncorrected multi-frequency maps are propagated through the foreground cleaning step, with the goal of computing clean CMB maps. We employ two component separation algorithms, the Bayesian-Separation of Components and Residuals Estimate Tool (B-SeCRET), and the Needlet Internal Linear Combination (NILC). We find that the recovered CMB maps obtained with algorithms that do not make any assumptions about the foreground properties, such as NILC, are only mildly affected by the angle miscalibration. However, polarization angle offsets strongly bias results obtained with the parametric fitting method. Once the miscalibration angles are corrected by EB nulling prior to the component separation, both component separation algorithms result in an unbiased estimation of the r parameter. While this work is motivated by the conceptual design study for LiteBIRD, its framework can be broadly applied to any CMB polarization experiment. In particular, the combination of simulation plus blind analysis provides a robust forecast by taking into account not only detector sensitivity but also systematic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Requirements on the gain calibration for LiteBIRD polarisation data with blind component separation

The detection of primordial B modes of the cosmic microwave background (CMB) could provide information about the early stages of the Universe's evolution. The faintness of this signal requires exquisite calibration accuracy and control of instrumental systematic effects which otherwise could bias the measurements. In this work, we study the impact of an imperfect relative polarisation gain calibration on the recovered value of the tensor-to-scalar ratio r for the LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We derive requirements on the relative calibration accuracy of the overall polarisation gain (Δg ν ) for each LiteBIRD frequency channel. Our results show that minimum variance techniques, as NILC, are less sensitive to systematic gain calibration uncertainties compared to a parametric approach, if the latter is not equipped with a proper modelling of these instrumental effects. In this study, the most stringent requirements are found in the channels where the CMB signal is relatively brighter, with the tightest constraints at 166 GHz (Δg ν ≈ 0.16%). This differs from the outcome of an analogous analysis performed with a parametric method, where the tightest requirements are obtained for the foreground-dominated channels. Gain calibration uncertainties, corresponding to the derived requirements, are then simultaneously propagated into all frequency channels. By doing so, we find that the overall impact on estimated r is lower than the total gain systematic budget for LiteBIRD approximately by a factor 5, due to the correlations of the impacts of gain calibration uncertainties in different frequency channels. In order to decouple the systematic effect from the specific choice of the model, we derive the requirements assuming constant spectral parameters for the foreground emission. To assess the robustness of the obtained results against more realistic scenarios, we repeat the analysis assuming sky models of intermediate and high complexity. In these further cases, we adopt an optimised NILC pipeline, called the Multi-Clustering NILC (MC-NILC). We find that the impact of gain calibration uncertainties on r is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor 1.8.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Two-dimensional plasma density evolution local to the inversion layer during sawtooth crash events using Beam Emission Spectroscopy

Here, we present methods for analyzing Beam Emission Spectroscopy (BES) data to obtain the plasma density evolution associated with rapid sawtooth crash events at the DIII-D tokamak. BES allows coverage over a 2D spatial plane, inherently local measurements, with fast time responses, and, therefore, provides a valuable new channel for data during sawtooth events. A method is developed to remove sawtooth-induced edge-light pulses contained in the BES data. The edge light pulses appear to be from the D α emission produced by edge recycling during sawtooth events, and are large enough that traditional spectroscopic filtering and data analysis techniques are insufficient to deduce physically meaningful quantities. A cross-calibration of 64 BES channels is performed by using a novel method to ensure accurate measurements. For the large-amplitude density oscillations observed, we discuss and use the non-linear relationship between the BES signal δI/ I 0 and the plasma density variation δn e / n e0 . The 2D BES images cover an 8 × 20 cm 2 region around the sawtooth inversion layer and show large-amplitude density oscillations, with additional significant spatial variations across the inversion layer that grows and peaks near the time of the temperature crash. The edge light removal technique and method of converting large-amplitude δI/ I 0 to δn e / n e0 presented here may help analyze other impulsive MHD phenomena in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spectral Fidelity of Earth's Terrestrial and Aquatic Ecosystems

The Surface Biology and Geology (SBG) investigation will create global maps of spectral surface reflectance and emissivity at a cadence of 16 days or better, with coverage to address global questions about Earth's geology, cryosphere and ecosystems. The revolutionary potential poses a commensurate challenge: creating contiguous maps free from regional biases induced by atmosphere, observation geometry, or inversion error. This will require an accurate calibration with precise knowledge of each channel's spectral response. In this work, we quantify the impact of spectral calibration on SBG's aquatic and terrestrial ecosystem objectives. We find that contemporary algorithms for ecosystem trait retrieval demand more accurate spectral calibration than historical missions. Errors due to drift or spatial nonuniformity in the wavelength calibration that have previously been considered acceptable can cause systematic errors larger than the instrument noise and of the same order as the variability SBG aims to measure. Moreover, their impact on atmospheric correction can induce climate-dependent systematic errors that thwart comparisons between ecosystems. These results underscore the importance of spectral response accuracy in SBG mission design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Calibration improvements expand filterscope diagnostic use

The filterscope diagnostic on DIII-D utilizes photomultiplier tubes to measure visible light emission from the plasma. The system has undergone a substantial upgrade since previous attempts to cross-calibrate the filterscope with other spectroscopic diagnostics were unsuccessful. The optics now utilize a dichroic mirror to initially split the light at nearly 99% transmission or reflectance for light below or above 550 nm. This allows the system to measure D α emission without degrading visible light emission from the plasma for wavelengths below 550 nm (to measure D β , D γ , W–I, C-III, etc.). Additional optimization of the optical components and calibration techniques reduce the error in the signal up to 10% in some channels compared to previous methods. Cross-calibration measurements with two other high resolution spectroscopic diagnostics now show excellent agreement for the first time. This expands the capabilities of the filterscope system allowing measurement of divertor detachment, emission profiles, edge-localized mode behavior, and plasma–wall interactions. Additionally, it enables direct comparisons against calculations from boundary plasma simulations. These were not possible before.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

A 2–4 keV multilayer mirrored channel for the NIF Dante system

During inertial confinement fusion experiments at the National Ignition Facility (NIF), a capsule filled with deuterium and tritium (DT) gas, surrounded by a DT ice layer and a high-density carbon ablator, is driven to the temperature and densities required to initiate fusion. In the indirect method, 2 MJ of NIF laser light heats the inside of a gold hohlraum to a radiation temperature of 300 eV; thermal x rays from the hohlraum interior couple to the capsule and create a central hotspot at tens of millions degrees Kelvin and a density of 100–200 g/cm3. During the laser interaction with the gold wall, m-band x rays are produced at ∼2.5 keV; these can penetrate into the capsule and preheat the ablator and DT fuel. Preheat can impact instability growth rates in the ablation front and at the fuel–ablator interface. Monitoring the hohlraum x-ray spectrum throughout the implosion is, therefore, critical; for this purpose, a Multilayer Mirror (MLM) with flat response in the 2–4 keV range has been installed in the NIF 37° Dante calorimeter. Precision engineering and x-ray calibration of components mean the channel will report 2–4 keV spectral power with an uncertainty of ±8.7%.

Rubery, M. S. (ORCID:0000000344152818)↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗