Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain rate”

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 487 records · Page 27

Experimental deformation of polycrystalline H2O ice at high pressure and low temperature - Preliminary results

A preliminary study is carried out of involving 70 constant strain deformation tests on pure polycrystalline H2O ice under conditions covering most of the stability field of ice I sub h. Brittle failure of I sub h is found to be promoted by lower P, lower T, and higher strain rates. Ductile flow is found to be promoted by higher P, higher T, and lower strain rates. The brittle failure of ice I sub h is found to be most unusual. The fracture strength is a positive function of P only below 50 MPa. At pressures greater than this, the fracture strength is independent of P, and the fracture plane lies approximately 45 deg from the load axis. It is believed that existing extrapolation based on existing experimental data to Ganymede and Callisto may be badly in error.

Durham, W. B.↗

Solid-propellant rocket motor internal ballistic performance variation analysis, phase 2

The Monte Carlo method was used to investigate thrust imbalance and its first time derivative throughtout the burning time of pairs of solid rocket motors firing in parallel. Results obtained compare favorably with Titan 3 C flight performance data. Statistical correlations of the thrust imbalance at various times with corresponding nominal trace slopes suggest several alternative methods of predicting thrust imbalance. The effect of circular-perforated grain deformation on internal ballistics is discussed, and a modified design analysis computer program which permits such an evaluation is presented. Comparisons with SRM firings indicate that grain deformation may account for a portion of the so-called scale factor on burning rate between large motors and strand burners or small ballistic test motors. Thermoelastic effects on burning rate are also investigated. Burning surface temperature is calculated by coupling the solid phase energy equation containing a strain rate term with a model of gas phase combustion zone using the Zeldovich-Novozhilov technique. Comparisons of solutions with and without the strain rate term indicate a small but possibly significant effect of the thermoelastic coupling.

Sforzini, R. H.↗

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

High-rate strength response of tantalum from dynamic hole closure experiments

The science and engineering communities have significant interest in experimental platforms to evaluate and improve models for dynamic material deformation. While well-developed platforms exist, there are still gaps to fill for strain and strain rate conditions accessed during impact and other high-rate loading scenarios. To fill one such gap for strength measurements, a platform was recently developed that accesses high strain rate (≥ 10 5 /s) and large strain (≥ 50%) conditions by measuring the transient closure of a cylindrical hole using in situ x-ray imaging. In the work reported here, further refinement of the platform is performed to reduce the potential effects of porosity and anelasticity on the measurement. This helps us to isolate the strength effects that are the focus of the experiment. The updated experimental configuration employs a two-layer flyer design and elongated target to reduce the magnitude of the tensile excursions associated with rarefaction wave interactions. This allows for a more direct assessment of strength models commonly used for dynamic simulations of metals. Here we apply the new technique to well-characterized tantalum material, allowing for a robust connection to other experimental techniques. Deformation localization can be a concern in large strain experiments, and to help inform future use of the experimental platform, we use simulations with a sub-zone treatment of shear banding to explore potential localization behavior. Here we develop and utilize an experimental configuration with improved isolation of strength effects that can be applied to an expanded range of materials.

42 ENGINEERING↗

Developing a robust strength model using physically-informed genetic programming

The strength of materials is influenced by a range of external conditions, such as temperature and deformation rate. Consequently, materials that demonstrate substantial variations in their mechanical behavior due to fluctuations in temperature and strain rate require complex strength models to accurately predict material performance in real-world applications. To predict such complex behavior, a robust and flexible strength model is necessary. In this work, we utilize genetic programming-based symbolic regression (GPSR) to develop data-driven strength models that accurately represent the measured stress–strain responses of tin across a wide range of strain, strain rate and temperature regimes. The GPSR models are constrained by physically-informed conditions, which leads to significant improvement in extrapolation. The best model is integrated into a multi-physics code to perform Taylor impact simulations, validating the model’s accuracy and robustness. In conclusion, the model predictions showed excellent agreement with experimental results, particularly when compared to predictions using traditional strength models.

Genetic programming↗

Observations and Modeling of Fiber-Optics Strain on Hydraulic Fracture Height Growth in HFTS-2

Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced dataset to allow us to observe and understand fracture geometries rigorously. Low frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. Furthermore, the simulation results improved the understating of the HFTS-2 LF-DAS data and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

58 GEOSCIENCES↗

Modeling creep behavior of fiber composites

A micromechanical model for the creep behavior of fiber composites is developed based on a typical cell consisting of a fiber and the surrounding matrix. The fiber is assumed to be linearly elastic and the matrix nonlinearly viscous. The creep strain rate in the matrix is assumed to be a function of stress. The nominal stress-strain relations are derived in the form of differential equations which are solved numerically for off-axis specimens under uniaxial loading. A potential function and the associated effective stress and effective creep strain rates are introduced to simplify the orthotropic relations.

Chen, J. L.↗

Solid Fuel Burning in Steady, Strained, Premixed Flow Fields: The Graphite/Air/Methane System

A detailed numerical investigation was conducted on the simultaneous burning of laminar premixed CH4/air flames and solid graphite in a stagnation flow configuration. The graphite and methane were chosen for this model, given that they are practical fuels and their chemical kinetics are considered as the most reliable ones among solid and hydrocarbon fuels, respectively. The simulation was performed by solving the quasi-one-dimensional equations of mass, momentum, energy, and species. The GRI 2.1 scheme was used for the gas-phase kinetics, while the heterogeneous kinetics were described by a six-step mechanism including stable and radical species. The effects of the graphite surface temperature, the gas-phase equivalence ratio, and the aerodynamic strain rate on the graphite burning rate and NO, production and destruction mechanisms were assessed. Results indicate that as the graphite temperature increases, its burning rate as well as the NO, concentration increase. Furthermore, it was found that by increasing the strain rate, the graphite burning rate increases as a result of the augmented supply of the gas-phase reactants towards the surface, while the NO, concentration decreases as a result of the reduced residence time. The effect of the equivalence ratio on both the graphite burning rate and NO, concentration was found to be non-monotonic and strongly dependent on the graphite temperature. Comparisons between results obtained for a graphite and a chemically inert surface revealed that the chemical activity of the graphite surface can result to the reduction of NO through reactions of the CH3, CH2, CH, and N radicals with NO.

Egolfopoulos, Fokion N.↗

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗

Observations and Modeling of Fiber Optic Strain on Hydraulic Fracture Height Growth in Hydraulic Fracturing Test Site 2 (HFTS-2)

Summary Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced data set to allow us to observe and understand fracture geometries rigorously. Low-frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1,000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. The simulation results improved the understating of the HFTS-2 LF-DAS data, and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

Engineering↗

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗

Characterization of elastic-viscoplastic properties of an AS4/PEEK thermoplastic composite

The elastic-viscoplastic properties of an AS4/PEEK (APC-2) thermoplastic composite were characterized at 24 C (75 F) and 121 C (250 F) by using a one-parameter viscoplasticity model. To determine the strain-rate effects, uniaxial tension tests were performed on unidirectional off-axis coupon specimens with different monotonic strain rates. A modified Bodner and Partom's model was also used to describe the viscoplasticity of the thermoplastic composite. The experimental results showed that viscoplastic behavior can be characterized quite well using the one-parameter overstress viscoplasticity model.

Yoon, K. J.↗

Effects of a temperature-dependent rheology on large scale continental extension

The effects of a temperature-dependent rheology on large-scale continental extension are investigated using a thin viscous sheet model. A vertically-averaged rheology is used that is consistent with laboratory experiments on power-law creep of olivine and that depends exponentially on temperature. Results of the calculations depend principally on two parameters: the Peclet number, which describes the relative rates of advection and diffusion of heat, and a dimensionless activation energy, which controls the temperature dependence of the rheology. At short times following the beginning of extension, deformation occurs with negligible change in temperature, so that only small changes in lithospheric strength occur due to attenuation of the lithosphere. However, after a certain critical time interval, thermal diffusion lowers temperatures in the lithosphere, strongly increasing lithospheric strength and slowing the rate of extension. This critical time depends principally on the Peclet number and is short compared with the thermal time constant of the lithosphere. The strength changes cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain. Results of the calculations are compared with observations from the Aegean, where maximum extensional strains are found in the south, near Crete, but maximum present-day strain rates are largest about 300 km further north.

Sonder, Leslie J.↗

Effects of a temperature-dependent rheology on large-scale continental extension

The effects of a temperature-dependent rheology on large-scale continental extension are investigated using a thin viscous sheet model. A vertically-averaged rheology is used that is consistent with laboratory experiments on power-law creep of olivine and that depends exponentially on temperature. Results of the calculations depend principally on two parameters: the Peclet number, which describes the relative rates of advection and diffusion of heat, and a dimensionless activation energy, which controls the temperature dependence of the rheology. At short times following the beginning of extension, deformation occurs with negligible change in temperature, so that only small changes in lithospheric strength occur due to attenuation of the lithosphere. However, after a certain critical time interval, thermal diffusion lowers temperatures in the lithosphere, strongly increasing lithospheric strength and slowing the rate of extension. This critical time depends principally on the Peclet number and is short compared with the thermal time constant of the lithosphere. The strength changes cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain. Results of the calculations are compared with observations from the Aegean, where maximum extensional strains are found in the south, near Crete, but maximum present-day strain rates are largest about 300 km further north.

Sonder, Leslie J.↗

Orthotropic Viscoelastic Modeling of Polymeric Battery Separator

The separator is one of the most important components for battery safety since its failure can cause short-circuit that may lead to thermal runaway. As such, understanding the mechanical behaviors of separators is critical to improve battery safety. Nevertheless, modeling the mechanical behaviors of separators is challenging since they usually depend on multiple factors such as orientation, temperature and strain rate due to the microstructure of separators. Here we develop an orthotropic linear viscoelastic model for separators, which takes both anisotropy and strain rate effects into account. This model is the first step of a future model that considers more factors such as temperature and large deformation. The model has been implemented as a user material model in a finite element package, where a discretization algorithm is developed to evaluate the stiffness-based hereditary integral with a kernel of Prony series. The user model has been applied to a polypropylene (PP) separator Celgard ® 2400. The simulation results agree well with experimental data, which include the rate dependent uniaxial tensile behavior of the separator in the two in-plane material directions and in shear, and the strain contour patterns of every strain component during biaxial tension.

25 ENERGY STORAGE↗

Distributed Fiber Optic Sensing for in-well hydraulic fracture monitoring

This study presents the results from in-well hydraulic fracture monitoring within a horizontal well in an unconventional reservoir utilizing Distributed Fiber Optic Sensing (DFOS). An in-house-developed Brillouin-based Distributed Strain Sensing (DSS) interrogator was deployed to obtain strain measurements, complemented by a commercial Raman-based Distributed Temperature Sensing (DTS) interrogator for temperature measurements and a commercial Rayleigh-based Low-Frequency Distributed Acoustic Sensing (LF-DAS) interrogator for strain-rate measurements. Examined over a ten-day period, the spatio-temporal distribution of temperature-compensated strain obtained from DSS and DTS revealed distinct signatures of the multi-stage hydraulic fracturing process. These signatures were analyzed with respect to fracture width growth and closure, residual strain effects, and fracture conductivity near the wellbore. Fracture widths within the fracture zone were estimated for individual stages. The findings were assessed with LF-DAS measurements for further evaluation. This work integrates DFOS-measured strain, temperature, and strain-rate data for monitoring in-well hydraulic fracturing, with the goal of supporting future studies in interpreting DFOS measurements for improved understanding of hydraulic fracturing in unconventional reservoirs.

58 GEOSCIENCES↗

Material Properties Measurements for Selected Materials

Hugoniot equation of state measurements were made on Coconino sandstone, Vacaville basalt, Kaibab limestone, Mono Crater, pumice and Zelux (a polycarbonate resin) for pressures to 2 Mb. A single data point was obtained for fused quartz at 1.6 Mb. In addition to the hugoniot studies, the uniaxial compressive stress behavior of Vacaville basalt and Zelux was investigated at strain rates from about 10(exp -5)/sec to 10(exp 3)/second. The data presented include the stress - strain relations as a function of strain rate for these two materials.

Green, S. J.↗

Rate Dependent Deformation and Strength Analysis of Polymer Matrix Composites

A research program is being undertaken to develop rate dependent deformation and failure models for the analysis of polymer matrix composite materials. In previous work in this program, strain-rate dependent inelastic constitutive equations used to analyze polymers have been implemented into a mechanics of materials based composite micromechanics method. In the current work, modifications to the micromechanics model have been implemented to improve the calculation of the effective inelastic strain. Additionally, modifications to the polymer constitutive model are discussed in which pressure dependence is incorporated into the equations in order to improve the calculation of constituent and composite shear stresses. The Hashin failure criterion is implemented into the analysis method to allow for the calculation of ply level failure stresses. The deformation response and failure stresses for two representative uniaxial polymer matrix composites, IM7/977-2 and AS4-PEEK, are predicted for varying strain rates and fiber orientations. The predicted results compare favorably to experimentally obtained values.

Goldberg, Robert K.↗