Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAGNETIC RESONANCE”

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 451 records · Page 25

Experimental evaluation of electrical conductivity imaging of anisotropic brain tissues using a combination of diffusion tensor imaging and magnetic resonance electrical impedance tomography

Anisotropy of biological tissues is a low-frequency phenomenon that is associated with the function and structure of cell membranes. Imaging of anisotropic conductivity has potential for the analysis of interactions between electromagnetic fields and biological systems, such as the prediction of current pathways in electrical stimulation therapy. To improve application to the clinical environment, precise approaches are required to understand the exact responses inside the human body subjected to the stimulated currents. In this study, we experimentally evaluate the anisotropic conductivity tensor distribution of canine brain tissues, using a recently developed diffusion tensor-magnetic resonance electrical impedance tomography method. At low frequency, electrical conductivity of the biological tissues can be expressed as a product of the mobility and concentration of ions in the extracellular space. From diffusion tensor images of the brain, we can obtain directional information on diffusive movements of water molecules, which correspond to the mobility of ions. The position dependent scale factor, which provides information on ion concentration, was successfully calculated from the magnetic flux density, to obtain the equivalent conductivity tensor. By combining the information from both techniques, we can finally reconstruct the anisotropic conductivity tensor images of brain tissues. The reconstructed conductivity images better demonstrate the enhanced signal intensity in strongly anisotropic brain regions, compared with those resulting from previous methods using a global scale factor.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Bulk and Nanocrystalline Cesium Lead-Halide Perovskites as Seen by Halide Magnetic Resonance

Lead-halide perovskites increasingly mesmerize researchers because they exhibit a high degree of structural defects and dynamics yet nonetheless offer an outstanding (opto)electronic performance on par with the best examples of structurally stable and defect-free semiconductors. This highly unusual feature necessitates the adoption of an experimental and theoretical mindset and the reexamination of techniques that may be uniquely suited to understand these materials. Surprisingly, the suite of methods for the structural characterization of these materials does not commonly include nuclear magnetic resonance (NMR) spectroscopy. The present study showcases both the utility and versatility of halide NMR and NQR (nuclear quadrupole resonance) for probing the structure and structural dynamics of CsPbX 3 (X = Cl, Br, I), in both bulk and nanocrystalline forms. The strong quadrupole couplings, which originate from the interaction between the large quadrupole moments of, e.g., the 35 Cl, 79 Br, and 127 I nuclei, and the local electric-field gradients, are highly sensitive to subtle structural variations, both static and dynamic. The quadrupole interaction can resolve structural changes with accuracies commensurate with synchrotron X-ray diffraction and scattering. It is shown that space-averaged site-disorder is greatly enhanced in the nanocrystals compared to the bulk, while the dynamics of nuclear spin relaxation indicates enhanced structural dynamics in the nanocrystals. The findings from NMR and NQR were corroborated by ab initio molecular dynamics, which point to the role of the surface in causing the radial strain distribution and disorder. These findings showcase a great synergy between solid-state NMR or NQR and molecular dynamics simulations in shedding light on the structure of soft lead-halide semiconductors.

36 MATERIALS SCIENCE↗

Surface Terminations of LaAlO 3 Perovskite Nanoparticles as Viewed by Solid-State Nuclear Magnetic Resonance

Nanocrystal surfaces generally undergo reconstructions that differentiate them from the bulk structures, often in nontrivial ways. Understanding these terminations is critical across diverse fields, from heterogeneous catalysis to the formation of topological states and the synthesis of semiconductor nanomaterials. Determining surface structures is currently an interdisciplinary task, most often involving high-resolution electron microscopy and surface electron diffraction. These methods, however, do not provide a global view of the ensemble of structures present in a sample. Here, we show how surface-sensitive solid-state nuclear magnetic resonance (SSNMR) spectroscopy methods can bridge this gap. In this context, we investigated the surface structure of lanthanum aluminate (LaAlO 3 ) perovskite nanoparticles. Four distinct surface terminations have previously been observed for this material, but their relative abundances were unknown. Using an array of double- and triple-resonance SSNMR methods probing the relative proximities of surface 1 H, 27 Al, 17 O, and 139 La nuclei, we conclude the surface to be majority terminated (80%) by AlO x with substantial (20%) LaO x terminated regions.

Materials↗

Cardiac Function After Multimodal Breast Cancer Therapy Assessed With Functional Magnetic Resonance Imaging and Echocardiography Imaging

Purpose: Breast intensity modulated radiation therapy (IMRT) reduces high-dose heart volumes but increases low-dose volumes. We prospectively assessed heart changes after 3D conformal RT (3DCRT) and IMRT for left-sided breast cancer. Heart dose was analyzed individually, 3DCRT patients were moderately exposed, and IMRT was performed only in patients with unacceptably high heart doses upon 3DCRT planning. Methods and Materials: In 49 patients (38 patients received 3DCRT; 11 patients received IMRT; and 20 patients received neoadjuvant or adjuvant chemotherapy) magnetic resonance imaging (MRI) and echocardiography were performed before and at 6, 12, and 24 months after treatment. Results: Mean heart dose for IMRT was 12.9 ± 3.9 Gy versus 4.5 ± 2.4 Gy for 3DCRT. Heart volumes receiving >40 Gy were 2.6% (3DCRT) versus 1.3% (IMRT); doses were >50 Gy only with 3DCRT. Temporary ejection fraction (EF) decrease was observed on MRI after 6 months (63%-59%, P=.005) resolving at 24 months. Only 3 patients had pronounced largely transient changes of EF and left ventricular enddiastolic diameter (LVEDD). Mitral (M) and tricuspid (T) annular plane systolic excursion (MAPSE and TAPSE) were reduced over the whole cohort (still within normal range). After 24 months left ventricular remodeling index decreased in patients receiving chemotherapy (0.80 vs 0.70, P=.028). Neither wall motion abnormalities nor late enhancements were found. On echocardiography, in addition to EF findings that were similar to those on MRI, global strain was unchanged over the whole cohort at 24 months after a transient decrease at 6 and 12 months. Longitudinal strain decreased in the whole cohort after 24 months in some segments, whereas it increased in others. Conclusions: Until 24 months after risk-adapted modern multimodal adjuvant therapy, only subclinical cardiac changes were observed in both 3DCRT patients with inclusion of small to moderate amounts of heart volume in RT tangents and in the patients treated with IMRT and reduced high-dose heart exposure.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Occurrence of Aberrant Lymph Node Spread on Magnetic Resonance Lymphography in Prostate Cancer Patients With a Biochemical Recurrence After Radical Prostatectomy

Purpose: To investigate the pattern of lymph node spread in prostate cancer patients with a biochemical recurrence after radical prostatectomy, eligible for salvage radiotherapy; and to determine whether the clinical target volume (CTV) for elective pelvic irradiation in the primary setting can be applied in the salvage setting for patients with (a high risk of) lymph node metastases. Methods and Materials: The charts of 47 prostate cancer patients with PSA recurrence after prostatectomy who had positive lymph nodes on magnetic resonance lymphography (MRL) were reviewed. Positive lymph nodes were assigned to a lymph node region according to the guidelines of the Radiation Therapy Oncology Group (RTOG) for delineation of the CTV for pelvic irradiation (RTOG-CTV). We defined four lymph node regions for positive nodes outside this RTOG-CTV: the para-aortal, proximal common iliac, pararectal, and paravesical regions. They were referred to as aberrant lymph node regions. For each patient, clinical and pathologic features were recorded, and their association with aberrant lymph drainage was investigated. The distribution of positive lymph nodes was analyzed separately for patients with a prostate-specific antigen (PSA) <1.0 ng/mL. Results: MRL detected positive aberrant lymph nodes in 37 patients (79%). In 20 patients (43%) a positive lymph node was found in the pararectal region. Higher PSA at the time of MRL was associated with the presence of positive lymph nodes in the para-aortic region (2.49 vs. 0.82 ng/mL; p = 0.007) and in the proximal common iliac region (1.95 vs. 0.59 ng/mL; p = 0.009). There were 18 patients with a PSA <1.0 ng/mL. Ten of these patients (61%) had at least one aberrant positive lymph node. Conclusion: Seventy-nine percent of the PSA-recurrent patients had at least one aberrant positive lymph node. Application of the standard RTOG-CTV for pelvic irradiation in the salvage setting therefore seems to be inappropriate.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Comparison of Macroscopic Pathology Measurements With Magnetic Resonance Imaging and Assessment of Microscopic Pathology Extension for Colorectal Liver Metastases

Purpose: To compare pathology macroscopic tumor dimensions with magnetic resonance imaging (MRI) measurements and to establish the microscopic tumor extension of colorectal liver metastases. Methods and Materials: In a prospective pilot study we included patients with colorectal liver metastases planned for surgery and eligible for MRI. A liver MRI was performed within 48 hours before surgery. Directly after surgery, an MRI of the specimen was acquired to measure the degree of tumor shrinkage. The specimen was fixed in formalin for 48 hours, and another MRI was performed to assess the specimen/tumor shrinkage. All MRI sequences were imported into our radiotherapy treatment planning system, where the tumor and the specimen were delineated. For the macroscopic pathology analyses, photographs of the sliced specimens were used to delineate and reconstruct the tumor and the specimen volumes. Microscopic pathology analyses were conducted to assess the infiltration depth of tumor cell nests. Results: Between February 2009 and January 2010 we included 13 patients for analysis with 21 colorectal liver metastases. Specimen and tumor shrinkage after resection and fixation was negligible. The best tumor volume correlations between MRI and pathology were found for T1-weighted (w) echo gradient sequence (r{sub s} = 0.99, slope = 1.06), and the T2-w fast spin echo (FSE) single-shot sequence (r{sub s} = 0.99, slope = 1.08), followed by the T2-w FSE fat saturation sequence (r{sub s} = 0.99, slope = 1.23), and the T1-w gadolinium-enhanced sequence (r{sub s} = 0.98, slope = 1.24). We observed 39 tumor cell nests beyond the tumor border in 12 metastases. Microscopic extension was found between 0.2 and 10 mm from the main tumor, with 90% of the cases within 6 mm. Conclusions: MRI tumor dimensions showed a good agreement with the macroscopic pathology suggesting that MRI can be used for accurate tumor delineation. However, microscopic extensions found beyond the tumor border indicate that caution is needed in selecting appropriate tumor margins.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Screening of resonant magnetic perturbation fields by poloidally varying toroidal plasma rotation

The angular frequency of the subsonic equilibrium toroidal flow in a tokamak plasma is typically assumed constant at magnetic flux surfaces, i.e., the rotation frequency does not vary along the poloidal angle of the flux surface. Yet, there are several cases where this symmetry is broken. An interesting situation is a recently observed complex flow pattern induced by magnetic field line ergodization, in the presence of the tri-dimensional (3-D) resonant magnetic perturbation (RMP) [Schmitz et al., Nucl. Fusion 56, 066008 (2016)]. A new flow model including poloidally varying rotation frequencies has been implemented in the full resistive linear MHD code MARS-F [Liu et al., Phys. Plasmas 7, 3681 (2000)], allowing poloidal variation of the angular frequency of the equilibrium toroidal rotation in a generic toroidal geometry. The effect of this asymmetric flow, on top of a poloidally symmetric toroidal flow, on the plasma response to RMP fields is numerically investigated. It is found that a poloidally varying toroidal flow component enhances the favourable average curvature induced plasma screening of the applied 3-D field, for low toroidal flow velocities. At faster flow, when the resistive-inertial response becomes important, the asymmetric toroidal flow reduces the plasma screening. The most significant effect is found to come from the m=1 component of the poloidal asymmetry in the toroidal rotation frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Introduction to Nuclear Magnetic Resonance

The purpose of this paper is to try to give a short overview of what the status is on nuclear magnetic resonance (NMR). It's a subject where one really has to spend some time to look at the physics in detail to develop a proper working understanding. I feel it's not appropriate to present to you density matrices, Hamiltonians of all sorts, and differential equations representing the motion of spins. I'm really going to present some history and status, and show a few very simple concepts involved in NMR. It is a form of radio frequency spectroscopy and there are a great number of nuclei that can be studied very usefully with the technique. NMR requires a magnet, a r.f. transmitter/receiver system, and a data acquisition system.

Manatt, Stanley L.↗

Characterization of the AT180 epitope of phosphorylated Tau protein by a combined nuclear magnetic resonance and fluorescence spectroscopy approach

Highlights: {yields} pThr231 of the Tau protein is necessary for the binding of the AT180 antibody. {yields} pSer235 of the Tau protein does not interfere with the AT180 recognition of pThr231. {yields} Epitope mapping is efficiently achieved by combining NMR and FRET spectroscopy. -- Abstract: We present here the characterization of the epitope recognized by the AT180 monoclonal antibody currently used to define an Alzheimer's disease (AD)-related pathological form of the phosphorylated Tau protein. Some ambiguity remains as to the exact phospho-residue(s) recognized by this monoclonal: pThr231 or both pThr231 and pSer235. To answer this question, we have used a combination of nuclear magnetic resonance (NMR) and fluorescence spectroscopy to characterize in a qualitative and quantitative manner the phospho-residue(s) essential for the epitope recognition. Data from the first step of NMR experiments are used to map the residues bound by the antibodies, which were found to be limited to a few residues. A fluorophore is then chemically attached to a cystein residue introduced close-by the mapped epitope, at arginine 221, by mutagenesis of the recombinant protein. The second step of Foerster resonance energy transfer (FRET) between the AT180 antibody tryptophanes and the phospho-Tau protein fluorophore allows to calculate a dissociation constant Kd of 30 nM. We show that the sole pThr231 is necessary for the AT180 recognition of phospho-Tau and that phosphorylation of Ser235 does not interfere with the binding.

60 APPLIED LIFE SCIENCES↗

Influence of the Magnetic Moment of the Probe of a Magnetic Resonance Force Microscope on the Spin-Wave Resonance Spectra

Abstract—The ferromagnetic resonance in an array of permalloy microstrips 3000 × 500 × 30 nm{sup 3} in size, which are ordered on a rectangular 3.5 × 6 μm{sup 2} lattice, is studied using magnetic resonance force microscopy with strong probe–sample interaction. This interaction induces intricate modifications of the observed spectra, which are manifested both in line splitting and in changes in shape. The dependences of the observed spectra on the sample–probe distance and the orientation of the magnetic moment of the probe are analyzed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Raman and nuclear magnetic resonance investigation of alkali metal vapor interaction with alkene-based anti-relaxation coating

The use of anti-relaxation coatings in alkali vapor cells yields substantial performance improvements compared to a bare glass surface by reducing the probability of spin relaxation in wall collisions by several orders of magnitude. Some of the most effective anti-relaxation coating materials are alpha-olefins, which (as in the case of more traditional paraffin coatings) must undergo a curing period after cell manufacturing in order to achieve the desired behavior. Until now, however, it has been unclear what physicochemical processes occur during cell curing, and how they may affect relevant cell properties. We present the results of nondestructive Raman-spectroscopy and magnetic-resonance investigations of the influence of alkali metal vapor (Cs or K) on an alpha-olefin, 1-nonadecene coating the inner surface of a glass cell. It was found that during the curing process, the alkali metal catalyzes migration of the carbon-carbon double bond, yielding a mixture of cis- and trans-2-nonadecene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of Phase Composition of Cobalt Nanoparticles Using {sup 59}Co Internal Field Nuclear Magnetic Resonance

It is well known that cobalt exhibits polymorphism, i.e., the co-existence of both the hcp and fcc phases. In particular, the method of synthesis and other thermodynamic conditions is known to play a crucial role in determining the particular phase of cobalt. In this work, we have compared the phase composition of the cobalt nanoparticles synthesized using two different solvents (water) and ethanol (Co@C). XRD measurements confirm the existence of fcc phase in commercial cobalt nanoparticles (Co@A), co-existence of fcc and hcp phases in Co@B, while the existence of the hcp phase in Co@C. We have studied these cobalt nanoparticles using {sup 59}Co internal field nuclear magnetic resonance (IFNMR) for verification of phase composition. Our studies reveal that the Co@A has fcc as a major phase with minor quantity hcp phase. Co@B exhibits approximately equal amount of fcc and hcp phase while Co@C exhibits hcp as a major phase with minor fcc phase. Our SEM micrograph studies confirm that the cobalt particles have spherical shape in the fcc phase. The cobalt particles exhibit both spherical and dendrite morphology confirming the co-existence of fcc and hcp phases, while the sample with pure hcp phase exhibits the dendrite morphology. Our studies also throw light on understanding the effect of solvent in the phase formation of the cobalt nanoparticles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Role of Early Postradiation Magnetic Resonance Imaging Scans in Children With Diffuse Intrinsic Pontine Glioma

Purpose: To determine optimal timing of assessing postradiation radiographic response on magnetic resonance imaging (MRI) scans in pediatric patients with diffuse intrinsic pontine glioma (DIPG). Methods and Materials: Patients were treated on a prospective study at the National Cancer Institute (Protocol no. 06-C-0219) evaluating the effects of radiotherapy (RT). Standard RT was administered in standard fractionation over 6 weeks. Postradiation MRI scans were performed at 2 and 6-8 weeks. Results: Eleven patients with DIPG were evaluated. Median age was 6 years (range, 4-13 years). Patients were treated with external-beam RT to 55.8 Gy (n = 10) or 54 Gy (n = 1), with a gross tumor volume to planning target volume expansion of 1.8-2.0 cm. All patients received prescribed dose and underwent posttreatment MRI scans at 2 and 6-8 weeks. Pretreatment imaging revealed compression of fourth ventricle (n = 11); basilar artery encasement (n = 9); tumor extension outside the pons (n = 11); and tumor hemorrhage (n = 2). At the 2-week scan, basilar artery encasement improved in 7 of 9 patients, and extent of tumor was reduced in 5 of 11 patients. Fourth ventricle compression improved in 6 of 11 patients but worsened in 3 of 11 patients. Presumed necrosis was observed in 5 of 11 patients at 2 weeks and in 1 additional patient at 6-8 weeks. There was no significant difference in mean anteroposterior and transverse diameters of tumor between the 2- and 6-8-week time points. Six of 11 patients had increasing ventricular size, with no evidence of obstruction. Conclusions: There is no significant difference in tumor size of DIPG patients who have received standard RT when measured at 2 weeks vs. 6-8 weeks after RT. The majority of patients had the largest change in tumor size at the 2-week post-RT scan, with evolving changes documented on the 6-8-week scan. Six of 11 patients had progressive ventriculomegaly without obstruction, suggestive of communicating hydrocephalus. To the best of our knowledge, this is the first documentation of this phenomenon in this cohort of patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A Technique for Generating Volumetric Cine-Magnetic Resonance Imaging

Purpose: The purpose of this study was to develop a techique to generate on-board volumetric cine-magnetic resonance imaging (VC-MRI) using patient prior images, motion modeling, and on-board 2-dimensional cine MRI. Methods and Materials: One phase of a 4-dimensional MRI acquired during patient simulation is used as patient prior images. Three major respiratory deformation patterns of the patient are extracted from 4-dimensional MRI based on principal-component analysis. The on-board VC-MRI at any instant is considered as a deformation of the prior MRI. The deformation field is represented as a linear combination of the 3 major deformation patterns. The coefficients of the deformation patterns are solved by the data fidelity constraint using the acquired on-board single 2-dimensional cine MRI. The method was evaluated using both digital extended-cardiac torso (XCAT) simulation of lung cancer patients and MRI data from 4 real liver cancer patients. The accuracy of the estimated VC-MRI was quantitatively evaluated using volume-percent-difference (VPD), center-of-mass-shift (COMS), and target tracking errors. Effects of acquisition orientation, region-of-interest (ROI) selection, patient breathing pattern change, and noise on the estimation accuracy were also evaluated. Results: Image subtraction of ground-truth with estimated on-board VC-MRI shows fewer differences than image subtraction of ground-truth with prior image. Agreement between normalized profiles in the estimated and ground-truth VC-MRI was achieved with less than 6% error for both XCAT and patient data. Among all XCAT scenarios, the VPD between ground-truth and estimated lesion volumes was, on average, 8.43 ± 1.52% and the COMS was, on average, 0.93 ± 0.58 mm across all time steps for estimation based on the ROI region in the sagittal cine images. Matching to ROI in the sagittal view achieved better accuracy when there was substantial breathing pattern change. The technique was robust against noise levels up to SNR = 20. For patient data, average tracking errors were less than 2 mm in all directions for all patients. Conclusions: Preliminary studies demonstrated the feasibility of generating real-time VC-MRI for on-board localization of moving targets in radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Evaluation of hydrogen storage in sandstone reservoirs using 1 H nuclear magnetic resonance spectroscopy

Evaluation of the hydrogen storage capacity of porous rocks is crucial for underground hydrogen storage. Using 1 H nuclear magnetic resonance (NMR) spectroscopy, we successfully characterized the hydrogen responses and identified storage mechanisms in Berea sandstone under varying water saturation. Here, the results indicate that the injected hydrogen behaves as a free gas phase and is capable of occupying the empty pore volume regardless of the saturation state. No hysteresis was observed during injection and production cycles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of elongation and triangularity on plasma response to resonant magnetic perturbations

Abstract Systematic toroidal modeling of the plasma response to the n = 1–4 ( n is the toroidal mode number) resonant magnetic perturbation (RMP) field is carried out in order to understand the plasma-shaping effect on controlling the type-I edge-localized modes (ELMs) in tokamak experiments. Considered are large variations of the plasma elongation and triangularity at a fixed edge safety factor q a , for limiter plasmas with both single-null (SN) and double-null (DN) divertor-like boundary shapes. Numerical results assuming conformal 3D RMP coils show that (i) the optimum coil phasing between the upper and lower rows for ELM control becomes increasingly sensitive to the plasma elongation with higher- n toroidal spectra, (ii) the optimum coil phasing is however essentially independent of the plasma triangularity for all n = 1–4 RMP fields, (iii) with the same coil current and the optimum coil phasing, high elongation generally favors ELM control but it may be more challenging for plasmas with intermediate elongation and with lower- n ( n = 1–2) RMPs, and (iv) higher triangularity is generally always better for ELM control with all n = 1–4 fields for both DN and SN divertor-like plasma boundary shapes.

Physics↗

Direct Numerical Simulation of Pore-Scale Flow in a Bead Pack: Comparison with Magnetic Resonance Imaging Observations

A significant body of current research is aimed at developing methods for numerical simulation of flow and transport in porous media that explicitly resolve complex pore and solid geometries, and at utilizing such models to study the relationships between fundamental pore-scale processes and macroscopic manifestations at larger (i.e., Darcy) scales. A number of different numerical methods for pore-scale simulation have been developed, and have been extensively tested and validated for simplified geometries. However, validation of pore-scale simulations of fluid velocity for complex, three-dimensional (3D) pore geometries that are representative of natural porous media is challenging due to our limited ability to measure pore-scale velocity in such systems. Recent advances in magnetic resonance imaging (MRI) offer the opportunity to measure not only the pore geometry, but also local fluid velocities under steady-state flow conditions in 3D and with high spatial resolution. In this paper, we present a 3D velocity field measured at sub-pore resolution (tens of micrometers) over a centimeter-scale 3D domain using MRI methods. We have utilized the measured pore geometry to perform 3D simulations of Navier-Stokes flow over the same domain using direct numerical simulation techniques. We present a comparison of the numerical simulation results with the measured velocity field. It is shown that the numerical results match the observed velocity patterns well overall except for a variance and small systematic scaling which can be attributed to the known experimental error in the MRI measurements. The comparisons presented here provide strong validation of the pore-scale simulation methods and new insights for interpretation of uncertainty in MRI measurements of pore-scale velocity. This study also provides a potential benchmark for future comparison of other pore-scale simulation methods.

computational fluid dynamics↗