Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

LES/RANS Simulation of a Supersonic Reacting Wall Jet

This work presents results from large-eddy / Reynolds-averaged Navier-Stokes (LES/RANS) simulations of the well-known Burrows-Kurkov supersonic reacting wall-jet experiment. Generally good agreement with experimental mole fraction, stagnation temperature, and Pitot pressure profiles is obtained for non-reactive mixing of the hydrogen jet with a non-vitiated air stream. A lifted flame, stabilized between 10 and 22 cm downstream of the hydrogen jet, is formed for hydrogen injected into a vitiated air stream. Flame stabilization occurs closer to the hydrogen injection location when a three-dimensional combustor geometry (with boundary layer development resolved on all walls) is considered. Volumetric expansion of the reactive shear layer is accompanied by the formation of large eddies which interact strongly with the reaction zone. Time averaged predictions of the reaction zone structure show an under-prediction of the peak water concentration and stagnation temperature, relative to experimental data and to results from a Reynolds-averaged Navier-Stokes calculation. If the experimental data can be considered as being accurate, this result indicates that the present LES/RANS method does not correctly capture the cascade of turbulence scales that should be resolvable on the present mesh. Instead, energy is concentrated in the very largest scales, which provide an over-mixing effect that excessively cools and strains the flame. Predictions improve with the use of a low-dissipation version of the baseline piecewise parabolic advection scheme, which captures the formation of smaller-scale structures superimposed on larger structures of the order of the shear-layer width.

Edwards, Jack R.↗
DOE OSTI2022

Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization

Abstract A fundamental understanding of the interaction between microstructure and underlying physical mechanisms is essential, especially for developing more accurate multi-physics models for heterogeneous materials. Effects of microstructure on the material response at the macroscale are modeled by using the generalized thermomechanics. In this study, strain gradient theory is employed as a higher-order theory on the macroscale with thermodynamics modeled as a first-order theory on the microscale. Hence, energy depends only on the temperature such that we circumvent an extension of Fourier’s law and analyze the “simplest” thermo-mechanical model in strain gradient elasticity. Developing multiphysics models for heterogeneous materials is indeed a challenge and even this “simplest” model in generalized thermomechanics creates dozens of parameters to be determined. We develop a thermo-mechanical framework, in which microstructure is modeled as a periodic structure and through asymptotic homogenization approach, higher-order parameters at macroscopic scale are calculated. To illustrate the importance of higher-order parameters in overall thermo-mechanical response of a heterogeneous materials, finite element method (FEM) is employed with the aid of open-source codes (FEniCS). Verification example of a bulk system and several case studies of porous structures demonstrate how such numerical framework can be beneficial in the design of materials with tailored microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Cascade electrocatalysis via AgCu single-atom alloy and Ag nanoparticles in CO 2 electroreduction toward multicarbon products

Electrocatalytic CO 2 reduction into value-added multicarbon products offers a means to close the anthropogenic carbon cycle using renewable electricity. However, the unsatisfactory catalytic selectivity for multicarbon products severely hinders the practical application of this technology. In this paper, we report a cascade AgCu single-atom and nanoparticle electrocatalyst, in which Ag nanoparticles produce CO and AgCu single-atom alloys promote C-C coupling kinetics. As a result, a Faradaic efficiency (FE) of 94 ± 4% toward multicarbon products is achieved with the as-prepared AgCu single-atom and nanoparticle catalyst under ~720 mA cm -2 working current density at -0.65 V in a flow cell with alkaline electrolyte. Density functional theory calculations further demonstrate that the high multicarbon product selectivity results from cooperation between AgCu single-atom alloys and Ag nanoparticles, wherein the Ag single-atom doping of Cu nanoparticles increases the adsorption energy of *CO on Cu sites due to the asymmetric bonding of the Cu atom to the adjacent Ag atom with a compressive strain.

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

An ILLIAC program for the numerical simulation of homogeneous incompressible turbulence

An algorithm and ILLIAC computer program, developed for the simulation of homogeneous incompressible turbulence in the presence of an applied mean strain, are described. The turbulence field is represented spatially by a truncated triple Fourier series (spectral method) and followed in time using a fourth-order Runge-Kutta algorithm. These include: (1) transformation of variables suggested by Taylor's sudden-distortion theory; (2) implicit viscous diffusion by use of an integrating factor; (3) implicit pressure calculation suggested by Taylor's sudden-distortion theory, and (4) inexpensive control of aliasing by random and phased coordinate shifts.

Rogallo, R. S.↗
NASA NTRS1986

Analysis of Delamination Growth from Matrix Cracks in Laminates Subjected to Bending Loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg. plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg. plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg. ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg. ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg. ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, G. B.↗
NASA NTRS1988

Analysis of delamination growth from matrix cracks in laminates subjected to bending loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, Gretchen Bostaph↗
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗
NASA NTRS2001

Monitoring of bacteria growth using a wireless, remote query resonant-circuit sensor: application to environmental sensing

A new technique is presented for in-vivo remote query measurement of the complex permittivity spectra of a biological culture solution. A sensor comprised of a printed inductor-capacitor resonant-circuit is placed within the culture solution of interest, with the impedance spectrum of the sensor measured using a remotely located loop antenna; the complex permittivity spectra of the culture is calculated from the measured impedance spectrum. The remote query nature of the sensor platform enables, for example, the in-vivo real-time monitoring of bacteria or yeast growth from within sealed opaque containers. The wireless monitoring technique does not require a specific alignment between sensor and antenna. Results are presented for studies conducted on laboratory strains of Bacillus subtilis, Escherichia coli JM109, Pseudomonas putida and Saccharomyces cerevisiae.

Non-NASA Center↗
NASA NTRSDate not supplied

Vapor Growth of ZnSe-based Compound Semiconductors in Low Gravity environment on International Space Station

Six low-gravity experiments of crystal growth of ZnSe-based compound semiconductors by physical vapor transport (PVT) are scheduled to be processed inside Low Gradient Furnace (LGF) on Materials Science Research Rack (MSRR) in Material Science Laboratory (MSL) on International Space Station (ISS) starting from February 27th, 2023. The flight samples consist of three unseeded growths of ZnSe, Cr-doped ZnSe and ZnSeTe as well as growth of ZnSe on top of three single crystal ZnSe seeds with orientation of (100), (110) and (111). The main scientific objective of the flight experiments is to establish the effects of gravity-driven fluid flows on Earth to (i) the non-uniform incorporation of impurities and defects and (ii) the deviation from stoichiometry observed in the grown crystals caused by growth interface fluctuations and evolution as the results of buoyancy-driven convection and irregular fluid-flows. The investigation will also evaluate the effects of gravity on the PVT process by examining (i) the growth kinetics on various seed orientations, (ii) the dopant segregation and distribution in the Cr doped ZnSe, and (iii) the compositional segregation and distribution in the ternary ZnSeTe crystal grown by PVT. Additionally, from X-ray diffraction/topography and by measuring the density and distribution of various crystalline defects it will assess the self-induced strain effects developed during processing at elevated temperatures caused by the weight of the crystals. During the ground-based preparations, the following tasks have been accomplished: • Establish the thermodynamic properties of the ZnSe binary and related ternary systems by measuring the partial pressures of the vapor phase coexisting with the condensed phases and perform phase diagram calculations for comparison. • Evaluate the effects of gravity-driven convection in the growth process by performing experiments under various vapor transport orientations relative to gravity vector, i.e., horizontal, vertically stabilized and destabilized configurations. • Establish a quantitative correlation between growth parameters and the characteristics and properties of the grown crystals by characterizing the crystals using a variety of techniques. • Evaluate the fundamentals of the current theories on vapor transport and crystal growth kinetics by performing in-situ and real-time monitoring techniques during growth. • Develop a one-dimensional diffusion limited theory of mass transport and two- and three-dimensional numerical simulations of both mass transport and heat transfer and thus establish the microgravity requirements for the flight experiments and to predict the characteristics of the grown crystals such as the solid-vapor interface shapes… etc. and establish a fundamental understanding of the crystal growth process from the comparison with experimental results.

ZnSe↗
DOE OSTI2026

High-Burnup LOCA Burst Susceptibility BISON Analysis in PWRs and BWRs

Accurately assessing high-burnup fuel behavior during loss-of-coolant accidents (LOCAs) is essential for understanding fuel fragmentation, relocation, and dispersal (FFRD) risks across the US light-water reactor fleet. This work updates previous Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program multiphysics LOCA analyses for a pressurized water reactor (PWR) and a boiling water reactor (BWR) by incorporating recent model and material property advancements in the BISON fuel performance code, including a high-burnup structure (HBS) model, revised cladding burst criteria, and updated thermal–mechanical correlations. This update was needed to support ongoing industry initiatives and upcoming regulatory changes. Full-core, rod-resolved operating histories generated using Virtual Environment for Reactor Analysis (VERA) and system-level LOCA conditions obtained from TRACE were applied to statistically representative rod samples in BISON to evaluate burst behavior and FFRD susceptibility. These calculations used two cladding burst correlations and three fuel pulverization models so that the predictions of these models could be compared. The updated PWR simulations show markedly improved numerical stability as the number of crashed simulations decreased by 95% compared to the previous study, and hence higher confidence in results. The updated PWR simulations predicted cladding bursts exclusively among once-burned, high-power rods, with two different cladding burst models identifying the same burst-susceptible population. Resulting FFRD susceptibility estimates are significantly reduced compared with earlier studies, driven by cooler predicted fuel and plenum temperatures, lower hoop strains, and reduced fission gas release in the updated models. In contrast, none of the BWR rods were predicted to burst under either burst criterion, reaffirming minimal BWR FFRD susceptibility even with updated HBS and material models. Comparisons between the PWR and BWR end-of-cycle predictions are made. Comparison with prior work highlights significant shifts in PWR fuel performance metrics and confirmation of earlier BWR conclusions. Overall, the updated results underscore the importance of having high-resolution detailed modeling capability and continuously integrating evolving material models and physics into high-resolution multiphysics simulations. The unified assessment presented here strengthens confidence in predicting high-burnup LOCA behavior by improving agreement between different cladding burst correlations. These results also provide an improved foundation for future BISON model development, FFRD susceptibility calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Ultrafast photo-induced phonon hardening due to Pauli blocking in MAPbI 3 single-crystal and polycrystalline perovskites

Metal-halide perovskite semiconductors have attracted intense interest over the past decade, particularly for applications in photovoltaics. Low-energy optical phonons combined with significant crystal anharmonicity play an important role in charge-carrier cooling and scattering in these materials, strongly affecting their optoelectronic properties. We have observed optical phonons associated with Pb–I stretching in both MAPbI 3 single crystals and polycrystalline thin films as a function of temperature by measuring their terahertz conductivity spectra with and without photoexcitation. An anomalous bond hardening was observed under above-bandgap illumination for both single-crystal and polycrystalline MAPbI 3 . First-principles calculations reproduced this photo-induced bond hardening and identified a related lattice contraction (photostriction), with the mechanism revealed as Pauli blocking. For single-crystal MAPbI 3 , phonon lifetimes were significantly longer and phonon frequencies shifted less with temperature, compared with polycrystalline MAPbI 3 . We attribute these differences to increased crystalline disorder, associated with grain boundaries and strain in the polycrystalline MAPbI 3 . Thus we provide fundamental insight into the photoexcitation and electron–phonon coupling in MAPbI 3 .

36 MATERIALS SCIENCE↗
NASA NTRS1980

Effect of material nonhomogeneity on crack propagation characteristics

The influence of material nonhomogeneity on the behavior of a moving crack is investigated. The model assumes a running crack in a material whose elastic properties may differ from those of the surrounding material. Theoretical calculations showed that the energy stored in elements ahead of the crack can be raised or lowered depending on the crack velocity, the crack length and the degree of material nonhomogeneity which is associated with the ratio of the shear moduli and the distance between the crack and the neighboring material with different elastic properties. Based on the strain energy density theory, predictions are made on how material nonhomogeneity can influence the initiation and/or arrest characteristics of cracks.

Sih, G. C.↗
NASA NTRS1987

Mechanisms of rolling contact spalling

The results of a study aimed at analyzing the mechanical material interactions responsible for rolling contact spalling of the 440 C steel, high pressure oxygen turbopump bearings are presented. A coupled temperature displacement finite element analysis of the effects of friction heating under the contact is presented. The contact is modelled as a stationary, heat generating, 2 dimensional indent in an elastic perfectly plastic half-space with heat fluxes up to 8.6 x 10000 KW/m sq comparable to those generated in the bearing. Local temperatures in excess of 1000 C are treated. The calculations reveal high levels of residual tension after the contact is unloaded and cools. Efforts to promote Mode 2/Mode 3 fatigue crack growth under cyclic torsion in hardened 440 C steel are described. Spalls produced on 440 C steel by a 3 ball/rod rolling contact testing machine were studied with scanning microscopy. The shapes of the cyclic, stress strain hysteresis loops displayed by hardened 440 C steel in cyclic torsion at room temperature are defined for the plastic strain amplitudes encountered in rolling/sliding contact. Results of these analyses are discussed in detail.

Kumar, A. M.↗
NASA NTRS1993

Temperature induced distortions in space telescope mirrors

In this paper, it is illustrated how measured instantaneous coefficients of thermal expansion (CTE) can be accurately taken into account when modeling the structural behavior of space based optical systems. In particular, the importance of including CTE spatial variations in the analysis of optical elements is emphasized. A comparison is made between the CTE's of three optical materials commonly used in the construction of space mirrors (ULE, Zerodur, and beryllium). The overall impact that selection of any one of these materials has on thermal distortions is briefly discussed. As an example of how temperature dependent spatial variations in thermal strain can be accurately incorporated in the thermo-structural analysis of a precision optical system, a finite element model is developed, which is used to estimate the thermally induced distortions in the Hubble Space Telescope's (HST) primary mirror. In addition to the structural analysis, the optical aberrations due to thermally induced distortions are also examined. These calculations indicate that thermal distortions in HST's primary mirror contribute mainly to defocus error with a relatively small contribution to spherical aberration.

Nied, H. F.↗
NASA NTRS1937

Strain measurements on small duralumin box beams in bending

Extensive strain-gage measurements were made chiefly on the tension side of five small rectangular box beams constructed of sheet duralumin. The main conclusion was that within the test range the tension cover may be considered as being fully effective but that at any given point on the beam there may be unaccountable variations of 5 percent from the calculated stresses on the thicker sheets tested (0.044 and 0.023 inch) and of 10 percent or more on the thinnest sheet tested (0.014 inch).

Kuhn, Paul↗
NASA NTRS2019

Sidewall Effects on Exact Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary Layer Interaction

Large-eddy simulations are performed using wall-resolved mesh for a Mach 2.29 impinging shock wave/boundary-layer interaction. Flow conditions are based on an experiment and therefore entire span was simulated, including the two sidewalls. Mean flow comparison with the experimental data showed that the interaction was larger in the simulation. Time-series analysis of a rake of pressure probes immediately downstream of the mean reflected shock position showed a peak in weighted power spectral density occurred about $St_{Lint}=0.01$, owing to a larger interaction length. Budgets of Reynolds-stress transport calculated across the span and along the corner bisector showed high degree of anisotropy. Merging of the secondary flows and separation along the corner gives rise to unstablecounter rotating vortices, which straddle the corner and grow in size. This also leads to a development of new behavior in the viscous sublayer along the corner bisector, where the pressure strain andmolecular diffusion mechanisms become prominent.

Vyas, Manan A.↗