SEARCH · Engineering Papers
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.
The data in this repository correspond to the figures of the paper "Impact of error fields and error field correction on heat fluxes in SPARC". They are created by simulations run with the M3D-C1 code, the MAFOT code and python postprocessing tools.
The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.
The 2x2 Demonstrator, a prototype for the Deep Underground Neutrino Experiment (DUNE) liquid argon (LAr) Near Detector, was exposed to the Neutrinos from the Main Injector (NuMI) neutrino beam at Fermi National Accelerator Laboratory (Fermilab). This detector is a prototype of a new modular design for a liquid argon time-projection chamber (LArTPC), comprising a two-by-two array of four modules, each further segmented into two optically isolated LArTPCs. The 2x2 Demonstrator features a number of pioneering technologies, including a low-profile resistive field shell to establish drift fields, native 3D ionization pixelated imaging, and a high-coverage dielectric light readout system. The 2.4-tonne active mass detector is flanked upstream and downstream by supplemental solid-scintillator tracking planes, repurposed from the MINERvA experiment, which track ionizing particles exiting the argon volume. The antineutrino beam data collected by the detector over a 4.5 day period in 2024 include over 30,000 neutrino interactions in the LAr active volume—the first neutrino interactions reported by a DUNE detector prototype. During its physics-quality run, the 2x2 Demonstrator operated at a nominal drift field of 500 V/cm and maintained good LAr purity, with a stable electron lifetime of approximately 1.25 ms. This paper describes the detector and supporting systems, summarizes the installation and commissioning, and presents the initial validation of collected NuMI beam and off-beam self-triggers. In addition, it highlights observed interactions in the detector volume, including candidate muon antineutrino events.
Radiofrequency quadrupole (RFQ) linear accelerators with variable aperture cell parameters offer an effective means of accelerating high-intensity ion beams exceeding 100 mA and have been implemented in various RFQ designs. However, in a four-rod RFQ, the use of dielectric perturbation methods to measure the voltage distribution behind the rod electrodes is not applicable when implementing variable aperture geometries, because changes in the rod shape influence the electric field at the location of the perturbation element. Therefore, it was necessary to design the rod profile for each cell such that the electric field distribution behind the rods remains unaffected by changes in the aperture. In this study, we present a method to adjust the voltage distribution by modifying the geometry of the rod electrodes according to each cell's aperture parameters. This design approach enables the integration of variable aperture cells into a four-rod RFQ linac, significantly enhancing the achievable beam current. To ensure uniform capacitance and maintain a constant rod-electrode voltage despite variations in the average aperture radius, we developed a systematic method for shaping the rod electrodes. This method was refined using two-dimensional (2D) and three-dimensional (3D) electromagnetic field simulations to optimize the electrode cross-sectional profiles. The fabricated rod electrodes were then installed in a four-rod RFQ, and the electric field distribution was measured using the perturbation method and tuned by adjusting the base plates. As a result, the measured variation of ±1.1 % confirms that this design methodology can produce a sufficiently uniform voltage distribution for high-current beam acceleration. These results demonstrate the practical feasibility of a variable-aperture four-rod RFQ linac.
We advance soft X-ray vector ptychographic tomography to map the 3D magnetization field in self-assembled superparamagnetic nanoparticles at a liquid–liquid interface, revealing how layered structures influence magnetic ordering. We observe that monolayers with low coordination numbers exhibit weak magnetic order, with magnetic vortices disrupting spin alignment. In contrast, bilayers and trilayers with higher coordination numbers display long-range magnetic order with strong spin correlations across larger distances and a suppression of magnetic vortices. We further quantify the average distance for vortex–antivortex pairs as 26.0 ± 2.0 nm, while vortex–vortex and antivortex–antivortex pairs exhibit larger separations, averaging 44.9 ± 5.2 and 54.1 ± 7.4 nm, respectively. These experimental results are supported by micromagnetic Monte Carlo simulations. Our findings illustrate how layered structures enhance magnetic order and spin correlation in superparamagnetic nanoparticle assemblies, providing a promising approach for tuning magnetic properties in applications such as data storage, microrobotics, and biomedicine.
We present an overview of a novel electromagnetic coil configuration for stellarators and its application to two near-term fusion systems. The novel coil configuration is the planar coil stellarator, able to implement precisely-quasisymmetric 3D magnetic fields using a set of planar, plasma-encircling coils and a set of planar, field-shaping coils situated on a surface surrounding the plasma. This configuration combines the stellarator’s advantages of steady-state operation, stability, low recirculating power fraction, and a mature physics basis, with the benefits of using simpler, planar coils which allow for a maintenance scheme leveraging large ports, and the ability to control magnets individually. The initial near-term use case considered is a steady state deuterium–deuterium stellarator neutron source, called Eos. The second near-term use case considered is a deuterium–tritium stellarator fusion pilot plant, called Helios, that would be approximately twice the linear dimension of the Eos design.
It is demonstrated that the peak toroidal impact phase of the post-disruption runaway electrons (REs) can be varied shot-to-shot by means of applied static non-axisymmetric (3D) magnetic fields, or resonant magnetic perturbations (RMPs). In the experiments, high-current (500 kA), purged RE plateaus (i.e. low impurity post-disruption plasmas with current carried completely by REs) were terminated on the DIII-D tokamak center post (CP). The global toroidal impact peak of the REs appears to be dominantly set by the toroidal phase of a wall-locked tearing mode, typically (m/n) = (2/1), but sometimes (m/n) = (3/2). In the absence of external RMP, the locked mode reliably locks close to a preferred toroidal phase Φ ≃ 75◦ probably set by intrinsic error fields. With applied RMP, the locked mode and RE impact phase is pushed to a new location. When the tearing mode is (m/n) = (2/1), the phase is well controlled and tracks the applied RMP phase, but in the case of an (m/n) = (3/2) tearing mode, the new phase location is not well controlled. These results could have application in ITER and other fusion power plant (FPP) - sized tokamaks as a method for spreading disruption first wall damage from REs
QUIC stands for the Quick Urban & Industrial Complex (QUIC) dispersion modeling system. QUIC is a fast response urban dispersion model that runs on a laptop. QUIC is comprised of a 3D wind field model called QUIC-URB, a transport and dispersion model called QUIC-PLUME, and graphical user interface called QUIC-GUI. QUIC also includes QUIC-PRESSURE to solve for pressure fields in and around buildings, a population exposure assessment tool called QUIC-POP, and an indoor infiltration calculator for computing indoor concentrations. Transport and dispersion for different types of airborne contaminants can be computed on building to neighborhood scales in tens of seconds to tens of minutes. QUIC will never give perfect answers, but it will account for the effects of buildings in an approximate way and provide more realism than non-building aware dispersion models.
In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.
We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.
In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.
Three-dimensional (3D) magnetic nanostructures are an emerging platform capable of creating complex topological magnetic fields. The control of localized nanoscale magnetic fields is seen to be of importance for diverse areas from bioapplications such as drug delivery, to particle trapping and controlling Majorana Fermions for quantum computing. Three-dimensional geometric confinement and proximity can create tailor-made spin textures not possible in two dimensions. The control of magnetization afforded here can allow the formation of unique stray field textures. Here, we report the creation of reconfigurable 3D topological magnetic field textures induced by an interwoven 3D nanostructure and applied field protocol. These field textures emerge due to distinct DWs formed in this structure and lead to the creation of an antivortex field, a hexapole cusp and a 3D skyrmion field tube of mixed chirality. In conclusion, our results therefore show a key step toward the design and control of topological magnetic fields on the nanoscale.
Abstract In this paper, a numerical dart leader model has been implemented to understand the leader's development and the corresponding electric field changes observed by the 3D broadband mapping and polarization (BIMAP‐3D) system. The model assumes the extending leader channel is equipotential and has a linear charge distribution induced by an ambient electric field. The charge distribution induced by the ambient field can be used to model the electric field change at the ground. We then find the ambient electric field which best fits the field change measurements at the two BIMAP stations. The estimated ambient electric field decreases in the direction of dart leader propagation. Our observations and modeling results are consistent with our earlier hypothesis that dart leader speed is proportional to the electric field at the leader tip. The model also supports our earlier analysis that leader speed variations near branch junctions were due to previous charge deposits near the junctions. The modeled tip electric field is generally lower than the breakdown field unless the pre‐dart‐leader channel has a significant temperature of 3,000 K. This is consistent with the fact that dart leaders typically do not form new branches into the virgin air. Furthermore, the tip field is generally close to the negative streamer stability field at ambient temperatures, explaining the nature of the narrow and well‐defined channel structure. In addition to the charge distribution and the ambient and tip electric field, the development of the channel potential and current distribution are also presented.
The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.
Dark field x-ray microscopy (DXFM) can visualize microstructural distortions in bulk crystals. Using the femtosecond x-ray pulses generated by x-ray free-electron lasers (XFELs), DFXM can achieve sub-μm spatial resolution and <100 fs time resolution simultaneously. In this paper, we demonstrate ultrafast DFXM measurements at the European XFEL to visualize an optically driven longitudinal strain wave propagating through a diamond single crystal. We also present two DFXM scanning modalities that are new to the XFEL sources: spatial 3D and 2D axial-strain scans with sub-μm spatial resolution. With this progress in XFEL-based DFXM, we discuss new opportunities to study multi-timescale spatiotemporal dynamics of microstructures.
Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.
The purpose of this data services request was to demonstrate integration of the Atmospheric Radiation Measurement (ARM) User Facility’s “met” datastreams with time series data from other earth science data sources using the BASIN-3D data synthesis software tool. BASIN-3D is an open-source Python library that enables researchers to integrate data across configured public and private data sources. It provides a common query language for researchers to request measurement locations and time series data based on specified locations, variables, time period, statistics, aggregation, and data quality. BASIN-3D acquires the data that match the query from each configured data source and translates the results into harmonized vocabularies, thus reducing researchers' data-wrangling effort. In addition, because the queries are executed on demand, researchers can easily regenerate their synthesized data sets as new data and/or data updates become available, eliminating one-off data products. BASIN-3D can output data using a variety of different data structures for end-user applications including Python pandas data frames and hdf5 output formats.
Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.