Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Science - Operating Systems”

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 919 records · Page 51

Optimization of an on-board imaging system for extremely rapid radiation therapy

Purpose: Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system. Methods: Three different CT detectors are proposed: low-resolution (all 4 × 4 mm pixels), medium-resolution (a combination of 4 × 4 mm pixels and 2 × 2 mm pixels), and high-resolution (all 1 × 1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called “realistic” images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a “gold standard” deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields—the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed. Results: By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm. Conclusions: The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

60 APPLIED LIFE SCIENCES↗

High-Throughput Data Processing at FRIB Using ESnet

Real-time or nearly real-time (nearline) data processing methods are critical tools as detector technologies and data acquisition (DAQ) systems allow for higher data rates and volumes. The introduction of the energy sciences network (ESnet), a U.S. Department of Energy (DOE) supported high-speed network for scientific research, creates opportunities to leverage the computing power of DOE facilities like the National Energy Research Scientific Computing Center (NERSC). As a first step toward realizing a DOE Office of Science Integrated Research Infrastructure (IRI) pattern, an automated workflow was developed to remotely process data obtained from a nuclear physics experiment at the Facility for Rare Isotope Beams (FRIB) at NERSC with data transferred between FRIB and NERSC over ESnet. The workflow demonstrated the ability to process one week’s worth of experimental data in approximately 90 min and was used successfully for nearline analysis during a recently completed FRIB experiment. Here, a summary of the workflow development and results of recent demonstrations will be presented.

Data processing↗

Experimental and Computational Investigation of Integrated Internal and Film Cooling Designs Incorporating a Thermal Barrier Coating

Abstract Few studies in the open literature have studied the effect of thermal barrier coatings (TBC) when used in combination with shaped hole film cooling and enhanced internal cooling techniques. The current study presents Reynolds-averaged Navier–Stokes (RANS) conjugate heat transfer simulations that identify trends in cooling design performance as well as experimental measurements of overall effectiveness using a flat-plate matched-Biot number model with a simulated TBC layer of 0.42D thickness, where D is the film cooling hole diameter. Coolant is fed to the film cooling holes in a co-flow configuration, and the results of both the smooth and rib-turbulated channels are compared. At a constant coolant flow rate, enhanced internal cooling was found to provide a 44% increase in spatially-averaged overall effectiveness, $\bar{\bar{\phi}}$, without a TBC. The results show that the addition of a TBC can raise $\bar{\bar{\phi}}$ on a film-cooled component surface by 47%. The optimum velocity ratio was found to decrease with the addition of enhanced cooling techniques and a TBC as the film provided minimal benefit at the expense of reduced internal cooling. While the computational results closely identified trends in overall system performance without a TBC, the model over-predicted effectiveness on the metal–TBC interface. The results of this study will inform turbine component design as material science advances increase the reliability of the TBC.

Engineering↗

Side readout of long scintillation crystal elements with digital SiPM for TOF-DOI PET

Purpose: Side readout of scintillation light from crystal elements in positron emission tomography (PET) is an alternative to conventional end-readout configurations, with the benefit of being able to provide accurate depth-of-interaction (DOI) information and good energy resolution while achieving excellent timing resolution required for time-of-flight PET. This paper explores different readout geometries of scintillation crystal elements with the goal of achieving a detector that simultaneously achieves excellent timing resolution, energy resolution, spatial resolution, and photon sensitivity. Methods: The performance of discrete LYSO scintillation elements of different lengths read out from the end/side with digital silicon photomultipliers (dSiPMs) has been assessed. Results: Compared to 3 × 3 × 20 mm{sup 3} LYSO crystals read out from their ends with a coincidence resolving time (CRT) of 162 ± 6 ps FWHM and saturated energy spectra, a side-readout configuration achieved an excellent CRT of 144 ± 2 ps FWHM after correcting for timing skews within the dSiPM and an energy resolution of 11.8% ± 0.2% without requiring energy saturation correction. Using a maximum likelihood estimation method on individual dSiPM pixel response that corresponds to different 511 keV photon interaction positions, the DOI resolution of this 3 × 3 × 20 mm{sup 3} crystal side-readout configuration was computed to be 0.8 mm FWHM with negligible artifacts at the crystal ends. On the other hand, with smaller 3 × 3 × 5 mm{sup 3} LYSO crystals that can also be tiled/stacked to provide DOI information, a timing resolution of 134 ± 6 ps was attained but produced highly saturated energy spectra. Conclusions: The energy, timing, and DOI resolution information extracted from the side of long scintillation crystal elements coupled to dSiPM have been acquired for the first time. The authors conclude in this proof of concept study that such detector configuration has the potential to enable outstanding detector performance in terms of timing, energy, and DOI resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

libEnsemble: A complete Python toolkit for dynamic ensembles of calculations

Almost all science and engineering applications eventually stop scaling: their runtime no longer decreases as available computational resources increase. Therefore, many applications will struggle to efficiently use emerging extreme-scale high-performance, parallel, and distributed systems. libEnsemble is a complete Python toolkit and workflow system for intelligently driving ensembles of experiments or simulations at massive scales. It enables and encourages multidisciplinary design, decision, and inference studies portably running on laptops, clusters, and supercomputers.

97 MATHEMATICS AND COMPUTING↗

Future Directions in Topological States of Matter: Beyond the Single Particle Picture

Correlations between electrons have brought about some of the most celebrated discoveries in quantum materials research, including all forms of superconductivity, the fractional quantum hall effect, and giant magnetoresistance. A more recent development has been the discovery that the topological aspects of electrons are important, producing new states of matter like topological insulators, and Dirac and Weyl fermions. This conference will bring together experts from both the "correlated" and "topological" communities to explore future directions that merge both fields, including emergent physics in the flat bands of twisted bilayer systems, and new routes to quantum computation using topological superconductors and spin liquids. We hope to come away with new research questions, and new material synthesis directions, that will strongly overlap with the Basic Energy Sciences mission to support the fundamental experimental and theoretical quantum materials research.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Angular Correlation Date Measurements with the GeRMAC system

Advanced modeling and simulation efforts have improved at Idaho National Laboratory in recent years with a solid foundation of experimental results. Current computational methods represent significant modeling capabilities but are limited by the accuracy and availability of nuclear data. The creation of pre- and post-processing software tools to address these limitations is fundamental to the improvement of nuclear science modeling capacities. One aspect of predictive modeling tools deals with gamma-rays emitted from radionuclides, including fissile or fissionable material, fission products, or activation products, produced in reactor experiments or other neutron environments. The resulting radionuclides decay in unique ways, providing complications upon measurement as a result of random and cascade, or true, coincidence summing. These effects are not easily quantified during modeling efforts of gamma-ray source terms., The germanium rotational measurements for angular correlation (GeRMAC) system was built to quantify the relative angles for gamma rays emitted by radionuclides of interest to investigate true coincidence, or cascade, summing as well as the nuclear energy levels of decay schemes of interest. Proof of concept studies utilize a series of laboratory check sources to provide validity, and it will soon be used to perform the same measurements for fission products of interest. The resulting data can be used to implement into a Monte Carlo code, such as Geant4, to provide more precise gamma-ray source terms following irradiations of materials.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

In silico evaluation and exploration of antibiotic tuberculosis treatment regimens

Improvement in tuberculosis treatment regimens requires selection of antibiotics and dosing schedules from a large design space of possibilities. Incomplete knowledge of antibiotic and host immune dynamics in tuberculosis granulomas impacts clinical trial design and success, and variations among clinical trials hamper side-by-side comparison of regimens. Our objective is to systematically evaluate the efficacy of isoniazid and rifampin regimens, and identify modifications to these antibiotics that improve treatment outcomes. We pair a spatio-temporal computational model of host immunity with pharmacokinetic and pharmacodynamic data on isoniazid and rifampin. The model is calibrated to plasma pharmacokinetic and granuloma bacterial load data from non-human primate models of tuberculosis and to tissue and granuloma measurements of isoniazid and rifampin in rabbit granulomas. We predict the efficacy of regimens containing different doses and frequencies of isoniazid and rifampin. We predict impacts of pharmacokinetic/pharmacodynamic modifications on antibiotic efficacy. We demonstrate that suboptimal antibiotic concentrations within granulomas lead to poor performance of intermittent regimens compared to daily regimens. Improvements from dose and frequency changes are limited by inherent antibiotic properties, and we propose that changes in intracellular accumulation ratios and antimicrobial activity would lead to the most significant improvements in treatment outcomes. Results suggest that an increased risk of drug resistance in fully intermittent as compared to daily regimens arises from higher bacterial population levels early during treatment. In conclusion, our systems pharmacology approach complements efforts to accelerate tuberculosis therapeutic development.

60 APPLIED LIFE SCIENCES↗

ASC L2 Milestone Review CTS-2 Production Readiness [Slides]

Open and secure CTS-2 are ready for production capability to include accessibility for users and integration with existing storage and network environment. Hardware testing, early user testing is accomplished, and the systems are generally available for ASC production users.

97 MATHEMATICS AND COMPUTING↗

Combining Novel Simulation Methods and Nucleation Theory to Uncover the Secrets of Gas Hydrates

Conventional computer simulation methods fail for some of the most important problems. With the design and application of innovative algorithms, this project achieved a breakthrough for the case of systems undergoing first-order phase transitions. We gave a complete simulation protocol based upon a well optimized version of our "generalized replica exchange method". The transition of primary interest was gas hydrate formation, a process of significance for climate science and natural gas retrieval. Since hydrates consist of guest molecules in the cages of a water matrix, β ice, the freezing and melting of water was also studied. New information was uncovered about the transition pathways and thermodynamics. Some highlights are 1. the finding that in a very dilute solution without deep supercooling, representative of real-world conditions and very challenging to conventional algorithms, methane can act as a catalyst to drive the formation of large amounts of β ice with empty cages as metastable intermediates, which might be filled by additional methane in a mechanism for hydrate formation, and 2. illumination of the role of metastable cubic ice in water freezing, with determination of the surface tensions of the cubic, hexagonal, and β ices, and the free energy difference of cubic vs hexagonal ice. Work was begun on lipid systems, bilayers and nanoreactors promising for energy-related photoreductions, and targets for future research. Our methods yielded what is arguably the most complete description of the composite lipid/water phases and the transition pathways among them.

36 MATERIALS SCIENCE↗

High-Strength, High-Ductility, High Entropy Alloys with High-Efficiency Native Oxide Solar Absorbers for Concentrating Solar Power Systems

This EPSCoR Project has been investigating the synergy between the excellent high-temperature mechanical behavior of FeMnNiAlCr high entropy alloys (HEA) and the high solar absorptance of their native oxides for high efficiency concentrated solar thermal power (CSP) systems working at >700°C. While HEAs have attracted substantial interest in recent years, most investigations have focused on their applications as structural materials rather than functional materials. This EPSCoR project discovered that FeMnNiAlCr HEAs can potentially be applied synergistically as both a structural and functional material for high-efficiency concentrating solar thermal power (CSP) systems working at >700°C. The HEA itself would be used in high-temperature tubing to carry molten salts or supercritical CO 2 , while its surface oxide would act as a high-efficiency solar thermal absorber. With Fe and Mn being the major components in these HEAs (adding up to ~70 at.% of the alloy), these materials are much more cost-effective than the Ni-based superalloys currently being investigated for high-temperature CSP systems. Through this research, these Fe-Mn based HEAs have demonstrated yield strengths 2-3x greater than that of stainless steel at 700°C and a creep lifetime >800 h at 700ºC under a typical CSP tubing mechanical load of 35 MPa. Their Mn-rich surface oxides maintain a high optical-to-thermal conversion efficiency of ~87% under 1000x solar concentration ratio for 20 simulated day-night thermal cycles between 750ºC and room temperature. In preliminary corrosion studies, these HEAs have sustained immersion in unpurified bromide molten salts for 14 days at 750°C with <2% weight loss, in contrast to 70% weight loss from a 316 stainless steel reference. The simultaneous achievement of promising mechanical, optical, and thermochemical properties in this FeMnNiAlCr system opens the door to new applications of HEAs in solar energy harvesting. Partnerships with Ames Laboratory and Oak Ridge National Laboratory (ORNL) also advanced our understanding of the fundamental structure-property relationships through atomic scale material characterization and first-principles computational modeling. The key research results in this project can potentially be extended to other HEAs and their native oxides. In terms of applications, the proposed FeMnNiAlCr HEA/native oxide system could potentially exceed the mechanical and the optical performance of existing tubing and solar coating materials under EERE’s CSP program at lower cost, which also aligns well with the EPSCoR Science and Technology strategies of New Hampshire in boosting the deployment of renewable energy.

14 SOLAR ENERGY↗

Infrastructure for Multiphysics Software Integration in High Performance Computing-Aided Science and Engineering

The project described in this report constructed and exercised an innovative multiphysics coupling toolkit called the Illinois Rocstar MultiPhysics Application Coupling Toolkit (IMPACT). IMPACT is an open source, flexible, natively parallel infrastructure for coupling multiple uniphysics simulation codes into multiphysics computational systems. IMPACT works with codes written in several high-performance-computing (HPC) programming languages, and is designed from the beginning for HPC multiphysics code development. It is designed to be minimally invasive to the individual physics codes being integrated, and has few requirements on those physics codes for integration. The goal of IMPACT is to provide the support needed to enable coupling existing tools together in unique and innovative ways to produce powerful new multiphysics technologies without extensive modification and rewrite of the physics packages being integrated. There are three major outcomes from this project: 1) construction, testing, application, and open-source release of the IMPACT infrastructure, 2) production of example open-source multiphysics tools using IMPACT, and 3) identification and engagement of interested organizations in the tools and applications resulting from the project. This last outcome represents the incipient development of a user community and application echosystem being built using IMPACT. Multiphysics coupling standardization can only come from organizations working together to define needs and processes that span the space of necessary multiphysics outcomes, which Illinois Rocstar plans to continue driving toward. The IMPACT system, including source code, documentation, and test problems are all now available through the public gitHUB.org system to anyone interested in multiphysics code coupling. Many of the basic documents explaining use and architecture of IMPACT are also attached as appendices to this document. Online HTML documentation is available through the gitHUB site. There are over 100 unit tests provided that run through the Illinois Rocstar Application Development (IRAD) lightweight testing infrastructure that is also supplied along with IMPACT. The package as a whole provides an excellent base for developing high-quality multiphysics applications using modern software development practices. To facilitate understanding how to utilize IMPACT effectively, two multiphysics systems have been developed and are available open-source through gitHUB. The simpler of the two systems, named ElmerFoamFSI in the repository, is a multiphysics, fluid-structure-interaction (FSI) coupling of the solid mechanics package Elmer with a fluid dynamics module from OpenFOAM. This coupling illustrates how to combine software packages that are unrelated by either author or architecture and combine them into a robust, parallel multiphysics system. A more complex multiphysics tool is the Illinois Rocstar Rocstar Multiphysics code that was rebuilt during the project around IMPACT. Rocstar Multiphysics was already an HPC multiphysics tool, but now that it has been rearchitected around IMPACT, it can be readily expanded to capture new and different physics in the future. In fact, during this project, the Elmer and OpenFOAM tools were also coupled into Rocstar Multiphysics and demonstrated. The full Rocstar Multiphysics codebase is also available on gitHUB, and licensed for any organization to use as they wish. Finally, the new IMPACT product is already being used in several multiphysics code coupling projects for the Air Force, NASA and the Missile Defense Agency, and initial work on expansion of the IMPACT-enabled Rocstar Multiphysics has begun in support of a commercial company. These initiatives promise to expand the interest and reach of IMPACT and Rocstar Multiphysics, ultimately leading to the envisioned standardization and consortium of users that was one of the goals of this project.

42 ENGINEERING↗

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Visualization for Insight and Data Analysis in Energy Research

This talk explores how advanced visualization technologies are transforming analytical reasoning and knowledge discovery in energy research, drawing on recent work at the National Laboratory of the Rockies' Computational Science Center. Through a series of scientific case studies, we demonstrate how immersive and high-resolution visualization environments enable scientists and engineers to identify previously unseen patterns and features - insights that often remain hidden in traditional desktop-based analysis. By embedding richer information into interactive analytics tools, these approaches support the exploration of complex, multivariate parameter spaces, where interaction itself catalyzes understanding. Beyond capability, we emphasize the critical role of visualization design grounded in perception and cognition, showing how visual encodings directly influence analytical outcomes. Spanning applications from materials science to integrated energy systems, these visualization approaches accelerate innovation and improve decision-making by enabling deeper, more reliable insight into increasingly complex energy data.

97 MATHEMATICS AND COMPUTING↗

Probing the thermal decomposition of plutonium (III) oxalate with IR and Raman spectroscopy, X-ray diffraction, and electron microscopy

The thermal decomposition of Pu(III) oxalate was analyzed by Raman microspectroscopy, infrared spectroscopy, scanning electron microscopy, and powder X-ray diffraction. These data show that crystalline Pu 2 (C 2 O 4 ) 3 •9H 2 O progressively loses water and oxalate ligands as it is heated, which leads to a decrease in long-range lattice ordering, though minimal changes are observed in gross crystalline morphology. The onset of PuO 2 formation was observed between 200 - 250 °C. Thermal decomposition of oxalate ligands leads to the formation of CO 2 and plutonium oxalate-carbonate moieties, which had not been observed in previously published thermogravimetric measurements of Pu(III) oxalate. Formation of plutonium oxalate-carbonate moieties is believed to be associated with a change in the plutonium oxidation state from 3+ to 4+, which occurs prior to PuO 2 formation. The data provided herein demonstrate the rich spectroscopic nature of a rather underexplored, and technologically relevant, plutonium system. Ideally these results will further future investigations into the Pu(III) oxalate system both experimentally and computationally.

36 MATERIALS SCIENCE↗

Collective dynamics and long-range order in thermal neuristor networks

Abstract In the pursuit of scalable and energy-efficient neuromorphic devices, recent research has unveiled a novel category of spiking oscillators, termed “thermal neuristors.” These devices function via thermal interactions among neighboring vanadium dioxide resistive memories, emulating biological neuronal behavior. Here, we show that the collective dynamical behavior of networks of these neurons showcases a rich phase structure, tunable by adjusting the thermal coupling and input voltage. Notably, we identify phases exhibiting long-range order that, however, does not arise from criticality, but rather from the time non-local response of the system. In addition, we show that these thermal neuristor arrays achieve high accuracy in image recognition and time series prediction through reservoir computing, without leveraging long-range order. Our findings highlight a crucial aspect of neuromorphic computing with possible implications on the functioning of the brain: criticality may not be necessary for the efficient performance of neuromorphic systems in certain computational tasks.

Science & Technology - Other Topics↗

Subtle metastability of the layered magnetic topological insulator MnBi 2 Te 4 from weak interactions

Abstract Layered quantum materials can host interesting properties, including magnetic and topological, for which enormous computational predictions have been done. Their thermodynamic stability is much less visited computationally, which however determines the existence of materials and can be used to guide experimental synthesis. MnBi 2 Te 4 is one of such layered quantum materials that was predicted to be an intrinsic antiferromagnetic topological insulator, and later experimentally realized but in a thermodynamically metastable state. Here, using a combined first-principles-based approach that considers lattice, charge, and spin degrees of freedom, we investigate the metastability of MnBi 2 Te 4 by calculating the Helmholtz free energy for the reaction Bi 2 Te 3 + MnTe → MnBi 2 Te 4 . We identify a temperature range (~500–873 K) in which the compound is stable with respect to the competing binary phases, consistent with experimental observation. We validate the predictions by comparing the calculated specific heats contributed from different degrees of freedom with experimental results. Our findings indicate that the degrees of freedom responsible for the van der Waals interaction, lattice vibration, magnetic coupling, and nontrivial band topology in MnBi 2 Te 4 not only enable emergent phenomena but also play a crucial role in determining its thermodynamic stability. This conclusion lays the foundation for the future computational material synthesis of novel layered systems.

36 MATERIALS SCIENCE↗

Pore-scale simulation of drying in porous media using a hybrid lattice Boltzmann: pore network model

In this work, a hybrid method coupling a pseudo-potential lattice Boltzmann model (LBM) and a pore network model (PNM) to simulate drying in porous media is proposed. Based on the watershed method, the porous medium is firstly decomposed into pore regions. According to the liquid–vapor phase distribution at a given time, the pore regions are further divided into four pore types, namely two-phase pores where a liquid–vapor interface exists, buffer pores next to the two-phase pores, single-liquid and single-vapor phase pores. The pseudo-potential LBM is used in the two-phase and buffer pores to simulate liquid drying and track the movement of the interfaces, while the single-phase PNM simulations are conducted in the buffer and single-phase pores to simulate vapor or liquid flow. LBM and PNM are coupled in the buffer pores through exchange of boundary information. The hybrid method is applied to simulate liquid drying in a porous medium. The whole-domain LBM simulation is considered as the reference solution to validate the hybrid method. Liquid saturation variation during the drying process and detailed phase and pressure distributions obtained by the two methods match quite well, demonstrating the accuracy of the hybrid method. For the specific case studied, the hybrid method saves more than 60% computational time compared to the whole-domain LBM simulation. In addition, the speedup of the hybrid method becomes more significant for a larger computational domain. In summary, the hybrid method developed in this work combines the accuracy of LBM and the efficiency of PNM to simulate drying in porous media at pore scale and can lead to significant reduction of computation time, thus allowing the pore-scale consideration of drying in larger porous systems.

36 MATERIALS SCIENCE↗