Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix completion”

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 613 records · Page 34

Interactive application of quadratic expansion of chi-square statistic to nonlinear curve fitting

This report contains a detailed theoretical description of an all-purpose, interactive curve-fitting routine that is based on P. R. Bevington's description of the quadratic expansion of the Chi-Square statistic. The method is implemented in the associated interactive, graphics-based computer program. Taylor's expansion of Chi-Square is first introduced, and justifications for retaining only the first term are presented. From the expansion, a set of n simultaneous linear equations is derived, then solved by matrix algebra. A brief description of the code is presented along with a limited number of changes that are required to customize the program of a particular task. To evaluate the performance of the method and the goodness of nonlinear curve fitting, two typical engineering problems are examined and the graphical and tabular output of each is discussed. A complete listing of the entire package is included as an appendix.

Badavi, F. F.↗

An Overview of the Null-Field Method. I. Formulation and Basic Results

In this paper we revisit the fundamentals of the null-field method with discrete sources. We prove the unique solvability of the null-field equations for the total field inside the particle and the internal field outside the particle, at all wavenumbers. For this purpose, we use the equivalence between the null-field and surface integral equations methods. Furthermore, we discuss the completeness property of different systems of discrete sources for a surface field approximation, and derive an infinite set of integral equations for the surface fields in a variety of discrete sources. Finally, we formulate the null-field scheme as an approach aiming to construct an approximate solution to the scattering problem. The way in which we introduce the matrix of a particle is different from the standard approach relying on the assumption that the system of regular vector spherical wave functions for the interior problem is a basis.

Null-field method↗

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations G(r,t), a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor S(k,ω). In this work, we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in S(k,ω) going from the noninteracting limit (U = 0) to strong coupling (U→∞), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find S(k,ω), from which QFI is then calculated. We find that QFI grows with U. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above U=2.5 (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth W(U) and the half-chain von Neumann entanglement entropy. We compute G(r,t) by Fourier transforming S(k,ω). The results indicate a crossover in the short-time short-distance dynamics at low U characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large U featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by U=3. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the G(r,t) dynamics found at low U.

1-dimensional systems↗

Notes on resonances and unitarity from celestial amplitudes

We study the celestial description of the O(N) sigma model in the large N limit as introduced by Coleman, Jackiw and Politzer. Focusing on three dimensions, we analyze the implications of a UV complete, all-loop order 4-point amplitude of pions in terms of correlation functions defined on the celestial circle. We find these retain many key features from the previously studied tree-level case, such as their relation to Generalized Free Field theories and crossing-symmetry, but also incorporate new properties such as IR/UV softness and S-matrix metastable states. In particular, to understand unitarity, we propose a form of the optical theorem that controls the imaginary part of the correlator based solely on the presence of these resonances. We also explicitly analyze the conformal block expansions and factorization of four-point functions into three-point functions. We find that summing over resonances is key for these factorization properties to hold. We end with some topics for future study.

1/N expansion↗

Temperature transformation of blended magnesium potassium phosphate cement binders

In this study, a multi-technique approach was utilised to determine the high temperature performance of magnesium potassium phosphate cement (MKPC) blended with fly ash (FA) or ground granulated blast furnace slag (GBFS) with respect to nuclear waste immobilisation applications. Conceptual fire conditions were employed (up to 1200 °C, 30 min) to simulate scenarios that could occur during interim storage, transportation or within a final geological disposal facility. After exposure up to 400 °C, the main crystalline phase, struvite-K (MgKPO{sub 4}·6H{sub 2}O), was dehydrated to poorly crystalline MgKPO{sub 4} (with corresponding volumetric and mass changes), with MgKPO{sub 4} recrystallisation achieved by 800 °C. XRD and SEM/EDX analysis revealed reaction occurred between the MgKPO{sub 4} and FA/GBFS components after exposure to 1000–1200 °C, with the formation of potassium aluminosilicate phases, leucite and kalsilite (KAlSi{sub 2}O{sub 6} and KAlSiO{sub 4}), commensurate with a reduced relative intensity (or complete elimination) of the dehydrated struvite-K phase, MgKPO{sub 4.} This was further supported by solid-state NMR ({sup 27}Al and {sup 29}Si MAS), where only residual features associated with the raw FA/GBFS components were observable at 1200 °C. The high temperature phase transformation of blended MKPC binders resulted in the development of a glass/ceramic matrix with all existing porosity infilled via sintering and the formation of a vitreous phase, whilst the physical integrity was retained (no cracking or spalling). This study demonstrates that, based on small-scaled specimens, blended MKPC binders should perform satisfactorily under fire performance parameters relevant to the operation of a geological disposal facility, up to at least 1200 °C.

36 MATERIALS SCIENCE↗

Generation and optimization of cross-sections for electron-C 4 F 7 N collisions

Abstract A complete and consistent set of electron-neutral collision cross-sections for the novel insulating gas C 4 F 7 N is reported. The set is composed of a combination of cross-sections previously reported in literature, optimized via a genetic algorithm in conjunction with a multi-term Boltzmann equation solver, and calculatedab initiousing the R-matrix code Quantemol-EC. The finalized set accurately reproduces reported macroscopic rate and transport coefficients as well as Townsend coefficients and critical electric field strengths in C 4 F 7 N and its mixtures with nitrogen, carbon dioxide, and argon.

Physics↗

Part B: Pattern control of horn antennas

During this period, the computations of the impedance elements were completed. These include interactions between the two electric current modes, the elecric current mode and the magnetic current mode, and the two magnetic current modes. An accurate and efficient formulation of computing interactions between electric current mode and magnetic current mode was accomplished. This, together with other subroutines allows for the fill-in of all the elements in the matrix. After the fill-in of the impedance elements in the matrix, the forward problem is accomplished. That is, given the specification of the horn and the excitating waveguide mode, the radiation pattern of the antenna based on the integral equation can be obtained. An example case was run for a standard X-band gain-horn (DBG-520). The H- and E-plane patterns of this horn antenna with perfectly conducting walls are compared with the gain pattern available from the manufacturer for up to the first side lobe. Good agreements are obtained although the cross polarization has not yet been accounted for. The effect of the lossy coating on the radiation pattern was also investigated. The resulting E-plane pattern shows about 3-dB improvement in the first sidelobe and 4-dB improvement in the second sidelobe.

Balanis, Constantine A.↗

Model correlation and damage location for large space truss structures: Secant method development and evaluation

On-orbit testing of a large space structure will be required to complete the certification of any mathematical model for the structure dynamic response. The process of establishing a mathematical model that matches measured structure response is referred to as model correlation. Most model correlation approaches have an identification technique to determine structural characteristics from the measurements of the structure response. This problem is approached with one particular class of identification techniques - matrix adjustment methods - which use measured data to produce an optimal update of the structure property matrix, often the stiffness matrix. New methods were developed for identification to handle problems of the size and complexity expected for large space structures. Further development and refinement of these secant-method identification algorithms were undertaken. Also, evaluation of these techniques is an approach for model correlation and damage location was initiated.

Smith, Suzanne Weaver↗

Multiscale Experiments and Models for Polymer Nanocomposites

This report summarizes our research activities at the University of Virginia conducted under the support of this NASA research grant from March 2004 to October 2004 and provides a summary of the research to be conducted in the remaining period of the project at the University of Illinois at Urbana-Champaign (UIUC). The original research objectives will remain the same after our transition to UIUC. In the beginning of this program we completed the analysis on the graphite polymer micro and nanocomposites that were studied at NASA Langley and UVa in the summer of 2003. A publication summarizing this work was submitted in the fall of 2004 to the journal of Experimental Mechanics and the reviewer comments are currently addressed. The research program aims at understanding the local interactions between the polymer matrix and hard reinforcement in polymer nanocomposites (PNCs). The materials selected for this study were PNCs with different fillers (14nm silica nanospheres or nanoclays) and volume fractions (0.5%, 1%, 2% and 5%).

Chasiotis, Ioannis↗

Reexamination of the Edge Crack Torsion Test for Determining the Mode III Delamination Toughness of Laminated Composites

The edge crack torsion (ECT) test is reexamined to address two primary concerns with this method. The first of these relates to inconsistencies in the data reduction procedure that was originally proposed and is most commonly used to compute the mode III interlaminar fracture toughness, G(sub IIIc), from the experimental data. The second relates to the observation that plies neighboring the preimplanted insert in the standard ECT layup are unable to contain initiation and growth of intralaminar cracks, and that growth of these cracks contributes to errors in the apparent values of G(sub IIIc). To address these issues, an alternative data reduction technique, which uses a polynomial expression to represent the ECT specimen’s compliance versus crack length relation, is proposed. The accuracy and potential drawbacks of this expression are investigated using FE simulations. Next, seven groups of ECT specimens are fabricated and tested. Each of these groups contains specimens with preimplanted inserts that are bounded by plies at different +q and -q orientations. The seven groups consist of values of q that vary from 0° to 90° in 15° increments. Tests on these specimens show that intralaminar cracking occurs in specimens from all batches except those for which q equals 15° or 30°. However, values of G(sub IIIc) measured using specimens with these two ply orientations were not equal and, for q = 30°, the values of G(sub IIIc) that were obtained were not completely independent of initial starter crack length. Potential issues accounting for these observations and other challenges associated with measuring G(sub IIIc) using the ECT test are discussed.

Polymer matrix composites↗

Cytoskeletal filament assembly and the control of cell spreading and function by extracellular matrix

This study was undertaken to analyze how cell binding to extracellular matrix produces changes in cell shape. We focused on the initial process of cell spreading that follows cell attachment to matrix and, thus, cell 'shape' changes are defined here in terms of alterations in projected cell areas, as determined by computerized image analysis. Cell spreading kinetics and changes in microtubule and actin microfilament mass were simultaneously quantitated in hepatocytes plated on different extracellular matrix substrata. The initial rate of cell spreading was highly dependent on the matrix coating density and decreased from 740 microns 2/h to 50 microns 2/h as the coating density was lowered from 1000 to 1 ng/cm2. At approximately 4 to 6 hours after plating, this initial rapid spreading rate slowed and became independent of the matrix density regardless of whether laminin, fibronectin, type I collagen or type IV collagen was used for cell attachment. Analysis of F-actin mass revealed that cell adhesion to extracellular matrix resulted in a 20-fold increase in polymerized actin within 30 minutes after plating, before any significant change in cell shape was observed. This was followed by a phase of actin microfilament disassembly which correlated with the most rapid phase of cell extension and ended at about 6 hours; F-actin mass remained relatively constant during the slow matrix-independent spreading phase. Microtubule mass increased more slowly in spreading cells, peaking at 4 hours, the time at which the transition between rapid and slow spreading rates was observed. However, inhibition of this early rise in microtubule mass using either nocodazole or cycloheximide did not prevent this transition. Use of cytochalasin D revealed that microfilament integrity was absolutely required for hepatocyte spreading whereas interference with microtubule assembly (using nocodazole or taxol) or protein synthesis (using cycloheximide) only partially suppressed cell extension. In contrast, cell spreading could be completely inhibited by combining suboptimal doses of cytochalasin D and nocodazole, suggesting that intact microtubules can stabilize cell form when the microfilament lattice is partially compromised. The physiological relevance of the cytoskeleton and cell shape in hepatocyte physiology was highlighted by the finding that a short exposure (6 hour) of cells to nocodazole resulted in production of smaller cells 42 hours later that exhibited enhanced production of a liver-specific product (albumin). These data demonstrate that spreading and flattening of the entire cell body is not driven directly by net polymerization of either microfilaments or microtubules.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Cell Biology↗

Modifications of system for elevated temperature tensile testing and stress-strain measurement of metal matrix composites

Composites consisting of tungsten alloy wires in superalloy matrices are being studied because they offer the potential for increased strength compared to current materials used at temperatures up to at least 1093 C (2000F). Previous research at the NASA Lewis Research Center and at other laboratories in the U.S., Europe, and Japan has demonstrated laboratory feasibility for fiber reinforced superalloys (FRS). The data for the mechanical and physical properties used to evaluate candidate materials is limited and a need exists for a more detailed and complete data base. The focus of this work is to develop a test procedure to provide a more complete FRS data base to quantitatively evaluate the composite's potential for component applications. This paper will describe and discuss the equipment and procedures under development to obtain elevated temperature tensile stress-strain, strength and modulus data for the first generation of tungsten fiber reinforced superalloy composite (TFRS) materials. Tensile stress-strain tests are conducted using a constant crosshead speed tensile testing machine and a modified load-strain measuring apparatus. Elevated temperature tensile tests are performed using a resistance wound commercial furnace capable of heating test specimens up to 1093 C (2000 F). Tensile stress-strain data are obtained for hollow tubular stainless steel specimens serving as a prototype for future composite specimens.

Diaz, J. O.↗

Inversion technique to obtain local ion temperature profiles for an axisymmetric plasma with toroidal and radial velocities

A matrix inversion technique is derived to calculate local ion temperature from line-integrated measurements of an extended emission source in an axisymmetric plasma which exactly corrects for both toroidal velocity and radial velocity components. Local emissivity and toroidal velocity can be directly recovered from line-integrated spectroscopic measurements, but an independent measurement of the radial velocity is necessary to complete the temperature inversion. The extension of this technique to handle the radial velocity is relevant for magnetic reconnection and merging compression devices where temperature inversion from spectroscopic measurements is desired. A simulation demonstrates the effects of radial velocity on the determination of ion temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inversion technique to obtain local ion temperature profiles for an axisymmetric plasma with toroidal and radial velocities

A matrix inversion technique is derived to calculate local ion temperature from line-integrated measurements of an extended emission source in an axisymmetric plasma which exactly corrects for both toroidal velocity and radial velocity components. Local emissivity and toroidal velocity can be directly recovered from line-integrated spectroscopic measurements, but an independent measurement of the radial velocity is necessary to complete the temperature inversion. The extension of this technique to handle the radial velocity is relevant for magnetic reconnection and merging compression devices where temperature inversion from spectroscopic measurements is desired. A simulation demonstrates the effects of radial velocity on the determination of ion temperature.

ion temperature↗

Specific Hardening Function Definition and Characterization of a Multimechanism Generalized Potential-based Viscoelastoplasticity Model

Given the previous complete-potential structure framework together with the notion of strain- and stress-partitioning in terms of separate contributions of several submechanisms (viscoelastic and viscoplastic) to the thermodynamic functions (stored energy and dissipation) a detailed viscoelastoplastic multimechanism characterization of a specific hardening functional form of the model is presented and discussed. TIMETAL 21S is the material of choice as a comprehensive test matrix, including creep, relaxation, constant strain-rate tension tests, etc. are available at various temperatures. Discussion of these correlations tests, together with comparisons to several other experimental results, are given to assess the performance and predictive capabilities of the present model particularly with regard to the notion of hardening saturation as well as the interaction of multiplicity of dissipative (reversible/irreversible) mechanisms.

Arnold, S. M.↗

Controlled release of entrapped nanoparticles from thermoresponsive hydrogels with tunable network characteristics

Thermoresponsive hydrogels have been studied intensively for creating smart drug carriers and controlled drug delivery. Understanding the drug release kinetics and corresponding transport mechanisms of nanoparticles (NPs) in a thermoresponsive hydrogel network is the key to the successful design of smart drug delivery systems. Here, we construct a mesoscopic model of rigid NPs entrapped in a hydrogel network in an aqueous solution, where the hydrogel network is formed by cross-linked semiflexible polymers of thermoresponsive poly(N-isopropylacrylamide) (PNIPAM). By varying the temperature crossing the lower critical solution temperature of PNIPAM, we can significantly change the hydrogel network characteristics. We systematically investigate how the matrix porosity and the nanoparticle size affect the transport kinetics of NPs at different temperatures. Quantitative results on the mean-squared displacement and the van Hove displacement distributions of NPs show that all NPs entrapped in the smart hydrogels undergo subdiffusion at both low and high temperatures. For a coil state, the transport of NPs in the hydrogels can be enhanced by decreasing the matrix porosity of the polymer network and NPs' size. However, when the solution temperature is increased above the critical temperature, the hydrogel network collapses following the coil-to-globule transition, with the NPs tightly trapped in some local regions inside the hydrogels. Consequently, the NP diffusion coefficient can be reduced by two orders of magnitude, or the diffusion processes can even be completely stopped. These findings provide new insights for designing controlled drug release from stimuli-responsive hydrogels, including autonomously switch on/off drug release in response to physical and chemical stimuli.

42 ENGINEERING↗

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Subsonic panel method for designing wing surfaces from pressure distribution

An iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical distribution of pressure. The calculations are initialized by using a surface panel method to analyze a baseline wing or wing-fuselage configuration. A first-order expansion to the baseline panel method equations is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter. In every iteration cycle, the matrix is used both to calculate the geometry perturbation and to analyze the perturbed geometry. The distribution of potential on the perturbed geometry is established by simple linear extrapolation from the baseline solution. The extrapolated potential is converted to pressure by Bernoulli's equation. Not only is the accuracy of the approach good for very large perturbations, but the computing cost of each complete iteration cycle is substantially less than one analysis solution by a conventional panel method.

Bristow, D. R.↗