Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Superhydrophobic Array Devices for the Enhanced Formation of 3D Cancer Models

During the metastatic cascade, cancer cells travel through the bloodstream as circulating tumor cells (CTCs) to a secondary site. Clustered CTCs have greater shear stress and treatment resistance, yet their biology remains poorly understood. We therefore engineered a tunable superhydrophobic array device (SHArD). The SHArD-C was applied to culture a clinically relevant model of CTC clusters. Using our device, we cultured a model of cancer cell aggregates of various sizes with immortalized cancer cell lines. These exhibited higher E-cadherin expression and are significantly more capable of surviving high fluid shear stress-related forces compared to single cells and model clusters grown using the control method, helping to explain why clustering may provide a metastatic advantage. Additionally, the SHArD-S, when compared with the AggreWell 800 method, provides a more consistent spheroid-forming device culturing reproducible sizes of spheroids for multiple cancer cell lines. Overall, we designed, fabricated, and validated an easily tunable engineered device which grows physiologically relevant three-dimensional (3D) cancer models containing tens to thousands of cells.

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

Experimental and Numerical Study of Coating Thickness Using Multi-slot Air Knives

Gas-jet wiping is a widely employed production technology for controlling the final zinc coating thickness on a moving substrate during continuous hot-dip galvanizing. This paper presents an experimental investigation and numerical analysis of a prototype multi-slot air knife, which offers an increase in wiping efficiency relative to the traditional single-slot jet geometry in the continuous galvanizing process. The applicability of the analytical coating weight model of Elsaadawy et al. (Metall Mater Trans B, 38:413-424, 2007) to predict the final coating weight was determined for the multi-slot geometry, where particular focus was devoted to the effect of geometric parameters. Experimental measurements under a variety of knife geometry and process conditions agreed with the coating weight predictions of the analytical model. It was also shown that the air-knife geometric parameters had a significant effect on the pressure profile and shear stress distribution applied by the air knives to the moving substrate. It was determined that the final coating thickness was significantly affected by the auxiliary jet width, D{sub a}, where lighter coating weights at higher strip velocities (up to 5.4 pct at V{sub s} = 1.5 m/s) could be achieved by using the multi-slot air-knives prototype vs. the conventional single-slot configuration.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Role of Fe in long-range ordered Ni 2 Cr precipitates in Ni-Cr-Fe model alloys during isothermal aging

The precipitation of new phases during long-term service at elevated temperatures is a concern for the thermal stability of engineering alloys. In Ni-Cr-based alloys, e.g., Alloy 690, the formation of long-range ordered Ni 2 Cr causes embrittlement and may impact the lifetime of nuclear power plant components. In this work, we quantify the formation and evolution of Ni 2 Cr precipitation in eleven Ni-Cr-Fe model alloys with 0, 5, 7, and 10wt % Fe contents, and with Ni/Cr atomic ratios of 1.8, 2.0, 2.2, 2.4. These alloys were isothermally aged up to 10,000h at temperatures between 330 and 475°C. The alloys were characterized by synchrotron-based x-ray diffraction and Vickers hardness testing to quantify Ni 2 Cr precipitate size, and the impact of precipitate size on the mechanical properties as a function of Fe content. After 10,000 h of aging at 475°C and 418°C, the formation of Ni 2 Cr was observed in all alloys with 0 and 5wt % Fe. After 10,000h of aging at 418°C, Ni 2 Cr precipitates were also observed in the 7wt % Fe containing Ni/Cr=2.0 sample. Further, no clear evidence of Ni 2 Cr was observed in any of the 10wt % Fe samples at any time and temperature combination. We find that the face-centered cubic matrix lattice contraction and Vickers hardness are correlated with the Ni 2 Cr formation. The greatest change in hardness and lattice contraction occurs in stoichiometric alloys (Ni/Cr=2.0) with 0wt % Fe at 475°C. The rate of change in the material properties for the 5wt % Fe alloys is reduced, however the magnitude of changes is similar to 0wt % Fe alloys. A precipitation hardening model developed for Ni-Cr alloys based on critical resolved shear stress with weakly coupled dislocations shows a clear link between Ni 2 Cr precipitate size and hardness. This trend held across all alloys with Ni 2 Cr formation regardless of Fe concentration. This important structure-property relationship can potentially help define Ni-Cr-Fe-based component lifetimes directly through an understanding of how Ni 2 Cr formation impacts mechanical properties as a function of Fe content.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Adaptive Determination of the Optimal Exchange Location in Wall-Modeled Large-Eddy Simulation

Wall-modeled large-eddy simulation introduces a modeling interface (or exchange location) separating the wall-modeled layer from the rest of the domain. The current state-of-the-art is to rely on user expertise when choosing where to place this modeling interface, whether this choice is tied to the grid or not. This paper presents a postprocessing algorithm that determines the exchange location systematically. The algorithm is based on a model for the error in the predicted wall shear stress and a model for the computational cost, and then finds the exchange location that minimizes a combination of the two. Here, the algorithm is tested both a priori and a posteriori using an equilibrium wall model for the flow over a wall-mounted hump, a boundary layer in an adverse pressure gradient, and a shock/boundary-layer interaction. The algorithm produces exchange locations that mostly agree with what an experienced user would suggest, with thinner wall-modeled layers in nonequilibrium flow regions and thicker wall-modeled layers where the boundary layer is closer to equilibrium. This suggests that the algorithm should be useful in simulations of realistic and highly complex geometries.

42 ENGINEERING↗
DOE OSTI2024

Assessment of Machine Learning Wall Modeling Approaches for Large Eddy Simulation of Gas Turbine Film Cooling Flows: An a Priori Study

Here, in this work, a priori analysis of machine learning (ML) strategies is carried out with the goal of data-driven wall modeling for large eddy simulation (LES) of gas turbine film cooling flows. High-fidelity flow datasets are extracted from wall-resolved LES (WRLES) of flow over a flat plate interacting with the coolant flow supplied by a single row of 7-7-7 shaped cooling holes inclined at 30 degrees with the flat plate at different blowing ratios (BR). The WRLES are performed using the high-order Nek5000 spectral element computational fluid dynamics (CFD) solver. Light gradient boosting machine (LightGBM) is employed as the ML algorithm for the data-driven wall model. Parametric tests are conducted to systematically assess the influence of a wide range of input flow features (velocity components, velocity gradients, pressure gradients, and fluid properties) on the accuracy of ML wall model with respect to prediction of wall shear stress. In addition, the use of spatial stencil and time delay is also explored within the ML wall modeling framework. It is shown that features associated with gradients of the streamwise and spanwise velocity components have a major impact on the prediction fidelity of wall model, while the effect of gradients of wall-normal velocity component is found to be negligible. Moreover, adding flow feature information from an x-y-z spatial stencil significantly improves the ML model accuracy and generalizability compared to just using local flow features from the matching location. Overall, highest prediction accuracy is achieved when both spatial stencil and time delay features are incorporated within the data-driven wall modeling paradigm.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Hydrogen effects on the deformation and slip localization in a single crystal austenitic stainless steel

Hydrogen is known to embrittle austenitic stainless steels, which are widely used in high-pressure hydrogen storage and delivery systems, but the mechanisms that lead to such material degradation are still being elucidated. The current work investigates the deformation behavior of single crystal austenitic stainless steel 316L through combined uniaxial tensile testing, characterization and atomistic simulations. Thermally precharged hydrogen is shown to increase the critical resolved shear stress (CRSS) without previously reported deviations from Schmid’s law. Molecular dynamics simulations further expose the statistical nature of the hydrogen and vacancy contributions to the CRSS in the presence of alloying. Slip distribution quantification over large in-plane distances (> 1 mm), achieved via atomic force microscopy (AFM), highlights the role of hydrogen increasing the degree of slip localization in both single and multiple slip configurations. The most active slip bands accumulate significantly more deformation in hydrogen precharged specimens, with potential implications for damage nucleation. For $\langle$110$\rangle$ tensile loading, slip localization further enhances the activity of secondary slip, increases the density of geometrically necessary dislocations and leads to a distinct lattice rotation behavior compared to hydrogen-free specimens, as evidenced by electron backscatter diffraction (EBSD) maps. Finally, the results of this study provide a more comprehensive picture of the deformation aspect of hydrogen embrittlement in austenitic stainless steels.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗
DOE OSTI2022

Evaluation of intracoronary hemodynamics identifies perturbations in vorticity

Coronary artery disease (CAD) is highly prevalent and associated with adverse events. Challenges have emerged in the treatment of intermediate coronary artery stenoses. These lesions are often interrogated with fractional flow reserve (FFR) testing to determine if a stenosis is likely to be causative for ischemia in a cardiac territory. This invasive test requires insertion of a pressure wire into a coronary vessel. Recently computational fluid dynamics (CFD) has been used to noninvasively assess fractional flow reserve in vessels reconstructed from medical imaging data. However, many of these simulations are unable to provide additional information about intravascular hemodynamics, including velocity, endothelial shear stress (ESS), and vorticity. We hypothesized that vorticity, which has demonstrated utility in the assessment of ventricular and aortic diseases, would also be an important hemodynamic factor in CAD. Three-dimensional (3D), patient-specific coronary artery geometries that included all vessels >1 mm in diameter were created from angiography data obtained from 10 patients who underwent diagnostic angiography and FFR testing (n = 9). A massively parallel CFD solver (HARVEY) was used to calculate coronary hemodynamic parameters including pressure, velocity, ESS, and vorticity. These simulations were validated by comparing velocity flow fields from simulation to both velocities derived from in vitro particle image velocimetry and to invasively acquired pressure wire-based data from clinical testing. There was strong agreement between findings from CFD simulations and particle image velocimetry experimental testing (p < 0.01). CFD-FFR was also highly correlated with invasively measured FFR (ρ = 0.77, p = 0.01) with an average error of 5.9 ± 0.1%. CFD-FFR also had a strong inverse correlation with the vorticity (ρ = -0.86, p = 0.001). Simulations to determine the effect of the coronary stenosis on intravascular hemodynamics demonstrated significant differences in velocity and vorticity (both p < 0.05). Further evaluation of an angiographically normal appearing non-FFR coronary vessel in patients with CAD also demonstrated differences in vorticity when compared with FFR vessels (p < 0.05). The use of highly accurate 3D CFD-derived intravascular hemodynamics provides additional information beyond pressure measurements that can be used to calculate FFR. Vorticity is one parameter that is modified by a coronary stenosis and appears to be abnormal in angiographically normal vessels in patients with CAD, highlighting a possible use-case in preventative screening for early coronary disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗
NASA NTRS1977

Edge effects in angle-ply composite laminates

This paper presents the results of a zeroth-order solution for edge effects in angle-ply composite laminates obtained using perturbation techniques and a limiting free body approach. The general solution for edge effects in laminates of arbitrary angle ply is applied to the special case of a (+ or - 45)s graphite/epoxy laminate. Interlaminar stress distributions are obtained as a function of the laminate thickness-to-width ratio and compared to finite difference results. The solution predicts stable, continuous stress distributions, determines finite maximum tensile interlaminar normal stress and provides mathematical evidence for singular interlaminar shear stresses in (+ or - 45) graphite/epoxy laminates.

Hsu, P. W.↗
NASA NTRS1981

Thermomechanical response of Gr/Pi composites

The effects of temperature changes upon the stresses and strains in composite laminates having carbon fibers in a polyimide matrix were evaluated. Composites having laminae in which there were nonlinear stress-strain relations for stresses transverse to the fibers and for axial shear stresses were treated. Material properties were considered to be temperature dependent. Separately, the effects of laminae viscoelastic response were also treated. The results suggest that, for this material, nonlinearities due either to stress or time dependent effects do not appear to be of major practical importance for conventional high temperature composite structures.

Rosen, B. W.↗
NASA NTRS1992

Improved assumed-stress hybrid shell element with drilling degrees of freedom for linear stress, buckling, and free vibration analyses

An improved four-node quadrilateral assumed-stress hybrid shell element with drilling degrees of freedom is presented. The formulation is based on Hellinger-Reissner variational principle and the shape functions are formulated directly for the four-node element. The element has 12 membrane degrees of freedom and 12 bending degrees of freedom. It has nine independent stress parameters to describe the membrane stress resultant field and 13 independent stress parameters to describe the moment and transverse shear stress resultant field. The formulation encompasses linear stress, linear buckling, and linear free vibration problems. The element is validated with standard tests cases and is shown to be robust. Numerical results are presented for linear stress, buckling, and free vibration analyses.

Rengarajan, Govind↗
NASA NTRS2015

Relation Between Residual and Hoop Stresses and Rolling Bearing Fatigue Life

Rolling-element bearings operated at high speed or high vibration may require a tight interference fit between the bore of the bearing and shaft to prevent rotation of the bearing bore around the shaft and fretting damage at the interfaces. Previous work showed that the hoop stresses resulting from tight interference fits can reduce bearing lives by as much as 65 percent. Where tight interference fits are required, case-carburized steel such as AISI 9310 or M50 NiL is often used because the compressive residual stresses inhibit subsurface crack formation and the ductile core inhibits inner-ring fracture. The presence of compressive residual stress and its combination with hoop stress also modifies the Hertz stress-life relation. This paper analyzes the beneficial effect of residual stresses on rolling-element bearing fatigue life in the presence of high hoop stresses for three bearing steels. These additional stresses were superimposed on Hertzian principal stresses to calculate the inner-race maximum shearing stress and the resulting fatigue life of the bearing. The load-life exponent p and Hertz stress-life exponent n increase in the presence of compressive residual stress, which yields increased life, particularly at lower stress levels. The Zaretsky life equation is described and is shown to predict longer bearing lives and greater load- and stress-life exponents, which better predicts observed life of bearings made from vacuum-processed steel.

Fatigue↗
NASA NTRS1971

Time-dependent polymer rheology under constant stress and under constant shear conditions.

A kinetic rate theory previously presented for describing non-Newtonian phenomena has been further modified to predict the flow behavior of viscoelastic materials under constant stress conditions. The thixotropic shear stress or shear rate is predicted by the kinetic theory, and the experimental stress or shear rate is obtained by modifying the thixotropic value by a stress or shear rate retardation term. The retardation term stems from a Maxwellian approach for stress retardation. In order to test the validity of this approach, transient and steady-state data were obtained for two solutions of polymethylmethacrylate in diethylphthalate. Both constant stress measurements and constant shear rate data were taken over a broad range.

Lee, K. H.↗
NASA NTRS1939

Experimental and analytical investigation of a monocoque wing model loaded in bending

Bending tests with transverse loads and with pure bending were undertaken on a double-web monocoque wing model in order to establish the relation between the state of stress and the results from the elementary bending theory. The longitudinal stresses in the stiffeners were measured with tensiometers and the shear stresses in the sheet were calculated from them. The measurements were made at both moderate loads with no buckles in the covering and at loads with which the critical stress in the individual panels was exceeded. For the comparison, the wing skin was considered as stiffened sheet according to the shear panel scheme. In this way, the statistically indeterminate calculation was confirmed by the test results.

Schapitz, E↗