Impact of non-linear elastic behavior on bentonite density evolution at the FEBEX
Not Available
Takayama, Yusuke↗
SEARCH · Search NASA
Results for “Non-linear elasticity”
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.
Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)
Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.
Constitutive model↗
Topology optimization of bistable elastic structures — An application to logic gates
Not Available
42 ENGINEERING↗
Tensile Modeling PVC Gels for Electrohydraulic Actuators
Polyvinyl chloride (PVC)-dibutyl adipate (DBA) gels are a fascinating dielectric elastomer actuator showing promise in soft robotics. When actuated with high voltages, the gel deforms towards the anode. A recent application of PVC gels in electrohydraulic actuators motivates elastic and hyperelastic constitutive relationships for tensile loading modes. PVC gels with plasticizer-to-polymer weight ratios of 2:1, 4:1, 6:1, and 8:1 w/w were evaluated. PVC gels exhibit a linear elastic region up to 25% strain. The elastic modulus decreased with increasing plasticizer content from 288.8 kPa, 56.1 kPa, 24.7 kPa, to 11 kPa. Poisson’s ratio also decreased with increasing plasticizer content from 0.42, 0.43, 0.39, to 0.35. We suggest that the decrease in polymer concentration facilitates a weakly interconnected polymer network susceptible to chain slippage that hinders the network response, thus lowering Poisson’s ratio. Our work suggests that PVC gels can be treated as isotropic and incompressible for large strains and hyperelastic modeling; however, highly plasticized gels tend to act less incompressible at small strains. The power scaling law between the elastic modulus and plasticizer weight ratio showed high agreement, making the elastic modulus deterministic for any plasticizer content. The Neo–Hookean, Mooney–Rivlin, Yeoh, Gent, Ogden, and extended tube hyperelastic constitutive models are investigated. The Yeoh model shows the highest feasibility when evaluated up to 3.5 stretch, showing a maximum normalized root-mean-square-error of 6.85%. Together, these findings establish a constitutive basis for PVC-DBA gels, incorporating small strain elasticity, large strain non-linear behavior, and network analysis while providing suggestive insight into the network structure required for accurately modeling the EPIC.
Polymer Science↗
Drag reduction and the Vogel exponent of a flexible beam in transient shear flows
Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Predicting Mechanical Properties from Microstructure Images in Fiber-Reinforced Polymers Using Convolutional Neural Networks
Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.
Sun, Yixuan (ORCID:0000000311093380)↗
Sensitivity of the U.S. economy to oil prices controlling for domestic production and imports
This paper investigates the sensitivity of U.S. economic performance to oil price changes, accounting for changes in the domestic petroleum supply-demand balance over the last decade. Herein, a non-linear (asymmetric) autoregressive distributed lag (NARDL) model is used to estimate the U.S. GDP elasticity with respect to the oil price, controlling for oil production, consumption and imports, and macroeconomic variables. The positive and negative components of the oil price both have statistically significant long-run impacts on the real U.S. GDP. The parameter estimates imply that a 1% positive and permanent oil price shock, all else the same, would have a long-run impact of –0.045% on the U.S. economy (an elasticity of –0.045), but short-run parameters on the positive oil price terms are not significant at the 10% level. The long-run parameter on the negative oil price component term implies an elasticity of –0.034. Thus, controlling for domestic U.S. oil production and oil trade, the long-run oil price elasticity of the U.S. GDP remains within the range of estimates from previous studies. The results also show that domestic oil production and consumption have short-run impacts on the U.S. GDP. The potential extent of interactions among these variables, and implications for the net economic impacts, under an oil price shock are subjects of future research.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
Effects of internal swelling on residual elasticity of a quasi-brittle material through a composite sphere model
This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.
36 MATERIALS SCIENCE↗
Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength
Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.
36 MATERIALS SCIENCE↗
TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells
This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.
42 ENGINEERING↗
A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity
Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Price response of waste resources under demand shocks: four case studies
To quantify potential price changes of waste feedstocks that may be caused by increased demand, this study explored four case studies: 1) recyclable materials at local waste collection and recycling facilities in seven US states, 2) recyclable materials on internationally traded markets, 3) a wood-powered plant in Reading, PA, and 4) an MSW-powered plant in Tulsa, OK. The four case studies provided a total of ninety-one waste material price observations yielding twenty-seven price-supply relationships expressed as price elasticity of supply (PES) and unit price change with change in demand. For the two case studies related to recyclable materials with a global market shock, price changes vary widely from $0.08-$7.84 per ton with each percent change in demand. For the two case studies related to local market shocks in biomass used for energy, observed price changes range from $0.08-$0.23 per ton per percent change in demand. Based on these results, a base range of $0.10-$0.20 per ton per percent change in demand is recommended for local market shocks. I.e., for a facility that collects organic wastes on the order of 200,000-400,000 tons per year at price near zero, a doubling of demand could increase procurement prices by about $10.00-$20.00 per ton. Actual price changes with respect to change in demand may be non-linear. Shocks with broader geographic extent may increase price response due to a reduced opportunity for surrounding markets to mitigate localized market changes.
09 BIOMASS FUELS↗
Proliferated Resilient Economical Half-Meter Aperture Space Telescopes (PREEMPT)
The PREEMPT LDRD was motivated by a major national need for better space-based imaging systems that are both high performing and affordable for the US Government. Current optical payloads for intelligence, surveillance and reconnaissance (ISR) and space domain awareness (SDA) missions can cost hundreds of millions of dollars and take years to develop, which makes it difficult to build the large constellations required for persistent, world-wide coverage. To address this, the team aimed to advance a different kind of large-aperture (>25 cm) telescope, called a monolithic telescope, in which key optical surfaces are built into a single piece of fused silica. This design greatly reduces payload and spacecraft complexity, the need for precision focus actuators, improves mechanical and thermal robustness, and lowers cost when compared with traditional Cassegrain telescopes that rely on many precisely aligned components. The project focused on five primary thrusts. The first thrust was to advance the concept of a V10, 25 cm, monolithic telescope forward from optical design to flight-ready stage. This was accomplished in partnership with Optimax, who delivered the first test unit in the early stages of the LDRD. The team developed several technologies necessary for this optic to be integrated into a flight demonstration. These include carbon fiber housings, highly detailed structural and thermal models and stress-reducing elastic averaging Hirth groove designs. These technologies resulted in a successful maturation of the optic, which is now slated to fly in late 2026/early 2027 for a demonstration mission. The second thrust was to advance the manufacturability of these optics. In collaboration with NIF’s optical manufacturing shop, we reduced polishing time from 480 hours to 65 hours through the implementation of optimized processes and new tools. The NIF team utilized a conceptual V8 (18 cm) optic to demonstrate this optimization, though it can be applied to the rest of the monolithic optic portfolio. Third, the team developed the first conceptual V20 (50 cm) payload, which is slated to be the next generation of LLNL optical payload systems. A set of structural, dynamic and thermal simulations were performed to identify potential challenges in the future development of this payload. Early-stage simulations suggest the payload is feasible, though thermal management will be the key focus area to maintain optimal performance. Fourth, a non-linear model of Viton was developed, to further enhance the reliability of our structural and dynamic models for future payloads. Viton acts as the primary interface material between the optic and its housing. Lastly, the team focused on successfully displaying the feasibility of using additively manufactured metal composites for optical space payloads. The team successfully demonstrated layer by layer deposition of Al-SiC composites, which have highly tunable structural and coefficient of thermal expansion (CTE) properties. These are crucial for optical payloads because CTE mismatch is one of the causes for degraded optical performance for telescopes in orbit. Overall, the work showed that monolithic telescopes could become a practical, lower-cost path to high-resolution space imaging for both national security and scientific missions.
42 ENGINEERING↗
A generalized machine learning workflow to visualize mechanical discontinuity
Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.
42 ENGINEERING↗
Across the Scales of the Nucleus: Understanding Short Range Correlations from Medium Modification to Probe Independence
The atomic nucleus presents an intricate system due to the non-linear forces described by Quantum Chromodynamics (QCD) that govern its structure. The range of scales involved is remarkable; the most massive nuclei weigh approximately five orders of magnitude more than the quarks that compose them. The nucleus can be analyzed at various levels, from quarks to hadrons to the nucleus as a whole. Short-Range Correlations (SRCs) within the nucleus play a significant role that spans these diverse scales. At the most fundamental level, SRCs influence the interaction between nucleons. The nucleon-nucleon (NN) interaction, arising from QCD, is crucial in determining nuclear properties. SRCs serve as valuable probes for measuring this NN interaction, as the nucleons within SRCs become effectively decoupled from the rest of the nucleus. Multiple experimental techniques, including electron scattering, have been employed to investigate the NN interaction through SRCs. However, our first project demonstrates that inclusive measurements alone are inadequate to constrain this interaction fully. Moving to the scale of the nucleus, SRCs contribute to the high-momentum tail of the nuclear spectral function. While the low-momentum region is characterized by nucleons exhibiting bulk properties, nucleons begin to pair into SRCs at higher momenta. Our research aims to bridge the understanding between the mean-field portion of the nucleus and its high-momentum SRC components. Additionally, SRCs affect the quark structure of protons, as evidenced by the EMC effect, which indicates that quarks behave differently when protons are embedded within a nucleus—an effect referred to as medium modification. This thesis explores the correlation between SRCs and medium modification across various experimental setups. Finally, we seek to establish an interpretation of the nuclear ground state. Accomplishing this requires demonstrating that our SRC observables are independent of the probe’s scale and scheme. The concluding project of this thesis illustrates how we utilize triple coincidence quasi-elastic scattering across a range of (Q2) values to develop a model-dependent framework for understanding SRC distributions within the nucleus’s ground-state wavefunction.
Denniston, A. W. Denniston [University of Tel Aviv↗
A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate
SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.
58 GEOSCIENCES↗
Bifurcation Buckling Analysis and Non-Linear Collapse Analysis of Teardrop Shaped Vacuum Chamber
The Spallation Neutron Source (SNS) accelerator is being upgraded to increase the beam power from 1.4MW at 1GeV to 2.8MW at 1.3GeV. The currents in the middle two injection chicane magnets cannot simply be scaled up to accommodate the increased injection energy of 1.3GeV due to potential excessive H- stripping; the magnets must be replaced with longer, lower-field magnets and the associated vacuum chambers need to be redesigned. A new teardrop-shaped vacuum chamber was initially designed to accommodate the new magnets and the updated beam paths and instrumentation. This paper focuses on the structural stability study of the teardrop shape vacuum chamber based on buckling analysis. Protection against collapse from buckling according to the ASME BPVC requirement has been evaluated in depth. First, a Type-1 bifurcation buckling analysis using a linear eigenvalue solution to determine the critical load factor was performed. Subsequently, a Type-3 nonlinear collapse analysis was conducted using the static Riks method with elastic-plastic material properties and imperfections explicitly considered in the model geometry. The critical buckling load for the teardrop shape vacuum chamber was confidently estimated based upon this two-stage approach.
Jiang, Hao↗
Semi-supervised Bayesian Low-shot Learning
Deep neural networks (NNs) typically outperform traditional machine learning (ML) approaches for complicated, non-linear tasks. It is expected that deep learning (DL) should offer superior performance for the important non-proliferation task of predicting explosive device configuration based upon observed optical signature, a task which human experts struggle with. However, supervised machine learning is difficult to apply in this mission space because most recorded signatures are not associated with the corresponding device description, or “truth labels.” This is challenging for NNs, which traditionally require many samples for strong performance. Semi-supervised learning (SSL), low-shot learning (LSL), and uncertainty quantification (UQ) for NNs are emerging approaches that could bridge the mission gaps of few labels and rare samples of importance. NN explainability techniques are important in gaining insight into the inferential feature importance of such a complex model. In this work, SSL, LSL, and UQ are merged into a single framework, a significant technical hurdle not previously demonstrated. Exponential Average Adversarial Training (EAAT) and Pairwise Neural Networks (PNNs) are chosen as the SSL and LSL methods of choice. Permutation feature importance (PFI) for functional data is used to provide explainability via the Variable importance Explainable Elastic Shape Analysis (VEESA) pipeline. A variety of uncertainty quantification approaches are explored: Bayesian Neural Networks (BNNs), ensemble methods, concrete dropout, and evidential deep learning. Two final approaches, one utilizing ensemble methods and one utilizing evidential learning, are constructed and compared using a well-quantified synthetic 2D dataset along with the DIRSIG Megascene.
97 MATHEMATICS AND COMPUTING↗