Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress field”

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 667 records · Page 37

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗

Wall Shear Stress, Wall Pressure and Near Wall Velocity Field Relationships in a Whirling Annular Seal

The mean and phase averaged pressure and wall shear stress distributions were measured on the stator wall of a 50% eccentric annular seal which was whirling in a circular orbit at the same speed as the shaft rotation. The shear stresses were measured using flush mounted hot-film probes. Four different operating conditions were considered consisting of Reynolds numbers of 12,000 and 24,000 and Taylor numbers of 3,300 and 6,600. At each of the operating conditions the axial distribution (from Z/L = -0.2 to 1.2) of the mean pressure, shear stress magnitude, and shear stress direction on the stator wall were measured. Also measured were the phase averaged pressure and shear stress. These data were combined to calculate the force distributions along the seal length. Integration of the force distributions result in the net forces and moments generated by the pressure and shear stresses. The flow field inside the seal operating at a Reynolds number of 24,000 and a Taylor number of 6,600 has been measured using a 3-D laser Doppler anemometer system. Phase averaged wall pressure and wall shear stress are presented along with phase averaged mean velocity and turbulence kinetic energy distributions located 0.16c from the stator wall where c is the seal clearance. The relationships between the velocity, turbulence, wall pressure and wall shear stress are very complex and do not follow simple bulk flow predictions.

Morrison, Gerald L.↗

Calculation of stress intensity factors in an isotropic multicracked plate. Part 1: Theoretical development

An essential part of describing the damage state and predicting the damage growth in a multicracked plate is the accurate calculation of stress intensity factors (SIF's). Here, a methodology and rigorous solution formulation for SIF's of a multicracked plate, with fully interacting cracks, subjected to a far-field arbitrary stress state is presented. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. This analytical derivation and numerical solution are obtained by using intelligent application of symbolic computations and automatic FORTRAN generation capabilities (described in the second part of this paper). As a result, a symbolic/FORTRAN package, named SYMFRAC, that is capable of providing accurate SIF's at each crack tip was developed and validated.

Binienda, W. K.↗

Elasto visco-plastic flow with special attention to boundary conditions

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Shimazaki, Y.↗

Predicting Two-Dimensional, Unsteady Turbulent Combustion

MIMOC calculates flow field and flame propagation in turbulent combustion tunnel. Program employs algorithm for turbulent combustion modeling described by Ghoniem, Chorin, and Oppenheim. Program used to study flow field in model combustor, formed by rearward-facing step in channel, in terms of vorticity field, velocity field, turbulent shear stresses, flame contours, and concentration field.

Marek, C. J.↗

SRM Internal Flow Tests and Computational Fluid Dynamic Analysis: Titan, ASRM, and Subscale Motor Analyses - Volume 3

A computational fluid dynamics (CFD) analysis has been performed on the aft slot region of the Titan 4 Solid Rocket Motor Upgrade (SRMU). This analysis was performed in conjunction with MSFC structural modeling of the propellant grain to determine if the flow field induced stresses would adversely alter the propellant geometry to the extent of causing motor failure. The results of the coupled CFD/stress analysis have shown that there is a continual increase of flow field resistance at the aft slot due to the aft segment propellant grain being progressively moved radially toward the centerline of the motor port. This 'bootstrapping' effect between grain radial movement and internal flow resistance is conducive to causing a rapid motor failure.

Source record↗

Driving Force Analysis in an Infinite Anisotropic Plate with Multiple Crack Interactions

The methodology and a rigorous solution formulation are presented for stress intensity factors (SIF's, k) and total strain energy release rates (SERR, G(sub T)) of a multicracked plate, that has fully interacting cracks and is subjected to a far-field arbitrary stress state. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. Parametric studies are conducted for two, three, and four crack problems. The sensitivity and characteristics of the model is demonstrated.

Binienda, Wieslaw K.↗

Altimetry data and the elastic stress tensor of subduction zones

The maximum shear stress (mss) field due to mass anomalies is estimated in the Apennines, the Kermadec-Tonga Trench, and the Rio Grande Rift areas and the results for each area are compared to observed seismicity. A maximum mss of 420 bar was calculated in the Kermadec-Tonga Trench region at a depth of 28 km. Two additional zones with more than 300 bar mss were also observed in the Kermadec-Tonga Trench study. Comparison of the calculated mss field with the observed seismicity in the Kermadec-Tonga showed two zones of well correlated activity. The Rio Grande Rift results showed a maximum mss of 700 bar occurring east of the rift and at a depth of 6 km. Recorded seismicity in the region was primarily constrained to a depth of approximately 5 km, correlating well to the results of the stress calculations. Two areas of high mss are found in the Apennine region: 120 bar at a depth of 55 km, and 149 bar at the surface. Seismic events observed in the Apennine area compare favorably with the mss field calculated, exhibiting two zones of activity. The case of loading by seamounts and icecaps are also simulated. Results for this study show that the mss reaches a maximum of about 1/3 that of the applied surface stress for both cases, and is located at a depth related to the diameter of the surface mass anomaly.

Caputo, Michele↗

Modeling lattice rotation fields from discrete crystallographic slip bands in superalloys

Here in this work, we investigate the relationship between an intense slip band (ISB) and the zone of large lattice rotations that forms ahead of the tip of the ISB. We develop a crystal plasticity finite element model of a discrete ISB lying within an oligocrystalline assembly and calculate the local crystalline stress and lattice rotation fields generated by the ISB. The calculations demonstrate that, first, a region of severe lattice rotations, commonly referred to as a microvolume, does not form without the ISB, and second, large amounts of accumulated slip in the ISB are required to enlarge the microvolume to sizes and rotation magnitudes observed experimentally. Ahead of the ISB tip, the quintessential plastic zone always forms, but the atypical microvolume forms when non-concentrated and spatially diffuse slip is activated by the ISB-induced stress field. This result suggests that the detrimental ISB/microvolume pair will likely appear in pairs of crystals in which transmission of the slip from the ISB is severely blocked by the grain boundary, a hypothesis that we verify with a few target cases.

36 MATERIALS SCIENCE↗

Computationally Efficient Beam Elements for Accurate Stresses in Sandwich Laminates and Laminated Composites with Delaminations

Laminated composites are prone to delamination failure due to the lack of reinforcement through the thickness. Therefore, during the design process the initiation and propagation of delaminations should be accounted for as early as possible. This paper presents computationally efficient nine degree-of-freedom (dof) and eight-dof shear locking-free beam elements using the mixed form of the refined zigzag theory (RZT(sup(m))). The corresponding nine-dof and eight-dof elements use the anisoparametric and constrained anisoparametric interpolation schemes, respectively, to eliminate shear locking in slender beams. The advantage of the present element over previous RZT beam elements is that no post-processing is required to accurately model the transverse shear stress while maintaining the computational efficiency of a low-order beam element. Comparisons with high-fidelity finite element models and three-dimensional elasticity solutions show that the elements can robustly and accurately predict the displacement field, axial stress and transverse shear stress through the thickness of a sandwich beam or a composite laminate with an embedded delamination. In fact, the accuracy and computational efficiency of predicting stresses in laminates with embedded delaminations make the present elements attractive choices for RZT-based delamination initiation and propagation methodologies available in the literature.

R. M. J. Groh↗

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗

Surface Engineering by Simultaneous Action of Multiple External Fields (Final Technical Report)

This research project aims at enabling surface engineering strategies based on the fundamental understanding and predictive modeling of the surface morphological response of solid materials subjected to the combined action of multiple external forces. Surface morphological response to mismatch strain in strained-layer heterostructures also is analyzed. Our research focuses on the surface morphological evolution and stability of electrically conducting (metallic, such as Cu) and semiconducting (such as Si) solids under the simultaneous application of mechanical stresses, electric fields, and temperature gradients. We emphasize on identifying the conditions under which the multiply driven surface morphology is stable and on exploring the complexity of the corresponding various morphologically stable surface patterns. We also emphasize on the driven evolution of nanoscale features aiming at stabilizing and controlling nanoscale patterns on surfaces by their manipulation through simultaneously applied multiple external forces. Specifically, we address systematically the morphological response to the combined action of multiple external forces of bulk solid surfaces, surfaces of thin films grown epitaxially on thick or thin substrates, as well as nanoscale surface features such as coherently strained islands grown epitaxially on substrate surfaces. The research is based on a modeling approach that combines theoretical analyses of surface morphological stability with self-consistent dynamical simulations of surface evolution based on properly parameterized continuum and multiscale surface transport models.

36 MATERIALS SCIENCE↗

Plasma bulk flow in Jupiter's dayside middle magnetosphere

Using the plasma data obtained during the Voyager 1 encounter and the full response function of the Plasma Science (PLS) experiment, convective plasma velocities have been determined in the dayside middle magnetosphere of Jupiter (r = 10-25 Jupiter radii). It is found that temperature anisotropies have very little effect on plasma velocity determination and that the plasma data are well approximated by convected, isotropic Maxwellian ion distribution functions. The insensitivity of the analysis to any thermal anisotropies which may exist allows a good determination of the bulk plasma flow velocity. In addition to the subcorotational azimuthal flow, there exists a substantial nonazimuthal component of plasma flow. This nonazimuthal flow is mostly aligned (antialigned) with the local magnetic field but also exhibits a cross-field component. The velocity pattern is inconsistent with enhanced plasma outflow in the active sector, as suggested by the corotating convection model of plasma transport. The contribution of field-aligned flow along the curved magnetic field lines to the stress on the magnetic field is evaluated. In the region studied, such flow contributes up to one half the stress produced by the azimuthal plasma flow.

Sands, Mark R.↗

Estimating the Relative Water Content of Leaves in a Cotton Canopy

Remotely sensing plant canopy water status remains a long-term goal of remote sensing research. Established approaches to estimating canopy water status the Crop Water Stress Index, the Water Deficit Index and the Equivalent Water Thickness involve measurements in the thermal or reflective infrared. Here we report plant water status estimates based upon analysis of polarized visible imagery of a cotton canopy measured by ground Multi-Spectral Polarization Imager (MSPI). Such estimators potentially provide access to the plant hydrological photochemistry that manifests scattering and absorption effects in the visible spectral region.Twice during one day, +- 3 hours from solar noon, we collected polarized imagery and relative water content data on a cotton test plot located at the Arid Land Agricultural Research Center, United States Department of Agriculture, Maricopa, AZ. The test plot, a small portion of a large cotton field, contained stressed plants ready for irrigation. The evening prior to data collection we irrigated several rows of plants within the test plot. Thus, ground MSPI imagery from both morning and afternoon included cotton plants with a range of water statuses. Data analysis includes classifying the polarized imagery into sunlit reflecting, sunlit transmitting, shaded foliage and bare soil. We estimate the leaf surface reflection and interior reflection based upon the per pixel polarization and sunview directions. We compare our cotton results with our prior polarization results for corn and soybean leaves measured in the lab and corn leaves measured in the field.

Cotton Canopy↗

Magnetism in curved geometries

Curvature impacts physical properties across multiple length scales, ranging from the macroscopic scale, where the shape and size vary drastically with the curvature, to the nanoscale at interfaces and inhomogeneities in materials with structural, chemical, electronic, and magnetic short-range order. In quantum materials, where correlations, entanglement, and topology dominate, the curvature opens the path to novel characteristics and phenomena that have recently emerged and could have a dramatic impact on future fundamental and applied studies of materials. Particularly, magnetic systems hosting non-collinear and topological states and 3D magnetic nanostructures strongly benefit from treating curvature as a new design parameter to explore prospective applications in the magnetic field and stress sensing, microrobotics, and information processing and storage. This Perspective gives an overview of recent progress in synthesis, theory, and characterization studies and discusses future directions, challenges, and application potential of the harnessing curvature for 3D nanomagnetism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shape and Stress Sensing of Multilayered Composite and Sandwich Structures Using an Inverse Finite Element Method

The marked increase in the use of composite and sandwich material systems in aerospace, civil, and marine structures leads to the need for integrated Structural Health Management systems. A key capability to enable such systems is the real-time reconstruction of structural deformations, stresses, and failure criteria that are inferred from in-situ, discrete-location strain measurements. This technology is commonly referred to as shape- and stress-sensing. Presented herein is a computationally efficient shape- and stress-sensing methodology that is ideally suited for applications to laminated composite and sandwich structures. The new approach employs the inverse Finite Element Method (iFEM) as a general framework and the Refined Zigzag Theory (RZT) as the underlying plate theory. A three-node inverse plate finite element is formulated. The element formulation enables robust and efficient modeling of plate structures instrumented with strain sensors that have arbitrary positions. The methodology leads to a set of linear algebraic equations that are solved efficiently for the unknown nodal displacements. These displacements are then used at the finite element level to compute full-field strains, stresses, and failure criteria that are in turn used to assess structural integrity. Numerical results for multilayered, highly heterogeneous laminates demonstrate the unique capability of this new formulation for shape- and stress-sensing.

Cerracchio, Priscilla↗

Notch-tensile behavior of Al0.1CrFeCoNi high entropy alloy

Notch-tensile behavior of an Al0.1CrFeCoNi high entropy alloy (Al0.1-HEA) was studied using V-notch geometry and with local strain mapping using digital image correlation (DIC). A notch-strength ratio of 1.51 indicated notch strengthening. Further analysis of stress-strain response supplemented with microstructural analysis revealed that, while the presence of notch results in strengthening due to work hardening by twinning induced plasticity, the notch also contributes strongly to geometrical softening in the non-uniform ductility regime, and accounts for the onset of failure. The strain localization behavior of Al0.1-HEA due to the presence of V-notch was compared with three conventional alloys: Inconel 625 nickel-based superalloy, 304 stainless steel and Ti–6Al–4V titanium alloy. The study revealed that the nature of notch widening with increasing strain was dependent on material characteristics. The extent of notch widening impacted the local strain field and stress distribution, thereby influencing the propensity for crack initiation and growth. The experimental results were verified by finite element analysis.

Sinha, Subhasis↗

Strain intensity factor approach for predicting the strength of continuously reinforced metal matrix composites

A method was previously developed to predict the fracture toughness (stress intensity factor at failure) of composites in terms of the elastic constants and the tensile failing strain of the fibers. The method was applied to boron/aluminum composites made with various proportions of 0 to + or - 45 deg plies. Predicted values of fracture toughness were in gross error because widespread yielding of the aluminum matrix made the compliance very nonlinear. An alternate method was developed to predict the strain intensity factor at failure rather than the stress intensity factor because the singular strain field was not affected by yielding as much as the stress field. Strengths of specimens containing crack-like slits were calculated from predicted failing strains using uniaxial stress-strain curves. Predicted strengths were in good agreement with experimental values, even for the very nonlinear laminates that contained only + or - 45 deg plies. This approach should be valid for other metal matrix composites that have continuous fibers.

Poe, C. C., Jr.↗