SEARCH · Engineering Papers
Results for “shielding effectiveness”
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.
Irradiation creep measurement and microstructural analysis of chromium nitride–coated zirconium alloy using pressurized tubes
Environmental barrier coatings for Zr-based materials are currently under development to reduce oxidation and embrittlement in light-water reactors. Chromium nitride is one such candidate for this application, particularly as accident-tolerant fuel cladding. However, quantifying the impact of coatings on the irradiation-induced creep of zircaloy (Zry) is critical as this mechanism often exceeds thermal creep rates under light-water reactor operating conditions and can be a limiting design characteristic. Additionally, examining irradiation effects in the microstructure at the coating interface is key to understanding the compatibility of the material system. Here, to accelerate the experimental measurement of irradiation creep and microstructure evolution in CrN-Zry, compact, pressurized creep tubes were fabricated and irradiated in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. Miniature, thin-walled rodlets fabricated from annealed Zr-Sn barstock were coated with CrN using physical vapor deposition (PVD) to nominal thicknesses of 4 and 8 μm. Coated and uncoated rodlet specimens were internally pressurized and welded, generating nominal circumferential hoop stresses of 0, 90, or 180 MPa under 300°C irradiation conditions. Twelve specimens were measured diametrically prior to irradiation using a low-cost, automated, contactless laser profilometer developed for this work. Specimens were irradiated in sealed capsules for one 25-day HFIR cycle, accumulating approximately 1.8 $\times$ 10 21 n/cm 2 fast fluence (E n > 10 MeV). The irradiated samples were retrieved and remeasured using the same profilometry system in a shielded hot cell facility. Irradiation creep between specimens was compared using standard statistical tests and showed that both thicknesses of CrN coating had a negligible effect on the irradiation creep strain of the Zry material. Microstructure characterization of pre- and post-irradiated CrN-Zry specimens showed minimal changes due to irradiation but did show a substantial O-rich region at the Zry-CrN interface.
Understanding asymmetries using integrated simulations of capsule implosions in low gas-fill hohlraums at the National Ignition Facility
Current capsule implosions at the National Ignition Facility (NIF) using high-density-carbon ablators and laser energies close to 2 MJ have shown neutron yields in excess of 50 kJ. Improving on this performance requires understanding of the different degradation mechanisms. For many NIF implosions, nuclear diagnostic signatures have inferred residual hot-spot velocities that correlate with fuel areal density variations consistent with a low-order mode-1 asymmetry. A current working hypothesis attributes these asymmetries to a combination of beam to beam variations in the laser delivery and possibly coupling to target features, such as the diagnostic holes needed for x-ray imaging. Recently, a new source has been identified, thickness variation in the ablator shell. To gain better understanding and eventually mitigate the causes of <ρR> asymmetries, 3D integrated simulations using the actual delivered laser powers are needed. To capture the effect of the diagnostic holes (DHs) using direct numerical simulation would require significantly large computational resources. Instead, our 3D simulations make use of a subgrid model developed using highly-resolved 2D simulations that include several details of the DH engineering complexity. Simulations of NIF shots using hohlraums without DHs, to isolate the effect of beam-to-beam variations, reproduce fairly well the observed nuclear diagnostic signatures. Similarly, reasonable agreement between data and simulations is also obtained in the presence of diagnostic holes. To account for the remaining discrepancies a sensitivity study of ablator thickness variation showed that 1% thickness asymmetries are comparable in effect to 1% peak drive mode-1 asymmetries. Furthermore, this study identified sensitivity to variations in the imbedded doped layer (needed to shield the DT ice from the high energy x-rays generated in the hohlraum) thickness even when the inner and outer surface of the ablator are perfectly spherical.
Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels
Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.
Substation Cable Layouts for EMP Coupling Analysis
Direct coupling of early-time high-altitude electromagnetic pulse (HEMP) to substation control cables is simulated for cable layouts based on surveys of seven electrical substations in the United States. An analytic transmission line modeling code is used to estimate worst-case coupled current at the terminations of cable segments in or near the control shack. Where applicable, an induced voltage due to cable shield grounding is also estimated. Various configurations are simulated, including cables with different elevations, lengths, radii, and terminations. Plots of the coupled HEMP effects are given, and general relationships between these effects and the substations geometric and material parameters are highlighted and discussed.
Strategies for mitigating residual magnetic field effect on pre-production PIP-II SSR2 cryomodule performance
This paper outlines the strategy aimed at mitigating the adverse effects of residual magnetic fields on the performance of pre-production SSR2 superconducting cavities within the context of the PIP-II project at Fermilab. Residual magnetic fields can significantly impact cavity performance, leading to reduced quality factor. To address this challenge, our strategy integrates various approaches including magnetic shielding, careful selection of materials, quality controls aimed at measuring magnetic permeability, magnetic hygiene to reduce residual magnetic field at the installation phase. Additionally, experimental studies are being planned to analyze the behavior of the cavities under different magnetic field conditions, and the effectiveness of advanced demagnetization procedures.
Effect of parallel flow on resonant layer responses in high beta plasmas
Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.
Characterization of a boron-loaded deuterated liquid scintillator for fast and thermal neutron detection
The preparation and characterization of novel boron-loaded deuterated liquid scintillators (BLDLS) are presented for the first time to investigate the performance of a detector sensitive to both fast and thermal neutrons intended for future neutron spectroscopy measurements. Three deuterated toluene-based scintillation cocktails were produced, one without boron, one loaded with ortho-carborane (natural boron isotopic abundance), and one with 96 wt% 10B enriched ortho-carborane. General optical and material properties were analyzed to include the composition, density, index of refraction, photoluminescence emission spectra, and absorbance spectra in addition to the relative light yield, resolution, relative efficiency, and pulse-shape-discrimination performance for all three scintillators and a hydrogen-based EJ-309 scintillator to evaluate the effects of boron loading in comparison with a common commercially-available scintillator and feasibility for future neutron spectroscopy measurements. Finally, the detection of shielded neutron sources was explored to capitalize on the 10B thermal neutron capture reaction. The results showed loading the scintillators with boron deteriorated the PSD performance slightly compared to the unloaded deuterated scintillator and resulted in a relative light yield loss of 26% and 41% for the 10B enriched and natural boron BLDLS, respectively. Lastly, the BLDLS have a higher relative detection efficiency than the unloaded deuterated scintillator for shielded neutron sources, based on total count rate, a desirable characteristic for missions related to finding and characterizing illicitly-trafficked special nuclear material.
Bayesian Inference with Latent Hamiltonian Neural Networks (L-HNNs)
When sampling for Bayesian inference, one popular approach is to use Hamiltonian Monte Carlo (HMC) and the No-U-Turn Sampler (NUTS). However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice. We propose Hamiltonian neural networks (HNNs) with HMC and NUTS for solving Bayesian inference problems [1, 2]. Once trained, HNNs do not require gradients of the target density while sampling. Moreover, they satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well suited for use within HMC and NUTS because stationarity can be shown. We also propose an HNN extension called latent HNNs (L-HNNs), which predict latent variable outputs. Compared to HNNs, L-HNNs offer improved expressivity and a reduction in integration errors. Finally, we propose employing L-HNNs in NUTS with an online error monitoring scheme to prevent degeneracy of the sampling in regions of low probability density. We demonstrate L-HNNs in NUTS with online error monitoring by using several example cases involving complex, heavy-tailed, and high local curvature probability densities. Overall, L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. Compared to traditional NUTS, L-HNNs in NUTS with online error monitoring improved the effective sample size (ESS) per gradient by an order of magnitude.
Demonstration of Super-X divertor exhaust control for transient heat load management in compact fusion reactors
Nuclear fusion could offer clean, abundant energy. However, managing the power exhausted from the core fusion plasma towards the reactor wall remains a major challenge. This is compounded in emerging compact reactor designs promising more cost-effective pathways towards commercial fusion energy. Alternative Divertor Configurations (ADCs) are a potential solution. In this work, we demonstrate exhaust control in ADCs, employing a novel method to diagnose the neutral gas buffer, which shields the target. Our work on the Mega Ampere Spherical Tokamak Upgrade shows that ADCs tackle key risks and uncertainties for fusion energy. Their highly reduced sensitivity to perturbations enables active exhaust control in otherwise unfeasible situations and facilitates an increased passive absorption of transients, which would otherwise damage the divertor. We observe a strong decoupling of each divertor from other reactor regions, enabling near-independent control of the divertors and core plasma. Our work showcases the real-world benefits of ADCs for effective heat load management in fusion power reactors.
Neural network reconstruction of the DIII-D tokamak plasma boundary using a reduced set of diagnostics
This study investigates the feasibility of reconstructing the last closed flux surface in the DIII-D tokamak using neural network models trained on reduced input feature sets, addressing an ill-posed task. Two models are compared: one trained solely on coil currents and another incorporating coil currents, plasma current and loop voltage. The model trained exclusively on coil currents achieved a mean point displacement of $0.04$ m on a held-out test set, while the inclusion of plasma current and loop voltage reduced the error to $0.03$ m. This comparison highlights the trade-offs between input feature complexity and reconstruction accuracy, demonstrating the potential of machine learning algorithms to perform effectively in data-limited environments, such as those expected in fusion power plants due to diagnostic constraints imposed by the presence of blankets and shielding.
Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges
The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.
Tunable Dopants with Intrinsic Counterion Separation Reveal the Effects of Electron Affinity on Dopant Intercalation and Free Carrier Production in Sequentially Doped Conjugated Polymer Films
Carrier mobility in doped conjugated polymers is limited by Coulomb interactions with dopant counterions. This complicates studying the effect of the dopant’s oxidation potential on carrier generation because different dopants have different Coulomb interactions with polarons on the polymer backbone. In this work, dodecaborane (DDB)-based dopants are used, which electrostatically shield counterions from carriers and have tunable redox potentials at constant size and shape. DDB dopants produce mobile carriers due to spatial separation of the counterion, and those with greater energetic offsets produce more carriers. Neutron reflectometry indicates that dopant infiltration into conjugated polymer films is redox-potential-driven. Remarkably, X-ray scattering shows that despite their large 2-nm size, DDBs intercalate into the crystalline polymer lamellae like small molecules, indicating that this is the preferred location for dopants of any size. These findings elucidate why doping conjugated polymers usually produces integer, rather than partial charge transfer: dopant counterions effectively intercalate into the lamellae, far from the polarons on the polymer backbone. Finally, it is shown that the IR spectrum provides a simple way to determine polaron mobility. Overall, higher oxidation potentials lead to higher doping efficiencies, with values reaching 100% for driving forces sufficient to dope poorly crystalline regions of the film.
Bubble formation in helium-implanted nanostructured ferritic alloys at elevated temperatures
Helium bubble formation was examined by scanning/transmission electron microscopy (S/TEM) in Fe-9/10Cr binary alloys and two dispersion strengthened nanostructured alloys (CNA3 and 14YWT containing 5–10 nm diameter carbide and oxide particles, respectively) after ex-situ and in-situ He implantation to ~10,000 appm at 500 to 900 °C. The combination of high-resolution STEM images and electron energy loss spectroscopy (EELS) revealed that the Y-Ti-O nanoparticles in 14YWT were uniformly distributed and exhibited a one-to-one relationship for bubble attachment to the nanoclusters. In the in-situ experiment at 900 °C, grain boundary cracking was severe in the Fe-10Cr model alloy, but not in the nanostructured alloys. From 500 to 900 °C, the bubble size generally increased with increasing irradiation temperature, while the bubble density decreased with increasing temperature. At the same temperatures, the bubble size in the implanted materials was in the order of Fe-9/10Cr > CNA3 > 14YWT, while the bubble density showed the opposite order. The observed bubble number densities for the nanostructured alloys are comparable to the nanoparticle density, suggesting that the nanoparticles in both alloys were effective in trapping He. Overall, our results indicate that very high He concentrations can be managed in nanostructured alloys by sequestering the helium into smaller bubbles (which leads to a lower volume swelling value) and to shield He from the grain boundaries. This can be attributed to the much higher sink strength associated with the nanoclusters or the He trapping ability between different types of nanoclusters.
Spin Rectification and Electrically Controlled Spin Transport in Molecular-Ferroelectrics-Based Spin Valves
We report a spin-rectification effect in a spin-valve structure consisting of ferroelectric croconic acid (C 5 H 2 O 5 ) sandwiched between ferromagnetic electrodes La 0.7 Sr 0.3 MnO 3 and Co, which can be switched between a high-resistance (off ) and a low-resistance (on) state by a poling voltage. In the off state, the magnetoresistance (MR) sign reverses with the measurement voltage with a 0.1-V offset, suggesting a spin-rectification behavior, while in the on state the MR remains negative. These observations can be understood in terms of electrically controlled interfacial energy-band alignment either from the electrostatic effect or from the interfacial redox process. In conclusion, the observed spin-rectification effect suggests thepossibility of diodelike devices for spin-polarized current.
VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2
The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.
Evaluation of silicon and indium gallium arsenide photodiodes as direct timing detectors for pulsed x-ray systems
Benchtop pulsed x-ray systems are commonly used to record dynamic material data on the order of nanoseconds, but pulse timing is often difficult to accurately determine. This study demonstrates that commercially available photodiodes can be used effectively for direct x-ray pulse detection without the need for visible light scintillators. X-ray pulses from four commercially available flash x-ray systems were quantified using one silicon and two indium gallium arsenide (InGaAs) photodiodes. The measured InGaAs pulse durations were strongly dependent on radiation dose in the non-linear operating regime, so the photodiodes were shielded to operate below the 2.5 V non-linear regime threshold. The average pulse duration and pulse arrival time jitter of the photodiodes for each x-ray source were within several nanoseconds with the exception of two sets of measurements that were affected by low instrument sensitivity and electrical noise. These results show that InGaAs photodiodes can be used as effective and repeatable stand-alone timing diagnostics for x-ray pulses as short as 20 ns or less.
High-resolution ultra-low field magnetic resonance imaging with a high-sensitivity sensing coil
We present high-resolution magnetic resonance imaging (MRI) at ultra-low field (ULF) with a proton Larmor frequency of around 120 kHz. The key element is a specially designed high-sensitivity sensing coil in the shape of a solenoid with a few millimeter gap between windings to decrease the proximity effect and, hence, increase the coil’s quality ([Formula: see text]) factor and sensitivity. External noise is strongly suppressed by enclosing the sensing coil in a copper cylindrical shield, large enough not to negatively affect the coil’s [Formula: see text] factor and sensitivity, measured to be 217 and 0.47 fT/Hz[Formula: see text], respectively. To enhance small polarization of proton spins at ULF, a strong pulsed 0.1 T prepolarization field is applied, making the signal-to-noise ratio (SNR) of ULF MRI sufficient for high-quality imaging in a short time. We demonstrate ULF MRI of a copper sulfate solution phantom with a resolution of [Formula: see text] and SNR of 10. The acquisition time is 6.3 min without averaging. The sensing coil size in the current realization can accommodate imaging objects of 9 cm in size, sufficient for hand, and it can be further increased for human head imaging in the future. Since the in-plane resolution of [Formula: see text] is typical in anatomical medical imaging, this ULF MRI method can be an alternative low-cost, rapid, portable method for anatomical medical imaging of the human body or animals. This ULF MRI method can supplement other MRI methods, especially when such methods are restricted due to high cost, portability requirement, imaging artifacts, and other factors.