Search NASA⌕ Search

SEARCH · Search NASA

Results for “Load Balancing”

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 1,765 records · Page 98

Coupled methane and NO x conversion on Pt + Pd/Al 2 O 3 monolith: Conversion enhancement through feed modulation and Mn 0.5 Fe 2.5 O 4 spinel addition

The performance of a dual-layer Pt-Pd/Al 2 O 3 washcoated monolith promoted by an Al 2 O 3 -supported Mn 0.5 Fe 2.5 O 4 spinel base layer was evaluated with the objective to develop a catalyst with a lower loading of Platinum group metal (PGM) for stoichiometric natural gas vehicle (NGV) emission control. Flow reactor results show that a combination of a lean/rich modulated feed and spinel addition gives enhanced methane and NO conversions compared to a steady-state feed having the same overall composition. Washcoated catalysts – PGM-only, spinel-only, PGM + spinel – were compared over a wide range of operating conditions. The results show a decrease in the methane conversion light-off temperature of up to 85 °C for an application-relevant feed (CH 4 + NO + H 2 + CO + O 2 + H 2 O + CO 2 ) while a moderate 40 °C decrease was observed for a simple feed (CH 4 + O 2 + H 2 O + CO 2 ). In the absence of spinel the modulation enhancement is negatively impacted at high methane conversion (> 80 %). The data show that the O 2 storage and release property of the spinel is needed to fully exploit the modulation enhancement. The mechanism responsible for the enhancement is linked to suppression of O 2 inhibition on the methane oxidation rate near the stoichiometric neutral point (λ = 1). Sufficiently fast cycling achieves a balance between a metallic and oxidic precious metal crystallites favorable for methane activation. The addition of spinel moves this balance closer to the λ = 1 feed, enabling a high oxidation efficiency with lower byproduct NH 3 selectivity. Enhancement in the NO conversion follows the methane conversion due to an easier NO reduction under the slightly rich of stoichiometric feed. Comparison of the PGM + spinel catalysts to PGM-only catalysts shows that the addition of spinel in combination with modulation can achieve the same activity with a lower PGM loading.

03 NATURAL GAS↗

Control-oriented core-SOL-divertor model to address integrated burn and divertor control challenges in ITER

The real-time regulation of a burning plasma’s temperature and density, or burn control, will be necessary to produce high fusion power in future tokamaks like ITER. This is made more challenging due to the plasma’s nonlinear characteristics and the interdependence between the core-plasma and edge-plasma regions. For example, a raising plasma temperature leads to increasing reactivity and therefore to more alpha-particle heating, which further increases temperature. Furthermore, a raise of the fusion power increases the heat flow through the scrape-off-layer (SOL), which can compromise the integrity of the divertor without proper safeguards. For control design, a model-based approach is attractive because it can directly incorporate the nonlinear, coupled, burning-plasma dynamics into the design. To facilitate this design approach, a control-oriented core-SOL-divertor (CSD) model is presented in this work. In this CSD model, a core-plasma model captures the nonlinear dynamics of the core’s density and temperature, and a SOL-divertor model defines the plasma conditions at the separatrix and divertor including the heat load on the target plates. The core-plasma and SOL-divertor models are coupled through the exchange of various variables. In particular, the SOL-divertor model yields the separatrix temperature and the influx of recycled particles into the core-plasma. Further, these variables influence the power and particle balances captured by the core-plasma model. In return, the core-plasma model determines the intensity of the heat and particles fluxes across the separatrix, and this outflow strongly impacts the SOL-divertor model. Therefore, the power and density of the core-plasma, which can be readily modulated through external heating systems and pellet injection, can be viewed as control knobs for the SOL-divertor region in addition to the gas puffing. In simulations of the CSD model, it is demonstrated how external actuation can be utilized to meet burn control and divertor control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Catalyzed Nano-Framework Stablized High Density Reversible Hydrogen Storage Systems

A wide range of high capacity on-board rechargeable material candidates have exhibited non-ideal behavior related to irreversible hydrogen discharge / recharge behavior, and kinetic instability or retardation. This project addresses these issues by incorporating solvated and other forms of complex metal hydrides, with an emphasis on borohydrides, into nano-scale frameworks of low density, high surface area skeleton materials to stabilize, catalyze, and control desorption product formation associated with such complex metal hydrides. A variety of framework chemistries and hydride / framework combinations were investigated to make a relatively broad assessment of the method's potential. In this project, the hydride / framework interactions were tuned to decrease desorption temperatures for highly stable compounds or increase desorption temperatures for unstable high capacity compounds, and to influence desorption product formation for improved reversibility. First principle modeling was used to explore heterogeneous catalysis of hydride reversibility by modeling H 2 dissociation, hydrogen migration, and rehydrogenation. Atomic modeling also demonstrated enhanced NaTi(BH 4 ) 4 stabilization at nano-framework surfaces modified with multi-functional agents. Amine multi-functional agents were found to have more balanced interactions with nano-framework and hydride clusters than other functional groups investigated. Experimentation demonstrated that incorporation of Ca(BH 4 ) 2 and Mg(BH 4 ) 2 in aerogels enhanced hydride desorption kinetics. Carbon aerogels were identified as the most suitable nano-frameworks for hydride kinetic enhancement and high hydride loading. High loading of NaTi(BH 4 ) 4 ligand complex in SiO 2 aerogel was achieved and hydride stability was improved with the aerogel. Although improvements of desorption kinetics was observed, the incorporation of Ca(BH 4 ) 2 and Mg(BH 4 ) 2 in nano-frameworks did not improve their H 2 absorption due to the formation of stable alkaline earth B 12 H 12 intermediates upon rehydrogenation. This project primarily investigated the effect of nano-framework surface chemistry on hydride properties, while the effect of pore size is the focus area of other efforts (e.g., HRL, Sandia National Laboratories (SNL) etc.) within the Metal Hydride Center of Excellence (MHCoE). The projects were complementary in gaining an overall understanding of the influence of nano-frameworks on hydride behavior.

08 HYDROGEN↗

Experimental Evaluation of PV Inverter Anti-Islanding with Grid Support Functions in Multi-Inverter Island Scenarios

As PV and other DER systems are connected to the grid at increased penetration levels, island detection may become more challenging for two reasons: 1.) In islands containing many DERs, active inverter-based anti-islanding methods may have more difficulty detecting islands because each individual inverter's efforts to detect the island may be interfered with by the other inverters in the island. 2.) The increasing numbers of DERs are leading to new requirements that DERs ride through grid disturbances and even actively try to regulate grid voltage and frequency back towards nominal operating conditions. These new grid support requirements may directly or indirectly interfere with anti-islanding controls. This report describes a series of tests designed to examine the impacts of both grid support functions and multi-inverter islands on anti-islanding effectiveness. Crucially, the multi-inverter anti-islanding tests described in this report examine scenarios with multiple inverters connected to multiple different points on the grid. While this so-called 'solar subdivision' scenario has been examined to some extent through simulation, this is the first known work to test it using hardware inverters. This was accomplished through the use of power hardware-in-the-loop (PHIL) simulation, which allows the hardware inverters to be connected to a real-time transient simulation of an electric power system that can be easily reconfigured to test various distribution circuit scenarios. The anti-islanding test design was a modified version of the unintentional islanding test in IEEE Standard 1547.1, which creates a balanced, resonant island with the intent of creating a highly challenging condition for island detection. Three common, commercially available single-phase PV inverters from three different manufacturers were tested. The first part of this work examined each inverter individually using a series of pure hardware resistive-inductive-capacitive (RLC) resonant load based anti-islanding tests to determine the worst-case configuration of grid support functions for each inverter. A grid support function is a function an inverter performs to help stabilize the grid or drive the grid back towards its nominal operating point. The four grid support functions examined here were voltage ride-through, frequency ride-through, Volt-VAr control, and frequency-Watt control. The worst-case grid support configuration was defined as the configuration that led to the maximum island duration (or run-on time, ROT) out of 50 tests of each inverter. For each of the three inverters, it was observed that maximum ROT increased when voltage and frequency ride-through were activated. No conclusive evidence was found that Volt-VAr control or frequency-Watt control increased maximum ROT. Over all single-inverter test cases, the maximum ROT was 711 ms, well below the two-second limit currently imposed by IEEE Standard 1547-2003. A subsequent series of 244 experiments tested all three inverters simultaneously in the same island. These tests again used a procedure based on the IEEE 1547.1 unintentional islanding test to create a difficult-to-detect island condition. For these tests, which used the two worst-case grid support function configurations from the single-inverter tests, the inverters were connected to a variety of island circuit topologies designed to represent the variety of multiple-inverter islands that may occur on real distribution circuits. The interconnecting circuits and the resonant island load itself were represented in the real-time PHIL model. PHIL techniques similar to those employed here have been previously used and validated for anti-islanding tests, and the PHIL resonant load model used in this test was successfully validated by comparing single-inverter PHIL tests to conventional tests using an RLC load bank.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thermal coupling apparatus for electroexplosive devices

This paper describes the equipment for nondestructive testing of electroexplosive devices using a thermal coupling technique. The technique utilizes a self-balancing bridge to apply a nominal amount of power to the bridge wire of the device. The bridge wire behaving as a thermal sensor detects heat flow into the device when heat is preferentially directed down the center of the explosive column toward the bridge wire. Thermal equilibrium is upset and the change in power from the self-balancing bridge is measured directly in watts. Faults such as omissions, incorrect materials, contamination, and so on, can be detected by observing the rate of heat flow to the bridge wire. The response is indicative of the thermal resistance and diffusivity of the system. The paper also presents data on electroexplosive devices loaded with inert and explosive materials.

Menichelli, V. J.↗

Coupled rotor-flexible fuselage vibration reduction using open loop higher harmonic control

A fundamental study of vibration prediction and vibration reduction in helicopters using active controls was performed. The nonlinear equations of motion for a coupled rotor/flexible fuselage system have been derived using computer algebra on a special purpose symbolic computer facility. The trim state and vibratory response of the helicopter are obtained in a single pass by applying the harmonic balance technique and simultaneously satisfying the trim and the vibratory response of the helicopter for all rotor and fuselage degrees of freedom. The influence of the fuselage flexibility on the vibratory response is studied. It is shown that the conventional single frequency higher harmonic control is capable of reducing either the hub loads or only the fuselage vibrations but not both simultaneously. It is demonstrated that for simultaneous reduction of hub shears and fuselae vibrations a new scheme called multiple higher harmonic control is required.

Papavassiliou, I.↗

A study of coupled rotor-fuselage vibration with higher harmonic control using a symbolic computing facility

A fundamental study of vibration prediction and vibration reduction in helicopters using active controls was performed. The nonlinear equations of motion for a coupled rotor/flexible fuselage system have been derived using computer algebra on a special purpose symbolic computing facility. The details of the derivation using the MACSYMA program are described. The trim state and vibratory response of the helicopter are obtained in a single pass by applying the harmonic balance technique and simultaneously satisfying the trim and the vibratory response of the helicopter for all rotor and fuselage degrees of freedom. The influence of the fuselage flexibility on the vibratory response is studied. It is shown that the conventional single frequency higher harmonic control (HHC) capable of reducing either the hub loads or only the fuselage vibrations but not both simultaneously. It is demonstrated that for simultaneous reduction of hub shears and fuselage vibrations a new scheme called multiple higher harmonic control (MHHC) is required. The fundamental aspects of this scheme and its uniqueness are described in detail, providing new insight on vibration reduction in helicopters using HHC.

Papavassiliou, I.↗

Lunar surface structural concepts and construction studies

The topics are presented in viewgraph form and include the following: lunar surface structures construction research areas; lunar crane related disciplines; shortcomings of typical mobile crane in lunar base applications; candidate crane cable suspension systems; NIST six-cable suspension crane; numerical example of natural frequency; the incorporation of two new features for improved performance of the counter-balanced actively-controlled lunar crane; lunar crane pendulum mechanics; simulation results; 1/6 scale lunar crane testbed using GE robot for global manipulation; basic deployable truss approaches; bi-pantograph elevator platform; comparison of elevator platforms; perspective of bi-pantograph beam; bi-pantograph synchronously deployable tower/beam; lunar module off-loading concept; module off-loader concept packaged; starburst deployable precision reflector; 3-ring reflector deployment scheme; cross-section of packaged starburst reflector; and focal point and thickness packaging considerations.

Mikulas, Martin↗

Aeroelasticity and structural optimization of composite helicopter rotor blades with swept tips

This report describes the development of an aeroelastic analysis capability for composite helicopter rotor blades with straight and swept tips, and its application to the simulation of helicopter vibration reduction through structural optimization. A new aeroelastic model is developed in this study which is suitable for composite rotor blades with swept tips in hover and in forward flight. The hingeless blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. Arbitrary cross-sectional shape, generally anisotropic material behavior, transverse shears and out-of-plane warping are included in the blade model. The nonlinear equations of motion, derived using Hamilton's principle, are based on a moderate deflection theory. Composite blade cross-sectbnal properties are calculated by a separate linear, two-dimensional cross section analysis. The aerodynamic loads are obtained from quasi-steady, incompressible aerodynamics, based on an implicit formulation. The trim and steady state blade aeroelastic response are solved in a fully coupled manner. In forward flight, where the blade equations of motion are periodic, the coupled trim-aeroelastic response solution is obtained from the harmonic balance method. Subsequently, the periodic system is linearized about the steady state response, and its stability is determined from Floquet theory.

Yuan, K. A.↗

Development of a solar-cell dust opacity measurement instrument for Mars Pathfinder

The atmosphere of Mars has a considerable load of suspended dust. Over time, this dust is deposited out of the atmosphere. The mechanism and the temporal and geographical variation of this deposition are not well characterized. Measurements of settling rates and dust properties are of considerable scientific interest. Atmospheric dust affects the atmospheric solar absorption and thus the heat balance of Mars, as well as serving as nucleation sites for water and CO2 frost. Knowledge of dust properties is of critical interest to design and prediction of the lifetime and power output of solar arrays, and also to design of mechanical mechanisms and radiators. An instrument has been designed and fabricated to measure the dust accumulation during the course of the Mars Pathfinder rover mission. The solar-cell coverglass transmission experiment will measure the change in optical opacity of a transparent coverglass as dust settles on the surface, and a quartz crystal monitor will measure the mass deposited.

Geoffrey A. Landis↗

Floating wind turbine system

A floating wind turbine system with a tower structure that includes at least one stability arm extending therefrom and that is anchored to the sea floor with a rotatable position retention device that facilitates deep water installations. Variable buoyancy for the wind turbine system is provided by buoyancy chambers that are integral to the tower itself as well as the stability arm. Pumps are included for adjusting the buoyancy as an aid in system transport, installation, repair and removal. The wind turbine rotor is located downwind of the tower structure to allow the wind turbine to follow the wind direction without an active yaw drive system. The support tower and stability arm structure is designed to balance tension in the tether with buoyancy, gravity and wind forces in such a way that the top of the support tower leans downwind, providing a large clearance between the support tower and the rotor blade tips. This large clearance facilitates the use of articulated rotor hubs to reduced damaging structural dynamic loads. Major components of the turbine can be assembled at the shore and transported to an offshore installation site.

Viterna, Larry A.↗

Simulation of dynamic fracture using peridynamics, finite element modeling, and contact.

Peridynamics is a nonlocal extension of classical solid mechanics that allows for the modeling of bodies in which discontinuities occur spontaneously. Because the peridynamic expression for the balance of linear momentum does not contain spatial derivatives and is instead based on an integral equation, it is well suited for modeling phenomena involving spatial discontinuities such as crack formation and fracture. In this study, both peridynamics and classical finite element analysis are applied to simulate material response under dynamic blast loading conditions. A combined approach is utilized in which the portion of the simulation modeled with peridynamics interacts with the finite element portion of the model via a contact algorithm. The peridynamic portion of the analysis utilizes an elastic-plastic constitutive model with linear hardening. The peridynamic interface to the constitutive model is based on the calculation of an approximate deformation gradient, requiring the suppression of possible zero-energy modes. The classical finite element portion of the model utilizes a Johnson-Cook constitutive model. Simulation results are validated by direct comparison to expanding tube experiments. The coupled modeling approach successfully captures material response at the surface of the tube and the emerging fracture pattern. The coupling of peridynamics and finite element analysis via a contact algorithm has been shown to be a viable means for simulating material fracture in a high-velocity impact experiment. A combined peridynamics/finite element approach was applied to model an expanding tube experiment performed by Vogler, et al., in which loading on the tube is a result of Lexan slugs impacting inside the tube. The Lexan portion of the simulation was modeled with finite elements and a Johnson-Cook elastic-plastic material model in conjunction with an equation-of-state law. The steel tube portion of the simulation was modeled with peridynamics, an elastic-plastic material model, and a critical stretch bond damage model. The application of peridynamics to the tube portion of the model allowed the capture of the formation of cracks and eventual fragmentation of the tube. The simulation results yielded good agreement with the experimental results published by Vogler, et al., for the velocity and displacement profiles on the surface of the tube and the resulting fragment distribution. Numerical difficulties were encountered that required removal of hexahedron elements from the Lexan portion of the model over the course of the simulation. The significant number of inverted and nearly-inverted elements appearing over the course of the simulation is believed to be a result of irregularities in the contact between the Lexan and AerMet portions of the model, and was likely exacerbated by the ultra-high strength of the AerMet tube. Future simulations are planned in which the Lexan portion of the simulation is modeled with peridynamics, or with an alternative method such as smoothed particle hydrodynamics, with the goal of reducing these numerical difficulties.

36 MATERIALS SCIENCE↗

An overlapping Domain Decomposition preconditioning method for monolithic solution of shear bands

Metals subjected to high strain rate impact often exhibit a sudden and profound drop in the material’s load bearing capability, a ductile failure phenomenon known as shear banding. Shear bands, characterized as material instabilities, are driven by shear heating and described as narrow regions that have sustained intense plastic deformation and high temperature rise. This coupled thermo-mechanical localization problem can be formulated as a nonlinear system with two balance equations, momentum and energy, and two constitutive laws for elasticity and plasticity. Here in our formulation, mixed finite elements are used to discretize the equations in space and an implicit finite difference scheme is used to advance the system in time, where at every time step, a Newton type method is used to solve the system monolithically. To that end, a block Jacobian matrix is formed analytically using Gâteaux derivatives. The resulting Jacobian is sparse, nonsymmetric and its sparsity pattern and eigenvalue content vary with the different stages of shear bands formation, consequently posing a significant challenge to iterative solvers. To address this issue, an overlapping Domain Decomposition preconditioner that takes into account the physics of shear bands and the domain in which it forms, is proposed. The key idea is to resolve the shear band domain accurately while only solving the remaining domain approximately, with selective updates when necessary. The preconditioner is formed on the basis of an additive Schwartz method that is applied to a Schur complement partition of the system. The proposed preconditioner is implemented in parallel and tested on three different benchmark problems as compared against off-the-shelf solvers. We study the effect of h- and k-refinement that are obtained from Isogeometric discretizations of the system, the overlapping strategy and its parallel performance. Excellent performance is demonstrated in serial and parallel on all benchmark examples.

42 ENGINEERING↗

Colloidal quantum dots for optoelectronics

Colloidal quantum dots (QDs) are semiconductor nanocrystals that have unique size-tunable optoelectronic properties and are suitable for wet processing. QD research aims to answer fundamental questions about the chemical and physical properties of nanoscale materials and use these tools for technological applications ranging from bio-imaging to quantum optics. At the core of this field is a set of synthetic, processing and analytical methods designed to produce QDs in uniform ensembles that meet the highest performance standards. Here, this Primer reviews QD fabrication methods with a focus on the applications of QDs in printed optoelectronics and quantum optics. After outlining the current state-of-the-art QD syntheses, the experimental and computational analysis of QDs is discussed. These topics are then connected to the methodologies, processes and concepts required for developing QD-based photodetectors, light-emitting devices and quantum optics applications. Special attention is paid to challenges in reproducibility and current limitations of the field, such as the need to balance non-restricted material composition with high performing technology while achieving long-term stability in QD devices under operating conditions. Finally, the ongoing advancement in QD synthesis, precise atomic-level analysis and computational methodologies are highlighted as key drivers towards rational QD design, particularly in understanding how structural changes under loading impact QD properties.

optical materials↗

Grid Forming Fast Frequency Response for PMSG-Based Wind Turbines

Electric power systems are undergoing a rapid transition from fuel-based generation using synchronous generators to renewable generation interfaced by power electronics. In this context, a key challenge is using renewable generation with limited controllability to contribute to power system stability. In this work, we investigate grid-forming control of permanent magnet synchronous generator (PMSG) wind turbines. The proposed control and curtailment strategy supports the entire spectrum of standard functions of grid-following (e.g. maximum power point tracking (MPPT) and grid-forming control (e.g., primary frequency control) without explicit mode switching. A detailed case study is used to compare the performance of the proposed control at operating points corresponding to MPPT and frequency control with standard grid-forming and grid-following controls. Further, the results demonstrate that (i) the proposed energy balancing grid forming control is self-synchronizing in MPPT mode, (ii) the (limited) energy storage and controllability of wind turbines can be adequately utilized to provide grid support, and (ii) the proposed control exhibits good performance under variable wind speeds and does not result in a significant increase of wind turbine fatigue loads.

17 WIND ENERGY↗

Development of Training Programs to Optimize Planetary Ambulation

Astronauts experience disturbances in functional mobility following their return to Earth due to adaptive responses that occur during exposure to the microgravity conditions of space flight. Despite significant time spent performing in-flight exercise routines, these training programs have not been able to mitigate postflight alterations in postural and locomotor function. Therefore, the goal of our two inter-related projects (NSBRI-ground based and ISS flight study, "Mobility") is to develop and test gait training programs that will serve to optimize functional mobility during the adaptation period immediately following space flight, thereby improving the safety and efficiency of planetary ambulation. The gait training program entails manipulating the sensory conditions of treadmill exercise to systematically challenge the balance and gait control system. This enhances the overall adaptability of locomotor function enabling rapid reorganization of gait control to respond to ambulation in different gravitational environments. To develop the training program, we are conducting a series of ground-based studies evaluating the training efficacy associated with variation in visual flow, body loading, and support surface stability during treadmill walking. We will also determine the optimal method to present training stimuli within and across training sessions to maximize both the efficacy and efficiency of the training procedure. Results indicate that variations in both visual flow and body unloading during treadmill walking leads to modification in locomotor control and can be used as effective training modalities. Additionally, the composition and timing of sensory challenges experienced during each training session has significant impact on the ability to rapidly reorganize locomotor function when exposed to a novel sensory environment. We have developed the capability of producing support surface variation during gait training by mounting a treadmill on a six-degree-of-freedom motion device. This hardware development will allow us to evaluate the efficacy of this type of training in conjunction with variation in visual flow and body unloading.

Bloomberg, J. J.↗

Unnatural hypoxic regimes

Abstract Coastal hypoxia is increasing worldwide in response to human‐caused changes in global climate and biogeochemical cycles. In this paper, we view anthropogenic trends in coastal hypoxia through the lens of disturbance ecology and complexity theory. Complexity theory provides a framework for describing how estuaries and other coastal aquatic ecosystems respond to hypoxia by understanding feedback loops. Can it also be valuable in understanding how these ecosystems behave under shifting (i.e., unnatural) disturbance regimes? When viewed as a disturbance regime, shifts in the spatial (areal extent and fragmentation) and temporal (frequency and duration of events) characteristics of coastal hypoxia can be used to track changes into a non‐stationary future. Here, we consider options for increasing the resilience of coastal aquatic ecosystems to future, unnatural hypoxic regimes. To start, we define desirable states as ecosystems with long trophic chains and slow nutrient and carbon dynamics that produce many ecosystem services. We then work backward to describe circumstances dominated by positive feedbacks that can lead ecosystems toward an undesirable state (i.e., depauperate communities and chemically reduced sediments). Processes of degradation and recovery can be understood in the context of island biogeography whereby species diversity in habitats fragmented by hypoxia is determined by the balance between rapid local extinction, slow recolonization from the edges of hypoxic patches, and opportunities for ecological succession during between disturbance events. We review potential future changes associated with changing global climate and highlight ways to enhance coastal resilience. In addition to efforts to slow climate change, potential interventions include reduced nutrient and carbon loadings from rivers, restoration of aquatic vegetation, and managing for key species, including those that promote sediment oxygenation, that enhance water clarity, or that promote grazing on epiphytic algae through top‐down control.

54 ENVIRONMENTAL SCIENCES↗

A mixed-mode phase field fracture model in anisotropic rocks with consistent kinematics

Under a pure tensile loading, cracks in brittle, isotropic, and homogeneous materials often propagate such that pure mode I kinematics are maintained at the crack tip. However, experiments performed on geo-materials, such as sedimentary rock, shale, mudstone, concrete and gypsum, often lead to the conclusion that the mode I and mode II critical fracture energies/surface energy release rates are distinctive. This distinction has great influence on the formation and propagation of wing cracks and secondary cracks from pre-existing flaws under a combination of shear and tensile or shear and compressive loadings. To capture the mixed-mode fracture propagation, a mixed-mode I/II fracture model that employs multiple critical energy release rates based on Shen and Stephansson, IJRMMS, 1993 is reformulated in a regularized phase field fracture framework. Here we obtain the mixed-mode driving force of the damage phase field by balancing the microforce. Meanwhile, the crack propagation direction and the corresponding kinematics modes are determined via a local fracture dissipation maximization problem. Several numerical examples that demonstrate mode II and mixed-mode crack propagation in brittle materials are presented. Possible extensions of the model capturing degradation related to shear/compressive damage, as commonly observed in sub-surface applications and triaxial compression tests, are also discussed.

42 ENGINEERING↗