Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLED STRESS”

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 217 records · Page 12

Determination of permeability for hydrocarbon release due to excavation-induced stress redistribution in rock salt

Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. Specifically, a multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure buildups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.

58 GEOSCIENCES↗

Cooling-induced reactivation of distant faults during long-term geothermal energy production in hot sedimentary aquifers

Deep geothermal energy (DGE) represents an opportunity for a sustainable and carbon-free energy supply. One of the main concerns of DGE is induced seismicity that may produce damaging earthquakes, challenging its widespread exploitation. It is widely believed that the seismicity risk can be controlled by using doublet systems circulating water to minimize the injection-induced pressure changes. However, cold water reinjection may also give rise to thermal stresses within and beyond the cooled region, whose potential impacts on fault reactivation are less well understood. Here, we investigate by coupled thermo-hydro-mechanical modeling the processes that may lead to fault reactivation in a hot sedimentary aquifer (HSA) in which water is circulated through a doublet. We show that thermal stresses are transmitted much ahead of the cooled region and are likely to destabilize faults located far away from the doublet. Meanwhile, the fault permeability mainly controls the fault reactivation timing, which entails the importance of employing appropriate characterization methods. This investigation is crucial for understanding the mechanisms controlling induced seismicity associated with DGE in a HSA and allows the success of future DGE projects.

15 GEOTHERMAL ENERGY↗

Sea surface warming and ocean-to-atmosphere feedback driven by large-scale offshore wind farms under seasonally stratified conditions

Offshore wind farms may induce changes in the upper ocean and near-surface atmosphere through coupled ocean-atmosphere feedbacks. Yet, the role of air-sea interactions mediated by offshore wind farms remains poorly understood. Using fully coupled ocean-atmosphere-wave model simulations for seasonally stratified conditions along the US East Coast, we show that simulated cumulative reductions in wind stress due to large-scale wind farm clusters lead to sea surface warming of 0.3° to 0.4°C and a shallower mixed layer. This warming drives upward heat fluxes, destabilizing the atmospheric boundary layer and enhancing wind stress, which partially offsets wake-induced wind deficits. These wake-ocean interactions influence near-surface meteorology and air-sea fluxes, suggesting that a coupled modeling approach may be necessary for assessing potential oceanographic impacts of offshore wind developments. However, ocean coupling exerts limited influence on winds at turbine-relevant heights or within downstream wakes, resulting in minimal impact on long-term energy. These findings suggest that models without ocean coupling may be adequate for wind energy applications.

17 WIND ENERGY↗

Polymorphism of garnet solid electrolytes and its implications for grain-level chemo-mechanics

Understanding and mitigating filament formation, short-circuit and solid electrolyte fracture is necessary for advanced all-solid-state batteries. Here, we employ a coupled far-field high-energy diffraction microscopy and tomography approach for assessing the chemo-mechanical behaviour for dense, polycrystalline garnet (Li 7 La 3 Zr 2 O 12 ) solid electrolytes with grain-level resolution. In situ monitoring of grain-level stress responses reveals that the failure mechanism is stochastic and affected by local microstructural heterogeneity. In this work, coupling high-energy X-ray diffraction and far-field high-energy diffraction microscopy measurements reveals the presence of phase heterogeneity that can alter local chemo-mechanics within the bulk solid electrolyte. These local regions are proposed to be regions with the presence of a cubic polymorph of LLZO, potentially arising from local dopant concentration variation. The coupled tomography and FF-HEDM experiments are combined with transport and mechanics modelling to illustrate the degradation of polycrystalline garnet solid electrolytes. The results showcase the pathways for processing high-performing solid-state batteries.

36 MATERIALS SCIENCE↗

A model to assess Zircaloy’s mechanical property changes following a transient beyond critical heat flux

Maintaining the integrity of nuclear fuel rods is essential for ensuring public health and safety in nuclear power generation. During reactor operation, this integrity is confirmed by demonstrating compliance with established regulatory acceptance criteria. For moderate-frequency events, such as limiting transients and anticipated operational occurrences (AOOs), the current fuel integrity criterion is based on preventing boiling transition. This criterion assumes that prevention of boiling transition will prevent excessive cladding heating and, thus, fuel failure during normal operations. While conservative, this approach places significant constraints on core design, fuel cycle economics, and a plant’s ability to perform major power uprates, leading to suboptimal fuel utilization and inefficient carbon-free energy production. A more efficient approach could be achieved by revising the failure criterion to a material-specific limit rather than strictly preventing the boiling transition, since boiling transition per se is not a cause of fuel cladding failure. Here, as a result, a new licensing framework based on material properties, termed time-at-temperature (t@T), is needed. This approach would allow for brief periods of post–critical heat flux operation during an AOO without compromising safety. Implementing the t@T licensing strategy requires a robust technical foundation in material properties, which must be established through comprehensive data collection on both unirradiated and irradiated fuel and cladding materials. This foundation would enable the development of a safety basis that ensures safe operation while providing greater flexibility and efficiency for reactor operation. This paper documents a thorough review of the available data to establish a baseline knowledge that can inform the development of cladding mechanical models, as well as identify experimental data gaps that need to be addressed in future research. Machine learning and data informatics were utilized to extract the importance of parameters on the t@T parameter. Industry tools were used to perform baseline analyses to define the relevant transient conditions for data analysis. The subsequent review successfully identified applicable experimental data, as well as sufficient data to evaluate changes in cladding mechanical properties following an AOO transient. Rather than developing new models, this work coupled existing irradiation annealing and recrystallization models to calculate changes in hardness, yield stress, and ultimate tensile stress following an AOO event. The findings from this review were summarized to highlight the experimental data needs required to fill remaining gaps and support the development of future t@T licensing methodologies.

Cladding performance↗

Scale–Bridging in Three–Dimensional Fracture Networks: Characterizing the Effects of Variable Fracture Apertures on Network–Scale Flow Channelization

Here, we incorporate observations of real fracture aperture variability observed in laboratory experiments into an ensemble of three-dimensional discrete fracture network (DFN) simulations to characterize how variations of this micro-scale feature can influence flow and transport behavior at the network scale. A shear fracture is created within a Marcellus shale sample, and the fracture aperture is measured using a triaxial direct-shear device coupled with real-time X-ray imaging at in-situ stress conditions. We construct an ensemble of fracture networks based on natural fractures in Marcellus shale and project regions of the experimental aperture field onto each fracture in the networks. Our calculations demonstrate that the degree of flow channelization, a network-scale flow field structure, is dramatically increased by local changes in the aperture field that in turn affects flow and transport properties.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗

Biogeochemical Changes in Estuaries

The biogeochemical composition of materials suspended and dissolved in estuarine waters is substantially transformed by biological and abiotic processes as water from rivers and subterranean aquifers travels through estuaries to the coastal ocean, influencing the rate of storage and/or exchange of carbon with the atmosphere. Biogeochemical signatures in the upper reaches of estuaries (e.g., low salinity) generally reflect characteristics of the rivers and groundwater draining into the estuary. The composition of organic matter (OM) where rivers meet estuaries is a mixture of terrigenous signatures dependent on watershed landscape features, aquatic signatures dependent on the prominence of freshwater primary productivity, and intermediate metabolites (i.e., byproducts of in situ transformation) of these terrigenous and aquatic OM sources. The prominence of aquatic primary production-derived OM and products of microbial decomposition generally increases towards the lower reaches of estuaries where they meet the coastal ocean (i.e.., continental shelf). The nature of biogeochemical transformations occurring along estuarine gradients is closely coupled with factors that are poised to change under multiple compounding disturbances and stress factors related to both climate change and other anthropogenic global changes. This chapter describes the fundamental mechanisms that drive the transport, transformation, and storage of OM in estuaries and how these dynamics have and may continue to change in response to key disturbances and stressors related to climate change and other anthropogenic global changes including warming, intensifying hydrological extremes (e.g., droughts and floods), sea level rise, and acidification.

Ward, Nicholas D.↗

The Brachypodium ENCODE Project – from sequence to function: Predicting physiological responses in grasses to facilitate engineering of biofuel crops

The overarching goal of this proposal was to gain deep knowledge of chromatin dynamics and gene networks related to stress responses in grasses. To accomplish this goal our Brachypodium ENCODE project (1) coupled drought-focused physiological analysis with multiple layers of ENCODE ‘omics’ profiling; (2) obtained fundamental information about the Brachypodium transcriptome, epigenome, chromatin dynamics, and gene regulatory networks; and (3) elucidated the regulatory dynamics of gene networks that can be manipulated to increase energy grass crop productivity.

09 BIOMASS FUELS↗

Conformational spread drives the evolution of the calcium–calmodulin protein kinase II

The calcium calmodulin (Ca 2+ /CaM) dependent protein kinase II (CaMKII) decodes Ca 2+ frequency oscillations. The CaMKIIα isoform is predominantly expressed in the brain and has a central role in learning. I matched residue and organismal evolution with collective motions deduced from the atomic structure of the human CaMKIIα holoenzyme to learn how its ring architecture abets function. Protein dynamic simulations showed its peripheral kinase domains (KDs) are conformationally coupled via lateral spread along the central hub. The underlying β-sheet motions in the hub or association domain (AD) were deconvolved into dynamic couplings based on mutual information. They mapped onto a coevolved residue network to partition the AD into two distinct sectors. A second, energetically stressed sector was added to ancient bacterial enzyme dimers for assembly of the ringed hub. The continued evolution of the holoenzyme after AD–KD fusion targeted the sector’s ring contacts coupled to the KD. Among isoforms, the α isoform emerged last and, it alone, mutated rapidly after the poikilotherm–homeotherm jump to match the evolution of memory. The correlation between dynamics and evolution of the CaMKII AD argues single residue substitutions fine-tune hub conformational spread. The fine-tuning could increase CaMKIIα Ca 2+ frequency response range for complex learning functions.

59 BASIC BIOLOGICAL SCIENCES↗

Atomistic Mechanisms of Stress-Dependent Molten Salt Corrosion in NiCr Alloys

Ni-based structural alloys in molten salt environments often experience simultaneous mechanical loading and corrosive attack, yet the mechanisms governing stress-corrosion interactions remain unclear. Prior studies largely emphasize tensile stress, while the role of compressive stress has received limited attention. Here, reactive molecular dynamics simulations are used to investigate the coupled effects of applied strain and corrosion in Ni 0.75 Cr 0.25 exposed to molten FLiNaK at 800 °C. A Σ5(210) grain boundary model is subjected to tensile (+4%) to compressive (−4%) uniaxial strains, and corrosion behavior is evaluated through fluorine adsorption, charge redistribution, and grain boundary evolution. Tensile strain accelerates intergranular corrosion susceptibility by reducing local atomic packing through elastic dilation and increasing excess free volume at the grain boundary, which enhances atomic mobility and salt infiltration. In contrast, compressive strain can suppress corrosion by promoting the formation of a ridge-like surface layer along the grain boundary, limiting salt access to the underlying alloy. These results provide atomistic insight into how stress states influence grain boundary corrosion in molten salts.

36 - MATERIALS SCIENCE↗

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗

Fluorescence-coupled micropipette aspiration assay to examine calcium mobilization caused by red blood cell mechanosensing

Abstract Mechanical stimuli such as tension, compression, and shear stress play critical roles in the physiological functions of red blood cells (RBCs) and their homeostasis, ATP release, and rheological properties. Intracellular calcium (Ca 2+ ) mobilization reflects RBC mechanosensing as they transverse the complex vasculature. Emerging studies have demonstrated the presence of mechanosensitive Ca 2+ permeable ion channels and their function has been implicated in the regulation of RBC volume and deformability. However, how these mechanoreceptors trigger Ca 2+ influx and subsequent cellular responses are still unclear. Here, we introduce a fluorescence-coupled micropipette aspiration assay to examine RBC mechanosensing at the single-cell level. To achieve a wide range of cell aspirations, we implemented and compared two negative pressure adjusting apparatuses: a homemade water manometer (− 2.94 to 0 mmH 2 O) and a pneumatic high-speed pressure clamp (− 25 to 0 mmHg). To visualize Ca 2+ influx, RBCs were pre-loaded with an intensiometric probe Cal-520 AM, then imaged under a confocal microscope with concurrent bright-field and fluorescent imaging at acquisition rates of 10 frames per second. Remarkably, we observed the related changes in intracellular Ca 2+ levels immediately after aspirating individual RBCs in a pressure-dependent manner. The RBC aspirated by the water manometer only displayed 1.1-fold increase in fluorescence intensity, whereas the RBC aspirated by the pneumatic clamp showed up to threefold increase. These results demonstrated the water manometer as a gentle tool for cell manipulation with minimal pre-activation, while the high-speed pneumatic clamp as a much stronger pressure actuator to examine cell mechanosensing directly. Together, this multimodal platform enables us to precisely control aspiration and membrane tension, and subsequently correlate this with intracellular calcium concentration dynamics in a robust and reproducible manner.

Wang, Haoqing↗

Multi-Physics Coupled Seismic Safety Analysis of Molten-Salt Reactors

Pool-type Molten Salt Reactors (MSRs) are a promising advanced nuclear reactor concept relying on passive physical behavior to provide increased safety characteristics. The high heat capacity, high boiling point, unpressurized liquid salt contains the nuclear fuel and as it passes through the vessel achieves criticality and produces nuclear power. The salt then carries away the heat produced along with the delayed neutron precursors towards heat exchangers and primary pumps before entering the core. This forms a highly coupled multiphysics problem between neutronics and fluid flow, which makes modeling these reactors challenging. The safety/licensing case heavily relies on modeling and simulation to prove their safety, and significant effort has been expanded by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop the appropriate simulation tools. Their safety/licensing will require a thorough analysis of their behavior during earthquake-based transients. Earthquakes are a transient condition commonly experienced in many regions, and are safely survived by dozens of reactors every year. Nonetheless, significant conservatisms were introduced with limited predictive modeling and simulation available. The first-of-a-kind capability developed by this research will allow for engineers to analyze ranges of design-basis and beyond-design basis accident scenarios to assess the integrity of the reactor vessel and surrounding system, possibly allowing for limiting the costly conservatisms in the design. This projects uses MASTODON, Griffin, and Pronghorn in the NEAMS ecosystem to model the coupled multiphysics problem posed by an earthquake in a molten salt reactor. It marches through a number of number of coupling schemes, from uncoupled simulations to tightly two-way coupled simulations. The appropriate level of coupling for these simulations will then be determined based on the accuracy in the quantities of interest (QoI) and computational effort involved in each scheme. The QoIs chosen for this paper include the mass flow rate across the heat exchanger (for Pronghorn simulations), power output of the reactor (for Griffin simulations), and the maximum Von Mises stress (for MASTODON simulations). This summary paper reports on the current progress of this project, with 2D tightly coupled multiphysics results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading

Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.

36 MATERIALS SCIENCE↗

Electrochemical-mechanical coupling failure mechanism of composite cathode in all-solid-state batteries

Composite cathode composed of active particles and solid electrolytes (SEs) can considerably enlarge the particle-SE contact areas and achieve high areal loadings in all-solid-state batteries (ASSBs). However, the challenging interfacial instability and particle damage problems remain unsolved. Herein, we establish a 3D electrochemical-mechanical coupled model to investigate the underlying failure mechanism by considering the governing electrochemical and physics processes. Micro-scale heterogeneous primary particles with random crystallographic orientation and size inside the LiNi 1/3 Co 1/3 Mn 1/3 O 2 (NCM111) secondary particle of the model result in the anisotropic Li diffusion and volume variation within the secondary particle, leading to significant nonuniformity of the Li concentration, and GPa-level stress distributions at primary particle boundaries, and finally causing the particle internal cracks. The particle volume shrinkage under the constraint of stiff Li 7 La 3 Zr 2 O 12 (LLZO) SE triggers the interface debonding (gap>50 nm) with increased interfacial impedance to degrade cell capacity. Higher C-rates result in larger residual stress (~100 MPa)/strain/debonding gap at dis-charging end, more likely to deteriorate the cell performance. Increasing the interfacial strength between the particle and SE can suppress the interface debonding but induces high stress (up to 10 GPa). In conclusion, results reveal the underlying mechanism of the electrochemical-mechanical coupling failure mechanism for composite cathode and provide promising guidance on the further improvement of a more robust composite cathode for ASSBs.

25 ENERGY STORAGE↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗