Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependence analysis”

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 883 records · Page 49

Analysis of the Hessian for Aerodynamic Optimization: Inviscid Flow

In this paper we analyze inviscid aerodynamic shape optimization problems governed by the full potential and the Euler equations in two and three dimensions. The analysis indicates that minimization of pressure dependent cost functions results in Hessians whose eigenvalue distributions are identical for the full potential and the Euler equations. However the optimization problems in two and three dimensions are inherently different. While the two dimensional optimization problems are well-posed the three dimensional ones are ill-posed. Oscillations in the shape up to the smallest scale allowed by the design space can develop in the direction perpendicular to the flow, implying that a regularization is required. A natural choice of such a regularization is derived. The analysis also gives an estimate of the Hessian's condition number which implies that the problems at hand are ill-conditioned. Infinite dimensional approximations for the Hessians are constructed and preconditioners for gradient based methods are derived from these approximate Hessians.

Arian, Eyal↗

The extracellular vesicle proteomes of Sorghum bicolor and Arabidopsis thaliana are partially conserved

Abstract Plant extracellular vesicles (EVs) are membrane-bound organelles involved mainly in intercellular communications and defense responses against pathogens. Recent studies have demonstrated the presence of proteins, nucleic acids including small RNAs, and lipids along with other metabolites in plant EVs. Here, we describe the isolation and characterization of EVs from sorghum (Sorghum bicolor). Nanoparticle tracking analysis, dynamic light scattering, and cryo-electron tomography showed the presence of a heterogeneous population of EVs isolated from the apoplastic wash of sorghum leaves. Cryo-electron microscopy revealed that EVs had a median size of 110 nm and distinct populations of vesicles with single or multiple lipid bilayers and low or high amounts of contents. The heterogeneity was further supported by data showing that only a subset of EVs that were stained with a membrane dye, Potomac Gold, were also stained with the membrane-permeant esterase-dependent dye, calcein acetoxymethyl ester. Proteomic analysis identified 437 proteins that were enriched in multiple EV isolations, with the majority of these also found in the EV proteome of Arabidopsis (Arabidopsis thaliana). These data suggest a partial conservation of EV contents and function between the monocot, sorghum, and a distantly related eudicot, Arabidopsis.

36 MATERIALS SCIENCE↗

A new waveform analysis technique to extract good energy and position resolution from a dual-axis duo-lateral position-sensitive detector

The dual-axis duo-lateral position-sensitive silicon detector was developed to detect charged particles with high quality position and energy resolution. When these detectors were used with conventional signal processing electronics, an empirically determined correction was used to improve energy resolution. In this work, the waveforms from the detector after preamplification are studied in detail to investigate position information contained in the waveforms. A 7.22 MeV/nucleon alpha particle beam was impinged directly on a masked dual-axis duo-lateral detector. Data obtained using a 228 Th alpha particle source was also used. Here, by studying the waveform characteristics that give rise to the position-dependent distortions, a new summed trigger analysis method has been developed to significantly improve linearity in position reconstruction without sacrificing energy resolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Fluctuations and Multifractally enhanced Superconductivity in Disordered Thin Films

The interplay between electron–electron interactions and weak localization (or anti-localization) phenomena in two-dimensional systems can significantly enhance the superconducting transition temperature. For this work, we develop the theory of quantum fluctuations within such multifractally enhanced superconducting states in thin films. In conditions of weak disorder, we employ the Finkel’stein nonlinear sigma model to derive an effective action for the superconducting order parameter and the quasiclassical Green’s function, meticulously accounting for the influence of quantum fluctuations. This effective action, applicable for interactions of any strength, reveals the critical role of well-known collective modes in a dirty superconductor, and its saddle-point analysis leads to modified Usadel and gap equations. These equations comprehensively incorporate the renormalizations stemming from the interplay between interactions and disorder, resulting in the non-trivial energy dependence of the gap function. Notably, our analysis establishes a direct relation between the self-consistent gap equation at the superconducting transition temperature and the known renormalization group equations for interaction parameters in the normal state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Facet-dependent electron transfer induces distinct arsenic reallocations on hematite

The interfacial electron transfer (ET) between electron shuttling compounds and iron (Fe) oxyhydroxides plays a crucial role in the reductive dissolution of Fe minerals and the fate of surface-bound arsenic (As). However, the impact of exposed facets of highly crystalline hematite on reductive dissolution and As immobilization is poorly understood. In this study, we systematically investigated the interfacial processes of the electron shuttling compound cysteine (Cys) on various facets of hematite and the reallocations of surface-bound As(III) or As(V) on the respective surfaces. Our results demonstrate that the ET process between Cys and hematite generates Fe(II) and leads to reductive dissolution, with more Fe(II) generated on {001} facets of exposed hematite nanoplates (HNPs). Reductive dissolution of hematite leads to significantly enhanced As(V) reallocations on hematite. Nevertheless, upon the addition of Cys, a raipd release of As(III) can be halted by its prompt re-adsorption, leaving the extent of As(III) immobilization on hematite unchanged throughout the course of reductive dissolution. This is due to that Fe(II) can form new precipitates with As(V), a process that is facet-dependent and influenced by water chemistry. Electrochemical analysis reveals that HNPs exhibit higher conductivity and ET ability, which is beneficial for reductive dissolution and As reallocations on hematite. In conclusion, these findings highlight the facet-dependent reallocations of As(III) and As(V) facilitated by electron shuttling compounds and have implications for the biogeochemical processes of As in soil and subsurface environments.

54 ENVIRONMENTAL SCIENCES↗

Experiments on the enhancement of compressible mixing via streamwise vorticity. I - Optical measurements

The present study was initiated in an attempt to increase the compressible mixing rate of a supersonic jet injected into a co-flowing supersonic stream by adding streamwise vorticity to the injected flow. This flowfield was imaged using the planar laser scattering technique, and mixing rates were derived entirely from captured images. The results of this image analysis indicate that increases in mixing rate of up to 60 percent are possible through the addition of streamwise vorticity. Furthermore, the degree of mixing enhancement appears to be dependent on the amount of swirl added. Analysis of the instantaneous turbulent structure suggests that this increase in mixing is due to a modification of the largest turbulent structures in the mixing layer.

Naughton, J. W.↗

Migrational Instabilities in Particle Suspensions

This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.

Goddard, Joe D.↗

Analysis Ready Data in Analytics Optimized Data Stores for Analysis of Big Earth Data in the Cloud

Cloud computing offers the possibility of making the analysis of Big Data approachable for a wider community due to affordable access to computing power, an ecosystem of usable tools for parallel processing, and migration of many large datasets to archives in the cloud, allowing data-proximal computing. Generally, data analysis acceleration in the cloud comes from running multiple nodes in a split-combine-apply strategy. Data systems such as the Earth Observing System Data and Information System are in a position to "pre-split" the data by storing them in a data store that is optimized for data parallel computing, i.e., an Analytics-Optimized Data Store (AODS). A variety of approaches to AODS are possible, from highly scalable databases to scalable filesystems to data formats optimized for cloud access (e.g., zarr and cloud-optimized datasets), with the optimal choice dependent on both the types of analysis and the geospatial structure of the data. A key question is how much preprocessing of the data to do, both before splitting and as the first part of the apply step. Again, the geospatial structure of the data and the analysis type influence the decision, with the added complexity of the user type. Trans-disciplinary users who are not well-versed in the nuances of quality-filtering and georeferencing of remote sensing orbit/swath/scene data tend to ask for more highly processed data, relying on the data provider to make sensible decisions on preprocessing parameters. (This accounts for the popularity of "Level 3" gridded data, despite the lower spatial resolution it provides.) In this case, data can be preprocessed before the split, resulting in higher performance in the rest of the "apply" step, which can be transformative for use cases such as interactive data exploration at scale. Discipline researchers who are experienced with remote sensing data often prefer more flexibility in customizing the preprocessing data into Analysis Ready Data, resulting in more need for on-the-fly preprocessing.

Lynnes, Christopher↗

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rapid AI-based Dissection of Ashes using Raman and XRF Spectroscopy (RADAR-X)

Waste-to-energy (WTE) facilities incinerate ~35 million tons of municipal solid waste annually in the United States. The incineration process reduces the mass and the volume of the waste fraction by over 75 and 95%, respectively. The fraction left after incineration remains as ash residues and is referred to as WTE ash, compromising of bottom and fly ash. In the United States, ~10 million tons of WTE ashes are generated annually and predominantly landfilled because of the lack of secondary end-use pathways. This incurs a significant financial burden (landfilling costs) on U.S. WTE facilities and also results in the loss of materials to the landfill. The primary objective of this research is to understand better the elemental and mineralogical composition of WTE ashes from diverse origins and find composition dependent upcycling pathways for diverting these ashes from landfills. This primary objective was addressed through three research tasks –(Task I) An AI-led Multi-Modal Approach for Compositional Analysis, (Task II) Developing a Dissolution-Based Test for Real Time Analysis, and (Task III) Establishing composition-dependent end uses. The chemical composition of WTE ash is dependent on two factors, i.e., the input waste composition and the operational parameters of a WTE facility (combustion conditions). Amongst these two factors, the input waste composition will likely show spatial and temporal variation. As a result, the chemical composition of WTE ash will also fluctuate. To understand the spatial and temporal variation in WTE ash composition, in Task I, we collected 128 ash samples (62 bottom ash and 66 fly ash samples) from 11 WTE facilities located in 11 U.S. states and characterized them via X-ray Fluorescence, powder X-ray Diffraction, and Raman Spectroscopy. The findings from this extensive characterization work indicated that the key elements in WTE fly ashes are Ca, Cl (greater than 10 wt. %), Si, S, K, Zn ( between 1 and 10 wt. %), Mg, Al, P, Ti, Fe, Cu, Br, and Pb (between 0.1 and 1 wt. %). Similarly, the key elements in WTE bottom ash fraction finer than 45μm are Ca (greater than 10 wt. %), Mg, Al, Si, S, Cl, K, Ti, Fe, Zn (between 1 and 10 wt. %), P, V, Cr, Mn, Cu, Br, and Pb (between 0.1 and 1 wt. %). Here, we note that the dominant fraction of WTE bottom ash is the coarse fraction. The coarse WTE bottom ash fraction (rich in silicon) was not characterized in this study because of excessive grinding requirements and their unsuitability as a supplementary cementitious material due to their coarse nature. The elements in WTE bottom ashes are present as calcite, anhydrite, vaterite, hydroxyapatite, quartz, bassanite, gehlenite, akermanite, hydrocalumite, and portlandite. Similarly, the mineralogical species present in WTE fly ashes are calcium chloride hydroxide, halite, calcite, anhydrite, sylvite, hydrocalumite, vaterite, hannebachite, quartz, and bassanite. Temporal variation in ash composition may also result in significant fluctuations in chemical compositions. Therefore, a WTE facility may need to monitor the ash composition (elemental and mineralogical composition) in real time. In Task I, we evaluated the possibility of using a portable X-ray fluorescence (XRF) spectrometer to monitor the elemental composition in real-time. Specifically, we collected XRF data on identical specimens via a portable XRF spectrometer (low-end) and a lab-based benchtop XRF spectrometer (high-end). The collected data was used to train an A.I. algorithm (portable XRF data as an input and benchtop XRF data as an output) to predict accurate elemental composition using portable XRF data. Finally, we developed a 2-minute photobleaching protocol to monitor the mineralogical characteristics of WTE ashes via Raman spectroscopy. Overall, the activities in Task I improved our understanding of ash composition and developed techniques to monitor elemental and mineralogical composition in near real-time. Based on the findings of Task I, we find that WTE ashes exhibit wide variability in mineralogy. For ICP-based elemental analysis, all the mineralogical species in WTE ashes must be brought into solution. This is traditionally accomplished with acid digestion using a combination of multiple acids. However, acid digestion with multiple acids is time-consuming and often fails to ensure complete digestion of the ash matrix. To address this limitation, in Task II, we developed an alkali-fusion-based digestion protocol using a combination of lithium tetraborate, lithium metaborate, and their combinations as possible alkali fluxes for digesting WTE ashes entirely and rapidly. The validity of the developed method was evaluated on two standard ash specimens, i.e., SRM 1633c coal fly ash and BCR-176R incineration fly ash specimen. The findings suggest that the developed protocol can ensure complete digestion of elements such as Al, Ba, Ca, Cr, Cu, Mg, Mn, P, Sr, V, Zn, Be, K, and rare earth elements. The recent changes in the energy market towards renewables and increased metal recycling have resulted in reduced supplies of supplementary cementitious materials (coal fly ash and slag). Therefore, in Task III, we evaluated the possibility of employing WTE ashes as SCMs. As the chemical composition of WTE ashes varies temporally (on an hourly basis), there was also a need to develop tests that can evaluate the suitability of material to act as supplementary cementitious material rapidly, i.e., in a few minutes. Therefore, in Task II, we also developed a rapid test to assess the suitability of a material to act as an SCM in 5 minutes. This represents a significant advance over the state-of-the-art R 3 test, which takes ~144 hours. This test was initially validated on amorphous aluminosilicates, such as calcined clays, and could be extended to evaluate other industrial by-products, such as WTE ashes. In Task III, we evaluated the possibility of employing WTE ashes for two applications, i.e., as an SCM and a lime substitute for clay stabilization. The findings from Task I indicated that WTE ashes are enriched in chlorine and, therefore, cannot be used directly as an SCM due to corrosion-related risks and altered hydration kinetics. Accordingly, we developed an ash treatment protocol to reduce the solubility of chlorine-containing species in WTE ashes. The developed treatment protocol also immobilized lead in certain mineral forms. As a result of the treatment, WTE ashes can be used as SCMs without any corrosion or heavy metal leaching concerns. The second application examined in this study was clay stabilization. WTE ashes are calcium-rich and can be an adequate lime replacement for clay stabilization. Our findings reveal that the sum of the concentrations of Ca(OH) 2 and CaClOH controls the clay stabilization capability of WTE ashes. In summary, in this work, we evaluated the elemental and mineralogical characteristics of U.S. WTE ashes from diverse origins and developed tests to evaluate the chemical characteristics of these ashes in real time through a portable XRF and a benchtop Raman spectrometer. Based on the chemical characteristics of these ashes, we developed an ash treatment process to enable the use of WTE ashes as an SCM and also evaluated the possibility of employing these ashes for clay stabilization. Overall, the findings from this work enables the diversion of WTE ashes from landfills for multiple end-uses, i.e., as an SCM or a lime substitute for clay stabilization.

36 MATERIALS SCIENCE↗

Time dependent inflow-outflow boundary conditions for 2D acoustic systems

An analysis of the number and form of the required inflow-outflow boundary conditions for the full two-dimensional time-dependent nonlinear acoustic system in subsonic mean flow is performed. The explicit predictor-corrector method of MacCormack (1969) is used. The methodology is tested on both uniform and sheared mean flows with plane and nonplanar sources. Results show that the acoustic system requires three physical boundary conditions on the inflow and one on the outflow boundary. The most natural choice for the inflow boundary conditions is judged to be a specification of the vorticity, the normal acoustic impedance, and a pressure gradient-density gradient relationship normal to the boundary. Specification of the acoustic pressure at the outflow boundary along with these inflow boundary conditions is found to give consistent reliable results. A set of boundary conditions developed earlier, which were intended to be nonreflecting is tested using the current method and is shown to yield unstable results for nonplanar acoustic waves.

Watson, Willie R.↗

Drosophila geotaxis as a tool for the study of aging

Age dependent changes in geotaxis profiles were examined in 27 wild-type populations of Drosophila, representing a diversity of species, semispecies and strains. In addition, four strains of D. melanogaster were tested. Tests were carried out at a minimum of three test ages, and involve the use of a calibrated, adjustable inclined plane that can be set at any angle between 0 and 85 deg. Among selected lines, decline in geotactic response occurs later in the long lived flies than in the controls. Longer lived flies continue to show an increase in negative geotactic response through age 14 days. These results suggest that common processes may be influencing the rate of decline in geotactic response and longevity. Further analysis of the mechanisms underlying age dependent changes in geotaxis may reveal factors which influence the aging process itself. The use of geotaxis aging markers in a broad range of Drosophila species reflecting varying degrees of genetic relatedness is proposed to test the universality vs. specificity of aging processes.

Schnebel, Edgar M.↗

An improved error assessment for the GEM-T1 gravitational model

Several tests were designed to determine the correct error variances for the Goddard Earth Model (GEM)-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment. These calibrations (utilizing each of the major data subsets within the solution) yield very stable calibration factors which vary by approximately 10 percent over the range of tests employed. Measurements of gravity anomalies obtained from altimetry were also used directly as observations to show that GEM-T1 is calibrated. The mathematical representation of the covariance error in the presence of unmodeled systematic error effects in the data is analyzed and an optimum weighting technique is developed for these conditions. This technique yields an internal self-calibration of the error model, a process which GEM-T1 is shown to approximate.

Lerch, F. J.↗

Analytical Modeling of the High Strain Rate Deformation of Polymer Matrix Composites

The results presented here are part of an ongoing research program to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric matrix materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical 5 plasticity theory definitions of effective stress and effective plastic strain are modified by applying variations of the Drucker-Prager yield criterion. To verify the revised formulation, the shear and tensile deformation of a representative toughened epoxy is analyzed across a wide range of strain rates (from quasi-static to high strain rates) and the results are compared to experimentally obtained values. For the analyzed polymers, both the tensile and shear stress-strain curves computed using the analytical model correlate well with values obtained through experimental tests. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. In the micromechanics, the unit cell is divided up into a number of independently analyzed slices, and laminate theory is then applied to obtain the effective deformation of the unit cell. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite (composed using the representative polymer analyzed for the correlation of the polymer constitutive equations) for several fiber orientation angles across a variety of strain rates. The computed values compare favorably to experimentally obtained results.

Goldberg, Robert K.↗

Stereoinversion via Alcohol Dehydrogenases Enables Complete Catabolism of β-1-Type Lignin-Derived Aromatic Isomers

Sphingobium sp. strain SYK-6 is an efficient aromatic catabolic bacterium that can consume all four stereoisomers of 1,2-diguaiacylpropane-1,3-diol (DGPD), which is a ring-opened β-1-type dimer. Recently, LdpA-mediated catabolism of erythro-DGPD was reported in SYK-6, but the catabolic pathway for threo-DGPD was as yet unknown. Here, in this study, we elucidated the catabolism of threo-DGPD, which proceeds through conversion to erythro-DGPD. When threo-DGPD was incubated with SYK-6, the Cα hydroxy groups of threo-DGPD (DGPD I and II) were initially oxidized to produce the Cα carbonyl form (DGPD-keto I and II). This initial oxidation step is catalyzed by Cα-dehydrogenases, which belong to the short-chain dehydrogenase/reductase (SDR) family and are involved in the catabolism of β-O-4-type dimers. Analysis of seven candidate genes revealed that NAD + -dependent LigD and LigL are mainly involved in the conversion of DGPD I and II, respectively. Next, we found that DGPD-keto I and II were reduced to erythro-DGPD (DGPD III and IV) in the presence of NADPH. Genes involved in this reduction were sought from Cα-dehydrogenase and ldpA-neighboring SDR genes. The gene products of SLG_12690 (ldpC) and SLG_12640 (ldpB) catalyzed the NADPH-dependent conversion of DGPD-keto I to DGPD III and DGPD-keto II to DGPD IV, respectively. Mutational analysis further indicated that ldpC and ldpB are predominantly involved in the reduction of DGPD-keto. Together, these results demonstrate that SYK-6 harbors a comprehensive catabolic enzyme system to utilize all four β-1-type stereoisomers through successive oxidation and reduction reactions of the Cα hydroxy group of threo-DGPD with a net stereoinversion using multiple dehydrogenases.

09 BIOMASS FUELS↗

Topology driven and soft phonon mode enabled Na-ion diffusion in quaternary chalcogenides, Na 3 ZnGaX 4 (X = S, and Se)

The compounds Na 3 ZnGaX 4 (X = S, Se) are potential solid electrolyte materials in sodium-based batteries, which have certain advantages over oxide materials and have shown significant ionic conductivity at ambient temperature. In this paper, we bring out atomic-level features of the diffusion process in these new materials using the microscopic techniques of inelastic neutron scattering (INS), quasielastic neutron scattering (QENS), and ab initio molecular dynamics (AIMD) simulations. The insights obtained from these techniques are unique and not available from other macroscopic experiments. Neutron scattering experiments have been performed at temperatures from 100 to 700 K. The simulations have been carried out up to 900 K. We have calculated the phonon spectra and the space–time correlation functions and found good agreement with the results of the neutron scattering experiments. The simulations enable detailed analysis of the atomic-site dependent dynamical information. We observe low-energy phonon modes of ~6 meV involving the vibrations of certain Na atoms in the lattice. This reveals that the Na at 32g Wyckoff sites (Na2) has sufficiently shallow potential among the two available crystallographic sites. This shallow potential facilitates diffusion. Furthermore, the specific structural topology of the network of interconnected zig-zag chains of the Na2 atomic sites provides the low-barrier energy pathways for diffusion. A small fraction of vacancy defects appears essential for diffusion. We further observe that the Na2 atoms undergo jump-like diffusion to the vacant next or the 2nd next neighbour sites at ~4 Å. While the QENS experiments reveal the jump-like diffusion and its time scale, detailed analysis of the AIMD simulations shows that the jumps appear mostly along zig-zag chains of the Na2 sites in the tetragonal ab-plane, as well as between the chains along the c-axis. In conclusion, the net diffusion is essentially 3-dimensional, with little anisotropy despite the anisotropy of the tetragonal crystal structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigations of correlation and coherence in turbulence from a large-eddy simulation

Abstract. Microscale flow descriptions are often given in terms of mean quantities, turbulent kinetic energy, and/or stresses. Those metrics, while valuable, give limited information about turbulent eddies and coherent turbulent structures. This work investigates the structure of an atmospheric boundary layer using coherence and correlation in space and time with a range of separation distances. We calculate spatial correlations over entire planes of velocity fluctuations, from which we can evaluate the correlation along different directions at different spacings. Similarly, coherence of the three velocity components over separations in the three directions is also investigated. We apply these analyses to a mesoscale–microscale coupled scenario with time-varying conditions and examine nuances in spatial correlations that are often overlooked. Through these analyses and results, this work highlights important differences observed in terms of coherence when comparing large-eddy simulation data to simpler models and suggests ways to improve these simpler models. We note that such differences are important for disciplines like wind energy structural dynamic analysis, in which blade loading and fatigue depend strongly on the structure of the turbulence. We emphasize the additional wealth of data that can be provided by typical atmospheric boundary layer large-eddy simulation when correlation and coherence analysis is included, and we also state the limitations of large-eddy simulation data, which inherently truncate the smaller scales of turbulence.

17 WIND ENERGY↗

Simulation of cryogenic turbopump annular seals

San Andres employed the NBS software package MIPROPS to account for density's dependence on pressure in the simulation of liquid annular seals. His example on a LH2 seal showed a significant change in the mass coefficient compared to a constant density model. San Andres, Yang, and Childs extended this analysis by including the pressure and temperature dependence of density, specific heat, viscosity, volumetric expansion, and thermal conductivity in a coupled solution of the energy, momentum, and continuity equations. Their example showed very significant changes in stiffness and inertia for a high speed (38,000 rpm), large L/D ratio (0.5) LOX seal, as compared to their constant temperature results. The current research rederived the San Andres-Yang-Childs (SYC) analysis and extended it to include not only the Moody friction model of SYC but also the Hir's friction model. The derivation begins with obtaining the local differential equations of continuity, momentum, and energy conservation in the seal. These equations are averaged across the film thickness to obtain the resulting 'bulk flow' differential equations. Shear stress and convective heat loss through the stator (seal) and rotor are related to the Moody and Hir's friction factor model. The Holman analogy is employed to relate heat conduction in or out of the fluid film's boundary layer to the friction induced shear stress.

Palazzolo, Alan B.↗