Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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

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

At least 937 records · Page 52

Machine Learning to Predict Joint Performance in Epoxy Composites Based on Process Parameters

Polymer matrix composites are gaining popularity in the aerospace industry due to their high specific strength, fatigue properties, and processability. However, based on current FAA certification guidelines, manufacturers utilizing current state-of-the art composites made with adhesive bonds commonly install redundant fasteners to guarantee the strength of these adhesively bonded composite parts.1,2 The number of fasteners in a single-aisle commercial transport aircraft is typically on the order of 105, which reduces manufacturing rate, increases cost tremendously, and reduces the advantage of the specific strength composites provide. Due to this, the Adhesive Free Bonding of Composites (AERoBOND) project at NASA Langley Research Center has developed a novel assembly process to manufacture complex composite parts without the use of adhesives and fasteners.1 However, optimization of the process is currently challenging due to the complex and interdependent process parameters. To assist with the optimization, four machine learning algorithms utilizing gradient boosting decision trees were created to provide predictions for the mechanical and characterization properties of the composite parts. Approximately 200 random states from each algorithm were tested, and the models from each state were isolated and analyzed based on their accuracy, a validation process, and their feature importance. This analysis concluded that the models created from the machine learning algorithms could accelerate a parametric study for the AERoBOND process by rapidly optimizing process parameters to achieve desired performance characteristics.

Brennen M Middleton↗

Improved Thermoplastic/Iron-Particle Transformer Cores

A method of fabricating improved transformer cores from composites of thermoplastic matrices and iron-particles has been invented. Relative to commercially available laminated-iron-alloy transformer cores, the cores fabricated by this method weigh less and are less expensive. Relative to prior polymer-matrix/ iron-particle composite-material transformer cores, the cores fabricated by this method can be made mechanically stronger and more magnetically permeable. In addition, whereas some prior cores have exhibited significant eddy-current losses, the cores fabricated by this method exhibit very small eddy-current losses. The cores made by this method can be expected to be attractive for use in diverse applications, including high-signal-to-noise transformers, stepping motors, and high-frequency ignition coils. The present method is a product of an experimental study of the relationships among fabrication conditions, final densities of iron particles, and mechanical and electromagnetic properties of fabricated cores. Among the fabrication conditions investigated were molding pressures (83, 104, and 131 MPa), and molding temperatures (250, 300, and 350 C). Each block of core material was made by uniaxial-compression molding, at the applicable pressure/temperature combination, of a mixture of 2 weight percent of LaRC (or equivalent high-temperature soluble thermoplastic adhesive) with 98 weight percent of approximately spherical iron particles having diameters in the micron range. Each molded block was cut into square cross-section rods that were used as core specimens in mechanical and electromagnetic tests. Some of the core specimens were annealed at 900 C and cooled slowly before testing. For comparison, a low-carbon-steel core was also tested. The results of the tests showed that density, hardness, and rupture strength generally increased with molding pressure and temperature, though the correlation was rather weak. The weakness of the correlation was attributed to the pores in the specimens. The maximum relative permeabilities of cores made without annealing ranged from 30 to 110, while those of cores made with annealing ranged from 900 to 1,400. However, the greater permeabilities of the annealed specimens were not associated with noticeably greater densities. The major practical result of the investigation was the discovery of an optimum distribution of iron-particle sizes: It was found that eddy-current losses in the molded cores were minimized by using 100 mesh (corresponding to particles with diameters less than or equal to 100 m) iron particles. The effect of optimization of particle sizes on eddy-current losses is depicted in the figure.

Wincheski, Russell A.↗

Modeling of Failure for Analysis of Triaxial Braided Carbon Fiber Composites

In the development of advanced aircraft-engine fan cases and containment systems, composite materials are beginning to be used due to their low weight and high strength. The design of these structures must include the capability of withstanding impact loads from a released fan blade. Relatively complex triaxially braided fiber architectures have been found to yield the best performance for the fan cases. To properly work with and design these structures, robust analytical tools are required that can be used in the design process. A new analytical approach models triaxially braided carbon fiber composite materials within the environment of a transient dynamic finite-element code, specifically the commercially available transient dynamic finite-element code LS-DYNA. The geometry of the braided composites is approximated by a series of parallel laminated composites. The composite is modeled by using shell finite elements. The material property data are computed by examining test data from static tests on braided composites, where optical strain measurement techniques are used to examine the local strain variations within the material. These local strain data from the braided composite tests are used along with a judicious application of composite micromechanics- based methods to compute the stiffness properties of an equivalent unidirectional laminated composite required for the shell elements. The local strain data from the braided composite tests are also applied to back out strength and failure properties of the equivalent unidirectional composite. The properties utilized are geared towards the application of a continuum damage mechanics-based composite constitutive model available within LS-DYNA. The developed model can be applied to conduct impact simulations of structures composed of triaxially braided composites. The advantage of this technology is that it facilitates the analysis of the deformation and damage response of a triaxially braided polymer matrix composite within the environment of a transient dynamic finite-element code such as LS-DYNA in a manner which accounts for the local physical mechanisms but is still computationally efficient. This methodology is tightly coupled to experimental tests on the braided composite, which ensures that the material properties have physical significance. Aerospace or automotive companies interested in using triaxially braided composites in their structures, particularly for impact or crash applications, would find the technology useful. By the development of improved design tools, the amount of very expensive impact testing that will need to be performed can be significantly reduced.

Goldberg, Robert K.↗

Mineralogy of Meteorite Groups

Approximately 275 mineral species have been identified in meteorites, reflecting diverse redox environments, and, in some cases, unusual nebular formation conditions. Anhydrous ordinary, carbonaceous and R chondrites contain major olivine, pyroxene and plagioclase; major opaque phases include metallic Fe-Ni, troilite and chromite. Primitive achondrites are mineralogically similar. The highly reduced enstatite chondrites and achondrites contain major enstatite, plagioclase, free silica and kamacite as well as nitrides, a silicide and Ca-, Mg-, Mn-, Na-, Cr-, K- and Ti-rich sulfides. Aqueously altered carbonaceous chondrites contain major amounts of hydrous phyllosilicates, complex organic compounds, magnetite, various sulfates and sulfides, and carbonates. In addition to kamacite and taenite, iron meteorites contain carbides, elemental C, nitrides, phosphates, phosphides, chromite and sulfides. Silicate inclusions in IAB/IIICD and lIE iron meteorites consist of mafic silicates, plagioclase and various sulfides, oxides and phosphates. Eucrites, howardites and diogenites have basaltic to orthopyroxenitic compositions and consist of major pyroxene and calcic plagioclase and several accessory oxides. Ureilttes .are made up mainly of calcic, chromian olivine and low-Ca clinopyroxene embedded in a carbonaceous matrix; accessory phases include the C polymorphs graphite, diamond, lonsdaleite and chaoite as well as metallic Fe-Ni, troilite and halides. Angrites are achondrites rich in fassaitic pyroxene (i.e. , AI-Ti diopside); minor olivine with included magnesian kirschsteinite is also present. Martian meteorites comprise basalts, Iherzolites, a dunite and an orthopyroxenite. Major phases include various pyroxenes and olivine; minor to accessory phases include various sulfides, magnetite, chromite and Ca-phosphates. Lunar meteorites comprise mare basalts with major augite and calcic plagioclase and anorthositic breccias with major calcic plagioclase. Several meteoritic phases were formed by shock metamorphism. Martensite (a2-fe,Ni) has a distorted body-centered-cubic structure and formed by a shear transformation from taenite during shock reheating and rapid cooling. The C polymorphs diamond, lonsdaleite and chaoite formed by shock from graphite. Suessite formed in the North Haig ureilite by reduction of Fe and Si (possibly from olivine) via reaction with carbonaceous matrix material. Ringwoodite, the spinel form of (Mg,Fe)2Si04, and majorite, a polymorph of (Mg,Fe)Si03 with the garnet structure, formed inside shock veins in highly shocked ordinary chondrites. Secondary minerals in meteorite finds that formed during terrestrial weathering include oxides and hy-. droxides formed directly from metallic Fe-Ni by oxidation, phosphates formed by the alteration of schreibersite, and sulfates formed by alteration of troilite.

Rubin, Alan E.↗

Mineral and organic matrix interaction in normally calcifying tendon visualized in three dimensions by high-voltage electron microscopic tomography and graphic image reconstruction

To define the ultrastructural accommodation of mineral crystals by collagen fibrils and other organic matrix components during vertebrate calcification, electron microscopic 3-D reconstructions were generated from the normally mineralizing leg tendons from the domestic turkey, Meleagris gallopavo. Embedded specimens containing initial collagen mineralizing sites were cut into 0.5-micron-thick sections and viewed and photographed at 1.0 MV in the Albany AEI-EM7 high-voltage electron microscope. Tomographic 3-D reconstructions were computed from a 2 degree tilt series of micrographs taken over a minimum angular range of +/- 60 degrees. Reconstructions of longitudinal tendon profiles confirm the presence of irregularly shaped mineral platelets, whose crystallographic c-axes are oriented generally parallel to one another and directed along the collagen long axes. The reconstructions also corroborate observations of a variable crystal length (up to 170 nm measured along crystallographic c-axes), the presence of crystals initially in either the hole or overlap zones of collagen, and crystal growth in the c-axis direction beyond these zones into adjacent overlap and other hole regions. Tomography shows for the first time that crystal width varies (30-45 nm) but crystal thickness is uniform (approximately 4-6 nm at the resolution limit of tomography); more crystals are located in the collagen hole zones than in the overlap regions at the earliest stages of tendon mineralization; the crystallographic c-axes of the platelets lie within +/- 15-20 degrees of one another rather than being perfectly parallel; adjacent platelets are spatially separated by a minimum of 4.2 +/- 1.0 nm; crystals apparently fuse in coplanar alignment to form larger platelets; development of crystals in width occurs to dimensions beyond single collagen hole zones; and a thin envelope of organic origin may be present along or just beneath the surfaces of individual mineral platelets. Implicit in the results is that the formation of crystals occurs at different sites and times by independent nucleation events in local regions of collagen. These data provide the first direct visual evidence from 3-D imaging describing the size, shape, orientation, and growth of mineral crystals in association with collagen of a normally mineralizing vertebrate tissue. They support concepts that c-axial crystal growth is unhindered by collage hole zone dimensions, that crystals are organized in the tendon in a series of generally parallel platelets, and that crystal growth in width across collagen fibrils may follow channels or grooves formed by adjacent hole zones in register.

Non-NASA Center↗

Comparison of the Computational Efficiency of the Original Versus Reformulated High-Fidelity Generalized Method of Cells

The High-Fidelity Generalized Method of Cells (HFGMC) micromechanics model has recently been reformulated by Bansal and Pindera (in the context of elastic phases with perfect bonding) to maximize its computational efficiency. This reformulated version of HFGMC has now been extended to include both inelastic phases and imperfect fiber-matrix bonding. The present paper presents an overview of the HFGMC theory in both its original and reformulated forms and a comparison of the results of the two implementations. The objective is to establish the correlation between the two HFGMC formulations and document the improved efficiency offered by the reformulation. The results compare the macro and micro scale predictions of the continuous reinforcement (doubly-periodic) and discontinuous reinforcement (triply-periodic) versions of both formulations into the inelastic regime, and, in the case of the discontinuous reinforcement version, with both perfect and weak interfacial bonding. The results demonstrate that identical predictions are obtained using either the original or reformulated implementations of HFGMC aside from small numerical differences in the inelastic regime due to the different implementation schemes used for the inelastic terms present in the two formulations. Finally, a direct comparison of execution times is presented for the original formulation and reformulation code implementations. It is shown that as the discretization employed in representing the composite repeating unit cell becomes increasingly refined (requiring a larger number of sub-volumes), the reformulated implementation becomes significantly (approximately an order of magnitude at best) more computationally efficient in both the continuous reinforcement (doubly-periodic) and discontinuous reinforcement (triply-periodic) cases.

Arnold, Steven M↗

Ballistics Model for Particles on a Horizontal Plane in a Vacuum Propelled by a Vertically Impinging Gas Jet

A simple trajectory model has been developed and is presented. The particle trajectory path is estimated by computing the vertical position as a function of the horizontal position using a constant horizontal velocity and a vertical acceleration approximated as a power law. The vertical particle position is then found by solving the differential equation of motion using a double integral of vertical acceleration divided by the square of the horizontal velocity, integrated over the horizontal position. The input parameters are: x(sub 0) and y(sub 0), the initial particle starting point; the derivative of the trajectory at x(sub 0) and y(sub 0), s(sub 0) = s(x(sub 0))= dx(y)/dy conditional expectation y = y((sub 0); and b where bx(sub 0)/y(sub 0) is the final trajectory angle before gravity pulls the particle down. The final parameter v(sub 0) is an approximation to a constant horizontal velocity. This model is time independent, providing vertical position x as a function of horizontal distance y: x(y) = (x(sub 0) + s(sub 0) (y-y(sub 0))) + bx(sub 0) -(s(sub 0)y(sub 0) ((y - y(sub 0)/y(sub 0) - ln((y/y(sub 0)))-((g(y-y(sub 0)(exp 2))/ 2((v(sub 0)(exp 2). The first term on the right in the above equation is due to simple ballistics and a spherically expanding gas so that the trajectory is a straight line intersecting (0,0), which is the point at the center of the gas impingement on the surface. The second term on the right is due to vertical acceleration, which may be positive or negative. The last term on the right is the gravity term, which for a particle with velocities less than escape velocity will eventually bring the particle back to the ground. The parameters b, s(sub 0), and in some cases v(sub 0), are taken from an interpolation of similar parameters determined from a CFD simulation matrix, coupled with complete particle trajectory simulations.

Lane, J. E.↗

Effects of Internal and External Hydrogen on Inconel 718

Internal hydrogen embrittlement (IHE) and hydrogen environment embrittlement (HEE) tensile and bend crack growth tests were performed on Inconel 718. For the IHE tests, the specimens were precharged to approximately 90 ppm hydrogen by exposure to 34.5 MPa H2 at 650 C. The HEE tests were performed in 34.5 MPa H2. Parameters evaluated were test temperature, strain rate for smooth and notch specimen geometries. The strain rate effect was very significant at ambient temperature for both IHE and HEE and decreased with increasing temperatures. For IHE, the strain rate effect was neglible at 260'C, and for HEE the strain rate effect was neglible at 400 C. At low temperatures, IHE was more severe than HEE, and at high temperatures HEE was more severe than IHE with a cross over temperature about 350 C. At 350 C, the equilibrium hydrogen concentration in Inconel 718 is about 50% lower than the hydrogen content of the precharged IHE specimens. Dislocation hydrogen sweeping of surface absorbed hydrogen was the likely transport mechanism for increasing the hydrogen concentration in the HEE tests sufficiently to produce the same degree of embrittlement as that of the more highly hydrogen charged IHE specimens. The main IHE fracture characteristic was formation of large, brittle flat facets, which decreased with increasing test temperature. The IHE fracture matrix surrounding the large facets ranged between brittle fine faceted to microvoid ductility depending upon strain rate, specimen geometry as well as temperature. The HEE fractures were characteristically fine featured, transgranular and brittle with a significant portion forming a "saw tooth" crystallographic pattern. Both IHE and HEE fractures were predominantly along the {1 1 1) slip and twin boundaries. With respect to embrittlement mechanism, it was postulated that dislocation hydrogen sweeping and hydrogen enhanced localized plasticity were active in HEE and IHE for concentrating hydrogen along (1 1 1) slip and twin planes. Final brittle failure occurred by hydrogen induced planer decohesion.

Walter, R. J.↗

Control of cytoskeletal mechanics by extracellular matrix, cell shape, and mechanical tension

We have investigated how extracellular matrix (ECM) alters the mechanical properties of the cytoskeleton (CSK). Mechanical stresses were applied to integrin receptors on the apical surfaces of adherent endothelial cells using RGD-coated ferromagnetic microbeads (5.5-microns diameter) in conjunction with a magnetic twisting device. Increasing the number of basal cell-ECM contacts by raising the fibronectin (FN) coating density from 10 to 500 ng/cm2 promoted cell spreading by fivefold and increased CSK stiffness, apparent viscosity, and permanent deformation all by more than twofold, as measured in response to maximal stress (40 dyne/cm2). When the applied stress was increased from 7 to 40 dyne/cm2, the stiffness and apparent viscosity of the CSK increased in parallel, although cell shape, ECM contacts, nor permanent deformation was altered. Application of the same stresses over a lower number ECM contacts using smaller beads (1.4-microns diameter) resulted in decreased CSK stiffness and apparent viscosity, confirming that this technique probes into the depth of the CSK and not just the cortical membrane. When magnetic measurements were carried out using cells whose membranes were disrupted and ATP stores depleted using saponin, CSK stiffness and apparent viscosity were found to rise by approximately 20%, whereas permanent deformation decreased by more than half. Addition of ATP (250 microM) under conditions that promote CSK tension generation in membrane-permeabilized cells resulted in decreases in CSK stiffness and apparent viscosity that could be detected within 2 min after ATP addition, before any measurable change in cell size.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗

Radar scene generation for tactical decision aids

The Mueller matrix and polarization covariance matrix for polarimetric radar systems was studied. The clutter is modeled by a layer of random permittivity, described by a three-dimensional correlation function, with variance, and horizontal and vertical correlation lengths. A general mixing formula was derived for discrete scatters immersed in a host medium. The results are applicable to general multiphase mixtures. The strong fluctuation theory was used to derive the backscattering cross sections, and was further extended to include higher order co-polarized and cross-polarized moments. A two-layer anisotropic random medium model was developed for the active and passive remote sensing of ice fields. A three-layer random medium model was adopted to study the volume scattering effects for the active and passive microwave remote sensing of snow-covered ice fields. The snow layer was simulated by an isotropic random medium and the ice layer by an anisotropic random medium. The vegetation canopy and snow-covered ice field were studied with a three-layer model, an isotropic random medium layer overlying an anisotropic random medium. The dyadic Green's function of the three-layer medium and the scattered electromagnetic intensities with Born approximation were calculated.

Kong, J. A.↗

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗

Extracellular matrix controls tubulin monomer levels in hepatocytes by regulating protein turnover

Cells have evolved an autoregulatory mechanism to dampen variations in the concentration of tubulin monomer that is available to polymerize into microtubules (MTs), a process that is known as tubulin autoregulation. However, thermodynamic analysis of MT polymerization predicts that the concentration of free tubulin monomer must vary if MTs are to remain stable under different mechanical loads that result from changes in cell adhesion to the extracellular matrix (ECM). To determine how these seemingly contradictory regulatory mechanisms coexist in cells, we measured changes in the masses of tubulin monomer and polymer that resulted from altering cell-ECM contacts. Primary rat hepatocytes were cultured in chemically defined medium on bacteriological petri dishes that were precoated with different densities of laminin (LM). Increasing the LM density from low to high (1-1000 ng/cm2), promoted cell spreading (average projected cell area increased from 1200 to 6000 microns2) and resulted in formation of a greatly extended MT network. Nevertheless, the steady-state mass of tubulin polymer was similar at 48 h, regardless of cell shape or ECM density. In contrast, round hepatocytes on low LM contained a threefold higher mass of tubulin monomer when compared with spread cells on high LM. Furthermore, similar results were obtained whether LM, fibronectin, or type I collagen were used for cell attachment. Tubulin autoregulation appeared to function normally in these cells because tubulin mRNA levels and protein synthetic rates were greatly depressed in round cells that contained the highest level of free tubulin monomer. However, the rate of tubulin protein degradation slowed, causing the tubulin half-life to increase from approximately 24 to 55 h as the LM density was lowered from high to low and cell rounding was promoted. These results indicate that the set-point for the tubulin monomer mass in hepatocytes can be regulated by altering the density of ECM contacts and changing cell shape. This finding is consistent with a mechanism of MT regulation in which the ECM stabilizes MTs by both accepting transfer of mechanical loads and altering tubulin degradation in cells that continue to autoregulate tubulin synthesis.

Non-NASA Center↗

Measurements of Erosion Wear Volume Loss on Bare and Coated Polymer Matrix Composites

An investigation was conducted to examine the erosion behavior of uncoated and coated polymer matrix composite (PMC) specimens subjected to solid particle impingement using air jets. The PMCs were carbon-Kevlar (DuPont, Wilmington, DE) fiber-epoxy resin composites with a temperature capability up to 393 K (248 F). Tungsten carbide-cobalt (WC-Co) was the primary topcoat constituent. Bondcoats were applied to the PMC substrates to improve coating adhesion; then, erosion testing was performed at the University of Cincinnati. All erosion tests were conducted with Arizona road-dust (ARD), impinging at angles of 20 and 90 on both uncoated and two-layer coated PMCs at a velocity of 229 m/s and at a temperature of 366 K (200 F). ARD contains primarily 10-m aluminum oxide powders. Vertically scanning interference microscopy (noncontact, optical profilometry) was used to evaluate surface characteristics, such as erosion wear volume loss and depth, surface topography, and surface roughness. The results indicate that noncontact, optical interferometry can be used to make an accurate determination of the erosion wear volume loss of PMCs with multilayered structures while preserving the specimens. The two-layered (WC-Co topcoat and metal bondcoat) coatings on PMCs remarkably reduced the erosion volume loss by a factor of approximately 10. The tenfold increase in erosion resistance will contribute to longer PMC component lives, lower air friction, reduced related breakdowns, decreased maintenance costs, and increased PMC reliability. The decrease in the surface roughness of the coated vanes will lead to lower air friction and will subsequently reduce energy consumption. Eventually, the coatings could lead to overall economic savings.

Miyoshi, Kazuhisa↗

Orientation of human optokinetic nystagmus to gravity: a model-based approach

Optokinetic nystagmus (OKN) was induced by having subjects watch a moving display in a binocular, head-fixed apparatus. The display was composed of 3.3 degrees stripes moving at 35 degrees/s for 45 s. It subtended 88 degrees horizontally by 72 degrees vertically of the central visual field and could be oriented to rotate about axes that were upright or tilted 45 degrees or 90 degrees. The head was held upright or was tilted 45 degrees left or right on the body during stimulation. Head-horizontal (yaw axis) and head-vertical (pitch axis) components of OKN were recorded with electro-oculography (EOG). Slow phase velocity vectors were determined and compared with the axis of stimulation and the spatial vertical (gravity axis). With the head upright, the axis of eye rotation during yaw axis OKN was coincident with the stimulus axis and the spatial vertical. With the head tilted, a significant vertical component of eye velocity appeared during yaw axis stimulation. As a result the axis of eye rotation shifted from the stimulus axis toward the spatial vertical. Vertical components developed within 1-2 s of stimulus onset and persisted until the end of stimulation. In the six subjects there was a mean shift of the axis of eye rotation during yaw axis stimulation of approximately 18 degrees with the head tilted 45 degrees on the body. Oblique optokinetic stimulation with the head upright was associated with a mean shift of the axis of eye rotation toward the spatial vertical of 9.2 degrees. When the head was tilted and the same oblique stimulation was given, the axis of eye rotation rotated to the other side of the spatial vertical by 5.4 degrees. This counterrotation of the axis of eye rotation is similar to the "Muller (E) effect," in which the perception of the upright is counterrotated to the opposite side of the spatial vertical when subjects are tilted in darkness. The data were simulated by a model of OKN with a "direct" and "indirect" pathway. It was assumed that the direct visual pathway is oriented in a body, not a spatial frame of reference. Despite the short optokinetic after-nystagmus time constants, strong horizontal to vertical cross-coupling could be produced if the horizontal and vertical time constants were in proper ratio and there were no suppression of nystagmus in directions orthogonal to the stimulus direction. The model demonstrates that the spatial orientation of OKN can be achieved by restructuring the system matrix of velocity storage. We conclude that an important function of velocity storage is to orient slow-phase velocity toward the spatial vertical during movement in a terrestrial environment.

Non-NASA Center↗

The effects of gravity darkening on the ultraviolet continuum polarization produced by circumstellar disks

We investigate the effects of gravity darkening on the UV continuum polarization produced by an axisymmetric disk that surrounds a rapidly rotating star. Although the model is a single scattering approximation, we do include the effects of attenuation (electron scattering plus hydrogen bound-free absorption) by the disk, using an approach similar to that of Sobolev (1963). Because of the gravity darkening of the star and the attenuation within the disk, the radiation field is not axially symmetric about the radius vector. This implies that the polarization source functions are no longer provided by the finite disk depolarization factors of Cassinelli, Nordsieck, & Murison (1987), which are functions of the intensity moments in a spherically symmetric atmosphere. We reformulate the polarization source functions using generalized intensity moment tensors (J, H(sub i), K(sub ij)) that are valid for an arbitrary radiation field and envelope geometry. We find that the polarization source functions are simplest when using intensity moments in the observer's reference frame. On the other hand, the intensity moments are most easily evaluated in the stellar reference frame. Using the rotation transformation properties of the generalized intensity moments, we relate the observer's moments to those evaluated in the stellar reference frame. Our procedure for determining the polarization source functions thus merely involves choosing a set of Euler angles for the coordinate rotations, and then evaluating the associated rotation matrix. The geometrical complications of polarization transfer are thus reduced to obtaining a coordinate rotation matrix.

Bjorkman, J. E.↗

Materials with periodic internal structure: Computation based on homogenization and comparison with experiment

The combination of thermal and mechanical loading expected in practice means that constitutive equations of metal matrix composites must be developed which deal with time-independent and time-dependent irreversible deformation. Also, the internal state of composites is extremely complicated which underlines the need to formulate macroscopic constitutive equations with a limited number of state variables which represent the internal state at the micro level. One available method for calculating the macro properties of composites in terms of the distribution and properties of the constituent materials is the method of homogenization whose formulation is based on the periodicity of the substructure of the composite. A homogenization procedure was developed which lends itself to the use of the finite element procedure. The efficiency of these procedures, to determine the macroscopic properties of a composite system from its constituent properties, was demonstrated utilizing an aluminum plate perforated by directionally oriented slits. The selection of this problem is based on the fact that, extensive experimental results exist, the macroscopic response is highly anisotropic, and that the slits provide very high stress gradients which severely test the effectiveness of the computational procedures. Furthermore, both elastic and plastic properties were investigated so that the application to practical systems with inelastic deformation should be able to proceed without difficulty. The effectiveness of the procedures was rigorously checked against experimental results and with the predictions of approximate calculations. Using the computational results it is illustrated how macroscopic constitutive equations can be expressed in forms of the elastic and limit load behavior.

Jansson, S.↗

Application and development of the Schwinger multichannel scattering theory and the partial differential equation theory of electron-molecule scattering

One version of the multichannel theory for electron-target scattering based on the Schwinger variational principle, the SMC method, requires the introduction of a projection parameter. The role of the projection parameter a is investigated and it is shown that the principal-value operator in the SMC equation is Hermitian regardless of the value of a as long as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix representation of this operator is also Hermitian. The use of such basis is consistent with the Schwinger variational principle because the Lippmann-Schwinger equation automatically builds in the correct boundary conditions. Otherwise, an auxiliary condition needs to be introduced, and Takatsuka and McKoy's original value of a is one of the three possible ways to achieve Hermiticity. In all cases but one, a can be uncoupled from the Hermiticity condition and becomes a free parameter. An equation for a based on the variational stability of the scattering amplitude is derived; its solution has an interesting property that the scattering amplitude from a converged SMC calculation is independent of the choice of a even though the SMC operator itself is a-dependent. This property provides a sensitive test of the convergence of the calculation. For a static-exchange calculation, the convergence requirement only depends on the completeness of the one-electron basis, but for a general multichannel case, the a-invariance in the scattering amplitude requires both the one-electron basis and the N plus 1-electron basis to be complete. The role of a in the SMC equation and the convergence property are illustrated using two examples: e-CO elastic scattering in the static-exchange approximation, and a two-state treatment of the e-H2 Chi(sup 1)Sigma(sub g)(+) yields b(sup 3)Sigma(sub u)(+) excitation.

Weatherford, Charles A.↗

Characterization of Meteorites by Focused Ion Beam Sectioning: Recent Applications to CAIs and Primitive Meteorite Matrices

Focused ion beam (FIB) sectioning has revolutionized preparation of meteorite samples for characterization by analytical transmission electron microscopy (TEM) and other techniques. Although FIB is not "non-destructive" in the purest sense, each extracted section amounts to no more than nanograms (approximately 500 cubic microns) removed intact from locations precisely controlled by SEM imaging and analysis. Physical alteration of surrounding material by ion damage, fracture or sputter contamination effects is localized to within a few micrometers around the lift-out point. This leaves adjacent material intact for coordinate geochemical analysis by SIMS, microdrill extraction/TIMS and other techniques. After lift out, FIB sections can be quantitatively analyzed by electron microprobe prior to final thinning, synchrotron x-ray techniques, and by the full range of state-of-the-art analytical field-emission scanning transmission electron microscope (FE-STEM) techniques once thinning is complete. Multiple meteorite studies supported by FIB/FE-STEM are currently underway at NASA-JSC, including coordinated analysis of refractory phase assemblages in CAIs and fine-grained matrices in carbonaceous chondrites. FIB sectioning of CAIs has uncovered epitaxial and other overgrowth relations between corundum-hibonite-spinel consistent with hibonite preceding corundum and/or spinel in non-equilibrium condensation sequences at combinations of higher gas pressures, dust-gas enrichments or significant nebular transport. For all of these cases, the ability of FIB to allow for coordination with spatially-associated isotopic data by SIMS provides immense value for constraining the formation scenarios of the particular CAI assemblage. For carbonaceous chondrites matrix material, FIB has allowed us to obtain intact continuous sections of the immediate outer surface of Murchison (CM2) after it has been experimentally ion processed to simulate solar wind space weathering. The surface amorphization and loss of OH produced by the irradiation provides important clues regarding space weathering on primitive asteroids such as the OSIRIS-Rex target 101955 Bennu.

Christoffersen, Roy↗