Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stress equilibrium”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
DOE OSTI2023

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↗
NASA NTRS1975

Analysis of bonded joints

A refined elastic analysis of bonded joints which accounts for transverse shear deformation and transverse normal stress was developed to obtain the stresses and displacements in the adherends and in the bond. The displacements were expanded in terms of polynomials in the thicknesswise coordinate; the coefficients of these polynomials were functions of the axial coordinate. The stress distribution was obtained in terms of these coefficients by using strain-displacement and stress-strain relations. The governing differential equations were obtained by integrating the equations of equilibrium, and were solved. The boundary conditions (interface or support) were satisfied to complete the analysis. Single-lap, flush, and double-lap joints were analyzed, along with the effects of adhesive properties, plate thicknesses, material properties, and plate taper on maximum peel and shear stresses in the bond. The results obtained by using the thin-beam analysis available in the literature were compared with the results obtained by using the refined analysis. In general, thin-beam analysis yielded reasonably accurate results, but in certain cases the errors were high. Numerical investigations showed that the maximum peel and shear stresses in the bond can be reduced by (1) using a combination of flexible and stiff bonds, (2) using stiffer lap plates, and (3) tapering the plates.

Srinivas, S.↗
NASA NTRS2001

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate, Epsilon(0) equations were considered. The emphasis of this paper is focused on the effects of the Epsilon(0)-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate (if change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters. Calculations show that the Epsilon(0)-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular Epsilon(0)-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the Epsilon(0)-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence.

Zhao, C. Y.↗
DOE OSTI2020

A phase-field model for non-isothermal phase transformation and plasticity in polycrystalline yttria-stabilized tetragonal zirconia

Here we propose an elastoplastic phase-field (PF) model to investigate the mechanics of tetragonal-to-monoclinic phase transformation (TMPT) and elastoplastic deformation of polycrystalline yttria-stabilized tetragonal zirconia (YSTZ). A Landau polynomial with non-vanishing chemical energy at the equilibrium temperature is introduced to account for the actual formation energies of the phases. The effects of different grain orientations, latent heat, and temperature on TMPT and deformation mechanisms are considered. The suppressive transformation effects of the grain boundaries (GBs) is modeled using an inhomogeneous kinetic coefficient in the bulk and GBs. The simulation results for single crystals demonstrate the capability of the model to reproduce the orientation-dependent compressive deformation of YSTZ similar to atomistic simulations and micropillar experiments. The single crystal with [100] crystallographic orientation along the loading direction (SC[100]) displays both TMPT and plasticity, SC[101] experiences only phase transformation, while SC[001] undergoes only plastic yielding. The TMPT induced by compressive loading exhibits shape memory effect (SME) below the equilibrium transformation temperature and pseudoelasticity above it, while the critical transformation stress increases with increasing loading temperature. The irrecoverable plastic strain is found to trap a part of the monoclinic phase, which prevents a complete reverse transformation. The polycrystalline cases also display SME and PE at low and high temperatures, respectively. Due to the orientation differences between grains and the stress concentrations at geometric nonlinearities, plastic deformation occurs in polycrystalline YSTZ for an applied load less than the yield stress. The results suggest a possible limitation of plasticity and an improvement of the shape recovery of YSTZ if one can control the orientation of the grains and/or increase the density of stacking faults at the GBs during material processing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Step-patterned survivorship curves: Mortality and loss of equilibrium responses to high temperature and food restriction in juvenile rainbow trout ( Oncorhynchus mykiss )

While survivorship curves typically exhibit smooth declines over time, step-patterned curves can occur with multiple stressors within a life stage. To explore this process, we examined the effects of heat (24°C) and food restriction on juvenile rainbow trout (Oncorhynchus mykiss Walbaum) in challenge experiments. We observed step-patterned survivorship curves determined by mortality and loss of equilibrium (LOE) endpoints. To examine the cause of heterogeneity in the stress responses from early to late mortality and LOE, we measured indices of energetic reserves. The step transition in the survivorship curves, the peak mortality rates, and start of when individuals reached a critical energetic threshold (14% dry mass; 4.0 kJ∙g -1 energy) all occurred at around days 10–15 of the challenge. The coherence in these temporal patterns suggest heterogeneity in the cohort stress responses, in which an early subgroup died from heat stress and a late subgroup died from starvation. Thus, their endpoint sensitivities resulted in step-patterned survivorship curves. We discuss the implications of the study for understanding effects of multiple stressors on population heterogeneity and note the possible significance of stress response selection under climate change in which heat stress and food limitations occur in concert.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Aeroelastic tailoring for oblique wing lateral trim

Composite material aeroelastic tailoring is presently explored as a means for the correction of the roll trim imbalance of oblique-wing aircraft configurations. The concept is demonstrated through the analysis of a realistic oblique wing by a static aeroelastic computational procedure encompassing the full potential transonic aerodynamic code FLO22 and a Ritz structural plate program that models the stiffness due to symmetrical-but-unbalanced composite wing skins. Results indicate that asymetric composite tailoring reduces the aileron deflection needed for roll equilibrium, and reduces control surface hinge moment and drag. Wing skin stresses are, however, very high.

Bohlmann, Jonathan D.↗
NASA NTRS1989

Evolution of coherent structures in the reattachment region of a separated flow

An experimental investigation of the evolution of coherent structures in a turbulent flow over a backward facing step has been made. Downstream of reattachment the large scale structures do not return rapidly to equilibrium boundary layer conditions and their coherent contribution to the Reynolds stresses continue to increase for some distance.

Jovic, S.↗
NASA NTRS2012

Global Patterns of Tectonism on Titan from Mountain Chains and Virgae

This research is based on the exploration of tectonic patterns on Titan from a global perspective. Several moons in the outer solar system display patterns of surface tectonic features that imply global stress fields driven or modified by global forces. Patterns such as these are seen in Europa's tidally induced fracture patterns, Enceladus's tiger stripes, and Ganymede's global expansion induced normal fault bands. Given its proximity to Saturn, as well as its eccentric orbit, tectonic features and global stresses may be present on Titan as well. Titan displays possible tectonic structures, such as mountain chains along its equator (Radebaugh et al. 2007), as well as the unexplored dark linear streaks termed virgae by the IAU. Imaged by Cassini with the RADAR instrument, mountain chains near the equator are observed with a predominante east-west orientation (Liu et al. 2012, Mitri et al. 2010). Orientations such as these can be explained by modifications in the global tidal stress field induced by global contraction followed by rotational spin-up. Also, due to Titan's eccentric orbit, its current rotation rate may be in an equilibrium between tidal spin-up near periapsis and spin-down near apoapsis (Barnes and Fortney 2003). Additional stress from rotational spin-up provides an asymmetry to the stress field. This, combined with an isotropic stress from radial contraction, favors the formation of equatorial mountain chains in an east-west direction. The virgae, which have been imaged by Cassini with both the Visual and Infrared Mapping Spectrometer (VIMS) and Imaging Science Subsystem (ISS) instruments, are located predominately near 30 degrees latitude in either hemisphere. Oriented with a pronounced elongation in the east-west direction, all observed virgae display similar characteristics: similar relative albedos as the surrounding terrain however darkened with an apparent neutral absorber, broken-linear or rounded sharp edges, and connected, angular elements with distinct, linear edges. Virgae imaged during northern latitude passes are oriented with their long dimensions toward Titan's antiSaturn point. If the virgae are of tectonic origin, for instance if the turn out to be i.e. grabens, they could serve as markers to Titan's global stress field. Using them in this way allows for a mapping of global tectonic patterns. These patterns will be tested for consistency against the various sources of global stress and orientations of mountain chains. By determining what drives Titan's tectonics globally, we will be able to place Titan within the context of the other outer planet icy satellites.

Cook, C.↗
DOE OSTI2020

On the inference of viscoelastic constants from stress relaxation experiments

Several constitutive theories have been proposed in the literature to model the viscoelastic response of materials, including widely used rheological constitutive models. These models are characterized by certain parameters ("time constants") that define the time scales over which the material relaxes. These parameters are primarily obtained from stress relaxation experiments using curve-fitting techniques. However, the question of how best to estimate these time constants remains open. As a step towards answering this question, we propose an optimal experimental design approach based on ideas from information geometry, namely Fisher information and Kullback-Leibler divergence. The material is modeled as a spring element in parallel with multiple Maxwell elements and described using a one- or two-term Prony series. Treating the time constants as unknowns, we develop expressions for the Fisher information and Kullback-Leibler divergence that allow us to maximize information gain from experimental data. Based on the results of this study, we propose that the largest time constant estimated from a stress relaxation experiment for a linear viscoelastic material should be at most one-fifth of the total time of the experiment in order to maximize information gain. Our results also provide confirmation that the equilibrium modulus of the material cannot be reliably determined from curve-fitting to data from a stress relaxation experiment.

Fisher information↗
NASA NTRS1987

An approximate semi-analytical method for prediction of interlaminar shear stresses in an arbitrarily laminated thick plate

An approximate semianalytical method for determination of interlaminar shear stress distribution through the thickness of an arbitrarily laminated thick plate has been presented. The method is based on the assumptions of transverse inextensibility and layerwise constant shear angle theory (LCST) and utilizes an assumed quadratic displacement potential energy based finite element method (FEM). Centroid of the triangular surface has been proved from a rigorous mathematical point of view (Aubin-Nitsche theory), to be the point of exceptional accuracy for the interlaminar shear stresses. Numerical results indicate close agreement with the available three-dimensional elasticity theory solutions. A comparison between the present theory and that due to an assumed stress hybrid FEM suggest that the (normal) traction-free-edge condition is not satisfied in the latter approach. Furthermore, the present paper is the first to present the results for interlaminar shear stresses in a two-layer thick square plate of balanced unsymmetric angle-ply construction. A comparison with the recently proposed Equilibrium Method (EM) indicates the superiority of the present method, because the latter assures faster convergence as well as simultaneous vanishing of the transverse shear stresses on both of the exposed surfaces of the laminate. Superiority of the present method over the EM, in the case of a symmetric laminate, is limited to faster convergence alone. It has also been demonstrated that the combination of the present method and the reduced (quadratic order) numerical integration scheme yields convergence of the interlaminar shear stresses almost as rapidly as that of the nodal displacements, in the case of a thin plate.

Chaudhuri, Reaz A.↗
NASA NTRS2000

Stress Recovery and Error Estimation for 3-D Shell Structures

The C1-continuous stress fields obtained from finite element analyses are in general lower- order accurate than are the corresponding displacement fields. Much effort has focussed on increasing their accuracy and/or their continuity, both for improved stress prediction and especially error estimation. A previous project developed a penalized, discrete least squares variational procedure that increases the accuracy and continuity of the stress field. The variational problem is solved by a post-processing, 'finite-element-type' analysis to recover a smooth, more accurate, C1-continuous stress field given the 'raw' finite element stresses. This analysis has been named the SEA/PDLS. The recovered stress field can be used in a posteriori error estimators, such as the Zienkiewicz-Zhu error estimator or equilibrium error estimators. The procedure was well-developed for the two-dimensional (plane) case involving low-order finite elements. It has been demonstrated that, if optimal finite element stresses are used for the post-processing, the recovered stress field is globally superconvergent. Extension of this work to three dimensional solids is straightforward. Attachment: Stress recovery and error estimation for shell structure (abstract only). A 4-node, shear-deformable flat shell element developed via explicit Kirchhoff constraints (abstract only). A novel four-node quadrilateral smoothing element for stress enhancement and error estimation (abstract only).

Riggs, H. R.↗
DOE OSTI2026

A critical review on additive manufacturing of refractory alloys from a data analytics perspective- beyond nickel-based superalloys

Refractory alloys (RAs) are promising materials due to their exceptional physicochemical properties, but most research remains at the laboratory scale. For broader adoption, advancements in manufacturing are essential. Because their high stability makes conventional methods like machining and casting difficult, additive manufacturing (AM) is emerging as an effective approach for fabricating refractory alloy components. However, AM's repeated non-equilibrium thermal cycles introduce undesired features (e.g. defects, anisotropic microstructures, and residual stresses), which are magnified due to RAs’ unique properties. This paper comprehensively reviews the state-of-the-art methods of AM for refractory alloys. It explores data analytics techniques to establish design rules based on multi-fidelity experimental and computational methods. Furthermore, it investigates integrated, collaborative efforts to harmonise standalone databases, information, knowledge, and predictive models at multi-physics, multi-stage, and multi-scale. Unlike the existing literature that focuses primarily on material systems or process fundamentals, this work provides an integrated perspective on AM of refractory alloys from a data analytics standpoint, highlighting the roles of integrated computational materials engineering (ICME), verification, validation, and uncertainty quantification (VV&UQ), and digital twin-driven qualification in overcoming data scarcity and accelerating rapid qualification.

Additive manufacturing↗
DOE OSTI2020

Mirror at the edge of the universe: Reflections on an accelerated boundary correspondence with de Sitter cosmology

An accelerated boundary correspondence (ABC) is solved for the de Sitter moving mirror cosmology. The beta Bogoliubov coefficients reveal the particle spectrum is a Planck distribution with temperature inversely proportional to horizon radius. The quantum stress tensor indicates a constant emission of energy flux consistent with eternal equilibrium, while the total energy carried by the particles remains finite. In conclusion, the curved spacetime transformation to flat spacetime with an accelerated boundary is illustrated, and also shown for anti-de Sitter (AdS) spacetime.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

The prediction and confirmation of critical epitaxial parameters

The coherency-incoherency transition in epitaxial crystals is said to take place at a critical misfit f(c) or, for a system in which a monolayer is subcritical, at a critical thickness h(c). In this paper, the physical principles and models used to predict critical parameters are analyzed and put into perspective. The dependence of the relevant principles on the equilibrium-nonequilibrium conditions under which the quantities are measured in practice is stressed. The advantages and disadvantages of the models used (essentially the Frenkel-Kontorowa and Volterra models) are highlighted.

Van Der Merwe, Jan↗
NASA NTRS1989

Effects of transverse shear on large deflection random response of symmetric composite laminates with mixed boundary conditions

Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived considering von Karman strains and transverse shear effects. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The effects of transverse shear on large deflection vibration of laminated plates with mixed boundary support conditions under random excitation are studied. Mean-square deflections and mean-square inplane strains are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for plates with two opposing edges simply supported and the remaining edges clamped.

Prasad, C. B.↗
NASA NTRS1992

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗