Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Material Science Image Analysis using Quant-CT in ImageJ

We introduce a computational analysis workflow to access properties of solid objects using nondestructive imaging techniques that rely on X-ray imaging. The goal is to process and quantify structures from material science sample cross sections. The algorithms can differentiate the porous media (high density material) from the void (background, low density media) using a Boolean classifier, so that we can extract features, such as volume, surface area, granularity spectrum, porosity, among others. Our workflow, Quant-CT, leverages several algorithms from ImageJ, such as statistical region merging and 3D object counter. It also includes schemes for bilateral filtering that use a 3D kernel, for parallel processing of sub-stacks, and for handling over-segmentation using histogram similarities. The Quant-CT supports fast user interaction, providing the ability for the user to train the algorithm via subsamples to feed its core algorithms with automated parameterization. Quant-CT plugin is currently available for testing by personnel at the Advanced Light Source and Earth Sciences Divisions and Energy Frontier Research Center (EFRC), LBNL, as part of their research on porous materials. The goal is to understand the processes in fluid-rock systems for the geologic sequestration of CO2, and to develop technology for the safe storage of CO2 in deep subsurface rock formations. We describe our implementation, and demonstrate our plugin on porous material images. This paper targets end-users, with relevant information for developers to extend its current capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Basin-wide gradients in dolomite d26Mg and 87Sr/86Sr identify the structural centre of the Williston Basin as a source of magnesium for dolomitization of the Middle Devonian Winnipegosis Formation

Dolomitized carbonates in the Williston Basin are common below the Middle Devonian Prairie Evaporite Formation, indicating a plausible role for brine-reflux as the source of Mg for dolomitization. However, the Lower Member Winnipegosis Formation, which directly underlies the evaporite sequence, exhibits dolomite 87Sr/86Sr ratios that are higher than expected for Middle Devonian seawater, instead suggesting that the dolomitizing fluids ascended from the deep basin where present-day formation waters are known to be more radiogenic. This interpretation is supported in this study by Mg isotope data showing a basin-scale gradient of increasing d26Mg values from the deep center (–2.03‰) to the edge (–1.34‰) of the Williston Basin along the lower Member Winnipegosis bed. Dolomite favors light isotopes of Mg when it forms, thus driving dolomitizing solutions, and the produced dolomite, towards heavier isotopic compositions in the direction of fluid flow. However, dolomite is also Mg-deficient when it forms, and its initial Mg isotope composition is susceptible to alteration by fluid-rock interactions occurring later in the depositional history of the deposit. We argue that the Mg isotope gradient in the Lower Member Winnipegosis was created by an ascending flow of isotopically light Mg-bearing fluid that differentially altered an earlier formed brine-reflux dolomite, with more alteration occurring deeper in the basin and less towards to edges. Alteration of pre-existing dolomite by an Mg isotope exchange mechanism requires large supplies of Mg. This fact, and the direction of fluid flow point to the underlying crystalline basement as an important Mg source for dolomitization in the Williston Basin. We speculate that ultramafic rocks supplied the Mg when they were hydrothermally altered during an anomalous heating event that affected the Williston Basin in the Late Devonian / Early Carboniferous. Seismic pumping drove these fluids upwards through fracture networks in the crust, and eventually into the deep center of the Williston Basin along vertically-oriented down-to-the-basement faults, pressuring deep aquifers like the Winnipegosis, which resulted in fluids migrating up the stratigraphic dip of the beds towards the edges of the basin. Plotting and contouring the dolomite d26Mg data on a map of the Williston Basin reveals the pathways taken by the migrating fluids. One prominent pathway extending outwards from the center of the basin towards the northeast bears a strong resemblance to petroleum migration fairways deduced for the underlying Late Ordovician Red River Formation. This study demonstrates that Mg isotopic mapping of large dolomite bodies is a useful tool for deciphering fluid flow histories in sedimentary basins and may assist in locating migrated petroleum and Mississippi Valley Type ore deposits.

Magnesium Isotopes, Strontium Isotopes, Dolomitiza↗
DOE OSTI2023

Material Plasma Exposure eXperiment High Heat Flux Microwave Absorber Design, Manufacture, and Articles Test

The Material Plasma Exposure eXperiment (MPEX) at Oak Ridge National Laboratory is in the final design phase. MPEX will be capable of exposing neutron-irradiated materials to plasmas for the study of plasma-material interaction. This facility will provide information about the complex effects of plasmas on materials and contribute to examining new materials that can withstand high heat fluxes and high ion fluences for future fusion devices. MPEX plasma is heated by 70-GHz or 105-GHz electron Bernstein wave/electron cyclotron heating (ECH), and the high-frequency microwaves are prone to scattering microwave power, which can have detrimental effects, especially on diagnostic components. A large portion of the injected ECH power is expected to be absorbed by plasma, but the remainder requires that microwave absorbers be placed immediately upstream and downstream of the ECH launcher to minimize stray microwaves leaving the ECH region. These microwaves can inadvertently heat components that cannot be shielded or otherwise protected. The microwave absorber design is based on an array of pyramid-shaped ceramic tiles brazed to a water-cooled explosion-bonded heat sink and a stainless steel plate to produce one tile module. Computational fluid dynamics and structural analyses were performed to optimize and validate the design. Multiple test coupons were produced to validate the process for brazing the two different tile materials to the Glidcop AL-15 baseplate. The articles were tested to evaluate the reliability and thermal performance through exposure to an electron beam with a heat flux of up to 1.5 MW/m2. Nondestructive testing was performed before and after testing to identify voids or separations that may have been introduced by the high heat flux. This paper discusses the details of high heat flux microwave absorber design, manufacturing details and associated challenges, and test results, demonstrating the effectiveness of the proposed design.

Hussain, Aftab↗
DOE OSTI2024

Studying Wind Loading on CSP Collectors to Improve Performance and Reliability

Electricity generation through Concentrating Solar-thermal Power (CSP) adds the advantage of thermal energy storage and heat production. CSP collector costs constitute nearly one-third of total plant costs. Managing wind loading, especially dynamic wind loading caused by the turbulent wind flow, is a significant design challenge. Traditional collector designs rely on wind tunnel experiments and numerical simulations, which do not fully capture dynamic effects. Here, we provide an overview of NREL's activities on wind loading on CSP collectors, aiming to improve their reliability and cost-efficiency. We present insights from field measurements at operational parabolic trough and power-tower CSP plants. Over two years, we measured turbulent wind fields and resulting structural loads in a parabolic trough plant, revealing wind, turbulence, and collector field interactions. Similar measurements are ongoing at the Crescent Dunes heliostat field. These campaigns show how atmospheric turbulence, wind direction, and collector orientation affect dynamic wind loading. Upwind collector structures can generate turbulent structures, causing fluctuating loads at downstream collectors, impacting fatigue lifetime and optical efficiency. Preliminary results indicate that interior collectors experience higher turning moments than those at the edges, suggesting higher drive wear. At Crescent Dunes, we use instrumented heliostats at the field's edge and interior to study the translation of turbulent wind to dynamic structural loads. The campaign is still ongoing, and we will present first findings. Another focus of our work is tying the experimental wind loading findings to optical performance of the collectors, using simulation tools. For example, we derive torsional errors of the parabolic troughs - the angular offset from the sun position along a collector row. This error is impacted by wind-induced mirror deformations. Future work will focus on continuing the measurement efforts and combining them with computational fluid dynamic simulations. By providing detailed field measurements and validated models, our efforts aim to enhance understanding of wind loading on CSP collectors, improving their structural integrity and optical performance.

CSP↗
NASA NTRSDate not supplied

Implicit Large-Eddy Simulation of Axisymmetric and Asymmetric Turbulent Shock Wave Boundary Layer Interaction at Mach 2.5

High-speed inlet flows feature turbulent shock wave boundary layer interactions at supersonic Mach numbers. A canonical Mach 2.5 shock wave boundary layer interaction experiment at the NASA Glenn Research Center was designed for axisymmetric and asymmetric interactions depending on the placement of the shock generator. Four different wall-resolved implicit large-eddy simulations of the experiment were carried out at one-tenth of the Reynolds number of the experiment: (1) Two axisymmetric interactions with shock generator radii of 9.2% and 14.7% of the test section diameter and (2) two asymmetric interactions with 14.7% shock generator radius and shock generator centerline offsets of one-third and one-fourth of the test section radius. For the two axisymmetric interactions, various sub-domains upstream of the interaction, within the interaction, and downstream of the interaction were analyzed with the proper orthogonal decomposition. For all three regions, the dominant modes are streamwise coherent structures with Strouhal numbers based on separation length of 0.15 and below. Fourier spectra of the wall-pressure coefficient for the interaction region have peaks at Strouhal numbers based on separation length of 0.15 and 0.04. For the two asymmetric interactions, the streamwise extent of the separation varies considerably in the azimuthal direction. For both shock generator offsets, the interaction is strongest for the azimuthal location that is closest to the shock generator. Compared to the axisymmetric interaction, the streamwise extent of the separated flow region is noticeably reduced and a substantial crossflow is observed. Fourier analyses of the unsteady flow indicate low-frequency content f or the upstream part of the interaction which is more two-dimensional. Away from the upstream region, with increasing sweep angle and crossflow, the low-frequency content is diminished.

Computational Fluid Dynamics↗
DOE OSTI2024

Pressure Spike in The LBNF Absorber Core s Gun Drilled Cooling Channel from an Accident Beam Pulse

he LBNF absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. The cooling water is supplied by the LBNF absorber radioactive water (RAW) cooling system which is designed as per ASME B31.3 normal fluid service. An uninhibited beam accident pulse striking the water channels was identified as a credible accident scenario. In this study, it is assumed that the beam pulse hits the absorber directly without interacting with any of the other upstream beamline components. The beam parameters used for the LBNF beam are 120 GeV, 2.4 MW with a 1.2 s cycle time. The accident pulse lasts for 10 µs. The maximum energy is deposited in the 3rd aluminum core block. For the sake of simplicity, it is assumed that the accident pulse strikes the 1 in. ID water channel directly. The analysis here simulates the pressure rise in the water during and after the beam pulse and its effects on the aluminum piping components that deliver water to the core blocks. The weld strengths as determined by the Load and Resistance Factor Design (LRDF) and the Allowable Strength Design (ASD) are compared to the forces generated in the weld owing to the pressure spike. A transient structural analysis was used to determine the equivalent membrane, peak, and bending stresses and they were compared to allowable limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Mixed Nanosphere Assemblies at a Liquid–Liquid Interface

The in-plane packing of gold (Au), polystyrene (PS), and silica (SiO 2 ) spherical nanoparticle (NP) mixtures at a water–oil interface is investigated in situ by UV–vis reflection spectroscopy. All NPs are functionalized with carboxylic acid such that they strongly interact with amine-functionalized ligands dissolved in an immiscible oil phase at the fluid interface. This interaction markedly increases the binding energy of these nanoparticle surfactants (NPSs). The separation distance between the Au NPSs and Au surface coverage are measured by the maximum plasmonic wavelength (λ max ) and integrated intensities as the assemblies saturate for different concentrations of non-plasmonic (PS/SiO 2 ) NPs. Further, as the PS/SiO 2 content increases, the time to reach intimate Au NP contact also increases, resulting from their hindered mobility. λ max changes within the first few minutes of adsorption due to weak attractive inter-NP forces. Additionally, a sharper peak in the reflection spectrum at NP saturation reveals tighter Au NP packing for assemblies with intermediate non-plasmonic NP content. Grazing incidence small angle X-ray scattering (GISAXS) and scanning electron microscopy (SEM) measurements confirm a decrease in Au NP domain size for mixtures with larger non-plasmonic NP content. The results demonstrate a simple means to probe interfacial phase separation behavior using in situ spectroscopy as interfacial structures densify into jammed, phase-separated NP films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Pressure Spike in The LBNF Absorber Core s Gun Drilled Cooling Channel from an Accident Beam Pulse

The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. The cooling water is supplied by the LBNF Absorber Radioactive Water (RAW) cooling system which is designed as per ASME B31.3 Normal Fluid Service [1]. An uninhibited beam accident pulse striking the water channels was identified as a credible accident scenario. In this study, it is assumed that the beam pulse hits the Absorber directly without interacting with any of the other upstream beamline components. The beam parameters used for the LBNF beam are 120 GeV, 2.4 MW with a 1.2 s cycle time. The accident pulse lasts for 10 μs. The maximum energy is deposited in the 3rd aluminum core block. For the sake of simplicity, it is assumed that the accident pulse strikes the 1 in. ID water channel directly. The analysis here simulates the pressure rise in the water during and after the beam pulse and its effects on the aluminum piping components that deliver water to the core blocks. The weld strengths as determined by the Load and Resistance Factor Design (LRDF) and the Allowable Strength Design (ASD) are compared to the forces generated in the weld owing to the pressure spike. A transient structural analysis was used to determine the equivalent membrane, peak, and bending stresses and they were com-pared to allowable limits.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

The Effect of Suprathermal Protons in the Heliosheath on the Global Structure of the Heliosphere and Heliotail

In the interaction between the solar wind (SW) and the local interstellar medium, various processes create ions with energies up to ~10 keV that are out of thermal equilibrium with the "core" population. Wave–particle interactions tend to isotropize the velocity distributions, but the collisionless nature of the SW precludes thermalization. Suprathermal protons can charge-exchange with interstellar hydrogen, producing energetic neutral atoms that are seen by the Interstellar Boundary EXplorer spacecraft. We have developed a model for the presence of several suprathermal populations in the SW downstream of the heliospheric termination shock. The model uses magnetohydrodynamics to satisfy the first three moments of the total ion distribution, and couples these through charge-exchange to neutral hydrogen, conserving mass, momentum, and energy in the combined system. The proton population is separated into a cool core and three suprathermal populations, and hydrogen atoms may charge-exchange with protons from any of those four populations. The phase-space properties of the pick-up ions are selected based on data and theoretical considerations. In this paper we quantify the impact of suprathermal protons on the global structure of the heliosphere by comparing our new model to a traditional Maxwellian fluid model, and a kappa-distribution model. We find that the differences in momentum and energy transfer rates from the protons onto neutral hydrogen between the models leads to different plasma properties in the heliotail, and also changes the size of the heliosphere. Including the energy-dependent charge-exchange cross section into the collision integrals reduces the magnitude of these differences.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Flowfield analysis of modern helicopter rotors in hover by Navier-Stokes method

The viscous, three-dimensional, flowfields of UH60 and BERP rotors are calculated for lifting hover configurations using a Navier-Stokes computational fluid dynamics method with a view to understand the importance of planform effects on the airloads. In this method, the induced effects of the wake, including the interaction of tip vortices with successive blades, are captured as a part of the overall flowfield solution without prescribing any wake models. Numerical results in the form of surface pressures, hover performance parameters, surface skin friction and tip vortex patterns, and vortex wake trajectory are presented at two thrust conditions for UH60 and BERP rotors. Comparison of results for the UH60 model rotor show good agreement with experiments at moderate thrust conditions. Comparison of results with equivalent rectangular UH60 blade and BERP blade indicates that the BERP blade, with an unconventional planform, gives more thrust at the cost of more power and a reduced figure of merit. The high thrust conditions considered produce severe shock-induced flow separation for UH60 blade, while the BERP blade develops more thrust and minimal separation. The BERP blade produces a tighter tip vortex structure compared with the UH60 blade. These results and the discussion presented bring out the similarities and differences between the two rotors.

Srinivasan, G. R.↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Selected OAST/OSSA space experiment activities in support of Space Station Freedom

The Space Experiments Division at NASA Lewis Research Center is developing technology and science space experiments for the Office of Aeronautics and Space Technology (OAST) and the Office of Space Sciences and Applications (OSSA). Selected precursor experiments and technology development activities supporting the Space Station Freedom (SSF) are presented. The Tank Pressure Control Experiment (TPCE) is an OAST-funded cryogenic fluid dynamics experiment, the objective of which is to determine the effectiveness of jet mixing as a means of equilibrating fluid temperatures and controlling tank pressures, thereby permitting the design of lighter cryogenic tanks. The information from experiments such as this will be utilized in the design and operation of on board cryogenic storage for programs such as SSF. The Thermal Energy Storage Flight Project (TES) is an OAST-funded thermal management experiment involving phase change materials for thermal energy storage. The objective of this project is to develop and fly in-space experiments to characterize void shape and location in phase change materials used in a thermal energy storage configuration representative of an advanced solar dynamic system design. The information from experiments such as this will be utilized in the design of future solar dynamic power systems. The Solar Array Module Plasma Interaction Experiment (SAMPIE) is an OAST-funded experiment to determine the environmental effects of the low earth orbit (LEO) space plasma environment on state-of-the-art solar cell modules biased to high potentials relative to the plasma. Future spacecraft designs and structures will push the operating limits of solar cell arrays and other high voltage systems. SAMPIE will provide key information necessary for optimum module design and construction. The Vibration Isolation Technology (VIT) Advanced Technology Development effort is funded by OSSA to provide technology necessary to maintain a stable microgravity environment for sensitive payloads on board spacecraft. The proof of concept will be demonstrated by laboratory tests and in low-gravity aircraft flights. VIT is expected to be utilized by many SSF microgravity science payloads. The Space Acceleration Measurement System (SAMS) is an OSSA-funded instrument to measure the microgravity acceleration environment for OSSA payloads on the shuttle and SSF.

Delombard, Richard↗
DOE OSTI2023

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

The Present and Future of QCD: QCD Town Meeting White Paper – An Input to the 2023 NSAC Long Range Plan

It is currently understood that there are four fundamental forces in nature: gravitational, electromagnetic, weak and strong forces. The strong force governs the interactions between quarks and gluons, elementary particles whose interactions give rise to the vast majority of visible mass in the universe. The mathematical description of the strong force is provided by the non-Abelian gauge theory Quantum Chromodynamics (QCD). While QCD is an exquisite theory, constructing the nucleons and nuclei from quarks, and furthermore explaining the behavior of quarks and gluons at all energies, remain to be complex and challenging problems. Such challenges, along with the desire to understand all visible matter at the most fundamental level, position the study of QCD as a central thrust of research in nuclear science. Experimental insight into the strong force can be gained using large particle accelerator facilities, which are necessary to probe the very short distance scales over which quarks and gluons interact. The Long Range Plans (LRPs) exercise of 1989 and 1996 led directly to the construction of two world-class facilities: the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) that is focused on studying how the structure of hadrons emerges from QCD (cold QCD research), and the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL) that aims at the discovery and study of a new state of matter, the quark-gluon plasma (QGP), at extremely high temperatures (hot QCD research). These past investments have produced major advances. Nucleons and nuclei are being studied with increasing precision with a unified description of the partonic structure utilizing multi-dimensional imaging. Significant progress has been made, paving the way towards a complete picture of how quarks and gluons give rise to the mass, spin, and momentum of the nucleon. In hot QCD, the QGP is created in the collisions of nuclei at RHIC and the Large Hadron Collider (LHC) and is observed to behave like a fluid with very low specific shear viscosity; the current goals are to understand how the fluid behavior emerges from QCD and to characterize the temperature (and chemical potential) dependence of the properties of the QGP. As this White Paper is written, current experimental programs at CEBAF, RHIC and the LHC continue to provide exciting near term opportunities to capitalize on the investments in experimental equipment and accelerator operations. Most importantly, the QCD community looks forward to the construction of the Electron Ion Collider (EIC) as a major new facility to push forward QCD research in the next decades, with significant focus on exploring the properties of gluons, the mediators of the strong force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Advanced System Thermal Fluids Solver Development for SAM

This work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key is the implementation of a high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. Leveraging the existing capabilities of the SAM code, significant code coverages were established in the finite volume method code. This in turn allows for a suite of test problems with different problem sizes and levels of complexity to be used to quantify the performance improvement of the finite volume method code. As evidently shown in this study, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE for the wide range of selected problems. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. In addition, for a complex reactor model, transient simulation was performed using the finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development. In this work, short-term priority development and testing items were identified, and long-term code adoption and integration plans were made for the eventual deployment of the finite volume method in the SAM code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Direct observation of a magnetic-field-induced Wigner crystal

Wigner predicted that when the Coulomb interactions between electrons become much stronger than their kinetic energy, electrons crystallize into a closely packed lattice. A variety of two-dimensional systems have shown evidence for Wigner crystals (WCs). However, a spontaneously formed classical or quantum WC has never been directly visualized. Neither the identification of the WC symmetry nor direct investigation of its melting has been accomplished. Here we use high-resolution scanning tunnelling microscopy measurements to directly image a magnetic-field-induced electron WC in Bernal-stacked bilayer graphene and examine its structural properties as a function of electron density, magnetic field and temperature. At high fields and the lowest temperature, we observe a triangular lattice electron WC in the lowest Landau level. The WC possesses the expected lattice constant and is robust between filling factor ν ≈ 0.13 and ν ≈ 0.38 except near fillings where it competes with fractional quantum Hall states. Increasing the density or temperature results in the melting of the WC into a liquid phase that is isotropic but has a modulated structure characterized by the Bragg wavevector of the WC. At low magnetic fields, the WC unexpectedly transitions into an anisotropic stripe phase, which has been commonly anticipated to form in higher Landau levels. Furthermore, analysis of individual lattice sites shows signatures that may be related to the quantum zero-point motion of electrons in the WC lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗