SEARCH · Search NASA

Results for “Flow 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2020

From pools to flow: The PROMISE framework for new insights on soil carbon cycling in a changing world

Soils represent the largest terrestrial reservoir of organic carbon, and the balance between soil organic carbon (SOC) formation and loss will drive powerful carbon-climate feedbacks over the coming century. To date, efforts to predict SOC dynamics have rested on pool-based models, which assume classes of SOC with internally homogenous physicochemical properties. However, emerging evidence suggests that soil carbon turnover is not dominantly controlled by the chemistry of carbon inputs, but rather by restrictions on microbial access to organic matter in the spatially heterogeneous soil environment. The dynamic processes that control the physicochemical protection of carbon translate poorly to pool-based SOC models; as a result, we are challenged to mechanistically predict how environmental change will impact movement of carbon between soils and the atmosphere. In this study, we propose a novel conceptual framework to explore controls on belowground carbon cycling: Probabilistic Representation of Organic Matter Interactions within the Soil Environment (PROMISE). In contrast to traditional model frameworks, PROMISE does not attempt to define carbon pools united by common thermodynamic or functional attributes. Rather, the PROMISE concept considers how SOC cycling rates are governed by the stochastic processes that influence the proximity between microbial decomposers and organic matter, with emphasis on their physical location in the soil matrix. We illustrate the applications of this framework with a new biogeochemical simulation model that traces the fate of individual carbon atoms as they interact with their environment, undergoing biochemical transformations and moving through the soil pore space. We also discuss how the PROMISE framework reshapes dialogue around issues related to SOC management in a changing world. We intend the PROMISE framework to spur the development of new hypotheses, analytical tools, and model structures across disciplines that will illuminate mechanistic controls on the flow of carbon between plant, soil, and atmospheric pools.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Active microrheology and dynamic phases for pattern-forming systems with competing interactions

We consider the driven dynamics of a probe particle moving through an assembly of particles with competing long-range repulsive and short-range attractive interactions, which form crystal, stripe, labyrinth, and bubble states as the ratio of attraction to repulsion is varied. We show that the probe particle exhibits a depinning-like threshold from an elastic regime, where the probe particle is trapped by interactions with the other particles, to a plastic flow regime, where the probe particle can break bonds in the surrounding medium. For a fixed particle density, the depinning threshold and sliding velocity of the probe particle vary nonmonotonically as the attraction term is increased. A velocity minimum appears near the crystal to stripe crossover, and there is a significant increase in the depinning threshold in the bubble regime when the probe particle is strongly confined inside the bubbles. For fixed attractive interaction but increasing particle density, the behavior is also nonmonotonic, and there are jumps and drops in the velocity and depinning threshold corresponding to points at which the system transitions between different structures. There are also several distinct flow states that can be characterized by the amount of plastic deformation induced by the probe particle in the surrounding medium. Each flow state generates a different amount of effective drag on the probe particle, and there can be jumps in the velocity-force curve at transitions between the states. In conclusion, we also find that when oriented stripes are present, the probe particle can move along the stripe in an edge transport state that has a finite Hall angle.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Response of Sulfonated Polystyrene Melts to Nonlinear Elongation Flows

Ionizable polymers form dynamic networks with domains controlled by two distinct energy scales, ionic interactions and van der Waals forces; both evolve under elongational flows during their processing into viable materials. A molecular level insight of their nonlinear response, paramount to controlling their structure, is attained by fully atomistic molecular dynamics simulations of a model ionizable polymer, polystyrene sulfonate. As a function of increasing elongational flow rate, the systems display an initial elastic response, followed by an ionic fraction-dependent strain hardening, stress overshoot, and eventually strain-thinning. As the sulfonation fraction increases, the chain elongation becomes more heterogeneous. Finally, flow-driven ionic assembly dynamics that continuously break and reform control the response of the system.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nonlinearity of the post-spinel transition and its expression in slabs and plumes worldwide

Phase transitions in the mantle control its internal dynamics and structure. The post-spinel transition marks the upper–lower mantle boundary, where ringwoodite dissociates into bridgmanite plus ferropericlase, and its Clapeyron slope regulates mantle flow across it. This interaction has previously been assumed to have no lateral spatial variations, based on the assumption of a linear post-spinel boundary in pressure and temperature. Here we present laser-heated diamond anvil cell experiments with synchrotron X-ray diffraction to better constrain this boundary, especially at higher temperatures. Combining our data with results from the literature, and using a global analysis based on machine learning, we find a pronounced nonlinearity in the post-spinel boundary, with its slope ranging from –4 MPa/K at 2100 K, to –2 MPa/K at 1950 K, and to 0 MPa/K at 1600 K. Changes in temperature over time and space can therefore cause the post-spinel transition to have variable effects on mantle convection and the movement of subducting slabs and upwelling plumes.

58 GEOSCIENCES↗

DOE OSTI2022

Temporal evolution of scalar modes in Richtmyer–Meshkov instability of inclined interface using high-speed PIV and PLIF measurements at 60 kHz

Abstract The current work presents simultaneous, high-speed measurements at 60,000 fields per second of velocity and mole fraction using particle image velocimetry (PIV) and planar laser induced acetone-fluorescence in a Richtmyer–Meshkov instability of an inclined interface (Atwood number, At = 0.22). Specifically, around 2 ms of temporal evolution of the vortex structures and their associated scalar modes immediately following the interface-reshock interaction is presented. Two initial interface conditions are discussed—(a) a sharp, inclined ‘single mode’ interface and (b) a ‘multi-mode’ interface where small perturbations are imposed on the single mode case. A 2D wavelet decomposition of the scalar flow field shows a highly intermittent distribution of small-scale variance throughout the interface even at late times. These are correlated strongly with the vortex structures and local turbulence intensity, where each small-scale scalar mode is sandwiched between two co-rotating vortex structures. This indicates that the interstitial regions between the vortices are significant hotspots of entrainment, which is then dispersed by the induced, counter-flow velocity fields. The multimode case demonstrates similar organization at large scales, while the scalar field is much more homogeneous at smaller scales. These observations highlight the importance of capturing the early time vortex evolution to accurately estimate any late time intermittency, especially where deposition of intense vorticity on sharp interfaces is present.

Engineering↗

DOE OSTI2023

Droplet-Gas Interactions in Nanoelectrospray Multiphase Flow

Electrospray (ES) is production of a charged droplet plume that transits from an electrically biased supply capillary to a counter electrode under the influence of an electric field. Drag interactions with the surrounding gas decelerate the droplets while transferring the momentum from droplets to accelerate the surrounding gas, resulting in an induced gas flow. This work seeks to characterize the structure of the resulting gas jet and identify the key mechanisms defining the flow structure. The phenomenon of gas jets produced by momentum transfer from nano-electrospray (nES) plumes is explored with schlieren visualization, thermal anemometry, and numerical simulations. Schlieren visualization experiments provide information on the flow structure in support of simulation predictions, and the hot thermistor anemometry measurements of gas velocities outside the spray demonstrate quantitatively validated simulation results. The study reveals the formation of a moderately high velocity coaxial gas jet within the nES plume and provides insight into the overall flow structure of the induced flow. The multiphase electrohydrodynamic simulations enable numerical experimentation to explore the fundamental physics of coupled droplet-gas transport and the resulting flow structure. The simulations, confirmed by the experiments, reveal gas jetting induced by nES with a narrow (sub hundred micrometers in diameter) core originating from the nES liquid-jet breakup region and a surrounding larger-in-extent zone (several hundred micrometers in diameter) of lower velocity gas flow within the electrospray plume. This behavior is due to nES ejecting a stream of droplets from the tip of a narrow liquid cone-jet, where the combined effect of many small droplets transferring momentum to a confined region of gas yields a narrow but high-velocity gas stream. Furthermore, the micro-jets produced by nES have practical utility for mass spectrometry, 3D printing and fabrication, and thermal management.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2020

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications. This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ~2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2 H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics. In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T g with increased confinement that is strongly influenced by the pore-liquid interface characteristics. This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. Finally, while the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗

DOE OSTI2022

A Macroalgal Cultivation Modeling System (MACMODS): Evaluating the Role of Physical-Biological Coupling on Nutrients and Farm Yield

Offshore aquaculture has the potential to expand the macroalgal industry. However, moving into deeper waters requires suspended structures that will present novel farm-environment interactions. Here, we present a computational modeling framework, the Macroalgal Cultivation Modeling System (MACMODS), to explore within-farm modifications to light, seawater flow, and nutrient fields across time and space scales relevant to macroalgae. A regional ocean model informs the site-specific setting, the Santa Barbara Channel in the Southern California Bight. A fine-scale hydrodynamic model predicts modified flows and turbulent mixing within the farm. A spatially resolved macroalgal growth model, parameterized for giant kelp, Macrocystis pyrifera , predicts kelp biomass. Key findings from model integration are that regional ocean conditions set overall farm performance, while fine-scale within-farm circulation and nutrient delivery are important to resolve variation in within-farm macroalgal performance. Therefore, we conclude that models resolving within-farm dynamics can provide benefit to farmers with insight on how farm design and regional ocean conditions interact to influence overall yield. Here, the presence of repeating longlines aligned with the mean current generate flow diversions around the farm as well as attached Langmuir circulations and increased turbulence intensity. These flow-induced phenomena lead to less biomass in the interior portion of the farm relative to the edges. We also find that there is an effluent “footprint” that extends as much as 20 km beyond the farm. In this regard, MACMODS can be used to not only evaluate farm design and cultivation practices that maximize yield but also explore interactions between the farm and ecosystem in order to minimize impacts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Evaluation of flow-induced plate deflection for University of Missouri research reactor low-enriched uranium fuel element

The University of Missouri Research Reactor (MURR), located on the campus of the University of Missouri in Columbia, Missouri, is one of the six United States (U.S.) High Performance Research Reactors (USHPRR), including one critical facility, that are actively collaborating with the U.S. Department of Energy (DOE) National Nuclear Security Administration (NNSA) Office of Material Management and Minimization (M3) Reactor Conversion Program to convert from highly enriched uranium (HEU, ≥20 wt% U-235) fuel to low-enriched uranium (LEU, <20 wt% U-235) fuel. A new type of very high-density LEU fuel based on a monolithic alloy of uranium and 10 wt% molybdenum (U-10Mo) is expected to allow conversion of some USHPRR, including MURR. In the design of its fuel elements, MURR is using thin parallel curved fuel plates separated by coolant channels. In this work, fluid-structure interaction (FSI) analysis of the MURR LEU fuel element is performed at the element level (as compared to the plate level analysis), which models all components of the LEU fuel element, including fuel plates and the supporting structures. Therefore, the effect of supporting structures on the flow distribution within the element and the fuel plate deflection are evaluated. In addition to the element nominal flow rate and dimensions, the tolerances in the geometry of the coolant channel and plate thickness, the effect of a comb on plate deflection, and the uncertainty of the flow rate per element are evaluated. For the LEU fuel plates, which are thinner than the current HEU plates, the predicted plate deflection is found to be small compared to the fabrication and assembly tolerances. Thus, the FSI-induced deflections are not expected to noticeably reduce the coolant flow rate or predicted safety margins in the limiting channels for the MURR LEU fuel element. In addition to the simulation work, a hydraulic performance test of the MURR LEU fuel element is currently being planned to support conversion to the use of LEU fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

X-ray imaging of liquid-liquid Mg cavitation

The coupling of dynamic loading platforms with synchrotron light sources have enabled experiments examining void formation in optically opaque materials. Much like spallation in a solid, cavitation occurs when rarefaction waves collide within a liquid. The resulting deformation leads to a complex flow field with numerous voids. It has been suggested that compaction of such a structure may be a significant source of ejecta (i.e. shallow bubble collapse). Experiments imaging supported and unsupported shocks interacting with a free surface in shock melted Mg with synchrotron radiation are presented. Here, the results are used to examine the structure of a cavitation plane that forms in the unsupported case and provide an estimate of the porosity within this region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Ratios of collective flow observables in high-energy isobar collisions are insensitive to final-state interactions

The ratios of bulk observables, such as harmonic flow v 2 and v 3 , between high-energy 96 Ru + 96 Ru and 96 Zr + 96 Zr collisions were recently argued to be a clean probe of the nuclear structure differences between 96 Ru and 96 Zr . Using a transport model simulation of isobar collisions, we quantify this claim from the dependence of the ratios v 2,Ru /v 2, Zr and v 3,Ru /v 3,Zr on various final-state effects, such as the shear viscosity, hadronization, and hadronic cascade. Although the v 2 and v 3 change by more than 50% when varying the final-state effects, the ratios are unchanged. In addition, these ratios are independent of the transverse momentum p T and hadron species, despite of up to a factor of 2 change in v n . Further, the ratio of mean transverse momentum $\langle$p T $\rangle$ is found to be controlled by the nuclear skin and nuclear radius but is only slightly impacted by the final-state effects. Therefore, these isobar ratios serve as a clean probe of the initial condition of the quark-gluon plasma, which, in turn, is controlled by the collective structure of the colliding nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A code-to-code benchmark for magneto-convection in a horizontal duct

Liquid metals and magnetic fields are used in many technical applications such as metallurgy, crystal growth and nuclear fusion reactors. When an electrically conducting fluid moves in a magnetic environment, electric currents and electromagnetic forces are generated that affect velocity and pressure losses in the flow. These magnetohydrodynamic (MHD) interactions have to be investigated to optimize the engineering processes. The characteristics of MHD flows depend on the geometrical configuration, the strength of the applied magnetic field, the electrical properties of fluid and structural materials and the thermal conditions. In the so-called blankets for fusion reactors, where liquid metals are used to breed the plasma fuel component tritium and to extract the generated heat, magneto-convective flows play a crucial role in determining heat and mass transfer. Therefore, the availability of numerical codes to simulate this type of flow is mandatory and their validation is a necessary step to guarantee the reliability of the results. For that reason, a benchmark problem has been defined to simulate liquid metal flows in a horizontal rectangular duct heated from below and exposed to a non-uniform magnetic field. Results obtained by five research groups using different codes are compared.

benchmark↗