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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

CFD modeling and simulation for corrosive wear of refractory in molten slag

The development of high-wear resistant refractories having minimal production costs is facilitated by characterizing the wear mechanisms associated with their corrosive wear. Static cup testing is a commonly used method for comparing the corrosion resistance performance of two or more refractory materials. Although the static cup test conditions are not as severe as dynamic tests, this study shows that the thermal gradient present within the system during heating and cooling stages serves to generate movement of the slag leading to mechanical wear. The thermal gradient within the refractory, and between the slag and the refractory, occurs during the ramp stage of the test and lasts until the soaking stage is reached bringing the system to a thermal equilibrium. Using computational fluid dynamics (CFD) capabilities embedded within ANSYS software, this study modelled and quantified the convection currents within the slag and associated shear stresses generated on the refractory walls due to the thermal gradient. A traditional ladle furnace was employed as a case study to verify the results of the studied CFD model. The corrosion rate of the refractory lining was found to depend on the mass transfer coefficient of the refractory dissolution into the slag, and a velocity term which governs the extent of corrosion at any given location. This velocity term is a function of slag viscosity, as well as the concentration gradient and/or temperature gradient at the triple points. In this study, wall shear stress was used as a reliable proxy for identifying high-velocity regions prone to excessive corrosive wear. Elevated wall shear stress near the slag/air and slag/molten steel interfaces align with observed corrosion grooves, which reflects the intensified corrosive wear at these locations.

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Time-reversal symmetry breaking in two-dimensional nonequilibrium viscous fluids

In this work, we study the rheological signatures of departure from equilibrium in two-dimensional viscous fluids with and without internal spin. Under the assumption of isotropy, we provide the most general linear constitutive relations for stress and couple stress in terms of the velocity and spin fields. Invoking Onsager's regression hypothesis for fluctuations about steady states, we derive the Green-Kubo formulas relating the transport coefficients to time-correlation functions of the fluctuating stress. In doing so, we show that one of the nonequilibrium transport coefficients, the odd viscosity, requires time-reversal symmetry breaking in the case of systems without internal spin. However, the Green-Kubo relations for systems with internal spin also show that there is a possibility for nonvanishing odd viscosity even when time-reversal symmetry is preserved. Furthermore, we find that breakdown of equipartition in nonequilibrium steady states results in the decoupling of the two rotational viscosities relating the vorticity and the internal spin.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1974

Effect of crosslinking density on biaxial relaxation of SBR by using reduced variables

The use of reduced variables to account for the effect of crosslinking density in a styrene-butadiene rubber (SBR) system is demonstrated for general biaxial stress states. Recently published results from stress relaxation tests on five SBR vulcanizates crosslinked to different degrees by tetramethylthiuram disulfide were superposed by using the crosslinking density as a reduction variable. The equilibrium shear modulus calculated from the master relaxation curve at long reduced times was in satisfactory agreement with other results for SBR. The time-axis shifts were related in a linear logarithmic manner to the crosslinking density but had a slope slightly less than values previously reported for elastomer systems.

Arenz, R. J.↗
NASA NTRS1992

Internally heated mantle convection and the thermal and degassing history of the earth

An internally heated model of parameterized whole mantle convection with viscosity dependent on temperature and volatile content is examined. The model is run for 4l6 Gyr, and temperature, heat flow, degassing and regassing rates, stress, and viscosity are calculated. A nominal case is established which shows good agreement with accepted mantle values. The effects of changing various parameters are also tested. All cases show rapid cooling early in the planet's history and strong self-regulation of viscosity due to the temperature and volatile-content dependence. The effects of weakly stress-dependent viscosity are examined within the bounds of this model and are found to be small. Mantle water is typically outgassed rapidly to reach an equilibrium concentration on a time scale of less than 200 Myr for almost all models, the main exception being for models which start out with temperatures well below the melting temperature.

Williams, David R.↗
NASA NTRS1979

Reynolds number effects on the turbulence field in compressible boundary layers

Detailed experiments were conducted in a zero pressure gradient, supersonic turbulent boundary layer, including measurements of the three components of velocity fluctuations and the turbulent shear stress, for Reynolds numbers ranging from 11.7 million to 105 million at a freestream Mach number of 2.3. The mean flow measurements established the fully developed and equilibrium nature of the boundary layer. Measurements of the turbulence field show that the vertical and transverse fluctuations are essentially equal throughout the boundary layer at all Reynolds numbers, a feature that is different from observations in incompressible flows. The data show that the boundary layer exhibits similarity in the turbulence profiles for the entire Reynolds number range and agrees with previous compressible and incompressible data using Morkovin's scaling to account for compressibility effects.

Acharya, M.↗
NASA NTRS1983

Restricted maximum principles for elastic bodies

A maximum principle for the equilibrium of an elastic material body which is free of body forces is described not all of the components of the displacement vector or of the principal stresses can simultaneously have a strict maximum or minimum at any point in the body which does not be either on the surface or on a material interface.

Rose, M. E.↗
NASA NTRS1989

A Monte Carlo-finite element model for strain energy controlled microstructural evolution - 'Rafting' in superalloys

This paper presents a specialized microstructural lattice model, MCFET (Monte Carlo finite element technique), which simulates microstructural evolution in materials in which strain energy has an important role in determining morphology. The model is capable of accounting for externally applied stress, surface tension, misfit, elastic inhomogeneity, elastic anisotropy, and arbitrary temperatures. The MCFET analysis was found to compare well with the results of analytical calculations of the equilibrium morphologies of isolated particles in an infinite matrix.

Gayda, J.↗
NASA NTRS1996

Emission Spectral Measurements in the Plenum of an Arc Jet Wind Tunnel

Arc jet wind tunnel facilities are used to evaluate thermal protection system materials for re-entry vehicles. The high speed, high temperature flowfield generated by the arc jet can simulate the extreme aerodynamic heating environment experienced during re-entry so that the survivability of heat shield materials and performance of various designs options can be tested. Although the re-entry heating environment can be approximated in the arc jet facility, the flowfield only partially simulates the actual re-entry flight conditions. Reynolds numbers are not matched so that surface shear stress distributions and mass transfer rates due to ablation or other mechanisms are not modeled correctly. Unlike flight conditions the arc freestream air is in non-equilibrium because of the rapid expansion that occurs in the supersonic nozzle. To properly study the actual re-entry flow environment, computational fluid dynamics, computational chemistry and radiation models must be used. Arc jet tunnel tests serve to validate these models. To perform accurate simulations inlet and boundary-conditions are needed, which come from measurements of the flowfield. The present study is concerned with measurements in the plenum region of an arc heater. In the past, conditions in the arc heater flowfield have been predicted using simulations since conventional measurement techniques could not be used in the harsh extremely high temperature environment. The present study is part of a recent push to utilize optical techniques to help better characterize the arc jet flowfields. Emission measurements have been made in the shock layer and the constrictor section of the arc heater to determine temperatures and species number densities. LIF measurements have been made in the free stream to determine temperature and velocity.

Donohue, Jim↗
NASA NTRS2021

Financial feasibility of water conservation in agriculture

Global water use for food production needs to be reduced to remain within planetary boundaries, yet the financial feasibility of crucial measures to reduce water use is poorly quantified. Here, we introduce a novel method to compare the costs of water conservation measures with the added value that reallocation of water savings might generate if used for expansion of irrigation. Based on detailed water accounting through the use of a high-resolution hydrology-crop model, we modify the traditional cost curve approach with an improved estimation of demand, increasing marginal cost per water conservation measure combination and add a correction to control for impacts on downstream water availability. We apply the method to three major river basins in the Indo-Gangetic plain (Indus, Ganges and Brahmaputra), a major global food producing region but increasingly water stressed. Our analysis shows that at basin level only about 10% (Brahmaputra) to just over 20% (Indus and Ganges) of potential water savings would be realised; the equilibrium price for water is too low to make the majority of water conservation measures cost effective. The associated expansion of irrigated area is moderate, about 7% in the Indus basin, 5% in the Ganges and negligible in the Brahmaputra, but farmers' gross profit increases more substantially, by 11%. Increasing the volumetric cost of irrigation water influences supply and demand in a similar way and has little influence on water reallocation. Controlling for the impact on return flows is important and more than halves the amount of water available for reallocation.

Geography / Land Utilization↗
NASA NTRS2004

Nonlinear Analysis of Bonded Composite Single-LAP Joints

This study presents a semi-analytical solution method to analyze the geometrically nonlinear response of bonded composite single-lap joints with tapered adherend edges under uniaxial tension. The solution method provides the transverse shear and normal stresses in the adhesive and in-plane stress resultants and bending moments in the adherends. The method utilizes the principle of virtual work in conjunction with von Karman s nonlinear plate theory to model the adherends and the shear lag model to represent the kinematics of the thin adhesive layer between the adherends. Furthermore, the method accounts for the bilinear elastic material behavior of the adhesive while maintaining a linear stress-strain relationship in the adherends. In order to account for the stiffness changes due to thickness variation of the adherends along the tapered edges, their in-plane and bending stiffness matrices are varied as a function of thickness along the tapered region. The combination of these complexities results in a system of nonlinear governing equilibrium equations. This approach represents a computationally efficient alternative to finite element method. Comparisons are made with corresponding results obtained from finite-element analysis. The results confirm the validity of the solution method. The numerical results present the effects of taper angle, adherend overlap length, and the bilinear adhesive material on the stress fields in the adherends, as well as the adhesive, of a single-lap joint

Oterkus, E.↗
NASA NTRS1989

'Hypersonic slip flows' and issues on extending continuum model beyond the Navier-Stokes level

Results of a direct simulation Monte Carlo method for a hypersonic flow about a flat plate at a 40 deg angle of attack have been compared with corresponding results from a theory for fully viscous shock layers (FVSLs). Using the 13-moment equations for a Maxwell gas, it is demonstrated that nonequilibrium and equilibrium FVSL flows can be correlated. With the exception of the pressure density, most of the flow properties along a streamline (including heat flux and shear and normal stresses) are correctly predicted to leading order by the Navier-Stokes model.

Cheng, H. K.↗
NASA NTRS1994

A dynamic localization model with stochastic backscatter

The modeling of subgrid scales in large-eddy simulation (LES) has been rationalized by the introduction of the dynamic localization procedure. This method allows one to compute rather than prescribe the unknown coefficients in the subgrid-scale model. Formally, the LES equations are supposed to be obtained by applying to the Navier-Stokes equations a 'grid filter' operation. Though the subgrid stress itself is unknown, an identity between subgrid stresses generated by different filters has been derived. Although preliminary tests of the Dynamic Localization Model (DLM) with k-equation have been satisfactory, the use of a negative eddy viscosity to describe backscatter is probably a crude representation of the physics of reverse transfer of energy. Indeed, the model is fully deterministic. Knowing the filtered velocity field and the subgrid-scale energy, the subgrid stress is automatically determined. We know that the LES equations cannot be fully deterministic since the small scales are not resolved. This stems from an important distinction between equilibrium hydrodynamics and turbulence. In equilibrium hydrodynamics, the molecular motions are also not resolved. However, there is a clear separation of scale between these unresolved motions and the relevant hydrodynamic scales. The result of molecular motions can then be separated into an average effect (the molecular viscosity) and some fluctuations. Due to the large number of molecules present in a box with size of the order of the hydrodynamic scale, the ratio between fluctuations and the average effect should be very small (as a result of the 'law of large numbers'). For that reason, the hydrodynamic balance equations are usually purely deterministic. In turbulence, however, there is no clear separation of scale between small and large eddies. In that case, the fluctuations around a deterministic eddy viscosity term could be significant. An eddy noise would then appear through a stochastic term in the subgrid-scale model and could be the source of backscatter.

Carati, Daniele↗
DOE OSTI2022

Local stress within a granular molecular solvent matrix, a mechanism for individual ion hydration

The activity of water is important in many natural systems and directed processes. To better understand water activity, the hydration of individual ions is approximated for a recently developed speciation-based solution model. Here a matrix of values for hydration of complete salts (anion and cation) was used as source data to calculate individual ion hydration values. Some electrolytes were excluded from calculations of individual ion hydration due to the prevalence of ion pairing and resulting liberation of water upon ion pairing. Evaluation of the periodic trends emerging from individual ion hydration values suggests that ionic volume, relative to the volume of molecular water, dictates both modeled hydration requirements (related to vapor-liquid-equilibrium-derived activity) as well as properties historically described by hydrodynamic radius (ion diffusion and solution viscosity). When water is conceptualized as a granular three-dimensional matrix, solutes may function as substitutional defects (dilation or constriction) inducing local compressive or tensile stress. These results have implications on the mechanism of the Hofmeister series and series reversals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

A new family of stable elements for the Stokes problem based on a mixed Galerkin/least-squares finite element formulation

Adding to the classical Hellinger-Reissner formulation, a residual form of the equilibrium equation, a new Galerkin/least-squares finite element method is derived. It fits within the framework of a mixed finite element method and is stable for rather general combinations of stress and velocity interpolations, including equal-order discontinuous stress and continuous velocity interpolations which are unstable within the Galerkin approach. Error estimates are presented based on a generalization of the Babuska-Brezzi theory. Numerical results (not presented herein) have confirmed these estimates as well as the good accuracy and stability of the method.

Franca, Leopoldo P.↗
NASA NTRS1978

Environmental effects on composites for aircraft

A number of ongoing, long-term environmental effects programs for composite materials are evaluated. The flight service experience was evaluated for 142 composite aircraft components after more than 5 years and 1 million successful component flight hours. Ground-based outdoor exposures of composite material coupons after 3 years of exposure at 5 sites have reached equilibrium levels of moisture pickup which are predictable. Solar ultraviolet-induced material loss is discussed for these same exposures. No significant degradation was observed in residual strength for either stressed or unstressed specimens, or for exposures to aviation fuels and fluids.

Pride, R. A.↗
NASA NTRS2004

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗
NASA NTRS2004

Combined In-Plane and Through-the-Thickness Analysis for Failure Prediction of Bolted Composite Joints

Although two-dimensional methods provide accurate predictions of contact stresses and bolt load distribution in bolted composite joints with multiple bolts, they fail to capture the effect of thickness on the strength prediction. Typically, the plies close to the interface of laminates are expected to be the most highly loaded, due to bolt deformation, and they are usually the first to fail. This study presents an analysis method to account for the variation of stresses in the thickness direction by augmenting a two-dimensional analysis with a one-dimensional through the thickness analysis. The two-dimensional in-plane solution method based on the combined complex potential and variational formulation satisfies the equilibrium equations exactly, and satisfies the boundary conditions and constraints by minimizing the total potential. Under general loading conditions, this method addresses multiple bolt configurations without requiring symmetry conditions while accounting for the contact phenomenon and the interaction among the bolts explicitly. The through-the-thickness analysis is based on the model utilizing a beam on an elastic foundation. The bolt, represented as a short beam while accounting for bending and shear deformations, rests on springs, where the spring coefficients represent the resistance of the composite laminate to bolt deformation. The combined in-plane and through-the-thickness analysis produces the bolt/hole displacement in the thickness direction, as well as the stress state in each ply. The initial ply failure predicted by applying the average stress criterion is followed by a simple progressive failure. Application of the model is demonstrated by considering single- and double-lap joints of metal plates bolted to composite laminates.

Kradinov, V.↗