Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain gradient modeling”

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 163 records · Page 9

Residual Stress in Cold Spray SS304L Measured Via Neutron Diffraction and Comparison of Analytical Models to Predict the Residual Stress

Here, this study employs neutron diffraction to investigate the relationship between residual stress and coating thickness in cold sprayed 304L austenitic stainless steel. Results show that shot peening predominantly impacts the residual stress profile, leading to substantial in-plane compressive force. The impact of laser heating, a widely used method to alter cold spray's microstructural properties, on the coating's residual stress is also analyzed. The findings indicate that the maximum compressive residual stress in the in-plane component is mainly independent of coating thickness, which suggests that the material properties determine the maximum residual stress. The cold sprayed deposits possessed compressive, nearly biaxial strain and stresses. After laser heating, these stresses were replaced by tensile residual stresses. Two analytical models, the Tsui and Clyne and the Boruah models, for predicting residual stresses are also evaluated, and both models provide reasonable fits to the experimental data. At this point, the deviations between the experimental results and the models are principally caused by the inability of the current models to address plastic deformation and relaxation, and the residual stresses generated by thermal gradients.

36 MATERIALS SCIENCE↗

High-throughput screening of multifunctional nanocoatings based on combinations of polyphenols and catecholamines

Biomimetic surface coatings based on plant polyphenols and catecholamines have been used broadly in a variety of applications. However, the lack of a rational cost-effective platform for screening these coatings and their properties limits the true potential of these functional materials to be unleashed. Here, we investigated the oxidation behavior and coating formation ability of a library consisting of 45 phenolic compounds and catecholamines. UV–vis spectroscopy demonstrated significant acceleration of oxidation and polymerization under UV irradiation. We discovered that several binary mixtures resulted in non-additive behavior (synergistic or antagonistic effect) yielding much thicker or thinner coatings than individual compounds measured by ellipsometry. To investigate the properties of coatings derived from new combinations, we used a miniaturized high-throughput strategy to screen 2,532 spots coated with single, binary, and ternary combinations of coating precursors in one run. We evaluated the use of machine learning models to learn the relation between the chemical structure of the precursors and the thickness of the nanocoatings. Formation and stability of nanocoatings were investigated in a high-throughput manner via discontinuous dewetting. 30 stable combinations (hits) were used to tune the surface wettability and to form water droplet microarray and spot size gradients of water droplets on the coated surface. No toxicity was observed against eukaryotic HeLa cells and Pseudomonas aeruginosa (strain PA30) bacteria after 24 h incubation at 37 C. The strategy introduced here for high-throughput screening of nanocoatings derived from combinations of coating precursors enables the discovery of new functional materials for various applications in science and technology in a cost-effective miniaturized manner.

36 MATERIALS SCIENCE↗

HOST structural analysis program overview

Hot-section components of aircraft gas turbine engines are subjected to severe thermal structural loading conditions, especially during the startup and takeoff portions of the engine cycle. The most severe and damaging stresses and strains are those induced by the steep thermal gradients induced during the startup transient. These transient stresses and strains are also the most difficult to predict, in part because the temperature gradients and distributions are not well known or readily predictable and, in part, because the cyclic elastic-viscoplastic behavior of the materials at these extremes of temperature and strain are not well known or readily predictable. A broad spectrum of structures related technology programs is underway to address these deficiencies at the basic as well as the applied level. The three key program elements in the HOST structural analysis program are computations, constitutive modeling, and experiments for each research activity. Also shown are tables summarizing each of the activities.

Thompson, Robert L.↗

Thermal-Structural Analysis of Sunshield Membranes

Future large infrared space telescopes, such as the James Webb Space Telescope (JWST), will require deployable sunshields to provide passive cooling for optics and instruments. Deployable sunshield structures for such applications typically consist of multiple thin-film membrane layers supported by deployable booms. The mechanical design of the sunshield must accommodate thermal strains due to layer-to-layer temperature differences as well as potentially large in-plane temperature gradients within individual film layers. This paper describes a thermal-structural analysis for predicting the stress state in a thin-film membrane subject to both mechanical thermal loads that could aid in the mechanical design of future sunshield structures. First the temperature field predicted by a thermal analysis is mapped to a structural finite element model, and then the structural response is predicted using a nonlinear static analysis. The structural model uses membrane elements in conjunction with a tension field material model to predict the response of the thin-film membrane layer. The tension field material model accounts for no-compression behavior associated with wrinkling and slackness. This approach was used to study the problem of a single membrane layer from the NASA reference concept for the JWST sunshield. Results from the analysis show that the membrane can experience a loss of tensile preload due to the presence of an in-plane temperature gradient representative of the cold-side layer temperature distribution predicted for the reference concept JWST.

Johnston, John↗

The response of a compressible turbulent boundary layer to short regions of concave surface curvature

An experimental investigation of the supersonic flow of a boundary layer over short regions of concave surface curvature, using radii of curvature equivalent to 10 and 50 (models I and II, respectively) initial boundary layer thicknesses and 8-deg turning angle, notes that the boundary layer experienced a strong adverse pressure gradient as the boundary layer passed through the curved region. Attention is given to the mean and turbulent field for each flow, with comparisons of results to those for an 8-deg ramp. In the ramp and model I, turbulence levels increased significantly, together with structural parameters; in model II, despite an increase in turbulence levels, structural parameters remained essentially unchanged. The former results suggest that the perturbation is rapid, and describable in terms of total rather than local strains.

Jayaram, Mohan↗

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗

A hereditary integral, transient network approach to modeling permanent set and viscoelastic response in polymers

An efficient numerical framework is presented for modeling viscoelasticity and permanent set of polymers. It is based on the hereditary integral form of transient network theory, in which polymer chains belong to distinct networks each with different natural equilibrium states. Chains continually detach from previously formed networks and reattach to new networks in a state of zero stress. The free energy of these networks is given in terms of the deformation gradient relative to the configuration at which the network was born . A decomposition of the kernel for various free energies allows for a recurrence relationship to be established, bypassing the need to integrate over all time history. The technique is established for both highly compressible and nearly incompressible materials through the use of neo-Hookean, Blatz–Ko, Yeoh, and Ogden-Hill material models. Multiple examples are presented showing the ability to handle rate-dependent response and residual strains under complex loading histories.

Finite element method↗

Determination of Creep Behavior of Thermal Barrier Coatings Under Laser Imposed High Thermal and Stress Gradient Conditions

A laser sintering/creep technique has been established to determine the creep behavior of thermal barrier coatings under steady-state high heat flux conditions. For a plasma sprayed zirconia-8 wt. % yttria coating, a significant primary creep strain and a low apparent creep activation energy were observed. Possible creep mechanisms involved include stress induced mechanical sliding and temperature and stress enhanced cation diffusion through the splat and grain boundaries. The elastic modulus evolution, stress response, and total accumulated creep strain variation across the ceramic coating are simulated using a finite difference approach. The modeled creep response is consistent with experimental observations.

Zhu, Dongming↗

Temporally continuous thermofluidic–thermomechanical modeling framework for metal additive manufacturing

Additive manufacturing (AM) is known to generate large magnitudes of residual stresses (RS) within builds due to steep and localized thermal gradients. In the current state of commercial AM technology, manufacturers generally perform heat treatments in effort to reduce the generated RS and its detrimental effects on part distortion and in-service failure. Computational models that effectively simulate the deposition process can provide valuable insights to improve RS distributions. Accordingly, it is common to employ Computational fluid dynamics (CFD) models or finite element (FE) models. While CFD can predict geometric and thermal-fluid behavior, it cannot predict the structural response (e.g., stress–strain) behavior. On the other hand, an FE model can predict mechanical behavior, but it lacks the ability to predict geometric and fluid behavior. Thus, an effectively integrated thermofluidic–thermomechanical modeling framework that exploits the benefits of both techniques while avoiding their respective limitations can offer valuable predictive capability for AM processes. In contrast to previously published efforts, the work herein describes a one-way coupled CFD-FEA framework that abandons major simplifying assumptions, such as geometric steady-state conditions, the absence of material plasticity, and the lack of detailed RS evolution/accumulation during deposition, as well as insufficient validation of results. Here, the presented framework is demonstrated for a directed energy deposition (DED) process, and experiments are performed to validate the predicted geometry and RS profile. Both single- and double-layer stainless steel 316L builds are considered. Geometric data is acquired via 3D optical surface scans and X-ray micro-computed tomography, and residual stress is measured using neutron diffraction (ND). Comparisons between the simulations and measurements reveal that the described CFD-FEA framework is effective in capturing the coupled thermomechanical and thermofluidic behaviors of the DED process. The methodology presented is extensible to other metal AM processes, including power bed fusion and wire-feed-based AM.

42 ENGINEERING↗

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗

The evolution of tectonic features on Ganymede

The bands of bright resurfaced terrain on Ganymede are probably broad grabens formed by global expansion and filled with deposits of ice. Grooves within the bands are thought to be extensional features formed during the same episode of expansion. The crust of Ganymede is modeled as a viscoelastic material subjected to extensional strain. With sufficiently high strain rates and stresses, deep normal faulting will occur, creating broad grabens that may then be filled. Continuing deformation at high strain rates and stresses will cause propagation of deep faults up into the flood deposits and normal faulting at the surface, while lower strain rates and stresses will cause formation of open extension fractures or, if the crustal strength is very low, grabens at the surface. The spacing between adjacent fractures may reflect the geothermal gradient at the time of deformation. Surface topography resulting from fracturing and normal faulting will decay with time as a result of viscous relaxation and mass-wasting.

Squyres, S. W.↗

Nonisothermal elasto-visco-plastic response of shell-type structures

A mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic or static thermomechanical loads is disussed. Among the system responses, which are associated with these load conditions, are thermal buckling and creep buckling. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and have been considered in the development of the mathematical model. Furthermore, this was accommodated in the solution procedures. A complete true ab-inito rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformation, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of the metric and curvature in two coordinate systems, a fixed (spatial) one and a convected (material) coordinate system. The relations between the time derivative and the covariant derivatives (gradients) have been developed for curved space and motion, so that the velocity components supply the connection between the equations of motion and the time rate of change of the metric and curvature tensors.

Simitses, G. J.↗

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗

Utah FORGE 2439: A Multi-Component Approach to Characterizing In-Situ Stress

Core-based in-situ stress estimation, Triaxial Ultrasonic Velocity (labTUV) data, and Deformation Rate Analysis (DRA) data for Utah FORGE well 16A(78)-32 using triaxial ultrasonic velocity and deformation rate analysis. Report documenting a multi-component approach to characterizing in-situ stress at the U.S. DOE FORGE EGS site: laboratory, modeling and field measurement. Core-based methods for in-situ stress estimation were applied using samples from 5 intervals within the Utah FORGE 16A(78)-32 well. At three of these locations, Triaxial Ultrasonic Velocity (labTUV) tests were performed, resulting in experimentally-determined relationships between wave velocities and stresses. Non-monotonic increase in the velocity-stress relationships are inferred provide evidence of stress history and are therefore used to estimate in-situ stress magnitudes. Additionally, Deformation Rate Analysis (DRA) tests were run on core plugs from various orientations at each of the 5 sampling locations. These, too, provide evidence of stress history based on stress-strain behavior. A novel Weight of Evidence (WoE) method was developed as a means of synthesizing in-situ stress evidence from these two types of tests. Results indicate the minimum horizontal stress gradient ranges from 0.58 psi/ft to 0.69 psi/ft, with 4 of the 5 values between 0.66 psi/ft and 0.69 psi/ft. The vertical stress gradient ranges from 1.05 psi/ft to 1.12 psi/ft, with 4 of the 5 zones given results between 1.09 psi/ft and 1.12 psi/ft. The maximum horizontal stress gradient ranges from 0.98 psi/ft to 1.34 psi/ft, with 4 of the 5 zones falling between 0.98 psi/ft and 1.24 psi/ft. The stress regime thus appears to be on the edge between normal faulting and strike-slip faulting, potentially flipping back and forth between the two regimes due to variability of rock properties, structures such as faults, and/or thermal anomalies.

15 GEOTHERMAL ENERGY↗

Thermal and Pressure Characterization of a Wind Tunnel Force Balance Using the Single Vector System

Wind tunnel research at NASA Langley Research Center s 31-inch Mach 10 hypersonic facility utilized a 5-component force balance, which provided a pressurized flow-thru capability to the test article. The goal of the research was to determine the interaction effects between the free-stream flow and the exit flow from the reaction control system on the Mars Science Laboratory aeroshell during planetary entry. In the wind tunnel, the balance was exposed to aerodynamic forces and moments, steady-state and transient thermal gradients, and various internal balance cavity pressures. Historically, these effects on force measurement accuracy have not been fully characterized due to limitations in the calibration apparatus. A statistically designed experiment was developed to adequately characterize the behavior of the balance over the expected wind tunnel operating ranges (forces/moments, temperatures, and pressures). The experimental design was based on a Taylor-series expansion in the seven factors for the mathematical models. Model inversion was required to calculate the aerodynamic forces and moments as a function of the strain-gage readings. Details regarding transducer on-board compensation techniques, experimental design development, mathematical modeling, and wind tunnel data reduction are included in this paper.

Lynn, Keith C.↗

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗