Search NASA⌕ Search

SEARCH · Search NASA

Results for “shear stress”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5

Orientation-selected micro-pillar compression of additively manufactured 316L stainless steels: Comparison of as-manufactured, annealed, and proton-irradiated variants

In this work, irradiation response and deformation mechanisms of additively manufactured (AM) 316L stainless steel were studied by atomic scale characterization and micro-pillar compression. The AM 316L stainless steels were fabricated by direct energy deposition, a laser-based additive manufacturing process. Irradiation with 2 MeV protons at 360°C was performed to create ~1.8 displacements-per-atom (dpa) damage in AM 316L. Deformation behaviors of the as-manufactured, annealed, and proton-irradiated variants were studied, focusing on the effects of manufacturing-induced pores, residual stress, and irradiation-introduced defects (dislocation loops and voids). Micro-pillars were prepared from grains of pre-selected orientation, avoiding contributions of grain boundaries and allowing determination of resolved shear stress on {111} glide planes. Transmission electron microscopy was used to characterize the pre- and post-deformation microstructure. It was found that in the as-manufactured alloy variant, moving dislocations were the major deformation carrier, with noticeable blocking by fabrication-induced pores, In the annealed variant, hardness was reduced, and deformation was also accomplished by dislocation gliding. In the proton-irradiated variant, significant twinning was observed. Comparing measured resolved shear stress and predicted critical stress for dislocation dissociation, we conclude that irradiation hardening became high enough to activate twinning. Therefore, the deformation mechanism changes from dislocation gliding to twinning. The study is important for both processing optimization and performance evaluation of AM alloys for reactor applications.

36 MATERIALS SCIENCE↗

Experimental investigation of internal two-phase flow structures and dynamics of quasi-stable sheet cavitation by fast synchrotron x-ray imaging

The quasi-stable sheet cavitation produced in a small Venturi channel is investigated using a fast synchrotron X-ray imaging technique aided with conventional high speed photography. The use of X-rays instead of visible light solves cavitation opacity related issues, and X-ray phase contrast-based edge enhancement enables high-definition visualization of the internal two-phase morphology. The simultaneous acquisition of time-resolved velocity and void fraction fields through post-processing the recorded X-ray images reveals, for the first time, the complex diphasic flow structures inside the sheet cavity, which is essentially divided into 6 characteristic parts. Distinct from the current mainstream view, the globallysteady sheet cavitation is found to be characterized by a weak but constantly-existing re-entrant flow that can penetrate the entire cavity. The turbulent velocity fluctuations inside the sheet cavity are also investigated. The turbulence level in the reverse flow region is observed to be as low as in the outer main flow demonstrating the relatively steady status of the reentrant flow. Unlike the streamwise and cross-stream fluctuations, the shear stress appears to be weakly correlated with the velocity gradient. The collapse of vapor phase and the vaporization at the upstream cavity interface are found to be the primary causes of shear stress intensification.

42 ENGINEERING↗

Unexpected Observation of Disorder and Multiple Phase-Transition Pathways in Shock-Compressed Zr

The response of materials under dynamic compression involves a complex interplay of various deformation mechanisms aimed at relieving shear stresses, yielding a remarkable diversity in material behavior. Here, in this Letter, we utilize femtosecond x-ray diffraction coupled with nanosecond laser compression to reveal an intricate competition between multiple shear-relieving mechanisms within an elemental metal. Our observations in shocked-compressed single-crystal Zr indicate a disorder-mediated shear relaxation at lower pressures. Above the phase-transition pressure, we observe the increasing contribution of structural phase transition in relieving shear stress. We detect not one but three concurrent pathways during the transition from the hcp to a hex-3 structure. These complex dynamics are partially corroborated through multimillion-atom molecular dynamics simulations employing a machine-learned interatomic potential. Our observation of multiple concurrent pathways and disorder during shock compression underscore the far greater intricacies in the dynamic response of metals than previously assumed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗

Acute Stent-Induced Endothelial Denudation: Biomechanical Predictors of Vascular Injury

Recent concern for local drug delivery and withdrawal of the first Food and Drug Administration-approved bioresorbable scaffold emphasizes the need to optimize the relationships between stent design and drug release with imposed arterial injury and observed pharmacodynamics. In this study, we examine the hypothesis that vascular injury is predictable from stent design and that the expanding force of stent deployment results in increased circumferential stress in the arterial tissue, which may explain acute injury poststent deployment. Using both numerical simulations and ex vivo experiments on three different stent designs (slotted tube, corrugated ring, and delta wing), arterial injury due to device deployment was examined. Furthermore, using numerical simulations, the consequence of changing stent strut radial thickness on arterial wall shear stress and arterial circumferential stress distributions was examined. Regions with predicted arterial circumferential stress exceeding a threshold of 49.5 kPa compared favorably with observed ex vivo endothelial denudation for the three considered stent designs. In addition, increasing strut thickness was predicted to result in more areas of denudation and larger areas exposed to low wall shear stress. We conclude that the acute arterial injury, observed immediately following stent expansion, is caused by high circumferential hoop stresses in the interstrut region, and denuded area profiles are dependent on unit cell geometric features. Such findings when coupled with where drugs move might explain the drug–device interactions.

Conway, Claire↗

Studying the Peak Force of Magnetic Linear Motion Devices Using Dimensional Neutrality

The merits of different sizing equations for linear motion magnetic force devices are assessed. A peak force scaling analysis for a linear magnetic coupling with respect to the volumetric force density, mass force density, magnetic shear stress and energy density is compared. By utilizing a 3-D analytic based model of a Halbach cylinder magnetic coupling it is shown that the maximizing of volumetric and mass force density is not an appropriate sizing metric for a linear magnetic actuator. Only energy density and magnetic shear stress are dimensionally neutral sizing metrics for linear actuators.

16 TIDAL AND WAVE POWER↗

Linear Reciprocating Tribometer for In Situ Neutron Reflectometry of Soft Matter

Neutron reflectometry is a technique for measuring structure near planar interfaces that has been previously used to non-destructively characterize the polymer density of hydrated, dilute, and soft materials. Previous investigations have conducted neutron reflectometry measurements of liquids, gels, emulsion, and polymer solutions at rest, in compression, and subject to shear stress. However, correlating structure with tribological properties of soft materials presents significant experimental challenges for prior instruments due to wall slip, sample thickness, and structural heterogeneity (e.g., depth-wise gradients). A linear reciprocating tribometer offers several advantages for in situ neutron reflectometry studies, including uniform velocity profiles, constant shear stress over large regions of interest, and independent control of normal force and sliding velocity during measurements. This work outlines basic considerations for the design of a custom linear reciprocating tribometer that operates in a neutron beamline and includes commissioning measurements. The tribometer is designed to compress soft and hydrated materials against linearly reciprocating silicon disks. The three key design considerations for this tribometer are (1) safety, (2) neutron transmission, and (3) sample positioning. This instrument design will enable in situ studies of soft matter and illuminate the role of interfacial structure on tribological phenomena.

42 ENGINEERING↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗

Liquid-phase turbulence measurements in air-water two-phase flows using particle image velocimetry

Liquid-phase turbulence plays a vital role in determining various gas-liquid two-phase flow parameters, such as void fraction distribution, bubble morphology, bubble-bubble interactions, and interfacial area concentration. In this study, a two-phase flow database including both the gas- and liquid-phases measurements was developed focusing on three bubbly flow conditions in an air-water two-phase flow loop with a vertical one-inch diameter circular pipe test section. A particle image velocimetry (PIV) system integrating an optical phase separation method, i.e., the planar laser-induced fluorescence (PLIF) technique using fluorescent particles and optical filtration, was applied to measure the liquid-phase turbulence information, including the time-averaged velocity, Reynolds stress, and turbulent kinetic energy for the liquid phase. The PIV measurements were taken at three ports along the test section at 14.5, 51.5, and 88.5 pipe inner diameters downstream of a bubble injector. In addition, a double-sensor conductivity probe was used to measure radial distributions of the local time-averaged void fraction and gas velocity. The measured liquid-phase turbulence was used to benchmark Sato's turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. Lastly, in the two bubbly flows tested that have low void fractions being less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction.

42 ENGINEERING↗

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗

Microstructural changes of proton irradiated Hastelloy-N and in situ micropillar compression testing of one single grain at different local damage levels

Here, in situ micropillar compression was used to study the deformation of proton-irradiated Hastelloy-N at different damage levels. Multiple pillars were prepared from a single domain along the cross-section of 2.5 MeV proton-irradiated Hastelloy-N. Depending on the location of micropillars, the critical resolved shear stress was obtained as a function of local damage levels. Such an approach eliminates the variation of yield stress due to the difference in the Schmid factor. Microstructural characterization showed that irradiation hardening is contributed by complicated defect structures, including (a) dislocation loops with many in corduroy-like alignments, (2) dislocations pile up, and (3) Cottrell atmosphere. Silicon atoms are found to segregate at dislocation lines and loops and form complicated Cottrell cloud patterns. The critical resolved shear stress, $t_{crss}$, and hardening exponents were both extracted as a function of displacements per atom values up to 2.3. There was a 60% increase in $t_{crss}$ at the highest damage level.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗

On the Morphodynamics of a Wide Class of Large‐Scale Meandering Rivers: Insights Gained by Coupling LES With Sediment‐Dynamics

Abstract In meandering rivers, interactions between flow, sediment transport, and bed topography affect diverse processes, including bedform development and channel migration. Predicting how these interactions affect the spatial patterns and magnitudes of bed deformation in meandering rivers is essential for various river engineering and geoscience problems. Computational fluid dynamics simulations can predict river morphodynamics at fine temporal and spatial scales but have traditionally been challenged by the large scale of natural rivers. We conducted coupled large‐eddy simulation and bed morphodynamics simulations to create a unique database of hydro‐morphodynamic data sets for 42 meandering rivers with a variety of planform shapes and large‐scale geometrical features that mimic natural meanders. For each simulated river, the database includes (a) bed morphology, (b) three‐dimensional mean velocity field, and (c) bed shear stress distribution under bankfull flow conditions. The calculated morphodynamics results at dynamic equilibrium revealed the formation of scour and deposition patterns near the outer and inner banks, respectively, while the location of point bars and scour regions around the apexes of the meander bends is found to vary as a function of the radius of curvature of the bends to the width ratio. A new mechanism is proposed that explains this seemingly paradoxical finding. The high‐fidelity simulation results generated in this work provide researchers and scientists with a rich numerical database for morphodynamics and bed shear stress distributions in large‐scale meandering rivers to enable systematic investigation of the underlying phenomena and support a range of river engineering applications.

54 ENVIRONMENTAL SCIENCES↗

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗