Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vertical axis”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Tomo2Mesh: Fast Reconstruction and Visualization of Tomography Data in Mesh Format

Tomo2Mesh is an open-source project targeted towards real-time reconstruction, segmentation, and visualization of computed tomography (CT) data in mesh format. The CT reconstruction scheme is based on filtered back-projection of voxel subsets. The segmentation scheme uses a 3D convolutional neural network. To allow for fast, real-time reconstruction, voxel subsets are first identified by coarse reconstruction. The detail in specific regions of interest is improved through full reconstruction of voxel subsets in that region. Data structures are implemented to store and process voxel subsets. Finally, voxel data is labeled using connected components to detect disconnected regions such as voids whose morphological attributes (e.g., Feret diameter, principal axis orientation, local number density, size, etc.) can be measured also in real-time. Finally, a fast marching cubes implementation processes labeled voxel data into a triangular face mesh (vertices and faces) in .ply format for visualization in Paraview or other mesh visualization tools. The code provides a simple programming interface for detecting, classifying, and visualizing regions of interest based on morphology. For example, detected voids can be classified as round pores or extended cracks. Highly porous neighborhoods can be identified based on local number density. The mesh texture (or color) is assigned based these morphological attributes to allow smart visualization scenarios in real-time (e.g., show only long cracks). At the time of first release (July 2022), extraction of face mesh for visualization for raw CT data from a 2 megapixel camera would take between 1-5 minutes for most scenarios.

TEKAWADE, ANIKET↗
DOE OSTI2025

Synergistic effects of Al, Ga, and In doping on ZnO nanorod arrays grown via citrate-assisted hydrothermal technique for highly efficient and fast scintillator screens

To be used as efficient alpha particle scintillator in the fields of nuclear security, nuclear nonproliferation and high-energy physics, scintillator screens must have high light output and fast decay properties. While there has been a great deal of progress in scintillation efficiency, achieving fast decay time properties are still a challenge. In this work, the near band edge (NBE) UV luminescence and alpha particle induced scintillation properties of vertically aligned densely packed ZnO nanorods (NRs) doped with Al, Ga, and In have been thoroughly investigated. The high crystalline hexagonal wurtzite structure with a strong orientation through the c -axis plane (002) and aspect ratios in the range 13–22 have been observed for all ZnO NRs. Electron paramagnetic resonance (EPR) analysis exhibited paramagnetic signals at g ≈ 1.96 for all ZnO NRs. A cost effective green hydrothermal synthesis technique was employed to grow well-aligned NRs. Using citrate as an additive acting as a strong reducing agent in the solution during the crystal growth, defects on the surface are significantly suppressed, thereby enhancing the NBE UV emission. Significantly higher NBE UV emission was observed from the top surface of ZnO NRs in cathodoluminescence (CL) microscopy. Results show that citrate assisted donor doping of ZnO NRs not only reduces the defect emission and NBE self-absorption, but also induces fast decay time (~ 600–700 ps), which makes ZnO NRs a good candidate for fast alpha particle scintillator screens used in associated particle imaging for time and direction tagging of individual neutrons generated in D–T and D–D neutron generators.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-mode THz quantum-cascade VECSELs based on disordered metasurfaces

Quantum-cascade vertical-external-cavity surface-emitting-lasers (VECSELs) based on disordered amplifying metasurfaces are demonstrated and explored as potential broadband, multi-mode THz sources. The disorder is introduced along one spatial axis of the metasurface by pseudo-randomly varying the width of its resonant ridge antennas. Compared to a quantum-cascade (QC) VECSEL based on a uniform metasurface, the disordered structure supports much more localized transverse modes with reduced spatial overlap within the QC gain material. This localization is hypothesized to facilitate the spatial hole burning of the gain material and, therefore, enable multi-mode lasing, particularly for short cavities on the order of a few wavelengths. Several devices have been fabricated and shown to differ from uniform QC-VECSELs in a few key ways, possessing highly nonlinear light-current characteristics, angle-dependent emission spectra and broadband multi-mode lasing. At most, 17 modes are simultaneously observed, spanning 680 GHz. The number of VECSEL modes is shown to have an inverse relationship with the cavity length, which is attributed to increased diffractive losses in the open plano-plano cavity. As the cavity length is tuned, the device emits over a quasi-continuous band from 3.15 to 3.97 THz or up to 4.24 THz with a longitudinal mode hop.

Morag, Eilam [Univ. of California, Los Angeles, CA↗
DOE OSTI2025

Hexagonal Boron Nitride: Physical Properties, Hydride Vapor‐Phase Epitaxy Growth of Large‐Diameter Quasi‐Bulk Wafers and Applications

Hexagonal boron nitride ( h ‐BN), with its ultrawide bandgap and 2D structure, holds an immense promise for advanced semiconductor applications. Scaling bulk crystals to large‐diameter wafers, crucial for complex device fabrication, remains a challenge with high temperature, high pressure, and metal flux solution methods. To address this, recent efforts have focused on hydride vapor‐phase epitaxy (HVPE) for producing large diameter thick h ‐BN quasi‐bulk wafers, reaching hundreds of micrometers. These HVPE‐grown quasi‐bulk crystals exhibit excellent c ‐axis long‐range order. Notably, the in‐plane mobility‐lifetime products for both electrons and holes surpass 10 −4 cm 2 V −1 , which are two orders of magnitude greater than the out‐of‐plane (vertical) values, highlighting the potential for high‐performance devices leveraging superior lateral transport. Lateral detectors fabricated from 100 μm thick B‐10 enriched h‐ BN wafers have achieved a record 60% thermal neutron detection efficiency. Based on its physical properties, h‐ BN appears to be an outstanding material of choice for light‐triggered electronic power switches capable of supporting high ‐ voltage and high ‐ power operations. These recent advancements in large‐diameter h ‐BN quasi‐bulk crystal growth, enabled by HVPE, pave the way for applications spanning deep UV photonics, high‐power electronics, high‐efficiency neutron detection, and quantum information technologies, establishing h ‐BN as both a versatile active semiconductor and an ideal substrate.

Jiang, Hongxing [Department of Electrical and Comp↗
DOE OSTI2024

La 4 Co 4 X ( X = Pb , Bi , Sb ) : A demonstration of antagonistic pairs as a route to quasi-low-dimensional ternary compounds

We outline how pairs of strongly immiscible elements, referred to here as antagonistic pairs, can be used to synthesize ternary compounds with low or quasi-reduced-dimensional motifs intrinsically built into their crystal structures. By identifying third elements that are mutually compatible with a given antagonistic pair, ternary compounds can be formed in which the third element segregates the immiscible atoms into spatially separated substructures. Quasi-low-dimensional structural units, such as sheets, chains, or clusters are a natural consequence of the immiscible atoms seeking to avoid close contact in the solid state. Further, as proof of principle, we present the discovery, crystal growth, and basic physical properties of La 4 ⁢Co 4 ⁢$\mathrm{X}$ (X = Pb, Bi, Sb), a family of intermetallic compounds based on the antagonistic pairs Co-Pb and Co-Bi. La 4 ⁢Co 4 ⁢$\mathrm{X}$ adopts an orthorhombic crystal structure (space group Pbam) containing quasi-two-dimensional Co slabs and La-X polyhedra that stack in an alternating manner along the α axis. Consistent with our proposal, the La atoms separate the Co and X substructures, ensuring there are no direct contacts between the members of the immiscible (antagonistic) pair. Within the Co slabs, the atoms occupy the vertices of corner sharing tetrahedra and triangles, and this bonding motif produces narrow electronic bands near the Fermi level that favor magnetism. The Co is moment bearing in each La 4 ⁢Co 4 $\mathrm{X}$ compound studied, and we show that whereas La 4 ⁢Co 4 ⁢Pb behaves as a three-dimensional antiferromagnet with T N =220K, La 4 ⁢Co 4 ⁢Bi and La 4⁢ Co 4 ⁢Sb have behavior consistent with low-dimensional magnetic coupling and ordering, with T N =153K and 143 K, respectively. In addition to the Pb-, Bi-, and Sb-based La 4 ⁢Co 4 ⁢$\mathrm{X}$ compounds, we also were likely able to produce an analogous La 4 ⁢Co 4 ⁢Sn in polycrystalline form, although we were unable to isolate single crystals. We anticipate that identifying and using mutually compatible third elements together with an antagonistic pair represents a generalizable design principle for discovering new materials and new structure types containing low-dimensional substructures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

X-ray wavefront sensor development at the Advanced Light Source

At Lawrence Berkeley National Laboratory's Advanced Light Source, we are developing x-ray wavefront sensors to support the creation and operation of beamlines with diffraction-limited quality. Our new approach to rapid, intermittent wavefront sensing operates in reflection at glancing incidence angles and is compatible with the high-power densities of modern beamlines. For soft x-ray applications especially, the wavefront sensor can operate upstream of the exit slit in a vertically dispersed beam. This single-shot technique supports lateral shearing interferometry and Hartmann wavefront sensing; it can be adapted to speckle-based techniques as well. The reflected beam is directed to an off-axis YAG crystal that produces scintillated visible light. A small mirror reflects the light to a microscope and camera, and the measured wavefront shape information can be used as feedback to adaptive x-ray mirror elements. A compact array of gratings enables measurement across a broad range of photon energies or wavefront curvatures. We describe recent demonstrations at soft x-ray and hard x-ray wavelengths measuring an adaptive x-ray mirror, and a toroidal focusing mirror.

Goldberg, Kenneth A↗
DOE OSTI2023

Capturing plume behavior in complex terrain: an overview of the Nevada National Security Site Meteorological Experiment (METEX21)

METEX21 was an atmospheric tracer release experiment executed at the Department of Energy’s Nevada National Security Site (NNSS) in the southwestern U.S to study terrain-induced wind and thermodynamic conditions that influence local-scale (<5-km) plume transport under varying atmospheric forcing conditions. Meteorological observations were collected using 10-m tall meteorological towers, 2-m tall tripods with 3-d sonic anemometers, a 3-m tall eddy covariance flux tower, Doppler profiling lidars, Doppler scanning lidars, weather-balloon launched radiosondes, and a tethered balloon equipped with wind, temperature, and aerosol sensors at heights up to 800 m. A smoke tracer was released along three transects in the horizontal and vertical directions and observed with video cameras, aerosol sensors and lidars (via aerosol backscatter). The observations showed evidence of large-scale/synoptic transience as well as local-scale upslope and downslope flows, along-axis valley flows, recirculation eddies on leeward slopes, and periods of strong shear and veer aloft. The release days were classified as either synoptically-driven or locally-driven, and a single case day is presented in detail for each. Synoptically-forced days show relatively narrow smoke plumes traveling down the valley from north to south (with the predominant wind direction), with little deviation in transport direction regardless of the elevation or ground locations of the smoke releases, except near the presence of leeside recirculation eddies. Locally-forced days exhibit a wider range of plume behavior due to the combination of thermally-induced valley and slope flows, which are often flowing in different cardinal directions, and wind shear found aloft at higher altitudes and elevations. We saw evidence of smoke lofting on top of the mesas due to strong upslope flows on these days. A major finding of this experiment was the effectiveness of scanning lidars to measure 2-dimensional plume transport out to a 2–3 km distance; much farther than could be visibly observed. METEX21 was the first of three planned tracer experiments at NNSS, and future experiments will incorporate multiple tracers to improve individual plume identification so that finer resolution flow details can be attained from these measurements, as well as deploy a larger suite of meteorological instrumentation, including more temperature profiling data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Determining tolerance requirements for spray-duct alignment in ducted fuel injection

Several ducted fuel injection (DFI) studies have highlighted the importance of accuracy in aligning the duct axis with that of its corresponding spray for optimal effectiveness, as misalignment adversely impacts the method’s performance. The need for accurate alignment could lead to added manufacturing complexity via tighter tolerances. This study systematically explores cases of horizontal, vertical, and rotational misalignment, analyzing their respective effects on DFI performance. Vertical and horizontal misalignments at the duct inlet plane were varied at magnitudes of 6.25%, 12.5%, and 25.0% of the duct diameter, corresponding to 0.125, 0.25, and 0.5 mm, respectively. Rotational misalignments were set at 1°, 2°, and 4°, corresponding to 3.65%, 7.30%, and 14.6%, respectively, of the duct diameter at its inlet plane. The investigation yields spray-duct alignment tolerance limits and highlights the influence of misalignment direction on emissions due to the interactions with swirl and squish inside the combustion chamber. The results indicate that the tolerance limits for the alignment are within 4° and 0.5 mm relative to the geometrically aligned position. If the misalignment exceeds 4°of rotation or 0.5 mm in the horizontal direction, the beneficial effects on soot reduction using this method are no longer observed. The findings contribute to an understanding that can be used to optimize DFI for cleaner and more efficient combustion in compression-ignition engines.

Şener, Ramazan↗
DOE OSTI2025

Utah FORGE: Well 16B(78)-32 Drill Core Fracture Analysis Images and Data

This dataset contains drilling core data from well 16B(78)-32, including PDF documents with flattened core images annotated by feature type and core interval, as well as spreadsheets detailing feature morphologies by depth, planar feature measurements, and planar feature orientations rotated to in situ conditions. Core was recovered from three intervals, one per stimulation stage, in the crystalline rocks affected by the stimulation of well 16A(78)-32. Seven core runs were conducted, yielding 135.8 feet of recovered core. Features in the core were categorized into planar fractures, semi-planar fractures, unbroken mineralized fractures, rough fractures, curviplanar fractures, concave-convex surfaces, and planar compositional features such as mylonite or dike-like structures. Planar features were measured while the core was positioned horizontally, with the core axis aligned to a downhole azimuth of 42 degrees. Planar core measurements from stimulations 2 and 3 that could be confidently correlated with FMI data were rotated to in situ orientations. This was done by rotating the planes along vertical and horizontal axes to match the azimuth and inclination data recorded in the directional survey of well 16B(78)-32, as well as applying an axial rotation to resemble the fracture orientations observed in the FMI log at corresponding depths. Coherent sets of planar fracture measurements were made by aligning the core within each 3-foot section of the dissected core barrel, and between adjacent 3-foot sections within a core run by matching rock fabrics, saw cuts and/or tool marks. Where coherent fracture measurements could not be made within a core run, data sets are denoted by a subscript (i.e. 2-Ta and 2-Tb both come from tangent core run number 2).

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Mass and momentum transport in the Tilted Rocket Rig experiment

Here, high-fidelity large eddy simulations of the inclined Rayleigh-Taylor instability in the Tilted Rocket Rig experimental configuration [Smeeton and Youngs, AWE Report No. O 35/87 (1987)] are performed using a tenth-order compact finite difference code. These simulations are analyzed for spatial distributions of turbulent kinetic energy, turbulent mass flux velocity, species mass fraction flux, species mass fraction variance, and Favre-averaged Reynolds stresses at two time instances, t = 37 ms and t = 52 ms. Additionally, the vertical distribution of the components of the unclosed budget equations over the center of the domain for these quantities are examined. The dominant terms of these budget equations are further decomposed to examine the principal contributions to these terms in each axis. Notably, the principal contribution to the horizontal turbulent mass flux velocity budget is found to be from a term which is commonly neglected in many Reynolds-averaged Navier-Stokes models.

42 ENGINEERING↗
DOE OSTI2021

Regolith, rock and fluid distributions at the Upper Colorado River Basin via a multicomponent seismic imaging approach

In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.

58 GEOSCIENCES↗
DOE OSTI2024

From Feast to Famine: A Systematic Study of Accretion onto Oblique Pulsars with 3D GRMHD Simulations

Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Spot size measurement of a deuterium–tritium dense plasma focus using neutron radiography

Neutron radiography is a technique uniquely suited to applications in nuclear diagnostics, non-destructive testing, and subcritical experiments. The spatial resolution of neutron radiographs is degraded by optical blur in the imaging system and the neutron source size, where the ideal source is point-like to optimize the point-spread function. A potential neutron source for radiography is the dense plasma focus (DPF), a coaxial Z-pinch that produces thermonuclear and beam-target neutrons. To assess if the source size is suitable for radiography, a neutron imaging system was used to measure the source size of the 4 MA Sodium DPF at the Nevada National Security Site operating with deuterium–tritium gas-fill. The source size was measured using the edge-spread function of tungsten objects, each having a rolled (convex) edge. The spot size was found to be 7–12 mm full-width at half-max (FWHM) assuming a Gaussian source, though comparison is presented for Lorentzian and Bennett distributions. The average FWHM was found to be 8.6 ± 1.2 mm vertically and 10.8 ± 1.2 mm horizontally with respect to the image plane, averaging over varied edges and alignments. The results were sensitive to source alignment and edge metrology, which introduced notable uncertainties. These results are consistent with separate experimental measurements as well as magnetohydrodynamics simulations of this DPF, which suggest that neutron production can originate from pinches ∼5–7 mm off-axis. These results suggest that the DPF should be used for radiography at low magnification (M < 1) where spot size does not dominate spatial blur.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Strain Recovery and Defect Characterization in Mg‐Implanted Homoepitaxial GaN on High‐Quality GaN Substrates

The evolution of defects due to high‐pressure annealing of magnesium ion‐implanted epitaxial GaN grown on high‐quality GaN substrates is investigated. Changes in the implant‐induced strain are quantified as a function of annealing temperature and time. After annealing at 1300 °C for 10 min, the implant‐induced strain is fully relieved and accompanied by the presence of extended defects such as basal plane stacking faults and prismatic loops. Approximately one‐third of the original implant‐induced strain remains after annealing at 700 °C, and 5% of the original strain remains at 1000 °C for 100 min. In all cases, nearly all of the recovered strain occurs within first few minutes of annealing. A prominent increase in the asymmetric (104) triple axis X‐ray rocking curve full width at 0.01 maximum (FW0.01M) is observed after annealing at 1300 °C for 10 min. After annealing at 1300 °C for 100 min, a subsequent decrease in FW0.01M is correlated with a reduction of the extended defect density from 4 × 10 8 to 3 × 10 7 cm −2 , determined through transmission electron microscope (TEM) measurements. Further reduction in the density of the extended defects by optimizing annealing temperature and time is expected to improve the performance of GaN‐based vertical power devices.

Wang, Yekan↗
DOE OSTI2020

Room temperature multiferroicity and magnetodielectric coupling in 0–3 composite thin films

Magnetoelectric (ME) composite thin films are promising candidates for novel applications in future nanoelectronics, spintronics, memory and other multifunctional devices as they exhibit much higher ME coupling and transition temperatures (Tc) than well-known single phase multiferroics discovered to date. Among the three types of multiferroic composite nanostructures, (2-2) layered and (1-3) vertically aligned composite nanostructures exhibit comparatively smaller ME coupling due to different shortcomings which restrict their use in many applications. Here we study the morphological, piezoresponse force microscopic (PFM), ferroelectric, magnetic, and magnetodielectric properties of 0-3 (magnetic nanoparticles (0) homogeneously distributed in ferroelectric matrices (3)) multiferroic composite thin films. The Pb(Fe 0.5 Nb 0.5 )O 3 (PFN) - Ni 0.65 Zn 0.35 Fe 2 O 4 (NZFO) particulate composite films were synthesized by pulsed laser deposition (PLD). These particulate composite thin films are completely c-axis oriented with very low surface roughness. We observed magnetic and ferroelectric T c above room temperature (RT) for all composite thin films. The PFN-NZFO 0-3 composites exhibit large polarization, high saturated magnetization with low coercive field, low dielectric loss along with magnetodielectric coupling at RT. These nanocomposites might be utilized in next generation nano/microelectronics and spintronic devices.

36 MATERIALS SCIENCE↗