Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 1,423 records · Page 79

Designing Ti-6Al-4V microstructure for strain delocalization using neural networks

Abstract The deformation behavior of Ti-6Al-4V titanium alloy is significantly influenced by slip localized within crystallographic slip bands. Experimental observations reveal that intense slip bands in Ti-6Al-4V form at strains well below the macroscopic yield strain and may serially propagate across grain boundaries, resulting in long-range localization that percolates through the microstructure. These connected, localized slip bands serve as potential sites for crack initiation. Although slip localization in Ti-6Al-4V is known to be influenced by various factors, an investigation of optimal microstructures that limit localization remains lacking. In this work, we develop a novel strategy that integrates an explicit slip band crystal plasticity technique, graph networks, and neural network models to identify Ti-6Al-4V microstructures that reduce the propensity for strain localization. Simulations are conducted on a dataset of 3D polycrystals, each represented as a graph to account for grain neighborhood and connectivity. The results are then used to train neural network surrogate models that accurately predict localization-based properties of a polycrystal, given its microstructure. These properties include the ratio of slip accumulated in the band to that in the matrix, fraction of total applied strain accommodated by slip bands, and spatial connectivity of slip bands throughout the microstructure. The initial dataset is enriched by synthetic data generated by the surrogate models, and a grid search optimization is subsequently performed to find optimal microstructures. Describing a 3D polycrystal with only a few features and a combination of graph and neural network models offer robustness compared to the alternative approaches without compromising accuracy. We show that while each material property is optimized through a unique microstructure solution, elongated grain shape emerges as a recurring feature among all optimal microstructures. This finding suggests that designing microstructures with elongated grains could potentially mitigate strain localization without compromising strength.

Ahmadikia, Behnam↗

O-GlcNAc transferase regulates collagen deposition and fibrosis resolution in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic pulmonary disease that is characterized by an excessive accumulation of extracellular matrix (ECM) proteins (e.g. collagens) in the parenchyma, which ultimately leads to respiratory failure and death. While current therapies exist to slow the progression, no therapies are available to resolve fibrosis. We characterized the O-linked N-Acetylglucosamine (O-GlcNAc) transferase (OGT)/O-GlcNAc axis in IPF using single-cell RNA-sequencing (scRNA-seq) data and human lung sections and isolated fibroblasts from IPF and non-IPF donors. The underlying mechanism(s) of IPF were further investigated using multiple experimental models to modulate collagen expression and accumulation by genetically and pharmacologically targeting OGT. Furthermore, we hone in on the transforming growth factor-beta (TGF-β) effector molecule, Smad3, by co-expressing it with OGT to determine if it is modified and its subsequent effect on Smad3 activation. We found that OGT and O-GlcNAc levels are upregulated in patients with IPF compared to non-IPF. We report that the OGT regulates collagen deposition and fibrosis resolution, which is an evolutionarily conserved process demonstrated across multiple species. Co-expression of OGT and Smad3 showed that Smad3 is O-GlcNAc modified. Blocking OGT activity resulted in decreased phosphorylation at Ser-423/425 of Smad3 attenuating the effects of TGF-β1 induced collagen expression/deposition. OGT inhibition or knockdown successfully blocked and reversed collagen expression and accumulation, respectively. Smad3 is discovered to be a substrate of OGT and its O-GlcNAc modification(s) directly affects its phosphorylation state. These data identify OGT as a potential target in pulmonary fibrosis resolution, as well as other diseases that might have aberrant ECM/collagen accumulation.

59 BASIC BIOLOGICAL SCIENCES↗

On polynomial preconditioning for indefinite Hermitian matrices

The minimal residual method is studied combined with polynomial preconditioning for solving large linear systems (Ax = b) with indefinite Hermitian coefficient matrices (A). The standard approach for choosing the polynomial preconditioners leads to preconditioned systems which are positive definite. Here, a different strategy is studied which leaves the preconditioned coefficient matrix indefinite. More precisely, the polynomial preconditioner is designed to cluster the positive, resp. negative eigenvalues of A around 1, resp. around some negative constant. In particular, it is shown that such indefinite polynomial preconditioners can be obtained as the optimal solutions of a certain two parameter family of Chebyshev approximation problems. Some basic results are established for these approximation problems and a Remez type algorithm is sketched for their numerical solution. The problem of selecting the parameters such that the resulting indefinite polynomial preconditioners speeds up the convergence of minimal residual method optimally is also addressed. An approach is proposed based on the concept of asymptotic convergence factors. Finally, some numerical examples of indefinite polynomial preconditioners are given.

Freund, Roland W.↗

Errors induced by the neglect of polarization in radiance calculations for Rayleigh-scattering atmospheres

Although neglecting polarization and replacing the rigorous vector radiative transfer equation by its approximate scalar counterpart has no physical background, it is a widely used simplification when the incident light is unpolarized and only the intensity of the reflected light is to be computed. We employ accurate vector and scalar multiple-scattering calculations to perform a systematic study of the errors induced by the neglect of polarization in radiance calculations for a homogeneous, plane-parallel Rayleigh-scattering atmosphere (with and without depolarization) above a Lambertian surface. Specifically, we calculate percent errors in the reflected intensity for various directions of light incidence and reflection, optical thicknesses of the atmosphere, single-scattering albedos, depolarization factors, and surface albedos. The numerical data displayed can be used to decide whether or not the scalar approximation may be employed depending on the parameters of the problem. We show that the errors decrease with increasing depolarization factor and/or increasing surface albedo. For conservative or nearly conservative scattering and small surface albedos, the errors are maximum at optical thicknesses of about 1. The calculated errors may be too large for some practical applications, and, therefore, rigorous vector calculations should be employed whenever possible. However, if approximate scalar calculations are used, we recommend to avoid geometries involving phase angles equal or close to 0 deg and 90 deg, where the errors are especially significant. We propose a theoretical explanation of the large vector/scalar differences in the case of Rayleigh scattering. According to this explanation, the differences are caused by the particular structure of the Rayleigh scattering matrix and come from lower-order (except first-order) light scattering paths involving right scattering angles and right-angle rotations of the scattering plane.

Mishchenko, M. I.↗

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗

Programmable display pushbuttons on the Space Station's telerobot control panel

The Man-Systems Telerobotics Laboratory at NASA's Johnson Space Center and supported by Lockheed, is working to ensure that the Flight Telerobotic Servicer (FTS) to be used on the Space Shuttle (Orbiter) and the Space Station has a well designed user interface from a Human Factors perspective. The FTS, which is a project led by NASA's Goddard Space Flight Center, will be a telerobot used for Space Station construction, maintenance, and satellite repair. It will be directly controlled from workstations on the Orbiter and the Space Station and monitored from a ground workstation. The FTS will eventually evolve into a more autonomous system, but in the short-term the system will be manually operated (teleoperated) for many tasks. This emphasizes the importance of the human/telerobot interface on this system. This phase of the FTS workstation evaluation covers a preliminary study of programmable display pushbuttons (PDP's). The PDP is constructed of a matrix of directly addressable electroluminescent (EL) pixels which can be used to form dot-matrix characters. PDP's can be used to display more than one message and to control more than one function. Since the PDP's have these features, then a single PDP may possibly replace the use of many single-function pushbuttons, rotary switches, and toggle switches, thus using less panel space. It is of interest to determine if PDP's can be used to adequately perform complex hierarchically structured task sequences. The objective of this investigation was to compare the performance of experienced and inexperienced Remote Manipulator System (RMS) operators while performing an RMS like task on simulated PDP and non-PDP computer prototypes so that guidelines governing the use of programmable display pushbuttons on the FTS workstation could be created. The functionality of the RMS on the Orbiter was used as a model for this evaluation since the functionality of the FTS at the time of this writing has not been solidified.

Stuart, Mark A.↗

Modulation of lens cell adhesion molecules by particle beams

Cell adhesion molecules (CAMs) are proteins which anchor cells to each other and to the extracellular matrix (ECM), but whose functions also include signal transduction, differentiation, and apoptosis. We are testing a hypothesis that particle radiations modulate CAM expression and this contributes to radiation-induced lens opacification. We observed dose-dependent changes in the expression of beta 1-integrin and ICAM-1 in exponentially-growing and confluent cells of a differentiating human lens epithelial cell model after exposure to particle beams. Human lens epithelial (HLE) cells, less than 10 passages after their initial culture from fetal tissue, were grown on bovine corneal endothelial cell-derived ECM in medium containing 15% fetal bovine serum and supplemented with 5 ng/ml basic fibroblast growth factor (FGF-2). Multiple cell populations at three different stages of differentiation were prepared for experiment: cells in exponential growth, and cells at 5 and 10 days post-confluence. The differentiation status of cells was characterized morphologically by digital image analysis, and biochemically by Western blotting using lens epithelial and fiber cell-specific markers. Cultures were irradiated with single doses (4, 8 or 12 Gy) of 55 MeV protons and, along with unirradiated control samples, were fixed using -20 degrees C methanol at 6 hours after exposure. Replicate experiments and similar experiments with helium ions are in progress. The intracellular localization of beta 1-integrin and ICAM-1 was detected by immunofluorescence using monoclonal antibodies specific for each CAM. Cells known to express each CAM were also processed as positive controls. Both exponentially-growing and confluent, differentiating cells demonstrated a dramatic proton-dose-dependent modulation (upregulation for exponential cells, downregulation for confluent cells) and a change in the intracellular distribution of the beta 1-integrin, compared to unirradiated controls. In contrast, there was a dose-dependent increase in ICAM-1 immunofluorescence in confluent, but not exponentially-growing cells. These results suggest that proton irradiation downregulates beta 1-integrin and upregulates ICAM-1, potentially contributing to cell death or to aberrant differentiation via modulation of anchorage and/or signal transduction functions. Quantification of the expression levels of the CAMs by Western analysis is in progress.

NASA Discipline Radiation Health↗

TURBOMAT: A Probabilistic Turbomachinery Aeroelastic Analysis Tool

An integration of aeroelastic analysis procedures with probabilistic analysis methods enables us to design safe reliable engines with quantified reliability. Towards this goal, a graphical user interface (GUI) based tool that integrates the codes Aeroelastic analysis of propfans (ASTROP2) and Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) is developed. The tool entitled TURBOMachinery Aeroelastic Analysis Tool (TURBOMAT), is developed utilizing the MATrix Laboratory (Matlab) Guide (Graphical User Interface Development) tool box. TURBOMAT provides a user friendly computational environment for rapid assessment of Turbomachinery blades flutter characteristics, subjected to uncertain loading conditions with variability in material and aerodynamic properties. The tool is seen as an education tool for new students and young engineers starting their careers in structural Aeroelasticity who want to learn and understand aeroelastic aspects of turbomachinery components, fans, compressors and turbines, including uncertainties in loading and material properties.A typical fan blade configuration geometry was chosen to demonstrate the tool. The results are presented in the form of probabilistic density function (PDF), the cumulative distribution function (CDF) and sensitivity factors. Both first order fast probability integration (FPI) and the Monto Carlo (MC) techniques are used in the analysis and compared. The tool enabled us to quantify blade flutter reliability as well as the ranking of uncertain variables and their importance to blade flutter response.

Aeroelastic Analysis↗

Hybrid analysis of radiative corrections to neutron decay with current algebra and effective field theory

We introduce a useful framework for high-precision studies of the neutron beta decay by merging the current algebra description and the fixed-order effective field theory calculation of the electroweak radiative corrections to the neutron axial form factor. We discuss the advantages of this hybrid method and show that it only requires a minimal amount of lattice QCD inputs to achieve a 10 –4 theory accuracy for the Standard Model prediction of the neutron lifetime and the axial-to-vector coupling ratio λ, both important to the search for physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Systematic Generation of Moment Invariant Bases for 2D and 3D Tensor Fields

Moment invariants have been successfully applied to pattern detection tasks in 2D and 3D scalar, vector, and matrix valued data. However so far no flexible basis of invariants exists, i.e., no set that is optimal in the sense that it is complete and independent for every input pattern. Here, we prove that a basis of moment invariants can be generated that consists of tensor contractions of not more than two different moment tensors each under the conjecture of the set of all possible tensor contractions to be complete. This result allows us to derive the first generator algorithm that produces flexible bases of moment invariants with respect to orthogonal transformations by selecting a single non-zero moment to pair with all others in these two-factor products. Since at least one non-zero moment can be found in every non-zero pattern, this approach always generates a complete set of descriptors.

97 MATHEMATICS AND COMPUTING↗

Creep Properties of NiAl-1Hf Single Crystals Re-Investigated

NiAl-1Hf single crystals have been shown to be quite strong at 1027 C, with strength levels approaching those of advanced Ni-based superalloys. Initial testing, however, indicated that the properties might not be reproducible. Study of the 1027 C creep behavior of four different NiAl-1Hf single-crystal ingots subjected to several different heat treatments indicated that strength lies in a narrow band. Thus, we concluded that the mechanical properties are reproducible. Recent investigations of the intermetallic NiAl have confirmed that minor alloying additions combined with single-crystal growth technology can produce elevated temperature strength levels approaching those of Ni-based superalloys. For example, General Electric alloy AFN 12 {Ni-48.5(at.%) Al-0.5Hf-1Ti-0.05Ga} has a creep rupture strength equivalent to Rene 80 combined with a approximately 30-percent lower density, a fourfold improvement in thermal conductivity, and the ability to form a self-protective alumina scale in aggressive environments. Although the compositions of strong NiAl single crystals are relatively simple, the microstructures are complex and vary with the heat treatment and with small ingot-toingot variations in the alloy chemistry. In addition, initial testing suggested a strong dependence between microstructure and creep strength. If these observations were true, the ability to utilize NiAl single-crystal rotating components in turbine machinery could be severely limited. To investigate the possible limitations in the creep response of high-strength NiAl single crystals, the NASA Glenn Research Center at Lewis Field initiated an in depth investigation of the effect of heat treatment on the microstructure and subsequent 1027 C creep behavior of [001]-oriented NiAl-1Hf with a nominal chemistry of Ni-47.5Al-1Hf-0.5Si. This alloy was selected since four ingots, grown over a number of years and possessing slightly different compositions, were available for study. Specimens taken from the ingots were subjected to several heat treatment schedules, examined by transmission electron microscopy, and tested in both compression and tension. An example of the microstructure found in a [001]-oriented NiAl-1Hf specimen after a solution treatment at 1317 C for 50 hr followed by air cooling is illustrated in the image on the left, where the NiAl matrix contains a uniform distribution of nanometer-scale Gphase (Ni16Hf6Si7) precipitates. Other heat treating schedules produced microstructures with nanometer-sized G-phase cubes and plates or, in an extreme case, produced a microstructure with all the G-phase converted to Heusler (Ni2AlHf) particles. The results of 1027 C creep strength and strain rate testing are illustrated which summarizes data from tensile and compressive testing of samples cut from all four NiAl-1Hf ingots and subjected to a variety of heat treatment schedules. With one exception, all the strength values lie in a narrow band that spans six orders of magnitude in strain rate. The only factor that produced results outside of this band was the heat treatment schedule that dissolved all the G-phase and replaced it with Heusler precipitates. The results portrayed in this figure lead to the important practical conclusion that the elevated-temperature creep properties of NiAl-1Hf single crystals are reproducible and are not affected by small variations in alloy chemistry from ingot to ingot or by different initial distributions of G-phase in the heat-treated alloy. The only variable in this study that produced a significant and delerious effect on mechanical strength was a post-solution heat treatment that lead to the complete disappearance of the G-phase in favor of Heusler precipitates.

Whittenberger, J. Daniel↗

The Myth of Fungible FTE: A Quantitative Assessment of Matrixed Resource Allocation

Matrix organizations allow scientific facilities to share specialized personnel across projects, operations, maintenance, and strategic initiatives. Nominal staffing allocations, however, may not capture the schedule consequences of fragmented individual commitments, limited access to specialist groups, and intermittent availability of key decision makers. We developed a stochas- tic, daily-time-step simulation of a hypothetical medium-sized accelerator-facility project com- prising sequential phases and parallel tasks. Each task requires role-specific work measured in FTE-days. Ordinary personnel may be unavailable because they contribute concurrently to other institutional activities, while designated key roles have independently specified daily un- availability probabilities. An organization-wide priority factor scales the number of people from each functional group who can effectively contribute to the project. It is interpreted as a composite proxy for project access and workforce fragmentation across competing commit- ments. We examined project completion time as a function of this factor and Project Lead unavailability using 100 Monte Carlo runs per condition. Increasing priority factor from 0.1 to 1.0 reduced median completion time from 1708.5 days (interquartile range 1681.5–1735.25) to 390 days (interquartile range 379–399). At priority factor = 0.1, increasing Project Lead unavailability from 0.5 to 0.9 increased median completion time from 1713.5 days (interquartile range 1691–1733.25) to 4,417 days (interquartile range 4271.75–4550.5). The model quantifies the commonly expected sensitivity of project schedules to fragmented resource commitments and limited coordination availability. Within this model, the results also indicate a possible threshold regime in which small increases in workforce availability yield only modest sched- ule improvements until sufficient capacity becomes accessible, after which project performance improves sharply. With further validation and calibration, this quantitative framework could support resource-allocation decisions during initial project planning and subsequent schedule rebaselining.

Bai, Mei [SLAC National Accelerator Laboratory (SL↗

Progress on the reevaluation and validation of the n+233U neutron cross sections

The set of 233 U resonance parameters of the ENDF/B-VIII.0 nuclear data library was adopted from the previous ENDF/B-VII.1 evaluation using the external levels to update the thermal values. Adoption of IAEA 2017 thermal standards ( σ f = 533.0 ± 2.2 b, σ c = 44.9 ± 0.9 b, and ν ‾ tot = 2.487 ± 0.011 ) and of the IAEA-recommended thermal-neutron induced prompt fission neutron spectrum (PFNS) with average PFNS energy of 2.030 ± 0.013 MeV requires a re-evaluation of 233 U neutron cross sections in the resolved resonance region. A newly produced evaluation is being tested on benchmarks carefully selected from the Handbook of International Criticality Safety Benchmark Experiments (ICSBEP) which are highly sensitive to 233 U data. An important goal of this work was to eliminate the strong negative gradient of the calculated effective multiplication factors with respect to the epithermal fission fraction observed in the validation of the ENDF/B-VIII.0 library for those assemblies. A significant improvement in integral performance of critical 233 U solutions is observed for the newly proposed evaluation. Further work addressing the fast neutron region is needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fermionic Partial Tomography via Classical Shadows

Here we propose a tomographic protocol for estimating any k-body reduced density matrix (k-RDM) of an n-mode fermionic state, a ubiquitous step in near-term quantum algorithms for simulating many-body physics, chemistry, and materials. Our approach extends the framework of classical shadows, a randomized approach to learning a collection of quantum-state properties, to the fermionic setting. Our sampling protocol uses randomized measurement settings generated by a discrete group of fermionic Gaussian unitaries, implementable with linear-depth circuits. We prove that estimating all k-RDM elements to additive precision ϵ requires on the order of ($^{n}_{k}$)k 3/2 log(n)/ϵ 2 repeated state preparations, which is optimal up to the logarithmic factor. Furthermore, numerical calculations show that our protocol offers a substantial improvement in constant overheads for k ≥ 2, as compared to prior deterministic strategies. We also adapt our method to particle-number symmetry, wherein the additional circuit depth may be halved at the cost of roughly 2–5 times more repetitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hot-wire calibration in subsonic/transonic flow regimes

A different approach for calibrating hot-wires, which simplifies the calibration procedure and reduces the tunnel run-time by an order of magnitude was sought. In general, it is accepted that the directly measurable quantities in any flow are velocity, density, and total temperature. Very few facilities have the capability of varying the total temperature over an adequate range. However, if the overheat temperature parameter, a(sub w), is used to calibrate the hot-wire then the directly measurable quantity, voltage, will be a function of the flow variables and the overheat parameter i.e., E = f(u,p,a(sub w), T(sub w)) where a(sub w) will contain the needed total temperature information. In this report, various methods of evaluating sensitivities with different dependent and independent variables to calibrate a 3-Wire hot-wire probe using a constant temperature anemometer (CTA) in subsonic/transonic flow regimes is presented. The advantage of using a(sub w) as the independent variable instead of total temperature, t(sub o), or overheat temperature parameter, tau, is that while running a calibration test it is not necessary to know the recovery factor, the coefficients in a wire resistance to temperature relationship for a given probe. It was deduced that the method employing the relationship E = f (u,p,a(sub w)) should result in the most accurate calibration of hot wire probes. Any other method would require additional measurements. Also this method will allow calibration and determination of accurate temperature fluctuation information even in atmospheric wind tunnels where there is no ability to obtain any temperature sensitivity information at present. This technique greatly simplifies the calibration process for hot-wires, provides the required calibration information needed in obtaining temperature fluctuations, and reduces both the tunnel run-time and the test matrix required to calibrate hotwires. Some of the results using the above techniques are presented in an appendix.

Nagabushana, K. A.↗

Lightweight Damage Tolerant, High-Temperature Radiators for Nuclear Power and Propulsion

NASA is increasingly emphasizing exploration to bodies beyond near-Earth orbit. New propulsion systems and new spacecraft are being built for these missions. As the target bodies get further out from Earth, high energy density systems, e.g., nuclear fusion, for propulsion and power will be advantageous. The mass and size of these systems, including supporting systems such as the heat exchange system, including thermal radiators, will need to be as small as possible. Conventional heat exchange systems are a significant portion of the total thermal management mass and size. Nuclear electric propulsion (NEP) is a promising option for high-speed, in-space travel due to the high energy density of nuclear fission power sources and efficient electric thrusters. Heat from the reactor is converted to power for use in propulsion or for system power. The heat not used in the power conversion is then radiated to space as shown in figure 1. Advanced power conversion technologies will require high operating temperatures and would benefit from lightweight radiator materials. Radiator performance dictates power output for nuclear electric propulsion systems. Pitch-based carbon fiber materials have the potential to offer significant improvements in operating temperature, thermal conductivity, and mass. These properties combine to allow significant decreases in the total mass of the radiators and significant increases in the operating temperature of the fins. A Center-funded project at NASA Marshall Space Flight Center has shown that high thermal conductivity, woven carbon fiber fins with no matrix material, can be used to dissipate waste heat from NEP systems and because of high specific power (kW/kg), will require less mass and possibly less total area than standard metal and composite radiator fins for radiating the same amount of heat. This project uses an innovative approach to reduce the mass and size required for the thermal radiators to the point that in-space NEP and power is enabled. High thermal conductivity carbon fibers are lightweight, damage tolerant, and can be heated to high temperature. Areal densities in the NASA set target range of 2 to 4 kg/m2 (for enabling NEP) are achieved and with specific powers (kW/kg) a factor of about 7 greater than conventional metal fins and about 1.5 greater than carbon composite fins. Figure 2 shows one fin under test. All tests were done under vacuum conditions.

Craven, Paul D.↗

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗

A random search algorithm for laboratory computers

The small laboratory computer is ideal for experimental control and data acquisition. Postexperimental data processing is often performed on large computers because of the availability of sophisticated programs, but costs and data compatibility are negative factors. Parameter optimization can be accomplished on the small computer, offering ease of programming, data compatibility, and low cost. A previously proposed random-search algorithm ('random creep') was found to be very slow in convergence. A method is proposed (the 'random leap' algorithm) which starts in a global search mode and automatically adjusts step size to speed convergence. A FORTRAN executive program for the random-leap algorithm is presented which calls a user-supplied function subroutine. An example of a function subroutine is given which calculates maximum-likelihood estimates of receiver operating-characteristic parameters from binary response data. Other applications in parameter estimation, generalized least squares, and matrix inversion are discussed.

Curry, R. E.↗