Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D velocity models”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Three-dimensional Hydrodynamic Simulations of Convective Nuclear Burning in Massive Stars Near Iron Core Collapse

Nonspherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multimessenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the postshock region. Here we report on a set of 4π 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final 10 minutes prior to and including iron core collapse. We consider initial 1D MESA models with masses of 14, 20, and 25 M ⊙ to survey a range of O/Si-shell density and compositional configurations. Further, we characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of α MLT = 1.5. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of ℓ = 1–3 near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of presupernova progenitor models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Magnetotelluric Investigations of the Kilauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. Data was acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kilauea. Melt fractions between 0.096 and 0.117 are calculated for the Kilauea and Pu'u'ō'ō low resistivity anomalies, respectively.

58 GEOSCIENCES↗
DOE OSTI2022

Development of a synthetic phase contrast imaging diagnostic for turbulence studies at Wendelstein 7-X

Abstract We present a synthetic phase contrast imaging (PCI) diagnostic for studying turbulence at the Wendelstein 7-X (W7-X) stellarator. We first describe the implemented instrument response model, which captures diffraction effects, detector noise, and the long-wavelength cutoff due to the phase plate of the PCI system. To verify the instrument response model, we show that it is capable of reproducing the PCI signal generated by the sound wave speaker used for calibration at W7-X. Next, we discuss the calculation of synthetic PCI signals based on the global, nonlinear gyrokinetic codes GENE-3D and EUTERPE, including results from some of the first stellarator simulations of this type with kinetic electrons (KEs) in GENE-3D. While the simulations used in this work lack a neoclassical radial electric field, which is crucial for reproducing experimental PCI signals, they do indicate that the dominant rotation direction and velocities of the turbulent fluctuations can be inferred from the wave number-frequency spectra of the PCI signals, as expected. The synthetic PCI wave number spectra are further shown to be similar to those of the line-integrated fluctuating electron density, with distinct differences between adiabatic and KE simulations, explainable by previously published turbulence models. For example, the wave number spectra of all adiabatic electron simulations analyzed here follow a power law with an exponent close to −5 for sufficiently large wave numbers. This indicates that universal features of electron density turbulence at W7-X may be studied using the PCI system.

EUTERPE↗
DOE OSTI2020

Three-dimensional microstructure-explicit and void-explicit mesoscale simulations of detonation of HMX at millimeter sample size scale

Fully three-dimensional (3D) microstructure-explicit and void-explicit mesoscale simulations of the shock-to-detonation (SDT) process of pressed granular HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) are performed. The overall size scale of the models is up to 3 × 3 × 15 mm3, with ∼30 000 grains and 206 265 voids. The models account for the heterogeneous material microstructure, constituent distribution, constituent morphology, and voids. Loading conditions considered involve piston velocities in the range of 600–1200 m/s or pressures in the range of 4–8 GPa. The focus is on analyzing the SDT process and the effects of microstructure and voids on the run-to-detonation distance (RDD). Companion two-dimensional (2D) simulations are also carried out to assess the differences between 2D and 3D. Statistically equivalent microstructure sample sets (SEMSSs) are generated and used for both 2D and 3D, allowing the prediction of the statistical and probabilistic Pop plots (PPs). The predictions are in general agreement with trends in available experimental data in the literature. It is found that both the microstructure (heterogeneous grain size, morphology, and size distribution) and voids significantly affect the RDD and the PPs. These effects are systematically delineated and quantified via the use of SEMSSs with different combinations of attributes. A recently developed probabilistic formulation for the PPs is used to characterize the results, allowing uncertainties in the relations between the shock pressure and RDD arising from material heterogeneities to be quantified. The probabilistic formulation is further used to quantify the confidence levels in the ranked order of influences of different combinations of microstructure and voids on the PPs.

Miller, Christopher (ORCID:0000000206844381)↗
DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Evaluation of pressure reconstruction techniques for Model Order Reduction in incompressible convective heat transfer

This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Attached Jet Velocity Profiles in Mixing Tanks

Free jets have been studied in detail over much of the last century, but the theory for offset and attached jets remains incomplete. Attached jets differ from free jets in that they lose momentum to nearby surfaces, attenuating their velocities. The velocity profiles of free circular jets are nearly Gaussian, with quantitative mathematical descriptions derived from first principles by Goertler and Tollmien (Rajaratnam, 1976). In contrast, mathematical descriptions of three-dimensional attached jets from circular nozzles remain much less mature. Agelin-Chaab and Tachie (2011) used particle imaging velocimetry of a three-dimensional attached jet to show that the scaled velocity decays with scaled distance from the nozzle with a power law exponent between -1.15 and -1.21, which is larger in magnitude than that of a free jet. However, quantitative analytical expressions for the velocity profiles of attached jets similar to those of free jets remain elusive. This paper addresses this critical gap. Here we evaluate the velocity profiles of three-dimensional offset jets emerging from circular nozzles that that become attached jets. These jets lose momentum due to interactions with nearby surfaces and are important to evaluating flows in mixing vessels and to suspending solids and trapped gases in radioactive waste tanks. Despite the importance of attached jets, prior insight has been purely experimental, limited to overly simplistic analytical models, or restricted to computationally expensive computational fluid dynamics case studies. We compare the expression of Verhoff (1963) to experimental results to find reasonable quantitative agreement. As stated by Agelin-Chaab and Tachie (2011), “detailed velocity measurements of 3D offset jets are rare.” Such remains the case. This study adds to the literature by providing information at two additional Reynolds numbers (1.43 106 and 1.87 106) and evaluating simple but accurate expressions for velocity profiles. These Reynolds numbers are an order of magnitude higher than prior reports. The semi-empirical stream wise velocity profile perpendicular to the surface proposed by Verhoff (1963) is in approximate agreement with these velocity profiles, which is surprising because these attached jets are three-dimensional instead of two-dimensional as evaluated by Verhoff. However, additional work is necessary to fully describe these profiles quantitatively.

jet, attached jet, offset jet↗
DOE OSTI2020

Antiferromagnetic transitions of Dirac fermions in three dimensions

We use determinant quantum Monte Carlo simulations to study the role of electron-electron interactions on three-dimensional (3D) Dirac fermions based on the π-flux model on a cubic lattice. We show that the Hubbard interaction drives the 3D Dirac semimetal to an antiferromagnetic (AF) insulator only above a finite critical interaction strength and the long-range AF order persists up to a finite temperature. We evaluate the critical interaction strength and temperatures using finite-size scaling of the spin structure factor. The critical behaviors are consistent with the (3+1)-dimensional Gross-Neveu universality class for the quantum critical point and 3D Heisenberg universality class for the thermal phase transitions. We further investigate correlation effects in the birefringent Dirac fermion system. It is found that the critical interaction strength U c is decreased by reducing the velocity of the Dirac cone, quantifying the effect of velocity on the critical interaction strength in 3D Dirac fermion systems. Our findings unambiguously uncover correlation effects in 3D Dirac fermions and may be observed using ultracold atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2021

Stellar Mass Black Hole Formation and Multimessenger Signals from Three-dimensional Rotating Core-collapse Supernova Simulations

We present self-consistent 3D core-collapse supernova simulations of a 40 ${M}_{\odot }$ progenitor model using the isotropic diffusion source approximation for neutrino transport and an effective general relativistic potential up to `0.9 s postbounce. We consider three different rotational speeds with initial angular velocities of ${{\rm{\Omega }}}_{0}=0$, 0.5, and 1 rad s -1 and investigate the impact of rotation on shock dynamics, black hole (BH) formation, and gravitational wave (GW) signals. The rapidly rotating model undergoes an early explosion at ~250 ms postbounce and shows signs of the low $T/| W| $ instability. We do not find BH formation in this model within ~460 ms postbounce. In contrast, we find BH formation at 776 ms postbounce and 936 ms postbounce for the nonrotating and slowly rotating models, respectively. The slowly rotating model explodes at ~650 ms postbounce, and the subsequent fallback accretion onto the proto–neutron star (PNS) results in BH formation. In addition, the standing accretion shock instability induces rotation of the PNS in the model that started with a nonrotating progenitor. Assuming conservation of specific angular momentum during BH formation, this corresponds to a BH spin parameter of a = J/M = 0.046. However, if no explosion sets in, all the angular momentum will eventually be accreted by the BH, resulting in a nonspinning BH. The successful explosion of the slowly rotating model drastically slows down the accretion onto the PNS, allowing continued cooling and contraction that results in an extremely high GW frequency (f ~ 3000 Hz) at BH formation, while the nonrotating model generates GW signals similar to our corresponding 2D simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗
DOE OSTI2026

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗
DOE OSTI2023

Magnetohydrodynamics in free surface liquid metal flow relevant to plasma-facing components

Abstract While flowing Liquid Metal (LM) Plasma-Facing Components (PFCs) represent a potentially transformative technology to enable long-pulse operation with high-power exhaust for fusion reactors, Magnetohydrodynamic (MHD) drag in the conducting LM will reduce the flow speed. Experiments have been completed in the linear open-channel LMX-U device [Hvasta et al 2018 Nucl. Fusion 58 01602] for validation of MHD drag calculations with either insulating or conducting walls, with codes similar to those used to design flowing LM PFCs for a Fusion Nuclear Science Facility [Kessel et al 2019 Fusion Sci. Technol . 75 886]. We observe that the average channel flow speed decreased with the use of conducting walls and the strength of the applied transverse magnetic field. The MHD drag from the retarding Lorentz force resulted in an increase of the LM depth in the channel that ‘piled up’ near the inlet, but not the outlet. As reproduced by OpenFOAM and ANSYS CFX calculations, the magnitude and characteristics of the pileup in the flow direction increased with the applied traverse magnetic field by up to 120%, as compared to the case without an applied magnetic field, corresponding to an average velocity reduction of ∼45%. Particle tracking measurements confirmed a predicted shear in the flow speed, with the surface velocity increasing by 300%, despite the 45% drop in the average bulk speed. The MHD effect makes the bulk flow laminarized but keeps surface waves aligned along the magnetic field lines due to the anisotropy of MHD drag. The 3D fringe field and high surface velocity generate ripples around the outlet region. It was also confirmed that the MHD drag strongly depends on the conductivity of the channel walls, magnetic field, and volumetric flow rate, in agreement with the simulations and a developed analytical model. These validated models are now available to begin to determine the conditions under which the ideal LM channel design of a constant flow speed and fluid depth could be attained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Novel Experimental Approach to Understand the Transport of Nanodrugs

Nanoparticle-based drugs offer attractive advantages like targeted delivery to the diseased site and size and shape-controlled properties. Therefore, understanding the particulate flow of the nanodrugs is important for effective delivery, accurate prediction of required dosage, and developing efficient drug delivery platforms for nanodrugs. In this study, the transport of nanodrugs including flow velocity and deposition is investigated using three model metal oxide nanodrugs of different sizes including iron oxide, zinc oxide, and combined Cu-Zn-Fe oxide synthesized via a modified polyol approach. The hydrodynamic size, size, morphology, chemical composition, crystal phase, and surface functional groups of the water-soluble nanodrugs were characterized via dynamic light scattering, transmission electron microscopy, scanning electron microscopy-energy dispersive X-ray, X-ray diffraction, and fourier transform infrared spectroscopy, respectively. Two different biomimetic flow channels with customized surfaces are developed via 3D printing to experimentally monitor the velocity and deposition of the different nanodrugs. A diffusion dominated mechanism of flow is seen in size ranges 92 nm to 110 nm of the nanodrugs, from the experimental velocity and mass loss profiles. The flow velocity analysis also shows that the transport of nanodrugs is controlled by sedimentation processes in the larger size ranges of 110–302 nm. However, the combined overview from experimental mass loss and velocity trends indicates presence of both diffusive and sedimentation forces in the 110–302 nm size ranges. It is also discovered that the nanodrugs with higher positive surface charges are transported faster through the two test channels, which also leads to lower deposition of these nanodrugs on the walls of the flow channels. The results from this study will be valuable in realizing reliable and cost-effective in vitro experimental approaches that can support in vivo methods to predict the flow of new nanodrugs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Astronomy and AstroPhysics↗