Search NASA⌕ Search

SEARCH · Search NASA

Results for “load pattern modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

Usage Pattern Analysis for the Summit Login Nodes

High performance computing (HPC) users interact with Summit through dedicated gateways, also known as login nodes. The performance and stability of these login nodes can have a significant impact on the user experience. In this study, the performance and stability of Summit’s five login nodes are evaluated by analyzing the log data from 2020 and 2021. The analysis focuses on the computing capability (CPU average load, users and tasks) and the storage performance, along with the associated job scheduler activity. The outcome of this study can serve as the foundation of a predictive modeling framework that enables the system admin of an HPC system to preemptively deploy countermeasures before the onset of a system failure.

Eiffert, Brett↗

Elastically-homogeneous lattice models of damage in geomaterials

This paper involves the development of the auxiliary stress approach for producing elastically-homogeneous lattice models of damage in geomaterials. The lattice models are based on random, three-dimensional assemblages of rigid-body-spring elements. Unlike conventional lattice or particle models, the elastic constants of a material (e.g., Young’s modulus and Poisson’s ratio) are represented properly in both global and local senses, without any need for calibration. The proposed approach is demonstrated and validated through analyses of homogeneous and heterogeneous systems under uni- and tri-axial loading conditions. Comparisons are made with analytical solutions and finite element results. Thereafter, the model is used to simulate a series of standard laboratory tests: (a) split-cylinder tests, and (b) uniaxial compressive tests of sedimentary rocks at the Horonobe Underground Research Laboratory in Hokkaido, Japan. Model inputs are based on physical quantities measured in the experiments. Finally, the simulation results agree well with the experimental results in terms of pre-peak stress-strain/displacement responses, strength measurements, and failure patterns.

3D irregular lattice↗

Vertical meteorological patterns and their impact on the energy demand of tall buildings

In accordance with the developing economy and growing population, an increasing number of tall buildings have been constructed over the last 20 years. In 2017, there were 144 new buildings worldwide that were 200 m or greater in height; 50% of these buildings were in China. Due to the vertical gradient impact of meteorological parameters, the energy performance of tall buildings differs from that of general buildings. Few studies exist on vertical meteorological changes using measured data at different heights. Most studies on dynamic energy simulation simulate meteorological parameters using models. This study explores vertical meteorological patterns using hourly dry-bulb temperature, relative humidity, and wind speed data from 2007 to 2017 for a 325 m meteorological tower in Beijing. The temperature decreased significantly with increasing altitude (~0.9 °C per 100 m), while the daily temperature difference decreased with increasing altitude. The absolute humidity did not change significantly with height. The wind speed increased with altitude at approximately 2 m/s per 100 m. The building simulation showed that the annual heating load at a height of 320 m increased by 85% from that at 8 m; the annual cooling load decreased by 20%. Such significant differences in thermal loads for 300-m-tall buildings confirm the need to consider vertical meteorological differences in building performance simulations for tall buildings. A greater number of measurement points at different heights improve the simulation accuracy. Finally, guidance on selecting the heights for vertical metrological measurements is provided based on the influences of building thermal loads.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Domain Structures of Ni and NiFe (Oxy)Hydroxide Oxygen-Evolution Catalysts from X-ray Pair Distribution Function Analysis

Ni-Fe (oxy)hydroxides, Ni( i-z )Fe z O x H y , are among the fastest-known water oxidation catalysts in alkaline media on a per-cation basis. At current densities relevant for electrolysis (e.g., >0.5 A/cm( -2 )), mass and electron transport through catalyst films with high mass loading are critical and depend substantially on the extended and intermediate catalyst architecture. We use X-ray pair distribution function (PDF) analysis to determine the intermediate nanostructures of electrodeposited Ni( i-z )Fe z O x H y films. We report the effects of electrodeposition technique (pulsed versus continuous), electrochemical cycling, and Fe content on the structure of the catalyst film. The PDF patterns for Ni( i-z )Fe z O x H y films are best simulated by model structures consisting of brucite-like β-Ni(OH) 2 fragments 1 to 3 layers in thickness. Only the oxidation state of the film significantly affects the intralayer scattering behavior (i.e., metal-oxygen bond distance). The interlayer interactions, however, are affected by Fe content and deposition conditions. The domain size of many of the systems are similar, extending to similar to ~5 nm, which are best modeled by sheets containing upward of similar to ~250 metal cations. Smaller domains were found for films deposited through a larger number of electrochemical cathodic current pulses. Films can be cycled between as-deposited, oxidized, and reduced states, with minimal loss of intrasheet coherence, indicating a degree of structural stability. We estimate heterogeneity in the domain structures by modeling the PDF data to linear combinations of oxyhydroxide fragments with different sizes and numbers of layers.

36 MATERIALS SCIENCE↗

Hardware-Based Emulator with Deep Learning Model for Building Energy Control and Prediction Based on Occupancy Sensors’ Data

Heating, ventilation, and air conditioning (HVAC) is the largest source of residential energy consumption. Occupancy sensors’ data can be used for HVAC control since it indicates the number of people in the building. HVAC and sensors form a typical cyber-physical system (CPS). In this paper, we aim to build a hardware-based emulation platform to study the occupancy data’s features, which can be further extracted by using machine learning models. In particular, we propose two hardware-based emulators to investigate the use of wired/wireless communication interfaces for occupancy sensor-based building CPS control, and the use of deep learning to predict the building energy consumption with the sensor data. We hypothesize is that the building energy consumption may be predicted by using the occupancy data collected by the sensors, and question what type of prediction model should be used to accurately predict the energy load. Another hypothesis is that an in-lab hardware/software platform could be built to emulate the occupancy sensing process. The machine learning algorithms can then be used to analyze the energy load based on the sensing data. To test the emulator, the occupancy data from the sensors is used to predict energy consumption. The synchronization scheme between sensors and the HVAC server will be discussed. We have built two hardware/software emulation platforms to investigate the sensor/HVAC integration strategies, and used an enhanced deep learning model—which has sequence-to-sequence long short-term memory (Seq2Seq LSTM)—with an attention model to predict the building energy consumption with the preservation of the intrinsic patterns. Because the long-range temporal dependencies are captured, the Seq2Seq models may provide a higher accuracy by using LSTM architectures with encoder and decoder. Meanwhile, LSTMs can capture the temporal and spatial patterns of time series data. The attention model can highlight the most relevant input information in the energy prediction by allocating the attention weights. The communication overhead between the sensors and the HVAC control server can also be alleviated via the attention mechanism, which can automatically ignore the irrelevant information and amplify the relevant information during CNN training. Our experiments and performance analysis show that, compared with the traditional LSTM neural network, the performance of the proposed method has a 30% higher prediction accuracy.

Ye, Zhijing↗

Space Shuttle ice prevention studies

Turbojet engine exhausts are being considered for the prevention of ice formation on the Space Shuttle External Tank (ET) when cryogenic propellants are loaded. The ice could dislodge and damage the Thermal Protection System tiles of the Shuttle Orbiter. An oil flow visualization technique was developed to quickly and inexpensively study the flowfields created by jet exhausts. The procedure can also determine the influence of geometry, jet exhausts, and wind conditions on the size and location of flow regions. A black tempra and kerosene mixture was applied to an ET model; the wind tunnel and jet flow were started; and when they stopped the flow patterns were photographed to determine which parameters had the greatest influence on the flowfields.

Porteiro, J. L. F.↗

Evidence for chaotic fault interactions in the seismicity of the San Andreas fault and Nankai trough

The dynamical behavior introduced by fault interactions is examined here using a simple spring-loaded, slider-block model with velocity-weakening friction. The model consists of two slider blocks coupled to each other and to a constant-velocity driver by elastic springs. For an asymmetric system in which the frictional forces on the two blocks are not equal, the solutions exhibit chaotic behavior. The system's behavior over a range of parameter values seems to be generally analogous to that of weakly coupled segments of an active fault. Similarities between the model simulations and observed patterns of seismicity on the south central San Andreas fault in California and in the Nankai trough along the coast of southwestern Japan.

Huang, Jie↗

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing during Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on numerical simulations with the NASA Global Modeling and Assimilation Office finite-volume general circulation model (fvGCM) with Microphyics of Clouds with the Relaxed Arakawa Schubert Scheme (McRAS), using aerosol forcing functions derived from the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The authors find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excite a planetary-scale teleconnection pattern in sea level pressure, temperature, and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, that is, South Asia, East Asia, and northern and western Africa. Significant atmospheric heating is found in regions with large loading of dust (over northern Africa and the Middle East) and black carbon (over Southeast Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east west dipole anomaly with strong cooling over the Caspian Sea and warming over central and northeastern Asia, where aerosol concentrations are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection pattern driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations. The surface temperature signature associated with the aerosol-induced teleconnection bears striking resemblance to the spatial pattern of observed long-term trend in surface temperature over Eurasia. Additionally, the boreal spring wave train pattern is similar to that reported by Fukutomi et al. associated with the boreal summer precipitation seesaw between eastern and western Siberia. The results of this study raise the possibility that global aerosol forcing during boreal spring may play an important role in spawning atmospheric teleconnections that affect regional and global climates.

Kim, Maeng-Ki↗

White Paper: Research & Development for the Time at Temperature Approach

Recent advancements in nuclear power research are greatly improving reactor safety and performance through the development of Accident Tolerant Fuel (ATF) and Low-Enriched Uranium Plus (LEU+). These innovations can address Departure from Nucleate Boiling (DNB) margins, which are vital for reactor safety. DNB happens when the coolant switches to film boiling, significantly decreasing heat transfer and posing a risk of fuel cladding failure. The U.S. Nuclear Regulatory Commission (NRC) employs conservative DNB criteria, which can potentially restrict the operational flexibility and efficiency of reactors. The Time at Temperature (TaT) approach could provide a more detailed and adaptable operational guideline by establishing acceptable time-temperature limits, accounting for the duration a material can withstand elevated temperatures without losing its integrity. This method allows reactors to operate more efficiently and safely, offering additional operational margins, faster power adjustments, and improved fuel cycle economics. TaT criteria allow for higher power levels and more flexible responses to operational transients, particularly applicable for anticipated operational occurrences (AOOs) that result in short durations of post-DNB conditions. It enhances plant operational flexibility, allows faster startup times, and enables quicker power level adjustments, optimizing fuel loading patterns and improving fuel cycle economics. Implementing TaT limits reduces core design constraints, lowers fuel usage, and reduces costs, essential for the long-term sustainability of Light Water Reactors (LWRs). TaT maximizes the use of advanced fuel technologies like ATF and LEU+, further enhancing their economic and environmental benefits. To apply the TaT approach in existing LWRs, collaborative research activities among various DOE-sponsored programs are essential. These efforts should incorporate fuel experiments, physics-based high-fidelity modeling, ML-based surrogate modeling, and optimization techniques. This whitepaper proposes four research and development areas: 1) Investigation of the feasibility of new operations of LWR with updated safety limits; 2) Assessment of reactor operation limits through uncertainty reduction; 3) Evaluation of power uprate in virtual environment; and 4) Lattice and reactor core design for power uprate. Each area includes why this research is in need and a suggested scope of work. These comprehensive research areas ensure practical and beneficial advancements for existing reactors, translating innovations in nuclear fuel and cladding technology into improved reactor performance and safety.

42 - ENGINEERING↗

Variational and phase response analysis for limit cycles with hard boundaries, with applications to neuromechanical control problems

Motor systems show an overall robustness, but because they are highly nonlinear, understanding how they achieve robustness is difficult. In many rhythmic systems, robustness against perturbations involves response of both the shape and the timing of the trajectory. This makes the study of robustness even more challenging. To understand how a motor system produces robust behaviors in a variable environment, we consider a neuromechanical model of motor patterns in the feeding apparatus of the marine mollusk Aplysia californica (Shaw et al. in J Comput Neurosci 38(1):25–51, 2015; Lyttle et al. in Biol Cybern 111(1):25–47, 2017). We established in (Wang et al. in SIAM J Appl Dyn Syst 20(2):701–744, 2021. https://doi.org/10.1137/20M1344974) the tools for studying combined shape and timing responses of limit cycle systems under sustained perturbations and here apply them to study robustness of the neuromechanical model against increased mechanical load during swallowing. Interestingly, we discover that nonlinear biomechanical properties confer resilience by immediately increasing resistance to applied loads. In contrast, the effect of changed sensory feedback signal is significantly delayed by the firing rates’ hard boundary properties. Our analysis suggests that sensory feedback contributes to robustness in swallowing primarily by shifting the timing of neural activation involved in the power stroke of the motor cycle (retraction). This effect enables the system to generate stronger retractor muscle forces to compensate for the increased load, and hence achieve strong robustness. The approaches that we are applying to understanding a neuromechanical model in Aplysia, and the results that we have obtained, are likely to provide insights into the function of other motor systems that encounter changing mechanical loads and hard boundaries, both due to mechanical and neuronal firing properties.

59 BASIC BIOLOGICAL SCIENCES↗

Acoustic measurements of three Prop-Fan models

Results of NASA sponsored acoustic tests of three 2 ft. diameter models of the Prop-Fan (a small diameter, highly loaded, many-bladed variable pitch advanced turboprop) are presented. The highly swept model designed for noise reduction produces substantially less near field noise at simulated 0.8 Mach number cruise conditions than the unswept or slightly swept models. It also produces less far field noise at conditions simulating takeoff and landing. The noise reduction mechanism is discussed. Correlation between harmonic noise measurements and theoretical predictions and between measured and predicted acoustic pressure pulses is good. Shadowgraph measurements which show the location of blade associated wave patterns were obtained. Predicted and measured wave locations show good general agreement. Full scale near and far field noise is predicted.

Brooks, B. M.↗

Numerical Study of Stratified Charge Combustion in Wave Rotors

A wave rotor may be used as a pressure-gain combustor effecting non-steady flow, and intermittent, confined combustion to enhance gas turbine engine performance. It will be more compact and probably lighter than an equivalent pressure-exchange wave rotor, yet will have similar thermodynamic and mechanical characteristics. Because the allowable turbine blade temperature limits overall fuel/air ratio to sub-flammable values, premixed stratification techniques are necessary to burn hydrocarbon fuels in small engines with compressor discharge temperature well below autoignition conditions. One-dimensional, unsteady numerical simulations of stratified-charge combustion are performed using an eddy-diffusivity turbulence model and a simple reaction model incorporating a flammability limit temperature. For good combustion efficiency, a stratification strategy is developed which concentrates fuel at the leading and trailing edges of the inlet port. Rotor and exhaust temperature profiles and performance predictions are presented at three representative operating conditions of the engine: full design load, 40% load, and idle. The results indicate that peak local gas temperatures will result in excessive temperatures within the rotor housing unless additional cooling methods are used. The rotor itself will have acceptable temperatures, but the pattern factor presented to the turbine may be of concern, depending on exhaust duct design and duct-rotor interaction.

Nalim, M. Razi↗

Simulation of Stochastic Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis–Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button–sized disk–shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300oC to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

residual tensile stress↗

Micromechanical combined stress analysis: MICSTRAN, a user manual

Composite materials are currently being used in aerospace and other applications. The ability to tailor the composite properties by the appropriate selection of its constituents, the fiber and matrix, is a major advantage of composite materials. The Micromechanical Combined Stress Analysis (MICSTRAN) code provides the materials engineer with a user-friendly personal computer (PC) based tool to calculate overall composite properties given the constituent fiber and matrix properties. To assess the ability of the composite to carry structural loads, the materials engineer also needs to calculate the internal stresses in the composite material. MICSTRAN is a simple tool to calculate such internal stresses with a composite ply under combined thermomechanical loading. It assumes that the fibers have a circular cross-section and are arranged either in a repeating square or diamond array pattern within a ply. It uses a classical elasticity solution technique that has been demonstrated to calculate accurate stress results. Input to the program consists of transversely isotropic fiber properties and isotropic matrix properties such as moduli, Poisson's ratios, coefficients of thermal expansion, and volume fraction. Output consists of overall thermoelastic constants and stresses. Stresses can be computed under the combined action of thermal, transverse, longitudinal, transverse shear, and longitudinal shear loadings. Stress output can be requested along the fiber-matrix interface, the model boundaries, circular arcs, or at user-specified points located anywhere in the model. The MICSTRAN program is Windows compatible and takes advantage of the Microsoft Windows graphical user interface which facilitates multitasking and extends memory access far beyond the limits imposed by the DOS operating system.

Naik, R. A.↗

A Climatology of Midlatitude Continental Clouds from the ARM SGP Site: Low-Level Cloud Macrophysical, Microphysical, and Radiative Properties - Part I

A record of single-layer and overcast low cloud (stratus) properties has been generated using approximately 4000 hours of data collected from January 1997 to December 2002 at the Atmospheric Radiation Measurement (ARM) Southern Great Plains Central Facility (SCF). The cloud properties include liquid-phase and liquid-dominant, mixed-phase, low cloud macrophysical, microphysical, and radiative properties including cloud-base and -top heights and temperatures, and cloud physical thickness derived from a ground-based radar and lidar pair, and rawinsonde sounding; cloud liquid water path (LWP) and content (LWC), and cloud-droplet effective radius (r(sub e)) and number concentration (N) derived from the macrophysical properties and radiometer data; and cloud optical depth (tau), effective solar transmission (gamma), and cloud/top-of-atmosphere albedos (R(sub cldy)/R(sub TOA)) derived from Eppley precision spectral pyranometer measurements. The cloud properties were analyzed in terms of their seasonal, monthly, and hourly variations. In general, more stratus clouds occur during winter and spring than in summer. Cloud-layer altitudes and physical thicknesses were higher and greater in summer than in winter with averaged physical thicknesses of 0.85 km and 0.73 km for day and night, respectively. The seasonal variations of LWP, LWC, N. tau, R(sub cldy), and R(sub TOA) basically follow the same pattern with maxima and minima during winter and summer, respectively. There is no significant variation in mean r(sub e), however, despite a summertime peak in aerosol loading, Although a considerable degree of variability exists, the 6-yr average values of LWP, LWC, r(sub e), N, tau, gamma, R(sub cldy) and R(sub TOA) are 150 gm(exp -2) (138), 0.245 gm(exp -3) (0.268), 8.7 micrometers (8.5), 213 cm(exp -3) (238), 26.8 (24.8), 0.331, 0.672, 0.563 for daytime (nighttime). A new conceptual model of midlatitude continental low clouds at the ARM SGP site has been developed from this study. The low stratus cloud amount monotonically increases from midnight to early morning (0930 LT), and remains large until around local noon, then declines until 1930 LT when it levels off for the remainder of the night. In the morning, the stratus cloud layer is low, warm, and thick with less LWC, while in the afternoon it is high, cold, and thin with more LWC. Future parts of this series will consider other cloud types and cloud radiative forcing at the ARM SCF.

Dong, Xiquan↗

Framework for Evaluating Grid Readiness for Integrating Electric Vehicles and Utility Scale Storage in India

Electric vehicles (EVs) and distributed solar are poised to grow substantially in India in the coming years following ambitious Government of India targets. The localized impact of these changes, both in terms of infrastructure investments and demand patterns, are not well understood, nor are tools for this type of analysis readily available. This paper describes the development of a framework that can analyze network readiness in terms of feeder impacts for distributed energy resources (DERs) including growing EV penetration, and the potential solutions introduced by utility-scale battery energy storage systems (BESSs). The building blocks of the feeder analysis require multiple data sets to be compiled into a usable network model. The work is driven by feeder head and distribution transformer (DT) loading data along with all the technical specifications and schematics of distribution feeders. These measured data sets are conditioned to remove any afflictions, and the cleaned load profiles are used to perform multi-year quasi-static time-series power flow analyses on detailed three-phase feeder models. EVs translated to grid-tied loads are represented with aggregated EV demand profile estimates based on possible charging scenarios in residential sector, and in workplace and/or public charging stations. Finally, BESSs are evaluated for their cost-effectiveness in the framework through an initial screening that includes the benefits of grid service applications to mitigate present and future overloading scenarios. The outcomes from this framework are expected to help utilities gauge the readiness of their distribution grid for integrating an increasing number of EVs as well as load growth during a projected multi-year time horizon, where BESSs can contribute in making the grid more reliable.

electric vehicles↗

A multi-level load shape clustering and disaggregation approach to characterize patterns of energy consumption behavior

This study presents representative electrical load shapes, disaggregated to the end-use level, for over 5000 customer clusters across California’s residential, commercial, industrial and agricultural sectors. We developed a novel, multi-level load shape clustering approach for residential and commercial sectors leveraging interval meter data for over 350,000 California utility customers collected as a part of the Phase 4 California Demand Response (DR) Potential Study. The clustering approach allowed us to identify typical consumption patterns and categorize customers based on their daily load shape displayed throughout the year. For example, we were able to identify customers with particular energy technologies such as electric vehicles and rooftop solar, as well as building occupancy types such as restaurants, grocery stores and even unoccupied buildings, based solely on whole-building interval data. We then combined the load shape-based clusters with other customer information including building type, climate, geographical area, total consumption and low-income status, to create a set of customer clusters based on both demographics and usage patterns. Total cluster electricity demand was then disaggregated into a wide variety of end-uses using weather normalization and other publicly available end-use load shape datasets. The resulting disaggregated cluster load shapes will be released in anonymized form as part of the Phase 4 DR Potential Study. They will have wide-ranging applications in energy research and policy analysis, including estimation of energy efficiency (EE) and DR potential on the end-use level, time-dependent valuation of EE savings, building stock modeling, and developing customer targeting strategies for EE and DR programs.

Murthy, Samanvitha↗

Skin-Stiffener Debond Prediction Based on Computational Fracture Analysis

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that is reinforced with stringers. Shear loading causes the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. Across the width of the stringer foot the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. For small applied loads the failure index is well below one across the entire width. With increasing load the failure index approaches one first near the edge of the stringer foot from which delamination is expected to grow. With increasing delamination lengths the buckling pattern of the panel changes and the failure index increases which suggests that rapid delamination growth from the initial defect is to be expected.

Krueger, Ronald↗