SEARCH · Engineering Papers

Results for “Full Matrix Capture”

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 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2022

Multiple Scattering Effects of Trees on L-band Microwave Using a Hybrid Method

The multiple scattering of electromagnetic waves in forests is studied using a two-step hybrid method. First, the T-matrix of a single tree is calculated based on its far-field computed with the full-wave simulations of the FEKO software. The T-matrix captures the multiple scattering caused by the tree structure. Second, the the interactions among different trees are considered using the T-matrices of the individual trees and the Foldy-Lax equations. The result of the two-step hybrid method is validated with FEKO by solving scattering from two trees directly. The multiple scattering effects are illustrated by the field solutions.

Yueh, Simon↗

NASA NTRS2021

Multi-Frequency NMM3D Simulations of Wave Propagation in Vegetation for Remote Sensing of Soil Moisture

To investigate the feasibility of using multi-frequencies to fur- ther improve the soil moisture retrieval, the recently devel- oped hybrid method is used to perform full-wave scattering simulations of a wheat field at three different frequencies of L-, S- and C- bands. In the hybrid method, the multiple scattering within a single wheat plant is first captured using the T-matrix of vec- tor cylindrical waves calculated from the full-wave solutions of HFSS. In the second step, the scattering among different plants is considered using the Foldy-Lax equations of multiple scattering theory (MST) and the extracted T-matrix. The transmission of microwaves through wheat field at the L-, S- and C-bands are calculated using the hybrid method. Results show that: (1) the transmission obtained from the full-wave simulations is much larger than those computed from the radiative transfer equations (RTE) model, (2) the hybrid method transmission has a weaker frequency dependence than those of the RTE, and (3) the attenuation caused by the vegetation layer would saturate with an increase in frequency.

Yueh, Simon↗

NASA NTRS2010

A Framework for Propagation of Uncertainties in the Kepler Data Analysis Pipeline

The Kepler space telescope is designed to detect Earth-like planets around Sun-like stars using transit photometry by simultaneously observing 100,000 stellar targets nearly continuously over a three and a half year period. The 96-megapixel focal plane consists of 42 charge-coupled devices (CCD) each containing two 1024 x 1100 pixel arrays. Cross-correlations between calibrated pixels are introduced by common calibrations performed on each CCD requiring downstream data products access to the calibrated pixel covariance matrix in order to properly estimate uncertainties. The prohibitively large covariance matrices corresponding to the ~75,000 calibrated pixels per CCD preclude calculating and storing the covariance in standard lock-step fashion. We present a novel framework used to implement standard propagation of uncertainties (POU) in the Kepler Science Operations Center (SOC) data processing pipeline. The POU framework captures the variance of the raw pixel data and the kernel of each subsequent calibration transformation allowing the full covariance matrix of any subset of calibrated pixels to be recalled on-the-fly at any step in the calibration process. Singular value decomposition (SVD) is used to compress and low-pass filter the raw uncertainty data as well as any data dependent kernels. The combination of POU framework and SVD compression provide downstream consumers of the calibrated pixel data access to the full covariance matrix of any subset of the calibrated pixels traceable to pixel level measurement uncertainties without having to store, retrieve and operate on prohibitively large covariance matrices. We describe the POU Framework and SVD compression scheme and its implementation in the Kepler SOC pipeline.

Clarke, Bruce D.↗

NASA NTRS2015

Fracture of Carbon Nanotube - Amorphous Carbon Composites: Molecular Modeling

Carbon nanotubes (CNTs) are promising candidates for use as reinforcements in next generation structural composite materials because of their extremely high specific stiffness and strength. They cannot, however, be viewed as simple replacements for carbon fibers because there are key differences between these materials in areas such as handling, processing, and matrix design. It is impossible to know for certain that CNT composites will represent a significant advance over carbon fiber composites before these various factors have been optimized, which is an extremely costly and time intensive process. This work attempts to place an upper bound on CNT composite mechanical properties by performing molecular dynamics simulations on idealized model systems with a reactive forcefield that permits modeling of both elastic deformations and fracture. Amorphous carbon (AC) was chosen for the matrix material in this work because of its structural simplicity and physical compatibility with the CNT fillers. It is also much stiffer and stronger than typical engineering polymer matrices. Three different arrangements of CNTs in the simulation cell have been investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. The effect of chemical crosslinking on composite properties is modeled by adding bonds between the CNTs and AC. The balance between weakening the CNTs and improving fiber-matrix load transfer is explored by systematically varying the extent of crosslinking. It is, of course, impossible to capture the full range of deformation and fracture processes that occur in real materials with even the largest atomistic molecular dynamics simulations. With this limitation in mind, the simulation results reported here provide a plausible upper limit on achievable CNT composite properties and yield some insight on the influence of processing conditions on the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

NASA NTRS2012

Wide Range SET Pulse Measurement

A method for measuring a wide range of SET pulses is demonstrated. Use of dynamic logic, faster than ordinary CMOS, allows capture of short pulses. A weighted binning of SET lengths allows measurement of a wide range of pulse lengths with compact circuitry. A pulse-length-conservative pulse combiner tree routes SETs from combinational logic to the measurement circuit, allowing SET measurements in circuits that cannot easily be arranged in long chains. The method is applied to add-multiplex combinational logic, and to an array of NFET routing switches, at .35 micron. Pulses are captured in a chain of Domino Logic AND gates. Propagation through the chain is frozen on the trailing edge by dropping low the second "enable" input to the AND gates. Capacitive loading is increased in the latter stages to create an approximately logarithmic weighted binning, so that a broad range of pulse lengths can be captured with a 10 stage capture chain. Simulations show pulses can be captured which are 1/5th the length of those typically captured with leading edge triggered latch methods, and less than the length of those captured with a trailing edge latch method. After capture, the pulse pattern is transferred to an SEU protected shift register for readout. 64 instances of each of two types of logic are used as targets. One is a full adder with a 4 to 1 mux on its inputs. The other is a 4 x 4 NFET routing matrix. The outputs are passed through buffered XNOR comparators to identify pulses, which are merged in a buffered not-nand (OR) tree designed to avoid pulse absorption as much as possible. The output from each of the two test circuits are input into separate pulse measurement circuits. Test inputs were provided so that the circuit could be bench tested and calibrated. A third SET measurement circuit with no inputs was used to judge the contribution from direct hits on the measurement circuit. Heavy ions were used with an LET range from 12 to 176. At LET of 21 and below, the very small number of SETs were not significantly higher in the test over the control circuits. At higher LET the test circuit SETs are one or two orders of magnitude greater than for the control circuit. The NFET circuit produces more and slightly longer SETs as expected. But the differences do not appear to be significant enough to modify strategies now used to avoid capture of SETs in chips such as FPGAs. Complete data and graphs will be in the full paper / presentation. In the summary figure below left, NOCL is the reference circuit without any input, and number of stages triggered is plotted. Simulation at right shows the smallest pulse captured (stage 2) at about 300 ps. Our conclusion is that the method is promising, but that improvements in the merge network are desirable before applying in a deep submicron process

Shuler, Robert L.↗

NASA NTRS2013

Improved Subcell Model for the Prediction of Braided Composite Response

In this work, the modeling of triaxially braided composites was explored through a semi-analytical discretization. Four unique subcells, each approximated by a "mosaic" stacking of unidirectional composite plies, were modeled through the use of layered-shell elements within the explicit finite element code LS-DYNA. Two subcell discretizations were investigated: a model explicitly capturing pure matrix regions, and a novel model which absorbed pure matrix pockets into neighboring tow plies. The in-plane stiffness properties of both models, computed using bottom-up micromechanics, correlated well to experimental data. The absorbed matrix model, however, was found to best capture out-of- plane flexural properties by comparing numerical simulations of the out-of-plane displacements from single-ply tension tests to experimental full field data. This strong correlation of out-of-plane characteristics supports the current modeling approach as a viable candidate for future work involving impact simulations.

Cater, Christopher R.↗

NASA NTRS1993

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

NASA NTRS2011

G-DYN Multibody Dynamics Engine

G-DYN is a multi-body dynamic simulation software engine that automatically assembles and integrates equations of motion for arbitrarily connected multibody dynamic systems. The algorithm behind G-DYN is based on a primal-dual formulation of the dynamics that captures the position and velocity vectors (primal variables) of each body and the interaction forces (dual variables) between bodies, which are particularly useful for control and estimation analysis and synthesis. It also takes full advantage of the spare matrix structure resulting from the system dynamics to numerically integrate the equations of motion efficiently. Furthermore, the dynamic model for each body can easily be replaced without re-deriving the overall equations of motion, and the assembly of the equations of motion is done automatically. G-DYN proved an essential software tool in the simulation of spacecraft systems used for small celestial body surface sampling, specifically in simulating touch-and-go (TAG) maneuvers of a robotic sampling system from a comet and asteroid. It is used extensively in validating mission concepts for small body sample return, such as Comet Odyssey and Galahad New Frontiers proposals.

Acikmese, Behcet↗

NASA NTRS2014

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗

NASA NTRS2003

Delaminations Investigated With Ultrasonic Spectroscopy

A previous study suggested that the ultrasonic spectroscopy technique identified possible disbonds or delaminations in polymer matrix composite (PMC) rings sectioned from flywheel rotors (ref. 1). These results went unsubstantiated by other nondestructive evaluation (NDE) methods. To explain the results, PMC rings were further investigated with ultrasonic spectroscopy (ref. 2). The ultrasonic spectroscopy system utilizes a continuous-swept sine waveform as the input. After the swept sine wave traverses the material, the captured waveform is subjected to two fast Fourier transforms. The second fast Fourier transform along with equalization of the frequency spectrum, allows for evaluation of the fundamental resonant frequency. The full-thickness resonance, the resonance corresponding to the location of the intentional disbond, and the frequency spectrum were examined in an effort to characterize the sensitivity of the NDE method to various delamination conditions.

Cosgriff, Laura M.↗

NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage models at the highest levels of assembly. The effort to generate/verify these models against their detailed counterparts is a challenge not yet conclusively solved. Higher level organizations often place a limit on the node count for delivered models, assuming smaller models result in faster computation. The burden of producing and verifying the accuracy of reduced models is primarily on the lower-level organizations, which consumes resources to produce these models. Limiting the number of allowable nodes may prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which require more nodes than centroid based models. A methodology using Thermal Desktop was described in 2010 using: (1) finite elements and edge nodes for the conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. While the methodology was clear, the implementation would have had to be done manually; with the inclusion of the OpenTD Application Programming Interface, this methodology can now be implemented programmatically and be a viable approach for the generation of reduced models. The approach was implemented at NASA-Goddard for the Capture-Containment-and-Return-System (CCRS) payload, resulting in the TCYEE tool. CCRS requires delivery with node limitations through the European Space Agency to its contractors. TCYEE generated the reduced model for delivery and the predictions compared favorably to the detailed model. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements. This paper describes the methodology, implementation, and the performance of reduced models compared to their detailed counterparts.

Computer Programming and Software↗

NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗

NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD API, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return (MSR) mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency (ESA) with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Hume L. Peabody↗

NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Thermal desktop↗

NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD Application Programming Interface, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Thermal Desktop↗

NASA NTRS2019

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗

NASA NTRS2014

The Tethered Moon

A reasonable initial condition on Earth after the Moonforming impact is that it begins as a hot global magma ocean1,2. We therefore begin our study with the mantle as a liquid ocean with a surface temperature on the order of 3000- 4000 K at a time some 100-1000 years after the impact, by which point we can hope that early transients have settled down. A 2nd initial condition is a substantial atmosphere, 100-1000 bars of H2O and CO2, supplemented by smaller amounts of CO, H2, N2, various sulfur-containing gases, and a suite of geochemical volatiles evaporated from the magma. Third, we start the Moon with its current mass at the relevant Roche limit. The 4th initial condition is the angular momentum of the Earth-Moon system. Canonical models hold this constant, whilst some recent models begin with considerably more angular momentum than is present today. Here we present a ruthlessly simplified model of Earth's cooling magmasphere based on a full-featured atmosphere and including tidal heating by the newborn Moon. Thermal blanketing by H2O-CO2 atmospheres slows cooling of a magma ocean. Geochemical volatiles - chiefly S, Na, and Cl - raise the opacity of the magma ocean's atmosphere and slow cooling still more. We assume a uniform mantle with a single internal (potential) temperature and a global viscosity. The important "freezing point" is the sharp rheological transition between a fluid carrying suspended crystals and a solid matrix through which fluids percolate. Most tidal heating takes place at this "freezing point" in a gel that is both pliable and viscous. Parameterized convection links the cooling rate to the temperature and heat generation inside the Earth. Tidal heating is a major effect. Tidal dissipation in the magma ocean is described by viscosity. The Moon is entwined with Earth by the negative feedback between thermal blanketing and tidal heating that comes from the temperature-dependent viscosity of the magma ocean. Because of this feedback, the rate that the Moon's orbit evolves is limited by the modest radiative cooling rate of Earth's atmosphere, which in effect tethers the Moon to the Earth. Consequently the Moon's orbit evolves orders of magnitude more slowly than in conventional models. Slow orbital evolution promotes capture by orbital resonances that may have been important in the Earth-Moon system

Earth-Moon system↗

NASA NTRS1995

Detection and Location of Transverse Matrix Cracks in Cross-Ply Gr/Ep Composites Using Acoustic Emission

Transverse matrix cracking in cross-ply gr/ep laminates was studied with advanced acoustic emission (AE) techniques. The primary goal of this research was to measure the load required to initiate the first transverse matrix crack in cross-ply laminates of different thicknesses. Other methods had been previously used for these measurements including penetrant enhanced radiography, optical microscopy, and audible acoustic microphone measurements. The former methods required that the mechanical test be paused for measurements at load intervals. This slowed the test procedure and did not provide the required resolution in load. With acoustic microphones, acoustic signals from cracks could not be clearly differentiated from other noise sources such as grip damage, specimen slippage, or test machine noise. A second goal for this work was to use the high resolution source location accuracy of the advanced acoustic emission techniques to determine whether the crack initiation site was at the specimen edge or in the interior of the specimen.In this research, advanced AE techniques using broad band sensors, high capture rate digital waveform acquisition, and plate wave propagation based analysis were applied to cross-ply composite coupons with different numbers of 0 and 90 degree plies. Noise signals, believed to be caused by grip damage or specimen slipping, were eliminated based on their plate wave characteristics. Such signals were always located outside the sensor gage length in the gripped region of the specimen. Cracks were confirmed post-test by microscopic analysis of a polished specimen edge, backscatter ultrasonic scans, and in limited cases, by penetrant enhanced radiography. For specimens with three or more 90 degree plies together, there was an exact 1-1 correlation between AE crack signals and observed cracks. The ultrasonic scans and some destructive sectioning analysis showed that the cracks extended across the full width of the specimen. Furthermore, the locations of the cracks from the AE data were in excellent agreement with the locations measured with the microscope. The high resolution source location capability of this technique, combined with an array of sensors, was able to determine that the cracks initiated at the specimen edges, rather than in the interior. For specimens with only one or two 90 degree plies, the crack-like signals were significantly smaller in amplitude and there was not a 1-1 correlation to observed cracks. This was similar to previous results. In this case, however, ultrasonic and destructive sectioning analysis revealed that the cracks did not extend across the specimen. They initiated at the edge, but did not propagate any appreciable distance into the specimen. This explains the much smaller AE signal amplitudes and the difficulty in correlating these signals to actual cracks in this, as well as in the previous study.

Prosser, W. H.↗