Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive behavior”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Design Curve Generation for 3D SiC Fiber Architecture

The design tool provides design curves that allow a simple and quick way to examine multiple factors that can influence the processing and key properties of the preforms and their final SiC-reinforced ceramic composites without over obligating financial capital for the fabricating of materials. Tool predictions for process and fiber fraction properties have been validated for a HNS 3D preform.The virtualization aspect of the tool will be used to provide a quick generation of solid models with actual fiber paths for finite element evaluation to predict mechanical and thermal properties of proposed composites as well as mechanical displacement behavior due to creep and stress relaxation to study load sharing characteristic between constitutes for better performance.Tool predictions for the fiber controlled properties of the SiCSiC CMC fabricated from the HNS preforms will be valuated and up-graded from the measurements on these CMC

matrix materials↗
DOE OSTI2024

Peridynamic neural operators: A data-driven nonlocal constitutive model for complex material responses

Neural operators, which can act as implicit solution operators of hidden governing equations, have recently become popular tools for learning the responses of complex real-world physical systems. Nevertheless, most neural operator applications have thus far been data-driven and neglect the intrinsic preservation of fundamental physical laws in data. Here, in this work, we introduce a novel integral neural operator architecture called the Peridynamic Neural Operator (PNO) that learns a nonlocal constitutive law from data. This neural operator provides a forward model in the form of state-based peridynamics, with objectivity and momentum balance laws automatically guaranteed. As applications, we demonstrate the expressivity and efficacy of our model in learning complex material behaviors from both synthetic and experimental data sets. We also compare the performances with baseline models that use predefined constitutive laws. We show that, owing to its ability to capture complex responses, our learned neural operator achieves improved accuracy and efficiency. Moreover, by preserving the essential physical laws within the neural network architecture, the PNO is robust in treating noisy data. The method shows generalizability to different domain configurations, external loadings, and discretizations.

42 ENGINEERING↗
NASA NTRS1977

Application of the endochronic theory of viscoplasticity to solid propellants and sandasphalt concrete

Solid propellants, sand-asphalt concrete and hard plastics showed rate sensitive mechanical behavior which, in addition, indicated that these materials have a permanent memory of the strain (or loading) path by which their present state was attained. A constitutive equation was formulated in general three dimensional tensorial form by means of irreversible thermodynamics. By using a very simple analytical form, it was shown that the mechanical behavior of solid propellants and sand-asphalt concrete can be readily described.

Peng, S. T. J.↗
NASA NTRSDate not supplied

Challenges in Scaling Analysis from Coupons to Elements in 3D Woven Composites for a Launch Vehicle End-Ring Structure

The limitations of structural design with 3D woven composites based on homogenized material models characterized using coupon test data, representative unit cell analysis, or a combination of the two are discussed in the context of a space launch system (SLS)-scale end ring. Selected building block test configurations at the coupon, element, and subcomponent scales are described in detail. Highly instrumented tests and X-ray computed tomography inspections are used to characterize the structural behavior and associated damage progression with increasing load. It is shown that a single homogenized material constitutive relationship, applied at coupon and element scale, is not sufficient to model the structural response observed in both tests. In fact, substantially different stress-strain laws are required to predict the response of the coupon and element scale specimens, despite similar loading conditions. A model with discrete representation of the binder tows and cohesive interfaces to model matrix cracking shows much better agreement with the test data at the two length scales. Thus, evidence is provided showing that predictive structural analysis of 3D woven composites requires discrete representation of the tows, matrix cracking morphology, and their interactions in the structural configuration. With accurate knowledge of the element scale behavior, engineering approximations readily enable analysis scale-up to the subcomponent level.

Andrew C Bergan↗
DOE OSTI2020

Distonic radical anion species in cysteine oxidation processes

Oxidation of cysteine residues constitutes an important regulatory mechanism in the function of biological systems. Much of this behavior is controlled by the specific chemical properties of the thiol side-chain group, where reac-tions with reactive oxygen species take place. Herein, we in-vestigated the entire cysteine oxidation cycle Cys-SH ? Cys-SOnH (n = 1, 2, 3) using cryogenic negative ion photoelectron spectroscopy and quantum-chemical calculations. The conven-tional view on the first reversible oxidation step (n = 1) is as-sociated with sulfinate species. Yet our results indicate that an alternative option exists in the form of novel distonic radical anion •OS-CH2CH(NH2)-COO- with an unpaired electron on the thiol group and excess negative charge on the carboxylate. Higher order oxidation states (n = 2, 3) are thought to be asso-ciated with irreversible oxidative damage, and our results show that excess negative charge in those cases migrates to -SOn- group. Furthermore, these species are stable towards 1e oxidation, as opposed to the n = 1 case that undergoes intra-molecular proton transfer. The molecular level insights report-ed in this work provide direct spectroscopic evidences towards unique chemical versatility of Cys-sulfenic acid (Cys-SOH) in post-translational modifications of protein systems. This research was supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences (X.-B. W. and M. V.) and was performed at the EMSL, a national scientific user facility sponsored by DOE's Office of Biological and Environmental Research and located at the Pacific Northwest National Laboratory. Z. B. Q. acknowledges additional support from the National Natural Science Foundation of China (Grant No. 21873003 and 21503003) and the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase) (Grant No. U1501501) and Super Computation of Shenzhen.

Qin, Zhengbo↗
DOE OSTI2021

Grid-Scale Ternary-Pumped Thermal Electricity Storage for Flexible Operation of Nuclear Power Generation under High Penetration of Renewable Energy Sources

In this work, the integration of a grid-scale ternary-Pumped Thermal Electricity Storage (t-PTES) with a nuclear power generation to enhance operation flexibility is assessed using physics-based models and digital real time simulation. A part of the electricity from the nuclear power generation is delivered to the grid, and the balance is used to power a heat pump that can be augmented by an auxiliary resistive load element to increase the charging rate of the thermal storage. This increases the thermal potential between hot and cold thermal stores (usually solid materials or molten salts inside large storage tanks). The thermal energy is transformed back into electricity by reversing the heat pump cycle. Different transient scenarios including startup, shutdown, and power change for grid-connected operation are simulated to determine the behavior of the hybrid nuclear-t-PTES system operating under variable loads that constitute a departure from conventional, baseload nuclear plant operation schemes. Ternary refers to the three modes operation: (i) heat pump (including heating coil), (ii) heat engine, and (iii) simultaneous operation of heat pump (including heating coil) and heat engine during changeover from pumping to generation or vice-versa. The controllability of t-PTES in the short timescales as a dynamic load is used to demonstrate operational flexibility of hybrid nuclear plants for flexible operation through advanced load management. The integration of t-PTES into nuclear power systems enhances the system flexibility and is an enabler for high penetration of renewable energy resources.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1983

An uncoupled viscoplastic constitutive model for metals at elevated temperature

An uncoupled constitutive model for predicting the transient response of thermal and rate dependent, inelastic material behavior is presented. The uncoupled model assumes that there is a temperature below which the total strain consists essentially of elastic and rate insensitive inelastic strains only. Above this temperature, the rate dependent inelastic strain (creep) dominates. The rate insensitive inelastic strain component is modeled in an incremental form with a yield function, flow rule and hardening law. Revisions to the hardening rule permit the model to predict temperature-dependent kinematic-isotropic hardening behavior, cyclic saturation, asymmetric stress-strain response upon stress reversal, and variable Bauschinger effect. The rate dependent inelastic strain component is modeled using a rate equation in terms of back stress, drag stress and exponent n as functions of temperature and strain. A sequence of hysteresis loops and relaxation tests are utilized to define the rate dependent inelastic strain rate. Evaluation of the model is performed by comparison with experiments involving various thermal and mechanical load histories on 5086 aluminum alloy, 304 stainless steel and Hastelloy-X.

Haisler, W. E.↗
NASA NTRS1987

Three-dimensional inelastic analysis methods for hot section components

The objective was to develop analytical methods capable of evaluating the cyclic time dependent inelasticity which occurs in hot section engine components. Because of the large excursions in temperature associated with hot section engine components, the techniques developed must be able to accommodate large variations in material behavior including plasticity and creep. To meet this objective, a matrix consisting of three constitutive mode-element formulations was developed. A separate program for each combination of constitutive mode-element model was written. The source codes of the nine programs range in size from 7300 lines for the Bodner/twenty code to 19,000 lines for the Haisler and Allen/nine code. All the codes were given a stand alone capability of performing cyclic nonlinear analysis.

Lackney, Joseph↗
NASA NTRS1991

Thermal barrier coating life prediction model development, phase 2

The objective of this program was to generate a life prediction model for electron-beam-physical vapor deposited (EB-PVD) zirconia thermal barrier coating (TBC) on gas turbine engine components. Specific activities involved in development of the EB-PVD life prediction model included measurement of EB-PVD ceramic physical and mechanical properties and adherence strength, measurement of the thermally grown oxide (TGO) growth kinetics, generation of quantitative cyclic thermal spallation life data, and development of a spallation life prediction model. Life data useful for model development was obtained by exposing instrumented, EB-PVD ceramic coated cylindrical specimens in a jet fueled burner rig. Monotonic compression and tensile mechanical tests and physical property tests were conducted to obtain the EB-PVD ceramic behavior required for burner rig specimen analysis. As part of that effort, a nonlinear constitutive model was developed for the EB-PVD ceramic. Spallation failure of the EB-PVD TBC system consistently occurred at the TGO-metal interface. Calculated out-of-plane stresses were a small fraction of that required to statically fail the TGO. Thus, EB-PVD spallation was attributed to the interfacial cracking caused by in-plane TGO strains. Since TGO mechanical properties were not measured in this program, calculation of the burner rig specimen TGO in-plane strains was performed by using alumina properties. A life model based on maximum in-plane TGO tensile mechanical strain and TGO thickness correlated the burner rig specimen EB-PVD ceramic spallation lives within a factor of about plus or minus 2X.

Meier, Susan Manning↗
NASA NTRS2003

Strength Differential Measured in Inconel 718: Effects of Hydrostatic Pressure Studied

Aeropropulsion components, such as disks, blades, and shafts, are commonly subjected to multiaxial stress states at elevated temperatures. Experimental results from loadings as complex as those experienced in service are needed to help guide the development of accurate viscoplastic, multiaxial deformation models that can be used to improve the design of these components. During a recent study on multiaxial deformation (ref. 1) on a common aerospace material, Inconel 718, it was shown that the material in the aged state exhibits a strength differential effect (SDE), whereby the uniaxial compressive yield and subsequent flow behavior are significantly higher than those in uniaxial tension. Thus, this material cannot be described by a standard von Mises yield formulation. There have been other formulations postulated (ref. 2) that involve other combinations of the stress invariants, including the effect of hydrostatic stress. The question remained as to which invariants are necessary in the flow model. To capture the physical mechanisms occurring during deformation and reflect them in the plasticity formulation, researchers examined the flow of Inconel 718 under various amounts of hydrostatic stress to determine whether or not hydrostatic stress is needed in the formulation. Under NASA Grant NCC3-464, monitored by the NASA Glenn Research Center, a series of tensile tests were conducted at Case Western Reserve University on aged (precipitation hardened) Inconel 718 at 650 C and with superimposed hydrostatic pressure. Dogbone shaped tensile specimens (3-mm-diameter gauge by 16-mm gauge length) and cylindrical compression specimens (3-mm-diameter gauge by 6-mm gauge length) were strain gauged and loaded in a high-pressure testing apparatus. Hydrostatic pressures were obtained with argon and ranged from 210 to 630 MPa. The aged Inconel 718 showed a pronounced difference in the tension and compression yield strength (i.e., an SDE), as previously observed. Also, there were no significant effects of hydrostatic pressure on either the tensile and compressive yield strength (see the graph) or on the magnitude of the SDE. This behavior is not consistent with the pressure-dependent theory of the SDE, which postulates that the SDE is associated with pressure-dependent and/or internal friction dependent deformation associated with non-Schmid effects at the crystal level (refs. 3 and 4). Flow in Inconel 718 appears to be independent of hydrostatic pressure, suggesting that this invariant may be removed from the phenomenological constitutive model. As part of an ongoing effort to develop advanced constitutive models, Glenn s Life Prediction Branch coordinated this work with that of research on the multiaxial deformation behavior of Inconel 718 being conducted at Pennsylvania State University under NASA Grant NCC597.

Lewandowski, John J.↗
NASA NTRS1990

Plasticity analysis of fibrous composite laminates under thermomechanical loads

Elastic-plastic behavior of fibrous composite laminates is analyzed for coupled in-plane mechanical loads and uniform thermal changes. Constitutive equations of the individual fibrous layers are derived from a vanishing fiber diameter model that represents the essential axial constraint between the phases. This permits derivation of closed form equations for the overall yield condition, stress concentration factors, and instantaneous compliance. Thermoelastic properties of the phases and yield stress of the matrix phase are functions of temperature. The effect of the model assumptions on the predicted behavior of composite laminates is examined by comparing the calculated response under cyclic thermal changes to available theoretical results and experimental measurements.

Bahei-El-din, Yehia A.↗
DOE OSTI2021

Formulation and calibration of two-dimensional constitutive models for composite structures based on panel tests

Steel Plate Concrete (SC) composite members have been widely adopted because of its cost-efficiency and enhanced structural behavior. While researchers have attempted to study its in-plane shear behavior in the past twenty years, very limited number of large-scale pure shear tests were performed due to the challenge of experimental set-up and the availability of facilities. In this paper, a series of uniaxial loading tests and two full-scale pure shear panel tests of SC members were reported, on which the “mechanics-based Membrane Model of SC elements (MM-SC)” is developed. The MM-SC model is based on the fixed-angle crack formulation and the smeared-crack formulation, in which the experimental-based uniaxial constitutive models are implemented, considering the local buckling of faceplate, the tension stiffening of steel plate, the strength degradation of cracked concrete and the confinement effect of concrete. The proposed MM-SC model is subsequently incorporated into the object-oriented software OpenSEES. Finally, the simulation results of proposed model well predict the SC test observations in terms of critical branch points and structural behaviors, including initial stiffness, cracking strength, post-crack stiffness, yield strength, maximum strength, and failure modes.

42 ENGINEERING↗
NASA NTRS1996

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗
NASA NTRS1987

A viscoplastic constitutive theory for metal matrix composites at high temperature

A viscoplastic theory is presented for representing the high-temperature deformation behavior of metal matrix composites. The point of view taken is a continuum one where the composite is considered a material in its own right, with its own properties that can be determined for the composite as a whole. It is presumed that a single preferential (fiber) direction is identifiable at each material point (continuum element) admitting the idealization of local transverse isotropy. A key ingredient in this work is the specification of an experimental program for the complete determination of the material functions and parameters for characterizing a particular metal matrix composite. The parameters relating to the strength of anisotropy can be determined through tension/torsion tests on longitudinally and circumferentially reinforced thin-walled tubes. Fundamental aspects of the theory are explored through a geometric interpretation of some basic features analogous to those of the classical theory of plasticity.

Robinson, D. N.↗
NASA NTRS1987

Deformation and fatigue behavior of SSME turbopump blade materials

Directionally solidified and single crystal superalloys which are intended for use as turbopump blade materials are anisotropic both elastically and plastically. Therefore, isotropic constitutive models must be modified. Several models which are now being developed are based on metallurgical theories of deformation in these types of alloys. However, these theories have not been fully justified, and the temperature and strain regimes over which they may be valid are poorly defined. The objective of this work is to study the deformation behavior of the alloys, in order to determine the validity of these models and to thereby support the ongoing research efforts in solid mechanics.

Milligan, Walter W.↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enhancing predictive understanding and accuracy in geological carbon dioxide storage monitoring: Simulation and history matching of tracer transport dynamics

Co-injection of conservative tracers with carbon dioxide (CO 2 ) is a viable tool for monitoring subsurface processes during geological CO 2 storage (GCS). This research investigates the simulation and history-matching of a gas tracer (sulfur hexafluoride, SF 6 ) during CO 2 flooding, employing a core flooding result in Berea sandstone. Four extensively used saturation functions are assessed for their efficacy in history matching of CO 2 /SF 6 injection at the core scale. The history-matching process incorporates particle swarm optimization (PSO) to fine-tune constitutive relationships parameters. Next, employing transport models at the aquifer scale, we interrogate the impact on tracer transport and mixing of saturation function uncertainties, arising from the non-uniqueness of constitutive relationships parameters and saturation function types. To assess the effects of geological heterogeneity on behavior of tracer breakthrough curves (BTCs), we employ two normalized parameters assessing the degree of mixing and SF 6 breakthrough time. The aquifer-scale investigation encompasses both homogeneous and heterogeneous systems with and without capillary heterogeneity effects. Our findings underscore the critical importance of addressing saturation function uncertainties, emphasizing the significance of auxiliary experiments and innovative methodologies to enhance predictive accuracy. The findings highlight significant disparities in arrival times, BTC peaks, tails, and mixing levels, even under optimal conditions. Heterogeneity, with or without capillary heterogeneity, plays a crucial role in shaping BTC variations, resulting in accelerated SF 6 breakthrough times and reduced BTC peaks. Evaluation of monitoring points distant from the injector reveals a dampening effect on the SF6 BTC peak, particularly in heterogeneous systems with capillary heterogeneity, where the peak is halved. These insights underscore the challenges associated with tracer monitoring and the necessity for enhanced methodologies to improve predictive accuracy in subsurface processes during GCS.

58 GEOSCIENCES↗
DOE OSTI2024

An experimental, theoretical, and kinetic modeling study of post-flame oxidation of ammonia

The post -flame oxidation rate of ammonia was investigated in a novel atmospheric pressure flow reactor at temperatures of 1280 ± 16 K and as a function of residence time and mixture composition (1-10% O 2 , dry and moist). The experimental results, as well as selected data from literature, were analyzed using an updated detailed chemical kinetic model. The medium temperature, very lean conditions enhance the importance of reactions of the nitroxyl (HNO) intermediate. High-level theory was used to calculate the rate constant for HNO + NH 2 , indicating that this step is significantly faster than values used in literature. Furthermore, a trajectory based approach was used to determine collision efficiencies for selected bath gases for HNO + M. The experimental results show that the NH 3 oxidation rate increases with temperature and O 2 concentration, while the presence of water vapor slightly inhibits reaction. Formation of NO and N 2 O was strongly promoted at higher levels of O 2 . Modeling results agreed well with the measurements, except at the lowest level of O 2 . The predicted oxidation rate of NH 3 was shown to result from a delicate balance between chain branching and terminating steps involving NH 2 , H 2 NO, and HNO. Recent theoretical work on reactions of these species by Klippenstein and coworkers and Stagni et al. was instrumental in improving modeling predictions. After initiation, NO reached a pseudo -steady-state level, where the pathways to NO were largely balanced by the NH 2 + NO reaction. Nitric oxide was partly oxidized to NO 2 , with the NH 2 + NO 2 reaction responsible for most of the N 2 O formation. Novelty and significance statement: This study provides the first detailed kinetic analysis of the lean postflame oxidation of ammonia, based on time -resolved flow reactor data in a novel reactor. In addition to the post -flame oxidation rate of ammonia, data for formation of NO and N 2 O were compared with modeling predictions. The medium temperature, very lean conditions enhance the importance of reactions of the HNO and H 2 NO intermediates. Inclusion in the model of results from recent high-level theoretical work, including present calculations for HNO + NH 2 and HNO + M, was crucial for capturing the observed behavior. It is argued that the post -initiation steady-state NH 3 oxidation rates constitute important data for model validation, along with ignition delays and laminar flame speeds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗