Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DYNAMICS”

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,657 records · Page 92

Decoding polymer self-dynamics using a two-step approach

The self-correlation function and corresponding self-intermediate scattering function in Fourier space are important quantities for describing the molecular motions of liquids. This work draws attention to a largely overlooked issue concerning the analysis of these space-time density-density correlation functions of polymers. Here, we show that the interpretation of non-Gaussian behavior of polymers is generally complicated by intrachain averaging of distinct self-dynamics of different segments. By the very nature of the mathematics involved, the averaging process not only conceals critical dynamical information, but also contributes to the observed non-Gaussian dynamics. To fully expose this issue and provide a thorough benchmark of polymer self-dynamics, we perform analyses of coarse-grained molecular dynamics simulations of linear and ring polymer melts as well as several theoretical models using a “two-step” approach, where interchain and intrachain averagings of segmental self-dynamics are separated. While past investigations primarily focused on the average behavior, our results indicate that a more nuanced approach to polymer self-dynamics is clearly required.

36 MATERIALS SCIENCE↗

Lyapunov exponent as a metric for assessing the dynamic content and predictability of large-eddy simulations

Metrics used to assess the quality of large-eddy simulations commonly rely on a statistical assessment of the solution. While these metrics are valuable, a dynamic measure is desirable to further characterize the ability of a numerical simulation for capturing dynamic processes inherent in turbulent flows. To address this issue, a dynamic metric based on the Lyapunov exponent is proposed which assesses the growth rate of the solution separation. Specifically, this metric is applied to two turbulent flow configurations: forced homogeneous isotropic turbulence and a turbulent jet diffusion flame. First, it is shown that, despite the direct numerical simulation (DNS) and large-eddy simulation (LES) being high-dimensional dynamical systems with O(107) degrees of freedom, the separation growth rate qualitatively behaves like a lower-dimensional dynamical system, in which the dimension of the Lyapunov system is substantially smaller than the discretized dynamical system. Second, a grid refinement analysis of each configuration demonstrates that as the LES filter width approaches the smallest scales of the system the Lyapunov exponent asymptotically approaches a plateau. Third, a small perturbation is superimposed onto the initial conditions of each configuration, and the Lyapunov exponent is used to estimate the time required for divergence, thereby providing a direct assessment of the predictability time of simulations. By comparing inert and reacting flows, it is shown that combustion increases the predictability of the turbulent simulation as a result of the dilatation and increased viscosity by heat release. The predictability time is found to scale with the integral time scale in both the reacting and inert jet flows. Fourth, an analysis of the local Lyapunov exponent is performed to demonstrate that this metric can also determine flow-dependent properties, such as regions that are sensitive to small perturbations or conditions of large turbulence within the flow field. Finally, it is demonstrated that the global Lyapunov exponent can be utilized as a metric to determine if the computational domain is large enough to adequately encompass the dynamic nature of the flow.

42 ENGINEERING↗

Dynamic entity formed by protein and its hydration water

The interaction between protein and water plays a pivotal role in shaping the structure, dynamics, and function of biomacromolecules. A comprehensive understanding of this intricate interplay necessitates a systematic evaluation of interaction strength and its consequential impact on the dynamics of proteins and water across diverse protein systems. Despite numerous works on understanding the dynamics of water and proteins and the coupling between them, there are still unanswered questions. Here, we combine neutron scattering and isotope labeling to probe the dynamics of proteins and their hydration water in a variety of protein systems. We consider proteins of different structures and varying thermostability as well as proteins within living cells with distinct growth temperatures. Simultaneous characterization of protein and hydration water dynamics across diverse systems was achieved. Moreover, we performed water sorption isothermal measurements on three representative proteins to correlate the observed dynamics with the strength of the interaction energies governing each system. The experimental results underscore that proteins manifesting stronger attractive interactions with water display diffusionlike dynamics with higher flexibility upon hydration, concomitant with a reduced mobility in hydration water. Significantly, our findings suggest that, in fact, it is the interaction between protein and its hydration water that facilitates the transfer of mobility from water to protein, with stronger interactions correlating to greater protein flexibility and slower hydration water diffusion. Published by the American Physical Society 2024

Ye, Yongfeng (ORCID:0009000792374198)↗

TH-AB-207A-11: Anatomical and Functional Assessment of Coronary Artery Disease Using Low-Dose Whole-Organ Dynamic Computed Tomography

Purpose: To evaluate the accuracy of a low-dose, first-pass-analysis (FPA) dynamic computed tomography angiography and perfusion (CTAP) technique, for whole-organ anatomical and functional assessment of coronary artery disease (CAD). Methods: An angioplasty balloon was advanced into the left anterior descending (LAD) coronary artery of five swine (35–45 kg) to induce several levels of stenosis at maximal hyperemia (intracoronary adenosine, 240 µg/min). Reference fluorescence microspheres and intravenous contrast (370 mg/mL iodine, 25 mL, 7 mL/s) were injected centrally and dynamic imaging was performed using a 320-slice CT scanner at 100 kVp and 200 mA. Twenty volume scans were acquired per stenosis level to capture complete aortic and myocardial enhancement curves, but only two volume scans were used for whole-organ dynamic FPA CTAP measurement. All CTAP measurements in the LAD were compared to the reference microsphere perfusion measurements using linear regression, concordance correlation, and Bland-Altman analysis. Results: The result of dynamic FPA CTAP measurement in the LAD was in good agreement with the reference microsphere perfusion measurement (P-CTAP = 1.01 P-MICRO + 0.16, R{sup 2} = 0.95). The root mean square error (RMSE) and difference (RMSD) of measurement were 0.51 mL/min/g and 0.47 mL/min/g, respectively. Bland-Altman analysis demonstrated negligible systematic measurement bias. Additionally, the concordance correlation coefficient (CCC) was found to be ρ = 0.97, indicating excellent agreement between dynamic FPA CTAP measurement and the reference microsphere perfusion measurement. Lastly, the effective dose of the proposed technique using the “simulated” two-volume scan CTAP acquisition protocol was 2.6 mSv; much lower than the ∼10 mSv effective dose of current dynamic CTP techniques alone. Conclusion: The results indicate the potential for significant improvements in CAD assessment through low-dose, quantitative dynamic FPA CTAP. Such improvements are afforded through whole-organ CT scanning technology, and have the potential to improve patient outcomes and lead to healthier patient lives. Conflict of Interest (only if applicable): Grant funding from Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

Emergent Subpopulation Behavior Uncovered with a Community Dynamic Metabolic Model of Escherichia coli Diauxic Growth

ABSTRACT Microbes have adapted to greatly variable environments in order to survive both short-term perturbations and permanent changes. A classical and yet still actively studied example of adaptation to dynamic environments is the diauxic shift of Escherichia coli , in which cells grow on glucose until its exhaustion and then transition to using previously secreted acetate. Here we tested different hypotheses concerning the nature of this transition by using dynamic metabolic modeling. To reach this goal, we developed an open source modeling framework integrating dynamic models (ordinary differential equation systems) with structural models (metabolic networks) which can take into account the behavior of multiple subpopulations and smooth flux transitions between time points. We used this framework to model the diauxic shift, first with a single E. coli model whose metabolic state represents the overall population average and then with a community of two subpopulations, each growing exclusively on one carbon source (glucose or acetate). After introduction of an environment-dependent transition function that determined the balance between subpopulations, our model generated predictions that are in strong agreement with published data. Our results thus support recent experimental evidence that diauxie, rather than a coordinated metabolic shift, would be the emergent pattern of individual cells differentiating for optimal growth on different substrates. This work offers a new perspective on the use of dynamic metabolic modeling to investigate population heterogeneity dynamics. The proposed approach can easily be applied to other biological systems composed of metabolically distinct, interconverting subpopulations and could be extended to include single-cell-level stochasticity. IMPORTANCE Escherichia coli diauxie is a fundamental example of metabolic adaptation, a phenomenon that is not yet completely understood. Further insight into this process can be achieved by integrating experimental and computational modeling methods. We present a dynamic metabolic modeling approach that captures diauxie as an emergent property of subpopulation dynamics in E. coli monocultures. Without fine-tuning the parameters of the E. coli core metabolic model, we achieved good agreement with published data. Our results suggest that single-organism metabolic models can only approximate the average metabolic state of a population, therefore offering a new perspective on the use of such modeling approaches. The open source modeling framework that we provide can be applied to model general subpopulation systems in more-complex environments and can be extended to include single-cell-level stochasticity.

59 BASIC BIOLOGICAL SCIENCES↗

Features of the dynamics of particles and fields at cyclotron resonances

Some features of the dynamics of particles and fields at cyclotron resonances have been discussed with the focus on chaotic dynamical regimes. It has been shown that the known criterion of the transition of the regular dynamics of particles to chaotic dynamics at cyclotron resonances sometimes describes this transition incorrectly. The reason for such a feature of the criterion has been revealed. The anomalous sensitivity of the dynamics of particles to external fluctuations at autoresonance has been analyzed. A theory of excitation of electromagnetic waves by a beam of phased oscillators under the conditions of isolated nonlinear cyclotron resonance has been developed. It has been shown that the chaotic dynamical regime is due to the periodic change in the structure of the phase portrait of particles in the wave field. It has been shown that higher moments can play a more significant role than lower moments in almost all chaotic dynamical regimes at cyclotron resonances. In this case, the known kinetic diffusion equations should be generalized with the inclusion of these higher moments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermodynamic versus Dynamic Controls on Extreme Precipitation in a Warming Climate from the Community Earth System Model Large Ensemble

The moisture budget is evaluated as a function of the probability distribution of precipitation for the end of the twentieth century and projected end of the twenty-first century in the Community Earth System Model Large Ensemble. For a given precipitation percentile, a conditional moisture budget equation relates precipitation minus evaporation ( P − E) to vertical moisture transport, horizontal moisture advection, and moisture storage. At high percentiles, moisture advection and moisture storage cancel and evaporation is negligible, so that precipitation is approximately equal to vertical moisture transport, and likewise for projected changes. Therefore, projected changes to extreme precipitation are approximately equal to the sum of thermodynamic and dynamic tendencies, representing changes to the vertical profiles of moisture content and mass convergence, respectively. The thermodynamic tendency is uniform across percentiles and regions as an intensification of the hydrological cycle, but the dynamic tendency is more complex. For extreme events, per degree of warming, in the mid-to-high latitudes the dynamic tendency is small, so that precipitation approximately scales by the Clausius–Clapeyron 7% K −1 increase. In the subtropics, a drying tendency originating from dynamics offsets the thermodynamic wetting tendency, with the net effect on precipitation varying among regions. The effect of this dynamic drying decreases with increasing percentile. In the deep tropics, a positive dynamic tendency occurs with magnitude similar to or greater than the positive thermodynamic tendency, resulting in generally a 10%–15% K −1 precipitation increase, and with a >25% K −1 increase over the tropical east Pacific. This reinforcing dynamical tendency increases rapidly for high percentiles.

54 ENVIRONMENTAL SCIENCES↗

Diagnosing the Dynamics of Observed and Simulated Ecosystem Gross Primary Productivity with Time Causal Information Theory Quantifiers

Data analysis and model-data comparisons in the environmental sciences require diagnostic measures that quantify time series dynamics and structure, and are robust to noise in observational data. This paper investigates the temporal dynamics of environmental time series using measures quantifying their information content and complexity. The measures are used to classify natural processes on one hand, and to compare models with observations on the other. The present analysis focuses on the global carbon cycle as an area of research in which model-data integration and comparisons are key to improving our understanding of natural phenomena. We investigate the dynamics of observed and simulated time series of Gross Primary Productivity (GPP), a key variable in terrestrial ecosystems that quantifies ecosystem carbon uptake. However, the dynamics, patterns and magnitudes of GPP time series, both observed and simulated, vary substantially on different temporal and spatial scales. Here we demonstrate that information content and complexity, or Information Theory Quantifiers (ITQ) for short, serve as robust and efficient data-analytical and model benchmarking tools for evaluating the temporal structure and dynamical properties of simulated or observed time series at various spatial scales. At continental scale, we compare GPP time series simulated with two models and an observations-based product. This analysis reveals qualitative differences between model evaluation based on ITQ compared to traditional model performance metrics, indicating that good model performance in terms of absolute or relative error does not imply that the dynamics of the observations is captured well. Furthermore, we show, using an ensemble of site-scale measurements obtained from the FLUXNET archive in the Mediterranean, that model-data or model-model mismatches as indicated by ITQ can be attributed to and interpreted as differences in the temporal structure of the respective ecological time series. At global scale, our understanding of C fluxes relies on the use of consistently applied land models. Here, we use ITQ to evaluate model structure: The measures are largely insensitive to climatic scenarios, land use and atmospheric gas concentrations used to drive them, but clearly separate the structure of 13 different land models taken from the CMIP5 archive and an observations-based product. In conclusion, diagnostic measures of this kind provide data-analytical tools that distinguish different types of natural processes based solely on their dynamics, and are thus highly suitable for environmental science applications such as model structural diagnostics.

54 ENVIRONMENTAL SCIENCES↗

Investigations of Beam Dynamics Issues at Current and Future Hadron Accelerators (Final Report)

There is a synergy between the fields of Beam Dynamics (BD) in modern particle accelerators and Applied Mathematics (AMa). We have formulated significant problems in BD and have developed and applied tools within the contexts of dynamical systems, topological methods, numerical analysis and scientific computing, probability and stochastic processes, and mathematical statistics. We summarize the three main areas of our AMa work since 2011. First, we continued our study of Vlasov-Maxwell systems. Previously, we developed a state of the art algorithm and code (VM3@A) to calculate coherent synchrotron radiation in single pass systems. In this cycle we carefully analyzed the major expense, namely the integral-over-history (IOH), and developed two approaches to speed up integration. The first strategy uses a representation of the Bessel function J0 in terms of exponentials. The second relies on “local sequences” developed recently for radiation boundary conditions, which are used to reduce computational domains. Although motivated by practicality, both strategies involve interesting and rather deep analysis and approximation theory. As an alternative to VM3@A, we are integrating Maxwell’s equations by a time-stepping method, bypass- ing the IOH, using a Discontinuous Galerkin (DG) method. DG is a generalization of Finite Element and Finite Volume methods. It is spectrally convergent, unlike the commonly used Finite Difference methods, and can handle complicated vacuum chamber geometries. We have applied this in several contexts and have obtained very nice results including an explanation of an experiment at the Canadian Light Source, where the geometry is quite complex. Second, we continued our study of spin dynamics in storage rings. There is much current and proposed activity where spin polarized beams are being used in testing the Standard Model and its modifications. Our work has focused on invariant spin fields (ISFs) and amplitude dependent spin tunes (ADSTs), which are essential for estimating beam polarization. Several algorithms have been developed since the 1980s for computing the ISF, among them the Heinemann- Hoffstaetter method of stroboscopic averaging (SA) which is implemented in the code SPRINT. SA, which computes the ISF by using spin tracking data, can find the ISF to computer precision if it exists and thus can give evidence for existence of the ISF. Central to our work is resolving the ISF conjecture, which says that, off orbital resonance, an ISF exists. Thus Heinemann developed, in his 2010 PhD thesis, a new framework which unifies and generalizes the notions of ISF and ADST by using bundle theory. This lead to a long paper which was a major collaborative effort during the recent cycle. In a nutshell, our bundle approach elegantly unifies the dynamics of spin-1/2 and spin-1 particles, e.g., protons and deuterons. In fact it is well known that these two kinds of dynamics are driven by the same spin transfer matrix and in our approach one sees the deeper reason for this: the spin transfer matrix carries the natural dynamics of a principal bundle whereas the difference between the spin-1/2 and spin-1 dynamics lies in their different geometrical situation, i.e., different underlying associated bundles. Thus one arrives at new results for polarized beams, among them the Invariant Reduction Theorem (IRT) and the Cross Section Theorem (CST). The IRT gives a necessary and sufficient condition for the ISF to exist. The SA technique revealed, 20 years ago, that each ISF can be viewed as a complex agglomerate of spin tracking data. The bundle approach goes one step further by using the IRT and the CST to glue spin tracking data into agglomerates which are candidates for ISFs. We gain insight because the “good” agglomerates, in the presence of an ISF, are very different from the “bad” ones. Finally we mention that the bundle approach has analogies to the approach used in geometrical Yang-Mills theory. Third, we studied X-Ray Free Electron Lasers (FELs), which are electron accelerators producing coherent undulator radiation over a wide range of frequencies from microwaves to x-rays. The photon beams produced in FEL undulators are used to study material samples in biology, material science etc. We developed a mathematical analysis, based on the 6D Lorentz system, of energetic electrons moving through a planar undulator excited by a Maxwell traveling wave field of wavelength λ. Our Method of Averaging perturbation analysis yields non-resonant and near-to-resonant normal form approximations as a function of λ, which we present in two averaging theorems. We prove the theorems in detail, error bounds and giving a tutorial on mathematically rigorous perturbation theory in a context where proofs are easily understood. To our knowledge the planar problem has not been analyzed with the generality here nor has the standard FEL pendulum system, which appears on resonance, been derived with error bounds. In addition to the domains of validity of the normal forms we obtain new insights that require further study, including a more general low gain theory. With a firm foundation in the non-collective case above we are studying the 3D collective case from start up from noise to high gain and saturation. We have formulated the noise start up as a problem of going from the microscopic Klimontovich-Maxwell to the macroscopic Vlasov-Maxwell with a Vlasov correction term. In the 1D setting, we seek an alternative to the phenomenological Vlasov approach which models shot noise by a perturbation on an initial “smooth” density. The 1D wave equation with a Klimontovich source is often the starting point for the 1D FEL high gain theory. We have a new representation of solutions which may lead to new insights.

43 PARTICLE ACCELERATORS↗

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗

Understanding Dynamics and Thermodynamics of ENSO and Its Complexity Simulated by E3SM and Other Climate Models

Despite the seeming success of most state-of-the-art climate models in simulating the El Niño-Southern Oscillation (ENSO), there is strong evidence that models achieve realistic levels of ENSO activity do so often for wrong reasons. This is owing to an often occurred near cancelation of large errors in terms of contributions to ENSO growth rate from coupled dynamic and thermodynamic feedback processes. Climate models remain deficient in simulating the observed ENSO’s spatial and temporal complexity that involves interplays of coupled dynamic and thermodynamic feedbacks, interactions across multiple scales, nonlinear processes in the tropical atmosphere and ocean system, biases in mean sate and physical processes, and influences external to equatorial Pacific coupled ENSO dynamics. Our proposed research aims at advancing predictive and process-level understandings of ENSO simulated in E3SM and other climate models under current and future climate conditions with two main objectives: (i) better understanding the aforementioned broad range interactive processes and sources that control fundamental properties of ENSO in E3SM and CIMP6 outputs as well as in observational (reanalysis) data sets, using a hierarchical of coupled dynamical frameworks consisting of theoretical analysis, intermediate complexity modeling, and coupled dynamic diagnostics; (ii) to use this understanding to explore pathways towards improving E3SM’s capability of simulating ENSO and its complexity. More specifically, we will focus on four main thrusts of research: (1) ENSO’s dynamic and thermodynamic feedbacks; (2) the across-scale interactions of ENSO with annual cycle and MJO/WWB/TIW (Madden Julian Oscillation/Westerly Wind Burst/Tropical Instability Wave) activity; (3) key nonlinear processes of ENSO involving atmospheric convective thresholds, nonlinear ocean dynamic heating, and thermocline outcropping; and (4) the impacts of climate mean-state biases/changes and perturbed physical processes on simulated ENSO and its complexity.

54 ENVIRONMENTAL SCIENCES↗

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Heterogeneous Polymer Dynamics Explored Using Static 1H NMR Spectra

NMR spectroscopy continues to provide important molecular level details of dynamics in different polymer materials, ranging from rubbers to highly crosslinked composites. It has been argued that thermoset polymers containing dynamic and chemical heterogeneities can be fully cured at temperatures well below the final glass transition temperature (Tg). In this paper, we described the use of static solid-state 1H NMR spectroscopy to measure the activation of different chain dynamics as a function of temperature. Near Tg, increasing polymer segmental chain fluctuations lead to dynamic averaging of the local homonuclear proton-proton (1H-1H) dipolar couplings, as reflected in the reduction of the NMR line shape second moment (M2) when motions are faster than the magnitude of the dipolar coupling. In general, for polymer systems, distributions in the dynamic correlation times are commonly expected. To help identify the limitations and pitfalls of M2 analyses, the impact of activation energy or, equivalently, correlation time distributions, on the analysis of 1H NMR M2 temperature variations is explored. It is shown by using normalized reference curves that the distributions in dynamic activation energies can be measured from the M2 temperature behavior. An example of the M2 analysis for a series of thermosetting polymers with systematically varied dynamic heterogeneity is presented and discussed.

36 MATERIALS SCIENCE↗

Gaseous Dynamical Friction in Presence of Black Hole Radiative Feedback

Dynamical friction is thought to be a principal mechanism responsible for orbital evolution of massive black holes (MBHs) in the aftermath of galactic mergers and an important channel for formation of gravitationally bound MBH binaries. We use 2D radiative hydrodynamic simulations to investigate the efficiency of dynamical friction in the presence of radiative feedback from an MBH moving through a uniform density gas. We find that ionizing radiation that emerges from the innermost parts of the MBH’s accretion flow strongly affects the dynamical friction wake and renders dynamical friction inefficient for a range of physical scenarios. MBHs in this regime tend to experience positive net acceleration, meaning that they speed up, contrary to the expectations for gaseous dynamical friction in absence of radiative feedback. The magnitude of this acceleration is however negligibly small and should not significantly alter the velocity of MBHs over relevant physical timescales. Our results suggest that suppression of dynamical friction is more severe at the lower mass end of the MBH spectrum which, compounded with inefficiency of the gas drag for lower mass objects in general, implies that <10{sup 7} M{sub ⊙} MBHs have fewer means to reach the centers of merged galaxies. These findings provide formulation for a sub-resolution model of dynamical friction in presence of MBH radiative feedback that can be easily implemented in large scale simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

The shock tube as a tool for dynamic testing of fabrics

Specimens of fabric and Mylar were loaded to failure both statically and dynamically. Static tests were accomplished by pressurizing the specimens and dynamic tests were conducted by firing a shock tube. The results of the study indicate that seam efficiency is degraded by dynamic loads and that a need exists to investigate a seaming procedure which gives an optimum dynamic efficiency. The study also illustrates the lack of correspondence between static and dynamic breaking strengths thereby indicating that fabrics which have dynamic applications should be evaluated using dynamic testing procedures such as the shock tube loading used in the present program.

Madden, R.↗

Dynamic fracture toughnesses of reaction-bonded silicon nitride

The room-temperature dynamic fracture response of reaction-bonded silicon nitride is investigated using a hybrid experimental-numerical procedure. In this procedure, experimentally determined crack velocities are utilized to drive a dynamic finite-element code or dynamic finite-difference code in its generation mode in order to extract numerically the dynamic stress intensity factor of the fracturing specimen. Results show that the dynamic fracture toughness vs crack velocity relations of the two reaction-bonded silicon nitrides do not follow the general trend in those relations of brittle polymers and steel. A definite slow crack velocity during the initial phase of dynamic crack propagation is observed in reaction-bonded silicon nitride, which results in a nonunique dynamic fracture toughness vs crack velocity relation. In addition, it is found that a propagating crack will continue to propagate under a static stress intensity factor substantially lower than K(IC).

Kobayashi, A. S.↗

The dynamic fission instability and the origin of the Moon

A theory for the formation of the Moon which involves the dynamic fission of a rapidly rotating protoplanet, which might then result in the formation of the Earth and the Moon is discussed. The fission hypothesis was originally based on analytic, linearized models of the growth of asymmetry in homogenous bodies. The fully nonlinear evolution of the dynamic instability in inviscid, compressible bodies was calculated by numerical techniques. It was found that the dynamic instability degenerates into the ejection of a ring of matter with a substantial fraction of the mass, leaving behind a central body with most of the mass. The linearized analytical approach and the numerical approach were used to show that dynamic fission probably does not occur in rocky protoplanets. The numerical calculations are performed with a fully three dimensional hydrodynamical code, which allows the nonlinear, time evolution of the instability to be followed. Sequences of uniformly rotating equilibria were constructed and are used as the initial models for the fission calculations. An initially imposed asymmetry consisting of a 10% binary perturbation in the density was found to disappear on the rotational period time scale. No dynamic instability occurred. This result are verified by including the velocity dissipation terms in the linearized analysis of the stability of a Maclaurin spheroid: the dynamic instability disappears when the simulated viscous dissipation terms are included. It is concluded that any rocky body, even with considerable partial melt or a molten core, should be stable to dynamic fission; any rotational instability that occurs can only result in equatorial mass loss.

Boss, A. P.↗

Predicted effect of dynamic load on pitting fatigue life for low-contact-ratio spur gears

How dynamic load affects the surface pitting fatigue life of external spur gears was predicted by using the NASA computer program TELSGE. Parametric studies were performed over a range of various gear parameters modeling low-contact-ratio involute spur gears. In general, gear life predictions based on dynamic loads differed significantly from those based on static loads, with the predictions being strongly influenced by the maximum dynamic load during contact. Gear mesh operating speed strongly affected predicted dynamic load and life. Meshes operating at a resonant speed or one-half the resonant speed had significantly shorter lives. Dynamic life factors for gear surface pitting fatigue were developed on the basis of the parametric studies. In general, meshes with higher contact ratios had higher dynamic life factors than meshes with lower contact ratios. A design chart was developed for hand calculations of dynamic life factors.

Lewicki, David G.↗