Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 559 records · Page 31

Development of a Tunable Electromechanical Acoustic Liner for Engine Nacelles

This report describes the development of a tunable electromechanical Helmholtz resonator (EMHR) for engine nacelles using smart materials technology. This effort addresses both near-term and long-term goals for tunable electromechanical acoustic liner technology for the Quiet Aircraft Technology (QAT) Program. Analytical models, i.e. lumped element model (LEM) and transfer matrix (TM) representation of the EMHR, have been developed to predict the acoustic behavior of the EMHR. The models have been implemented in a MATLAB program and used to compare with measurement results. Moreover, the prediction performance of models is further improved with the aid of parameter extraction of the piezoelectric backplate. The EMHR has been experimentally investigated using standard two-microphone method (TMM). The measurement results validated both the LEM and TM models of the EMHR. Good agreement between predicted and measured impedance is obtained. Short- and open circuit loads define the limits of the tuning range using resistive and capacitive loads. There is approximately a 9% tuning limit under these conditions for the non-optimized resonator configuration studied. Inductive shunt loads result in a 3 degree-of-freedom DOF) system and an enhanced tuning range of over 20% that is not restricted by the short- and open-circuit limits. Damping coefficient ' measurements for piezoelectric backplates in a vacuum chamber are also performed and indicate that the damping is dominated by the structural damping losses, such as compliant boundaries, and other intrinsic loss mechanisms. Based on models of the EMHR, a Pareto optimization design of the EMHR has been performed for the EMHR with non-inductive loads. The EMHR with non-inductive loads is a 2DOF system with two resonant fiequencies. The tuning ranges of the two resonant frequencies of the EMHR with non-inductive loads cannot be optimized simultaneously; a trade-off (i.e., a Pareto solution) must be reached. The Pareto solution provides the information for a designer that shows how design trade-offs can be used to satisfy specific design requirements. The optimization design of the EMHR with inductive loads aims at optimal tuning of these three resonant fiequencies. The results indicate that it is possible to keep the acoustic reactance of the resonator close to a constant over a given frequency range. An effort to mimic the second layer of the NASA 2DOF liner using a piezoelectric composite diaphragm has been made. The optimal acoustic reactance of the second layer of the NASA 2DOF liner is achieved using a thin PVDF composite diaphragm, but matching the acoustic resistance requires further investigation. Acoustic energy harvesting is achieved by connecting the EMHR to an energy reclamation circuit that converts the ac voltage signal across the piezoceramic to a conditioned dc signal. Energy harvesting experiment yields 16 m W continuous power for an incident SPL of 153 dB. Such a level is sufficient to power a variety of low power electronic devices. Finally, technology transfer has been achieved by converting the original NASA ZKTL FORTRAN code to a MATLAB code while incorporating the models of the EMHR. Initial studies indicate that the EMHR is a promising technology that may enable lowpower, light weight, tunable engine nacelle liners. This technology, however, is very immature, and additional developments are required. Recommendations for future work include testing of sample EMHR liner designs in NASA Langley s normal incidence dual-waveguide and the grazing-incidence flow facility to evaluating both the impedance characteristics as well as the energy reclamation abilities. Additional design work is required for more complex tuning circuits with greater performance. Poor electromechanical coupling limited the electromechanical tuning capabilities of the proof of concept EMHR. Different materials than those studies and perhaps novel composite material systems may dramatically improvehe electromechanical coupling. Such improvements are essential to improved mimicking of existing double layer liners.

Liu, Fei↗

Scattering and Extinction by Spherical Particles Immersed in an Absorbing Host Medium

Many applications of electromagnetic scattering involve particles immersed in an absorbing rather than lossless medium, thereby making the conventional scattering theory potentially inapplicable. To analyze this issue quantitatively, we employ the FORTRAN program developed recently on the basis of the first-principles electromagnetic theory to study far-field scattering by spherical particles embedded in an absorbing infinite host medium. We further examine the phenomenon of negative extinction identified recently for monodisperse spheres and uncover additional evidence in favor of its interference origin. We identify the main effects of increasing the width of the size distribution on the ensemble-averaged extinction efficiency factor and show that negative extinction can be eradicated by averaging over a very narrow size distribution. We also analyze, for the first time, the effects of absorption inside the host medium and ensemble averaging on the phase function and other elements of the Stokes scattering matrix. It is shown in particular that increasing absorption significantly suppresses the interference structure and can result in a dramatic expansion of the areas of positive polarization. Furthermore, the phase functions computed for larger effective size parameters can develop a very deep minimum at side-scattering angles bracketed by a strong diffraction peak in the forward direction and a pronounced backscattering maximum.

Far-field electromagnetic scattering; Absorbing ho↗

Unraveling Solar Wind Space Weathering of Carbon-Rich Asterods: Low-Flux VS. High-Flux Ion Irradiation of Murchison

Introduction: Space weathering processesalter the microstructural, compositional, and optical properties of airless planetary surfaces.Solar wind space weatheringforms ion-damaged rims within the outermost ~100 nm of regolith grains through atomic-displacements and implantation of energetic H+and He+ions[1]. Transmission electron microscope(TEM)studies of space weathered Itokawa and lunar olivine particles reveal predominantly nanocrystalline solar wind-damaged rims with only minor regions of amorphization[2-5]. This finding differs from laboratoryion irradiation experiments in which samples become amorphous at fluences 2-3 orders of magnitude lower than the independentlyconstrained fluences for still-crystalline natural samples [6-7].Ion flux is generally considered to have a secondaryeffectcompared to fluence in these irradiation experiments. However, solar wind irradiation of natural samples occurs at a flux that is ~4-5 orders of magnitude lower than whatisused inlaboratory simulations (e.g., [6,8-9]). Constraining how the very low ion flux of solar wind influences the microstructural, chemical, and spectral alteration of silicate minerals will improve our understanding of solar wind space weathering and aid in the analysis of returned samples from the Hayabusa2 and OSIRIS-REx missions. Here we report results from low-and high-flux ion irradiation experiments on the Murchison (CM2)meteoritewithcomparable total fluences. Methods: We performed low-flux and high-flux He+(4 keV) and H+(1 keV) irradiation experiments on separate, dry-cut Murchison slabs. Low-flux He+irradiation reached a total fluence of 2.1 x 1016ions/cm2(~400 years of exposure at Bennu) using a flux of 3.6 x 1011ions/cm2/s; low-flux H+irradiation reached a total fluence of 3.9 x 1016ions/cm2(~30 years) using a flux of 6.6 x 1011ions/cm2/s. High-flux He+irradiation reached a total fluence of 2.0 x 1016ions/cm2(~400 years)using a flux of 9.1 x 1012ions/cm2/s; high-flux H+irradiation reached a total fluence of 5.8 x 1016ions/cm2(~50 years)using flux of 8.4 x 1012ions/cm2/s. By irradiating to similar total fluences in both the low-flux and high-flux experiments, we can better isolate and evaluatethe effects of ion flux. Differences between unirradiatedand irradiated surfaces are determined using threeanalytical techniques:(1) in situX-ray photoelectron spectroscopy(XPS) reveals changes in surface chemistry,(2) visible to near-infrared spectroscopy (VNIR; 0.35-2.5 μm) highlights changes in spectraltrends, such as albedo and slope, and (3) TEM and energy-dispersive X-ray spectroscopy (EDS) identifieschanges in nanoscale structure and elemental composition within olivine and matrix electron-transparent thin sections prepared using focused ion beam scanning electron microscopy (FIB-SEM). Results and Discussion: XPS analyses of the low-fluxand high-fluxHe+-and H+-irradiated samples showa minorreduction in surface carbon contentand chemicalreduction ofFe3+to Fe2+. The latter result is consistent with the higher Fe2+/Fe3+ratios observed in phyllosilicate-rich space weathered rims of Ryugu grains [10].VNIR analyses suggest that the low-flux H+-irradiated spectrum has a slightly bluer slope (over 0.65-2.5 μm) compared to the unirradiated VNIR spectrum from [9] while the low-flux He+-irradiated spectral slope is comparable to the unirradiated surface. Both the high-flux He+-and H+-irradiated VNIR spectra exhibitslightlyhigheralbedos at shorter wavelengths relativeto their unirradiated counterparts, with He+irradiation yielding a greater change. Thus far, the low-flux He+-irradiated matrix and olivine FIBsections as well as the low-flux H+-irradiated olivine FIBsection have been examined with TEM. The surfaces of all three samplesshowless vesiculation than high-flux/high-fluenceirradiated matrix and olivine samples from [9]. The depth of phyllosilicate amorphization in thelow-fluxHe+-irradiated matrix FIBsection varieslaterally, with crystalline domains occurring as shallow as ~25 nm in some regions. Thelow-flux He+-irradiated olivine FIBsection has acontinuous,~40-90 nm thick ion-affected layer with polycrystalline and amorphous microstructures. The low-flux H+-irradiated olivine sampleexhibits a distinct, ~10-40 nm ion-affected layer whose uppermost ~20 nmin thicker regions is predominantly amorphous.Continuing to compare results from this work to those from [9] will provide a comprehensive understanding of how ion flux impacts the microstructural, chemical, and optical modification of various mineral phases and, in turn,help identify space weathering features from parent body alteration in Hayabusa2 and OSIRIS-REx returned samples.

D L Laczniak↗

On the equinoctial orbit elements.

This paper investigates the equinoctial orbit elements for the two-body problem, showing that the associated matrices are free from singularities for zero eccentricities and zero and ninety degree inclinations. The matrix of the partial derivatives of the position and velocity vectors with respect to the orbit elements is given explicitly, together with the matrix of inverse partial derivatives, in order to facilitate construction of the matrizant (state transition matrix) corresponding to these elements. The Lagrange and Poisson bracket matrices are also given. The application of the equinoctial orbit elements to general and special perturbations is discussed.

Broucke, R. A.↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Vibration transmission through rolling element bearings in geared rotor systems

A new mathematical model is proposed to examine the vibration transmission through rolling element bearings in geared rotor systems. Current bearing models, based on either ideal boundary conditions for the shaft or purely translational stiffness element description, cannot explain how the vibratory motion may be transmitted from the rotating shaft to the casing. Experimental results have shown that the casing plate motion is primarily flexural. Here, this issue is clarified qualitatively and quantitatively by developing a comprehensive bearing stiffness matrix of dimension 6 to model precision rolling element bearings using basic principles. The proposed bearing stiffness matrix is partially verified using available analytical and experimental data, and is completely characterized. The study extends the proposed bearing formulation to analyze the overall geared rotor system dynamics including casing and mounts. The bearing stiffness matrix is included in discrete system models using lumped parameter and/or dynamic finite element techniques.

Singh, Rajendra↗

Vibration transmission through rolling element bearings. III - Geared rotor system studies

The bearing matrix formulations proposed by Lim and Singh (1990) are extended to analyze the overall dynamics of a geared rotor system which includes a spur gear pair, shafts, rolling-element bearings, a motor, a load, a casing, and flexible or rigid mounts. For this purpose, discrete vibration models are developed and used to predict vibration transmission through the bearings and to investigate the effects of the bearing, casing, and mount dynamics on the dynamic characteristics of the internal rotating system. Analytical predictions show that the theory is capable of predicting the bearing and mount moment transmissibilities in addition to the force transmissibilities. The predicted flexural vibrations of the casing plate are in good agreement with measurements conducted on an experimental set-up that consisted of a high-precision beam and pinion, and four identical rolling element bearings contained in a flexible casing mounted rigidly on a massive foundation.

Lim, T. C.↗

On the Behavior of Phosphorus During the Aqueous Alteration of CM2 Carbonaceous Chondrites

During the earliest period of solar system formation, water played an important role in the evolution of primitive dust, both after accretion of planetesimals and possible before accretion within the protoplanetary disk. Many chondrites show evidence of variable degrees of aqueous alteration, the CM2 chondrites being among the most studied [1]. This group of chondrites is characterized by mineral assemblages of both primary and secondary alteration phases. Hence, these meteorites retain a particularly important record of the reactions that occurred between primary high temperature nebular phases and water. Studies of these chondrites can provide information on the conditions and environments of aqueous alteration and the mobility of elements during alteration. This latter question is at the core of a debate concerning the location of aqueous alteration, i.e. whether alteration occurred predominantly within a closed system after accretion (parent body alteration) or whether some degree of alteration occurred within the solar nebula or on ephemeral protoplanetary bodies prior to accretion. At the core of the parent body alteration model is the hypothesis that elemental exchange between different components, principally chondrules and matrix, must have occurred. chondrules and matrix, must have occurred. In this study, we focus on the behavior of the minor element, phosphorus. This study was stimulated by observations of the behavior of P during the earliest stages of alteration in glassy mesostasis in type II chondrules in CR chondrites and extends the preliminary observations of on Y791198 to other CM chondrites.

Brearley, Adrian J.↗

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗

Modeling Progressive Failure of Bonded Joints Using a Single Joint Finite Element

Enhanced finite elements are elements with an embedded analytical solution which can capture detailed local fields, enabling more efficient, mesh-independent finite element analysis. In the present study, an enhanced finite element is applied to generate a general framework capable of modeling an array of joint types. The joint field equations are derived using the principle of minimum potential energy, and the resulting solutions for the displacement fields are used to generate shape functions and a stiffness matrix for a single joint finite element. This single finite element thus captures the detailed stress and strain fields within the bonded joint, but it can function within a broader structural finite element model. The costs associated with a fine mesh of the joint can thus be avoided while still obtaining a detailed solution for the joint. Additionally, the capability to model non-linear adhesive constitutive behavior has been included within the method, and progressive failure of the adhesive can be modeled by using a strain-based failure criteria and re-sizing the joint as the adhesive fails. Results of the model compare favorably with experimental and finite element results.

Stapleton, Scott E.↗

Telescoping Mechanics: A New Paradigm for Composite Behavior Simulation

This report reviews the application of telescoping mechanics to composites using recursive laminate theory. The elemental scale is the fiber-matrix slice, the behavior of which propagates to laminate. The results from using applications for typical, hybrid, and smart composites and composite-enhanced reinforced concrete structures illustrate the versatility and generality of telescoping scale mechanics. Comparisons with approximate, single-cell, and two- and three-dimensional finite-element methods demonstrate the accuracy and computational effectiveness of telescoping scale mechanics for predicting complex composite behavior.

Chamis, C. C.↗

Concurrent Cholesky factorization of positive definite banded Hermitian matrices

First, the Cholesky factorization is extended to cover uniformly partitioned banded positive definite matrices of rank n which may be real symmetric or Hermitian. Then, two stratagems are given for the use of the algorithm in concurrent machines where the number of processing elements is less than required to factor the matrix in as few serial steps as possible, and where uniformly high efficiency is expected from all processing elements. Expressions are given for the efficiency factor e appearing in the speed-up expression q = eN, and these are specialized for the N node hypercube machine as a function of partition size s, the number N of processing elements of the hypercube machine, and the cost mu of interelement transmission relative to computation. It is shown that the efficiency factor e is inversely proportional to mu/s, and that e is almost independent of N when N is large and mu/s = 0. The task is completed in n/s serial steps with no limit on n. The half bandwidth b of the matrix is 2 Ns.

Utku, S.↗

Metal matrix composites microfracture: Computational simulation

Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.

Mital, Subodh K.↗

Metal matrix composites microfracture - Computational simulation

Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.

Mital, S. K.↗

Analysis of notched metal matrix composites under tensile loading

This paper presents techniques based on a three-dimensional finite-element analysis for the analysis of continuous fiber reinforced metal matrix composite. Examples are shown for specific metal matrix composites such as boron/aluminum and silicon-carbide/aluminum. Specimen stress-strain behavior and stress at first fiber failure were predicted for boron/aluminum laminates containing circular holes and crack-like slits. The predictions compared very well with test data for laminates containing 0 deg fibers and reasonably well for (+/-45)2s laminates. Mesh configuration was shown to have an effect on the calculation of stresses local to the notch. The presence of thin interface layers of matrix material had a significant influence on the slit tip stress state, causing sharper stress gradients near the notch. Interface layers reduced the slit-tip fibers stress in a (+/-45)2s silicon-carbide/aluminum laminate but increased them in a (0/90)2s laminate.

Bigelow, C. A.↗

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗

Probabilisitc Geobiological Classification Using Elemental Abundance Distributions and Lossless Image Compression in Recent and Modern Organisms

Last year we presented techniques for the detection of fossils during robotic missions to Mars using both structural and chemical signatures[Storrie-Lombardi and Hoover, 2004]. Analyses included lossless compression of photographic images to estimate the relative complexity of a putative fossil compared to the rock matrix [Corsetti and Storrie-Lombardi, 2003] and elemental abundance distributions to provide mineralogical classification of the rock matrix [Storrie-Lombardi and Fisk, 2004]. We presented a classification strategy employing two exploratory classification algorithms (Principal Component Analysis and Hierarchical Cluster Analysis) and non-linear stochastic neural network to produce a Bayesian estimate of classification accuracy. We now present an extension of our previous experiments exploring putative fossil forms morphologically resembling cyanobacteria discovered in the Orgueil meteorite. Elemental abundances (C6, N7, O8, Na11, Mg12, Ai13, Si14, P15, S16, Cl17, K19, Ca20, Fe26) obtained for both extant cyanobacteria and fossil trilobites produce signatures readily distinguishing them from meteorite targets. When compared to elemental abundance signatures for extant cyanobacteria Orgueil structures exhibit decreased abundances for C6, N7, Na11, All3, P15, Cl17, K19, Ca20 and increases in Mg12, S16, Fe26. Diatoms and silicified portions of cyanobacterial sheaths exhibiting high levels of silicon and correspondingly low levels of carbon cluster more closely with terrestrial fossils than with extant cyanobacteria. Compression indices verify that variations in random and redundant textural patterns between perceived forms and the background matrix contribute significantly to morphological visual identification. The results provide a quantitative probabilistic methodology for discriminating putatitive fossils from the surrounding rock matrix and &om extant organisms using both structural and chemical information. The techniques described appear applicable to the geobiological analysis of meteoritic samples or in situ exploration of the Mars regolith. Keywords: cyanobacteria, microfossils, Mars, elemental abundances, complexity analysis, multifactor analysis, principal component analysis, hierarchical cluster analysis, artificial neural networks, paleo-biosignatures

Storrie-Lombardi, Michael C.↗