Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closure problem”

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.

288 records · Page 16

A Method for Calculating Strain Energy Release Rates in Preliminary Design of Composite Skin/Stringer Debonding Under Multi-Axial Loading

Three simple procedures were developed to determine strain energy release rates, G, in composite skin/stringer specimens for various combinations of unaxial and biaxial (in-plane/out-of-plane) loading conditions. These procedures may be used for parametric design studies in such a way that only a few finite element computations will be necessary for a study of many load combinations. The results were compared with mixed mode strain energy release rates calculated directly from nonlinear two-dimensional plane-strain finite element analyses using the virtual crack closure technique. The first procedure involved solving three unknown parameters needed to determine the energy release rates. Good agreement was obtained when the external loads were used in the expression derived. This superposition technique was only applicable if the structure exhibits a linear load/deflection behavior. Consequently, a second technique was derived which was applicable in the case of nonlinear load/deformation behavior. The technique involved calculating six unknown parameters from a set of six simultaneous linear equations with data from six nonlinear analyses to determine the energy release rates. This procedure was not time efficient, and hence, less appealing. A third procedure was developed to calculate mixed mode energy release rates as a function of delamination lengths. This procedure required only one nonlinear finite element analysis of the specimen with a single delamination length to obtain a reference solution for the energy release rates and the scale factors. The delamination was extended in three separate linear models of the local area in the vicinity of the delamination subjected to unit loads to obtain the distribution of G with delamination lengths. This set of sub-problems was Although additional modeling effort is required to create the sub- models, this local technique is efficient for parametric studies.

Krueger, Ronald↗

Theoretical Studies of Some HEDM Species: Cyclic O4, Cyclic O3 and Cubane

Calculations have been carried out for the HEDM species (cyclic O4, cyclic O3, and cubane) using CASSCF/derivative and CASSCF/ICCI methods. Cyclic O4 is of interest both as a potential HEDM species and because of its possible role in the ozone deficit problem in atmospheric chemistry. We have studied the pathway for decomposition from the D(2d) minimum and also have found the approximate location of the singlet triplet crossing. The barrier to decomposition is found to be about 9 kcal/mol and is not limited by the singlet triplet crossing. For cyclic O3 we have focused on the crossings between the lowest five surfaces (X(1)A(1), s(1)A(1), (1)A(2), (1)B(1), and (1)B(2)) to provide some insight into ways to form cyclic O3 photochemically. The crossing region between the X(1)A(1) and 2(1)A(1) surfaces is in agreement with the work of Xantheas et al. The calculations show that vertical excitation from the ground state to the (1)A(2) state leads to a crossing with the (1)A(1) manifold near the crossing region of the X(1)A(1) and 2(1)A(1) surfaces. We have studied the decomposition pathways for cubane to benzene plus acetylene and to cyclooctatetraene. We have also studied the ground and excited states for the photochemical ring closure step. The state which closes to cubane can be described as a double triplet pi to pi* excitation with respect to the ground state. Thus, this state has only a small oscillator strength with respect to the ground state. However, there is a singlet pi to pi* state at nearly the same energy and excitation to this state followed by intersystem crossing could lead to the triplet pi to pi* state.

Walch, Stephen P.↗

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated temperatures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple)trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A Baseline CFD model is defined and its accuracy established to be comparable to(and which usually exceeds) the accuracy of existing HTC (heat-transfer coefficient) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled),generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.1

99 GENERAL AND MISCELLANEOUS↗

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated tempera-tures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple) trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A baseline CFD model is defined and its accuracy established to be comparable to (and which usually exceeds) the accuracy of existing HTC (heat-transfer correlation) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled), generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.

36 MATERIALS SCIENCE↗

Miniature Latching Valve

A miniature latching valve has been invented to satisfy a need for an electrically controllable on/off pneumatic valve that is lightweight and compact and remains in the most recently commanded open or closed state when power is not supplied. The valve includes a poppet that is moved into or out of contact with a seat to effect closure or opening, respectively, of the flow path. Motion of the poppet is initiated by electrical heating of one of two opposing pairs of nickel/titanium shape-memory alloy (SMA) wires above their transition temperature: heated wires contract to their remembered length, applying tension to pull the poppet toward or away from the seat. A latch consisting mainly of a bistable Belleville washer (a conical spring) made of a hardened stainless steel operates between two stable positions corresponding to the fully closed or fully open state, holding the poppet in one of these positions when power is not applied to either pair of SMA wires. To obtain maximum actuation force and displacement, the SMA wires must be kept in tension. The mounting fixtures at the ends of the wires must support large tensile stresses without creating stress concentrations that would limit the fatigue lives of the wires. An earlier design provided for each wire to be crimped in a conical opening with a conical steel ferrule that was swaged into the opening to produce a large, uniformly distributed holding force. In a subsequent design, the conical ferrule was replaced with a larger crimped cylindrical ferrule depicted in the figure. A major problem in designing the valve was to protect the SMA wires from a bake-out temperature of 300 C. The problem was solved by incorporating the SMA wires into an actuator module that is inserted into a barrel of the valve body and is held in place by miniature clip rings.

Johnson, A. David↗

Effect of time delay on flying qualities: An update

Flying qualities problems of modern, full-authority electronic flight control systems are most often related to the introduction of additional time delay in aircraft response to a pilot input. These delays can have a significant effect on the flying qualities of the aircraft. Time delay effects are reexamined in light of recent flight test experience with aircraft incorporating new technology. Data from the X-29A forward-swept-wing demonstrator, a related preliminary in-flight experiment, and other flight observations are presented. These data suggest that the present MIL-F-8785C allowable-control system time delay specifications are inadequate or, at least, incomplete. Allowable time delay appears to be a function of the shape of the aircraft response following the initial delay. The cockpit feel system is discussed as a dynamic element in the flight control system. Data presented indicate that the time delay associated with a significant low-frequency feel system does not result in the predicted degradation in aircraft flying qualities. The impact of the feel system is discussed from two viewpoints: as a filter in the control system which can alter the initial response shape and, therefore, the allowable time delay, and as a unique dynamic element whose delay contribution can potentially be discounted by special pilot loop closures.

Smith, R. E.↗

Effect of time delay on flying qualities - An update

Flying qualities problems of modern, full-authority electronic flight control systems are most often related to the introduction of additional time delay in aircraft response to a pilot input. These delays can have a significant effect on the flying qualities of the aircraft. This paper reexamines time delay effects in light of recent flight test experience with aircraft incorporating new technology. Data from the X-29A forward-swept-wing demonstrator, a related preliminary in-flight experiment, and other flight observations are presented. These data suggest that the present MIL-F-8785C allowable-control system time delay specifications are inadequate or, at least, incomplete. Allowable time delay appears to be a function of the shape of the aircraft response following the initial delay. The cockpit feel system is discussed as a dynamic element in the flight control system. Data presented indicate that the time delay associated with a significant low-frequency feel system does not result in the predicted degradation in aircraft flying qualities. The impact of the feel system is discussed from two viewpoints: as a filter in the control system which can alter the initial response shape and, therefore, the allowable time delay, and as a unique dynamic element whose delay contribution can potentially be discounted by special pilot loop closures.

Smith, R. E.↗

Separation of crack extension modes in orthotropic delamination models

In the analysis of an interface crack between dissimilar elastic materials, the mode of crack extension is typically not unique, due to oscillatory behavior of near-tip stresses and displacements. This behavior currently limits the applicability of interfacial fracture mechanics as a means to predict composite delamination. The Virtual Crack Closure Technique (VCCT) is a method used to extract mode 1 and mode 2 energy release rates from numerical fracture solutions. The mode of crack extension extracted from an oscillatory solution using the VCCT is not unique due to the dependence of mode on the virtual crack extension length, Delta. In this work, a method is presented for using the VCCT to extract Delta-independent crack extension modes for the case of an interface crack between two in-plane orthotropic materials. The method does not involve altering the analysis to eliminate its oscillatory behavior. Instead, it is argued that physically reasonable, Delta-independent modes of crack extension can be extracted from oscillatory solutions. Knowledge of near-tip fields is used to determine the explicit Delta dependence of energy release rate parameters. Energy release rates are then defined that are separated from the oscillatory dependence on Delta. A modified VCCT using these energy release rate definitions is applied to results from finite element analyses, showing that Delta-independent modes of crack extension result. The modified technique has potential as a consistent method for extracting crack extension modes from numerical solutions. The Delta-independent modes extracted using this technique can also serve as guides for testing the convergence of finite element models. Direct applications of this work include the analysis of planar composite delamination problems, where plies or debonded laminates are modeled as in-plane orthotropic materials.

Beuth, Jack L.↗

Growing crops for space explorers on the moon, Mars, or in space

An option in the long-duration exploration of space, whether on the Moon or Mars or in a spacecraft on its way to Mars or the asteroids, is to utilize a bioregenerative life-support system in addition to the physicochemical systems that will always be necessary. Green plants can use the energy of light to remove carbon dioxide from the atmosphere and add oxygen to it while at the same time synthesizing food for the space travelers. The water that crop plants transpire can be condensed in pure form, contributing to the water purification system. An added bonus is that green plants provide a familiar environment for humans far from their home planet. The down side is that such a bioregenerative life-support system--called a controlled environment life-support system (CELSS) in this paper--must be highly complex and relatively massive to maintain a proper composition of the atmosphere while also providing food. Thus, launch costs will be high. Except for resupply and removal of nonrecycleable substances, such a system is nearly closed with respect to matter but open with respect to energy. Although a CELSS facility is small compared to the Earth's biosphere, it must be large enough to feed humans and provide a suitable atmosphere for them. A functioning CELSS can only be created with the help of today's advanced technology, especially computerized controls. Needed are energy for light, possibly from a nuclear power plant, and equipment to provide a suitable environment for plant growth, including a way to supply plants with the necessary mineral nutrients. All this constitutes the biomass production unit. There must also be food preparation facilities and a means to recycle or dispose of waste materials and there must be control equipment to keep the facility running. Humans are part of the system as well as plants and possibly animals. Human brain power will often be needed to keep the system functional in spite of the best computer-driven controls. The particulars of a CELSS facility depend strongly on where it is to be located. The presence of gravity on the Moon and Mars simplifies the design for a facility on those bodies, but a spacecraft in microgravity is a much more challenging environment. One problem is that plants, which are very sensitive to gravity, might not grow and produce food in the virtual absence of gravity. However, the experience with growing super-dwarf wheat in the Russian space station Mir, while not entirely successful because of the sterile wheat heads, was highly encouraging. The plants grew well for 123 days, producing more biomass than had been produced in space before. This was due to the high photon flux available to the plants and the careful control of substrate moisture. The sterile heads were probably due to the failure to remove the gaseous plant hormone, ethylene, from the Mir atmosphere. Since ethylene can easily be removed, it should be possible to grow wheat and other crops in microgravity with the production of viable seeds. On the ground Biosphere-2 taught us several lessons about the design and construction of a CELSS facility, but Bios-3 came much closer to achieving the goals of such a facility. Although stability was never completely reached, Bios-3 was much more stable than Biosphere-2 apparently because every effort was made to keep the system simple and to use the best technology available to maintain control. Wastes were not recycled in Bios-3 except for urine, and inedible plant materials were incinerated to restore CO2 to the atmosphere. Since much meat (about 20% of calories) was imported, closure in the Bios-3 experiments was well below 100%. But then, a practical CELSS on the Moon might also depend on regular resupply from Earth. Several important lessons have been learned from the CELSS research described in this review.

Review↗

The NASA Ames Closed Environmental Research Chamber: Present Status

The Closed Environmental Research Chamber (CERC) at the NASA Ames Research Center was created to investigate both components and complete systems for life support of advanced space exploration missions. This facility includes a Main Chamber, an Airlock, a Sample Transfer Lock, a Vacuum System, an Air Recompression System, a dedicated control room and a pit area for housing supporting and environmental control systems. The Main Chamber provides 310 sq ft of internal working/living space on two levels. It is planned that the CERC will be a human-rated facility for habitation simulation under mass balance closure conditions. The internal pressure will be variable over the range of 14.7 psia to 5 psia with accompanying capability for variation in atmosphere composition to maintain the oxygen partial pressure at 160 mm Hg. The CERC will be provided with a core set of primary life support subsystems for temperature and humidity control, C02 removal and trace contaminant control. Interfacing with external life support technology test b~ds with be provided, along with connection to centralized, microprocessor-based data acquisition and control systems. This paper will discuss the current status of the CERC facility and show how it is being used to address the advanced technology requirements necessary to implement an integrated working and living environment for a planetary habitat. In particular, it will be shown how the CERC, along with a human-powered centrifuge, a planetary terrain simulator and advanced displays and a virtual reality capability will work together to develop and demonstration applicable technologies for future planetary habitats. Artificial intelligence and expert system programming techniques will be used extensively to provide an automated environment for a 4-person crew. There will be several robotic mechanisms performing exploration tasks external to the habitat that will be controlled through the virtual environment to provide representative workloads for the crew. Finally, there will be a discussion of how effective are innovative new multidisciplinary test facilities to the investigation of the wide range of human and machine problems inherent in exploration missions.

Gross, Anthony R.↗

Analysis of Impacts of Terminal Restrictions on Departures in D10 TRACON

The D10 – Dallas-Fort Worth Terminal Radar Approach CONtrol (TRACON) is an air traffic facility that controls and manages aircraft and airspace that consists of multiple airports located in the North Texas area. From an air traffic management standpoint, one of the main challenges in the D10 TRACON is the departure fix capacity as multiple airports compete for resources. This problem is magnified when en route demand/capacity imbalance and inclement weather around the TRACON reduce the capacity at the terminal fixes. This leads to multiple, dynamic Traffic Management Initiative (TMI) restrictions being issued by the Air Traffic Control (ATC) on departing flights. This, in turn, propagates delay to the surface of each airport within the metroplex. The NASA Airspace Technology Demonstration 2 (ATD-2) Phase 3 is deployed in D10 to demonstrate new technologies developed to manage the Integrated Arrival, Departure, and Surface (IADS) traffic in a metroplex environment where multiple airports are interacting and sharing resources at the terminal boundary. This paper uses the ATD-2 terminal restriction data collected in the D10 TRACON to quantify the impact of restrictions on the demand, analyze the relation between terminal restrictions and departure taxi time on airport surface, and establish relationships between restrictions and surface delay. We found that the restrictions on departure flights have a direct adverse effect on departure excess taxi time on the airport surface.

ATD-2↗

Analysis of Impacts of Terminal Restrictions on Departures in D10 TRACON

The D10 – Dallas-Fort Worth Terminal Radar Approach CONtrol (TRACON) is an air traffic facility that controls and manages aircraft and airspace that consists of multiple airports located in the North Texas area. From an air traffic management standpoint, one of the main challenges in the D10 TRACON is the departure fix capacity as multiple airports compete for resources. This problem is magnified when en route demand/capacity imbalance and inclement weather around the TRACON reduce the capacity at the terminal fixes. This leads to multiple, dynamic Traffic Management Initiative (TMI) restrictions being issued by the Air Traffic Control (ATC) on departing flights. This, in turn, propagates delay to the surface of each airport within the metroplex. The NASA Airspace Technology Demonstration 2 (ATD-2) Phase 3 is deployed in D10 to demonstrate new technologies developed to manage the Integrated Arrival, Departure, and Surface (IADS) traffic in a metroplex environment where multiple airports are interacting and sharing resources at the terminal boundary. This paper uses the ATD-2 terminal restriction data collected in the D10 TRACON to quantify the impact of restrictions on the demand, analyze the relation between terminal restrictions and departure taxi time on airport surface, and establish relationships between restrictions and surface delay. We found that the restrictions on departure flights have a direct adverse effect on departure excess taxi time on the airport surface.

ATD-2↗

A Multi-Scale Sampling Strategy for Detecting Physiologically Significant Signals in AVIRIS Imagery

Models of photosynthetic production at ecosystem and global scales require multiple input parameters specifying physical and physiological surface features. While certain physical parameters (e.g., absorbed photosynthetically active radiation) can be derived from current satellite sensors, other physiologically relevant measures (e.g., vegetation type, water status, carboxylation capacity, or photosynthetic light-use efficiency), are not generally directly available from current satellite sensors at the appropriate geographic scale. Consequently, many model parameters must be assumed or derived from independent sources, often at an inappropriate scale. An abundance of ecophysiological studies at the leaf and canopy scales suggests strong physiological control of vegetation-atmosphere CO2 and water vapor fluxes, particularly in evergreen vegetation subjected to diurnal or seasonal stresses. For example hot, dry conditions can lead to stomatal closure, and associated "downregulation" of photosynthetic biochemical processes, a phenomenon often manifested as a "midday photosynthetic depression". A recent study with the revised simple biosphere (SiB2) model demonstrated that photosynthetic downregulation can significantly impact global climate. However, at the global scale, the exact significance of downregulation remains unclear, largely because appropriate physiological measures are generally unavailable at this scale. Clearly, there is a need to develop reliable ways of extracting physiologically relevant information from remote sensing. Narrow-band spectrometers offer many opportunities for deriving physiological parameters needed for ecosystem and global scale photosynthetic models. Experimental studies on the ground at the leaf- to stand-scale have indicated that several narrow-band features can be used to detect plant physiological status. One physiological signal is caused by xanthophyll cycle pigment activity, and is often expressed as the Photochemical Reflectance Index (PRI). Because the xanthophyll cycle pigments are photoregulatory pigments closely linked to photosynthetic function, this index can be used to derive relative photosynthetic rates. An additional signal with physiological significance is the 970 nm water absorption band, which provides a measure of liquid water content. This feature has been quantified both using a simple 2-band ratio (900/970 nm, here referred to as the "Water Band Index" or WBI;), and using the "continuum removal" method. Current atmospheric correction methods for AVIRIS imagery also obtain quantitative expressions of surface liquid water absorption based on the 970 nm water band and may be comparable to ground-based estimates of water content using this feature. However, physiological interpretations of both the PRI and the WBI are best understood at the leaf and canopy scales, where complications of atmospheric interference and complex stand and landscape features can be minimized, and where experimental manipulations can be readily applied. Currently it is not known whether these physiological indices can be used to derive meaningful physiological information from AVIRIS imagery. In addition to the problem of atmospheric interference, another challenge is that any simple physiological index can be confounded by multiple factors unrelated to physiology, and this problem can become more severe at progressively larger spatial scales. For example, previous work has suggested that both the PRI and the WBI, are strongly correlated with other optical measures of canopy structure (e.g., the Normalized Difference Vegetation Index or green vegetation fraction), indicating a confounding effect of structure on physiological signals at the larger, landscape scale. Furthermore, the normal operating mode of most imaging spectrometers does not allow simultaneous, ground truthing at a level of detail needed for physiological sampling. Additionally, manipulative experiments of physiology are difficult to apply at a geographic scale suitable for comparison with remote imagery, which often works at spatial scales that are several orders of magnitude larger than those typically used for physiological studies. These limitations require the consideration of alternative approaches to validating physiological information derived from AVIRIS data. In this report, we present a multi-scale sampling approach to detecting physiologically significant signals in narrow-band spectra. This approach explores the multi-dimensional data space provided by narrow-band spectrometry, and combines AVIRIS imagery at a large scale, with ground spectral sampling at an intermediate scale, and detailed ecophysiological measurements at a fine scale, to examine seasonally and spatially changing relationships between multiple structural and physiological variables. Examples of this approach are provided by simultaneous sampling of the Normalized Difference Vegetation Index (NDVI), an index of fractional PAR interception and green vegetation cover, the Water Band Index (WBI, an index of liquid water absorption, and the Photochemical Reflectance Index (PRI, an index of xanthophyll cycle pigment activity and photosynthetic light-use efficiency. By directly linking changing optical properties sampled on the ground with measurable physiological states, we hope to develop a basis for interpreting similar signals in AVIRIS imagery.

Gamon, John A.↗

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗

Leveraging Improvements in Precipitation Measuring from GPM Mission to Achieve Prediction Improvements in Climate, Weather and Hydrometeorology

The main scientific goal of the GPM mission, currently planned for start in the 2007 time frame, is to investigate important scientific problems arising within the context of global and regional water cycles. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood/draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like "core" satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing Operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe. The constellation's orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. GPM is organized internationally, currently involving a partnership between NASA in the US and the National Space Development Agency in Japan. Additionally, the program is actively pursuing agreements with other international partners and domestic scientific agencies and institutions, as well as participation by individual scientists from academia, government, and the private sector to fulfill mission goals and to pave the way for what ultimately is expected to become an internationally-organized operational global precipitation observing system. Notably, the broad societal applications of GPM are reflected in the United Nation s identification of this mission as a foremost candidate for its Peaceful Uses of Space Program. In this presentation, an overview of the GPM mission design will be presented, followed by an explanation of its scientific agenda as an outgrowth of making improvements in rain retrieval accuracy, microphysics dexterity, sampling frequency, and global coverage. All of these improvements offer new means to observe variability in precipitation and water cycle fluxes and to achieve improved predictability of weather, climate, and hydrometeorology. Specifically, the scientific agenda of GPM has been designed to leverage the measurement improvements to improve prognostic model performance, particularly quantitative precipitation forecasting and its linked phenomena at short, intermediate, and extended time scales. The talk will address how GPM measurements will enable better detection of accelerations and decelerations in regional and global water cycle processes and their relationship to climate variability, better impacts of precipitation data assimilation on numerical weather prediction and global climate reanalysis, and better performance from basin scale hydrometeorological models for short and long term flood-drought forecasting and seasonal fresh water resource assessment. Improved hydrometeorological forecasting will be possible by using continuous global precipitation observations to obtain better closure in water budgets and to generate more realistic forcing of the models themselves to achieve more accurate estimates of interception, infiltration, evaporation/transpiration fluxes, storage, and runoff.

Smith, Eric A.↗

4.32 Billion Year Old Impact Melts at Apollo 14: Dating the Procellarum Basin?

Introduction: Recent work has improved our understanding of lunar crustal structure and basin evolution on the Moon, but its early impact history and implications for solar system dynamics remain poorly established [1,2]. In particular, there is a derth of absolute ages, especially for the period >4.0 Ga, that could be related to basin-forming events. Here we present U-Pb isotopic data for Zr-rich minerals found in impact-melt fragments from Apollo 14 soil sample 14163 that, when combined with previously published data, dates an impact event at 4324±15 Ma. The compositions and ages of these impact-melt fragments provide unique information about the timing of early impact events and the composition of the lunar crust. Results: Fourteen rocklets ranging in size from 1 to 3 mm were extracted from soil 14163. Most are impactmelt rocks with 10-30% clasts of mostly pyroxene and plagioclase (100 to 500 μm) in a crystalline matrix. Textures of the matrix vary from subophitic, formed by intergrowths of 10-20 μm plagioclase and pyroxene crystals, to poikilitic, with plagioclase and pyroxene reaching 20-50 μm size. Some fragments that contain no visible clasts have similar textures are also interpreted as impact-melt rocks. All fragments contain notably large proportions of euhedral to subhedral ilmenite, zircon, apatite/merrillite and less abundant zirconolite and baddeleyite in the melt matrices. These grains often form intricate intergrowths with each other and rock-forming minerals, indicating their crystallization from the melt (Fig. 1). However, some slightly larger (~50 μm) zircon and phosphate grains can be interpreted as relict clasts based on their granular textures and relationships with the surrounding phases (Fig. 1). This implies that zircon and phosphate minerals were present in the target rocks. Some of these grains were profoundly remelted during the impact, which resulted in an oversaturation of the melt in Zr and P and crystallization of new grains of Zrrich minerals and phosphates during melt solidification. The presence of zircon and phosphates indicates that the melts were significantly enriched in KREEP components and ties their provenance to the Procellarum- KREEP Terrane (PKT) [3]. Combined U-Pb zircon data obtained for all fragments form two main clusters on a concordia diagram (Fig. 2), between about 4.3 and 3.9 Ga. Phosphate data mostly concentrate on the younger end of this age range but some analyses are almost as old as the older zircon grains (Fig. 2). Figure 1: Zircon (Zr) and phosphate (Apt; Mer) grains in impact melt fragments from Apollo 14 soil 14161. A-zircon grains crystallized from impact melt; B-phosphate grain inherited from the target; C- granular zircon grains Our interpretation of these data is that all U-rich minerals experienced variable resetting of the U-Pb system, with phosphate, where closure temperature is significantly lower than that in zircon, affected more profoundly by Pb loss than zircon. Our best estimate of the time of formation of zircon and phosphate is based on statistically valid analysis of grains interpreted texturally as grown from the impact melt (Fig. 2). Ten of these analyses define an age of 4324±15 Ma (MSWD=3.0, probability of fit P=0.002). Our best estimate for the time of resetting is obtained by combining data from phosphates that are statistically indistinguishable from 3.9 Ga within the analytical uncertainties. This group is represented by 18 analyses of 15 phosphate grains from different fragments (Fig. 2) and defines an age of 3922±6 Ma (MSWD=1.2, P=0.23). Conversely, a minimum age of the target lithologies, remelted in the impact that produced the rocklets studied here (Fig. 2), can be determined from the five oldest analyses of zircon clasts at 4338±13 Ma (MSWD=1.5, P=0.2), which is indistinguishable from the age of the impact melt within the uncertainties. The obtained ages, combined with textural evidence, imply that the impact melt was formed at 4324±15 Ma, and that it occurred in a zircon-rich target with a minimum age of 4338±13 Ma. Further reworking occurred during a second impact event at 3922±6 Ma. Discussion: The younger age of 3922±6 Ma can be interpreted as the time of the Imbrium impact. The older age of 4324±15 Ma would then be the time of formation of the impact melt, which was then caught in the Imbrium ejecta either at the Apollo 14 landing site or within the target rocks of the Imbrium impact. The presence of abundant accessory phases such as zircons and phosphates is consistent with a substantial KREEP component in the analyzed particles. The current distribution of KREEP on the lunar surface appears to be strongly influenced by Imbrum ejecta [4] so interpretation of the 4.32 Ga age depends in part on assumed structure of the pre-impact crust and distribution of KREEP within the crust. If KREEP was present only in the deep crust at 4.32 Ga, then a basin-scale impact possibly analogous to Imbrium or larger may be necessary to excavate a KREEPy impact melt at this time. However, if KREEPy materials were present closer to the surface perhaps due to redistribution related to Mg-suite magmatism, then smaller impacts might be able to rework KREEPy compositions at shallower depths. The coherence of the data on the particles analysed here suggests a large volume of melt that has been preserved since 4.32 Ga, consistent with a large impact event. Relict zircons and mineral clasts suggests that the igneous crust in the vicinity of this impact was well developed by at least 4.34 Ga, similar to the model age of KREEP and older than the isochron ages of many Mg-suite cumulates [5]. A problem that confronts all lunar sample studies using the current collection is that the pre-Imbrium geology of the PKT (the source of these 4.32 Ga impact melts) is not well constrained. The South Pole-Aitken basin contains regions that are moderately enriched in Th, but its ejecta is Th-poor [6]; therefore these A14 fragments are probably not SPA ejecta. Alternatively, these fragments may represent formation of a hypothesized Procellarum basin [7] although the lack of a clearly defined basin ejecta signature is a potential problem with that interpretation. In any case, they provide a younger limit on the age of lunar differentiation and formation of KREEP within the lunar crust. Figure 2: U-Pb data for zircon and phosphate grains from impact melt fragments. A-all data; B-data used for age calculations. References: [1] Orgel C., Michael G., Fassett C. I., van der Bogert C. H., Riedel C., Kneissl T., and Hiesinger H. (2018) J. Geophys. Res. Planets 123, 748- 762. [2] Evans A. J., Andrews-Hanna J. C., Head J. W., Soderblom J. M., Solomon S. C., and Zuber M. T. (2018) J. Geophys. Planets. 123, 1596-1617. [3] Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L. and Wieczorek, M.A. (2000) J. Geophys. Res: Planets 105, 4197-4216. [4] Haskin L. A. (1998) J. Geophys. Res. Planets 103, 1679-1689. [5] Borg L.E., Gaffney A.M., and Shearer C.K. (2015) MAPS 50, 715-732. [6] Moriarty, D.P., Watkins, R.N., Valencia, S.N., Kendall, J.D., Evans, A.J., Dygert, N. and Petro, N.E. (2021) J. Geophys. Res. Planets 126. [7] Zhu, M.H., Wünnemann, K., Potter, R.W., Kleine, T. and Morbidelli, A. (2019) J. Geophys. Res: Planets 124, 2117-2140

M D Norman↗

Anticipated Improvements in Precipitation Physics and Understanding of Water Cycle from GPM Mission

The GPM mission is currently planned for start in the late-2007 to early-2008 time frame. Its main scientific goal is to help answer pressing scientific problems arising within the context of global and regional water cycles. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood/draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like core satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe -- continuously. The constellation s orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. GPM is organized internationally, involving existing, pending, projected, and under-study partnerships which will link NASA and NOAA in the US, NASDA in Japan, ESA in Europe, ISRO in India, CNES in France, and possibly AS1 in Italy, KARI in South Korea, CSA in Canada, and AEB in Brazil. Additionally, the program is actively pursuing agreements with other international collaborators and domestic scientific agencies and institutions, as well as participation by individual scientists from academia, government, and the private sector to fulfill mission goals and to pave the way for what is expected to become an internationally-organized operational global precipitation observing system. Notably, the broad societal applications of GPM are reflected in the United Nation s identification of GPM as a foremost candidate for its Peaceful Uses of Space Program. An overview of the GPM mission design is given, followed by an explanation of its scientific agenda as an outgrowth of making improvements in rain retrieval accuracy, microphysics dexterity, sampling frequency, and global coverage. All of these improvements offer new means to observe variability in precipitation and water cycle fluxes and to achieve improved predictability of weather, climate, and hydrometeorology. Specifically, the scientific agenda of GPM has been designed to leverage the measurement improvements to improve prognostic model performance, particularly quantitative precipitation forecasting and its linked phenomena at short, intermediate, and extended time scales. The talk addresses how GPM measurements will enable better detection of accelerations and decelerations in regional and global water cycle processes and their relationship to climate variability, better impacts of precipitation data assimilation on numerical weather prediction and global climate reanalysis, and better performance from basin scale hydrometeorological models for short and long term flood-drought forecasting and seasonal fresh water resource assessment. These improvements become possible by using more accurate, more microphysically-centric, more frequent, and fully global precipitation observations to achieve better water budget closure and to provide more realistic forcing and assessment of prediction models.

Smith, Eric A.↗