SEARCH · Engineering Papers

Results for “Closure”

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.

SourceResearch connectionFrom yearThrough year

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

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Investigations of turbulent scalar fields using probability density function approach

Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Gao, Feng↗

NASA NTRS1993

Cooled spool piston compressor

A hydraulically powered gas compressor receives low pressure gas and outputs a high pressure gas. The housing of the compressor defines a cylinder with a center chamber having a cross-sectional area less than the cross-sectional area of a left end chamber and a right end chamber, and a spool-type piston assembly is movable within the cylinder and includes a left end closure, a right end closure, and a center body that are in sealing engagement with the respective cylinder walls as the piston reciprocates. First and second annual compression chambers are provided between the piston enclosures and center housing portion of the compressor, thereby minimizing the spacing between the core gas and a cooled surface of the compressor. Restricted flow passageways are provided in the piston closure members and a path is provided in the central body of the piston assembly, such that hydraulic fluid flows through the piston assembly to cool the piston assembly during its operation. The compressor of the present invention may be easily adapted for a particular application, and is capable of generating high gas pressures while maintaining both the compressed gas and the compressor components within acceptable temperature limits.

Morris, Brian G.↗

NASA NTRS1993

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗

NASA NTRS1993

Flow interaction experiment. Volume 1: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

NASA NTRS1993

Flow interaction experiment. Volume 2: Aerothermal modeling, phase 2

Nikjooy, M.↗

NASA NTRS1995

Analytical and phenomenological studies of rotating turbulence

A framework, which combines mathematical analysis, closure theory, and phenomenological treatment, is developed to study the spectral transfer process and reduction of dimensionality in turbulent flows that are subject to rotation. First, we outline a mathematical procedure that is particularly appropriate for problems with two disparate time scales. The approach which is based on the Green's method leads to the Poincare velocity variables and the Poincare transformation when applied to rotating turbulence. The effects of the rotation are now reflected in the modifications to the convolution of a nonlinear term. The Poincare transformed equations are used to obtain a time-dependent analog of the Taylor-Proudman theorem valid in the asymptotic limit when the non-dimensional parameter mu is identical to Omega(t) approaches infinity (Omega is the rotation rate and t is the time). The 'split' of the energy transfer in both direct and inverse directions is established. Secondly, we apply the Eddy-Damped-Quasinormal-Markovian (EDQNM) closure to the Poincare transformed Euler/Navier-Stokes equations. This closure leads to expressions for the spectral energy transfer. In particular, an unique triple velocity decorrelation time is derived with an explicit dependence on the rotation rate. This provides an important input for applying the phenomenological treatment of Zhou. In order to characterize the relative strength of rotation, another non-dimensional number, a spectral Rossby number, which is defined as the ratio of rotation and turbulence time scales, is introduced. Finally, the energy spectrum and the spectral eddy viscosity are deduced.

Mahalov, Alex↗

NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗

NASA NTRS1994

A Review Of Modelling Small-Crack Behavior And Fatigue-Life Predictions For Aluminum Alloys

The small-crack effect, where small fatigue cracks grow faster and at lower stress-intensity factors than large cracks, has been found to be significant for many materials and loading conditions. In this paper, plasticity effects and crack-closure modelling of small fatigue cracks are reviewed. A crack-closure model with a cyclic-plastic zone-corrected effective stress-intensity factor range (related to the cyclic J-integral) and microstructural data on crack-initiation sites were used to calculate small-crack growth rates and fatigue lives for unnotched and notched specimens made of two aluminum alloys. The crack-closure transient from the plastic wake was shown to be the dominant cause of the small-crack effect and plasticity effects on the cyclic-plastic zone-corrected stress-intensity factor range were negligible except at extremely high stress levels. Small-crack growth rates and fatigue lives under both constant-amplitude and spectrum loading from tests and analyses agreed well.

Newman, J. C., Jr.↗

NASA NTRS2000

North Atlantic Aerosol Radiative Impacts Based on Satellite Measurements and Aerosol Intensive Properties from TARFOX and ACE-2

We estimate the impact of North Atlantic aerosols on the net short-wave flux at the tropopause by combining satellite-derived aerosol optical depth (AOD) maps with model aerosol properties determined via closure analyses in TARFOX and ACE 2. We exclude African dust, primarily by restricting latitudes to 25-60 N. The analyses use in situ aerosol composition measurements and air- and ship-borne sun-photometer measurements of AOD spectra. The aerosol model yields computed flux sensitivities (dFlux/dAOD) that agree with measurements by airborne flux radiometers in TARFOX. Its midvisible single-scattering albedo is 0.9. which is in the range obtained from in situ measurements of scattering and absorption in both TARFOX and ACE 2. Combining satellite-derived AOD maps with the aerosol model yields maps of 24-hour average net radiative flux changes. For simultaneous AVHRR, radiance measurements exceeded the sunphotometer AODs by about 0.04. However. shipboard sunphotometer and AVHRR AODs agreed Within 0.02 for data acquired during satellite overflights on two other days. We discuss attempts to demonstrate column closure within the MBL by comparing shipboard sunphotometer AODs and values calculated from simultaneous shipboard in-situ aerosol size distribution measurements. These comparisons were mostly unsuccessful, but they illustrate the difficulties inherent in this type of closure analysis. Specifically, AODs derived from near-surface in-situ size distribution measurements are extremely sensitive to the assumed hygroscopic growth model that itself requires an assumption of particle composition as a function of height and size, to the radiosonde-measured relative humidity, and to the vertical profile of particle number. We investigate further the effects of hygroscopic particle growth within the MBL by using shipboard lidar aerosol backscatter profiles together with the sunphotometer AOD.

Russell, P. B.↗

NASA NTRS1997

Direct Aerosol Radiative Forcing: Calculations and Measurements from the Tropospheric

Radiative forcing is defined as the change in the net (downwelling minus upwelling) radiative flux at a given level in the atmosphere. This net flux is the radiative power density available to drive climatic processes in the earth-atmosphere system below that level. Recent research shows that radiative forcing by aerosol particles is a major source of uncertainty in climate predictions. To reduce those uncertainties, TARFOX was designed to determine direct (cloud-free) radiative forcing by the aerosols in one of the world's major industrial pollution plumes--that flowing from the east coast of the US over the Atlantic Ocean. TARFOX measured a variety of aerosol radiative effects (including direct forcing) while simultaneously measuring the chemical, physical, and optical properties of the aerosol particles causing those effects. The resulting data sets permit a wide variety of tests of the consistency, or closure, among the measurements and the models that link them. Because climate predictions use the same or similar model components, closure tests help to assess and reduce prediction uncertainties. In this work we use the TARFOX-determined aerosol, gas, and surface properties to compute radiative forcing for a variety of aerosol episodes, with inadvisable optical depths ranging from 0.07 to 0.6. We calculate forcing by several techniques with varying degrees of sophistication, in part to test the range of applicability of simplified techniques--which are often the only ones feasible in climate predictions by general circulation models (GCMs). We then compare computed forcing to that determined from: (1) Upwelling and downwelling fluxes (0.3-0.7 mm and 0.7-3.0 mm) measured by radiometers on the UK MRF C-130. and (2) Daily average cloud-free absorbed solar and emitted thermal radiative flux at the top of the atmosphere derived from the AVHRR radiometer on the NOAA- 14 satellite. The calculations and measurements all yield aerosol direct radiative forcing in the range -50 to -190 W sq m per unit inadvisable optical depth. The magnitudes are about 15 to 100 times larger than the global-average direct forcing expected for the global-average sulfate aerosol optical depth of 0.04. The reasons for the larger forcing in TARFOX include the relatively large optical depths and the focus on cloud-free, daytime conditions over the dark ocean surface. These are the conditions that produce the actual major radiative forcing events that contribute to any global-average climate effect. Detailed comparisons of calculated and measured forcings for specific events are used for more refined tests of closure.

Russell, P. B.↗

NASA NTRS2002

Fluid Mechanics of Capillary-Elastic Instabilities in Microgravity Environment

The aim of this project is to investigate the closure and reopening of lung airways due to surface tension forces, coupled with airway elasticity. Airways are liquid-lined, flexible tubes and closure of airways can occur by a Rayleigh instability of the liquid lining, or an instability of the elastic support for the airway as the surface tension of the air-liquid interface pulls the tube shut, or both. Regardless of the mechanism, the airway is closed because the liquid lining has created a plug that prevents axial gas exchange. In the microgravity environment, surface tension forces dominate lung mechanics and would lead to more prevalent, and more uniformly distributed air-way closure, thereby creating a potential for respiratory problems for astronauts. Once closed the primary option for reopening an airway is by deep inspiration. This maneuver will pull the flexible airways open and force the liquid plug to flow distally by the incoming air stream. Airway reopening depends to a large extent on this plug flow and how it may lead to plug rupture to regain the continuity of gas between the environment and the alveoli. In addition to mathematical modeling of plug flows in liquid-lined, flexible tubes, this work has involved benchtop studies of propagating liquid plugs down tube networks that mimic the human airway tree. We have extended the work to involve animal models of liquid plug propagation in rat lungs. The liquid is radio-opaque and x-ray video imaging is used to ascertain the movement and distribution of the liquid plugs so that comparisons to theory may be made. This research has other uses, such as the delivery of liquids or drugs into the lung that may be used for surfactant replacement therapy or for liquid ventilation.

Grotberg, James B.↗

NASA NTRS2004

Apparent I-Xe Cooling Rates of Chondrules Compared with Silicates from the Colomera Iron Meteorite

In I-Xe dating, a regular pattern of increasing Xe-129/128Xe-128 ratio with increasing extraction temperature is often observed. If one makes the crude assumption that the temperatures at which the Xe is extracted in the laboratory is approximately the same as the temperature at which those sites closed 4.6 Ga ago, a (zeroth order) model cooling rate can be found. In order to test and refine this model we can apply the cooling theory of Dodson to those extraction steps approaching the I-Xe isochrons. Using an Arrhenius plot for these temperature fractions, and assuming that an only single phase is involved, the effective diffusion parameters can be estimated (frequency factor and activation energy). From the apparent (zeroth order) cooling rate, the closure temperature can be estimated from the Dodson equation. This model closure temperature can then be compared with the actual laboratory temperature at which the isochron begins. The ratio of the closure temperature and the temperature corresponding to the start of the isochron provides the ratio of the two temperature scales, incorporation and extraction. The actual cooling rate is then given by the apparent (zeroth order) cooling rate times the temperature scale factor. Figure 1 shows Arrhenius plots for I-Xe data

Hohenberg, C. M.↗

NASA NTRS2007

Three Dimensional Constraint Effects on the Estimated (Delta)CTOD during the Numerical Simulation of Different Fatigue Threshold Testing Techniques

Variation in constraint through the thickness of a specimen effects the cyclic crack-tip-opening displacement (DELTA CTOD). DELTA CTOD is a valuable measure of crack growth behavior, indicating closure development, constraint variations and load history effects. Fatigue loading with a continual load reduction was used to simulate the load history associated with fatigue crack growth threshold measurements. The constraint effect on the estimated DELTA CTOD is studied by carrying out three-dimensional elastic-plastic finite element simulations. The analysis involves numerical simulation of different standard fatigue threshold test schemes to determine how each test scheme affects DELTA CTOD. The American Society for Testing and Materials (ASTM) prescribes standard load reduction procedures for threshold testing using either the constant stress ratio (R) or constant maximum stress intensity (K(sub max)) methods. Different specimen types defined in the standard, namely the compact tension, C(T), and middle cracked tension, M(T), specimens were used in this simulation. The threshold simulations were conducted with different initial K(sub max) values to study its effect on estimated DELTA CTOD. During each simulation, the DELTA CTOD was estimated at every load increment during the load reduction procedure. Previous numerical simulation results indicate that the constant R load reduction method generates a plastic wake resulting in remote crack closure during unloading. Upon reloading, this remote contact location was observed to remain in contact well after the crack tip was fully open. The final region to open is located at the point at which the load reduction was initiated and at the free surface of the specimen. However, simulations carried out using the constant Kmax load reduction procedure did not indicate remote crack closure. Previous analysis results using various starting K(sub max) values and different load reduction rates have indicated DELTA CTOD is independent of specimen size. A study of the effect of specimen thickness and geometry on the measured DELTA CTOD for various load reduction procedures and its implication in the estimation of fatigue crack growth threshold values is discussed.

Seshadri, Banavara R.↗

NASA NTRS2004

Single-column Model Intercomparison for a Stably Stratified Atmospheric Boundary Layer

The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, which has a large impact on the medium-range weather forecasts and on climate integrations. A non-strongly stratified arctic case is simulated by nineteen single-column turbulence schemes. The statistics from the Large-eddy simulation (LES) intercomparison made for the same case by eight different models are used as a guiding reference. The single-column parameterizations include research schemes and operational schemes from major forecast and climate research centres. First order schemes, a large number of turbulence kinetic energy closures, and other proposals have submitted results. There is a large spread in the results; in general, the operational schemes mix more efficiently than the research ones, and the TKE and other higher order closures give results closer to the LES statistics. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.

Cuxart, J.↗

NASA NTRS2008

Energetic Consistency and Coupling of the Mean and Covariance Dynamics

The dynamical state of the ocean and atmosphere is taken to be a large dimensional random vector in a range of large-scale computational applications, including data assimilation, ensemble prediction, sensitivity analysis, and predictability studies. In each of these applications, numerical evolution of the covariance matrix of the random state plays a central role, because this matrix is used to quantify uncertainty in the state of the dynamical system. Since atmospheric and ocean dynamics are nonlinear, there is no closed evolution equation for the covariance matrix, nor for the mean state. Therefore approximate evolution equations must be used. This article studies theoretical properties of the evolution equations for the mean state and covariance matrix that arise in the second-moment closure approximation (third- and higher-order moment discard). This approximation was introduced by EPSTEIN [1969] in an early effort to introduce a stochastic element into deterministic weather forecasting, and was studied further by FLEMING [1971a,b], EPSTEIN and PITCHER [1972], and PITCHER [1977], also in the context of atmospheric predictability. It has since fallen into disuse, with a simpler one being used in current large-scale applications. The theoretical results of this article make a case that this approximation should be reconsidered for use in large-scale applications, however, because the second moment closure equations possess a property of energetic consistency that the approximate equations now in common use do not possess. A number of properties of solutions of the second-moment closure equations that result from this energetic consistency will be established.

Cohn, Stephen E.↗

NASA NTRS2006

Submicrosecond Power-Switching Test Circuit

A circuit that changes an electrical load in a switching time shorter than 0.3 microsecond has been devised. This circuit can be used in testing the regulation characteristics of power-supply circuits . especially switching power-converter circuits that are supposed to be able to provide acceptably high degrees of regulation in response to rapid load transients. The combination of this power-switching circuit and a known passive constant load could be an attractive alternative to a typical commercially available load-bank circuit that can be made to operate in nominal constant-voltage, constant-current, and constant-resistance modes. The switching provided by a typical commercial load-bank circuit in the constant-resistance mode is not fast enough for testing of regulation in response to load transients. Moreover, some test engineers do not trust the test results obtained when using commercial load-bank circuits because the dynamic responses of those circuits are, variously, partly unknown and/or excessively complex. In contrast, the combination of this circuit and a passive constant load offers both rapid switching and known (or at least better known) load dynamics. The power-switching circuit (see figure) includes a signal-input section, a wide-hysteresis Schmitt trigger that prevents false triggering in the event of switch-contact bounce, a dual-bipolar-transistor power stage that drives the gate of a metal oxide semiconductor field-effect transistor (MOSFET), and the MOSFET, which is the output device that performs the switching of the load. The MOSFET in the specific version of the circuit shown in the figure is rated to stand off a potential of 100 V in the "off" state and to pass a current of 20 A in the "on" state. The switching time of this circuit (the characteristic time of rise or fall of the potential at the drain of the MOSFET) is .300 ns. The circuit can accept any of three control inputs . which one depending on the test that one seeks to perform: a repetitive waveform from a signal generator, momentary closure of a push-button switch, or closure or opening of a manually operated on/off switch. In the case of a signal generator, one can adjust the frequency and duty cycle as needed to obtain the desired AC power-supply response, which one could display on an oscilloscope. Momentary switch closure could be useful for obtaining (and, if desired, displaying on an oscilloscope set to trigger on an event) the response of a power supply to a single load transient. The on/off switch can be used to switch between load states in which static-load regulation measurements are performed.

Folk, Eric N.↗

NASA NTRS2010

Symbolic Computation of Strongly Connected Components Using Saturation

Finding strongly connected components (SCCs) in the state-space of discrete-state models is a critical task in formal verification of LTL and fair CTL properties, but the potentially huge number of reachable states and SCCs constitutes a formidable challenge. This paper is concerned with computing the sets of states in SCCs or terminal SCCs of asynchronous systems. Because of its advantages in many applications, we employ saturation on two previously proposed approaches: the Xie-Beerel algorithm and transitive closure. First, saturation speeds up state-space exploration when computing each SCC in the Xie-Beerel algorithm. Then, our main contribution is a novel algorithm to compute the transitive closure using saturation. Experimental results indicate that our improved algorithms achieve a clear speedup over previous algorithms in some cases. With the help of the new transitive closure computation algorithm, up to 10(exp 150) SCCs can be explored within a few seconds.

Zhao, Yang↗

NASA NTRS2014

Crack Growth Behavior in the Threshold Region for High Cycle Fatigue Loading

The research results described in this paper presents a new understanding of the behavior of fatigue crack growth in the threshold region. It is believed by some crack growth experts that the ASTM load shedding test method does not produce true or valid threshold properties. The concern involves the observed fanning of threshold region da/dN data plots for some materials in which the low R-ratio data fans out or away from the high R-ratio data. This data fanning or elevation of threshold values is obviously caused by an increase in crack closure in the low R-ratio tested specimens. This increase in crack closure is assumed by some investigators to be caused by a plastic wake on the crack surfaces that was created during the load shedding test phase. This study shows that the increase in crack closure is the result of an extensive occurrence of crack bifurcation behavior in some materials, particularly in aluminum alloys, when the crack tip cyclic yield zone size becomes less than the grain size of the alloy. This behavior is related to the high stacking fault energy (SFE) property of aluminum alloys which results in easier slip characteristics. Therefore, the particular fanning behavior in aluminum alloys is a function of intrinsic dislocation property of the materials and that the fanned data represents valid material properties. However, for corrosion sensitive steel alloys used in this study the fanning was caused by a build-up of iron oxide at the crack tip from fretting corrosion.

Forman, Royce G.↗