Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic profiling”

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,279 records · Page 71

SU-E-T-358: Monte Carlo Dose Calculation of Small Field Electron Beams

Purpose: Dynamic radiotherapy involving electron beams such as Dynamic Electron Arc Radiotherapy (DEAR) requires accurate dose modelling of small field sizes, similar to the requirement of IMRT field on the small photon field. The current commercial electron Monte Carlo algorithms such as eMC v11 in Eclipse were developed for standard field sizes and do not support the planning of dynamic therapy yet. The purpose of this study is to develop a method to accurately model small field electron beam dosimetry using Monte Carlo simulations. Methods: Comparison between eMC, phantom measurements (diode), and Monte Carlo (MC) simulations (BEAMnrc/DOSYZnrc) were performed for a Varian TrueBeam linac. MC simulations utilized Varian TrueBeam phase space files which had been validated in another study. Static single small field was assessed by comparing dose distributions in water for a 16 MeV beam for circular (2 cm diameter) and rectangular (1×10 cm{sup 2}) cut-out. MC was performed with a resolution of 2.5×2.5×2 mm{sup 2} and statistical uncertainty < 4%. The dose distribution was averaged over adjacent bins to improve precision. Depth dose and orthogonal profiles were evaluated. Results: Small field PDDs differ from those with standard cones. For both circular and rectangular cutouts, the difference in range R8 0-R1 0 is less than 2 mm and in dose within 2%. For the orthogonal profiles, field size and penumbra differences were within 1 mm at depth of maximum dose. The eMC displayed a distinctive “step” in the out-field dose profile in disagreement with both measurement and MC results and needs further investigation. Conclusion: MC was able to characterize the small field dosimetry with good agreement with the measurement data, and thus offers the opportunity for treatment planning of dynamic radiotherapy. Analyses for all other electron energies and cut-out sizes are under way and results will be included in the presentation.

07 ISOTOPE AND RADIATION SOURCES↗

Design of a novel carbon/carbon composite microvascular solar receiver

Solar thermal power tower systems are the primary technology being proposed for solar electricity from thermal energy. Operational limits on these towers are often driven by mechanical properties under significant thermal loads, particularly at the receiver where incoming flux is converted to thermal energy. While the solar receiver’s efficiency is largely driven by its optical properties, thermomechanical stresses on the receiver limits the operational envelope. One pathway to higher efficiency is greater allowable solar fluxes on the receiver but novel materials are required. The present study uses computational fluid dynamics to describe a parametric design space for a microvascular carbon/carbon composite solar receiver as a new material option for high flux solar receivers. Simulations are conducted for different microvascular geometries considering the role of material properties and heat transfer fluids, for the impact on thermal efficiency, and allowable strain. Results show that microscale receiver modules made of the proposed carbon/carbon composite could achieve thermal efficiencies over 90% and full-scale receivers can achieve up to 85% thermal efficiency for the design explored considering realistic strain limits, flux levels, and material properties. These values are highly dependent on the heat transfer fluid pairing, the through plane thermal conductivity of the carbon/carbon composite, the path architecture of the microscale receiver, and the incident solar flux profiles.

14 SOLAR ENERGY↗

Application of dynamic metabolomics to examine in vivo skeletal muscle glucose metabolism in the chronically high-fat fed mouse

Rationale: Defects in muscle glucose metabolism are linked to type 2 diabetes. Mechanistic studies examining these defects rely on the use of high fat-fed rodent models and typically involve the determination of muscle glucose uptake under insulin-stimulated conditions. While insightful, they do not necessarily reflect the physiology of the postprandial state. In addition, most studies do not examine aspects of glucose metabolism beyond the uptake process. Here we present an approach to study rodent muscle glucose and intermediary metabolism under the dynamic and physiologically relevant setting of the oral glucose tolerance test (OGTT). Methods and results: In vivo muscle glucose and intermediary metabolism was investigated following oral administration of [U-{sup 13}C] glucose. Quadriceps muscles were collected 15 and 60 min after glucose administration and metabolite flux profiling was determined by measuring {sup 13}C mass isotopomers in glycolytic and tricarboxylic acid (TCA) cycle intermediates via gas chromatography–mass spectrometry. While no dietary effects were noted in the glycolytic pathway, muscle from mice fed a high fat diet (HFD) exhibited a reduction in labelling in TCA intermediates. Interestingly, this appeared to be independent of alterations in flux through pyruvate dehydrogenase. In addition, our findings suggest that TCA cycle anaplerosis is negligible in muscle during an OGTT. Conclusions: Under the dynamic physiologically relevant conditions of the OGTT, skeletal muscle from HFD fed mice exhibits alterations in glucose metabolism at the level of the TCA cycle. - Highlights: • Dynamic metabolomics was used to investigate muscle glucose metabolism in vivo. • Mitochondrial TCA cycle metabolism is altered in muscle of HFD mice. • This defect was not pyruvate dehydrogenase mediated, as has been previously thought. • Mitochondrial TCA cycle anaplerosis in muscle is virtually absent during the OGTT.

60 APPLIED LIFE SCIENCES↗

Structure and dynamics of cold water super-Earths: the case of occluded CH{sub 4} and its outgassing

In this work, we study the transport of methane in the external water envelopes surrounding water-rich super-Earths. We investigate the influence of methane on the thermodynamics and mechanics of the water mantle. We find that including methane in the water matrix introduces a new phase (filled ice), resulting in hotter planetary interiors. This effect renders the super-ionic and reticulating phases accessible to the lower ice mantle of relatively low-mass planets (∼5 M{sub E} ) lacking a H/He atmosphere. We model the thermal and structural profile of the planetary crust and discuss five possible crustal regimes which depend on the surface temperature and heat flux. We demonstrate that the planetary crust can be conductive throughout or partly confined to the dissociation curve of methane clathrate hydrate. The formation of methane clathrate in the subsurface is shown to inhibit the formation of a subterranean ocean. This effect results in increased stresses on the lithosphere, making modes of ice plate tectonics possible. The dynamic character of the tectonic plates is analyzed and the ability of this tectonic mode to cool the planet is estimated. The icy tectonic plates are found to be faster than those on a silicate super-Earth. A mid-layer of low viscosity is found to exist between the lithosphere and the lower mantle. Its existence results in a large difference between ice mantle overturn timescales and resurfacing timescales. Resurfacing timescales are found to be 1 Ma for fast plates and 100 Ma for sluggish plates, depending on the viscosity profile and ice mass fraction. Melting beneath spreading centers is required in order to account for the planetary radiogenic heating. The melt fraction is quantified for the various tectonic solutions explored, ranging from a few percent for the fast and thin plates to total melting of the upwelled material for the thick and sluggish plates. Ice mantle dynamics is found to be important for assessing the composition of the atmosphere. We propose a mechanism for methane release into the atmosphere, where freshly exposed reservoirs of methane clathrate hydrate at the ridge dissociate under surface conditions. We formulate the relation between the outgassing flux and the tectonic mode dynamical characteristics. We give numerical estimates for the global outgassing rate of methane into the atmosphere. We find, for example, that for a 2 M{sub E} planet outgassing can release 10{sup 27}-10{sup 29} molecules s{sup –1} of methane to the atmosphere. We suggest a qualitative explanation for how the same outgassing mechanism may result in either a stable or a runaway volatile release, depending on the specifics of a given planet. Finally, we integrate the global outgassing rate for a few cases and quantify how the surface atmospheric pressure of methane evolves over time. We find that methane is likely an important constituent of water planets' atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus’ genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

59 BASIC BIOLOGICAL SCIENCES↗

Asynchronous aging and turnover of human circulating and tissue-resident memory T cells across sites

Memory T cells are maintained in tissues as circulating effector-memory (T EM ) and tissue-resident (T RM ) populations for protective immunity, though the role of site and subset in memory persistence remains undefined. Here, in this work, we investigated age-associated dynamics of human T cells in lymphoid organs, mucosal sites, and blood over 10 decades of life using retrospective radiocarbon ( 14 C) birth dating, along with cellular, transcriptome, and epigenetic profiling. Memory T cells across peripheral sites exhibited continuous turnover with mean lifespans of 1–2 years, while the spleen contained longer-lived T cells. Over age, T EM cells expressed senescent markers and a GZMK transcriptional signature, while T RM cells maintained site-specific resident phenotypes without exhibiting features of senescence. Both T EM and T RM cells showed age-associated DNA hypomethylation, though T RM cells exhibited more epigenetically regulated genes. Together, our findings reveal asynchronous aging of human memory T cells by subset and site, as well as persistence of T RM cells without immunosenescence.

T cells↗

Silicon-On-Sapphire Metasurfaces Generate Arrays of Dark and Bright Traps for Neutral Atoms

We demonstrated crystalline silicon-on-sapphire (c-SOS) metasurfaces that convert a Gaussian beam into arrays of complex optical traps, including arrays of optical bottle beams that trap atoms in dark regions interleaved with bright tweezer arrays. The high refractive index and indirect band gap of crystalline silicon make it possible to design high-resolution near-infrared (λ > 700 nm) metasurfaces that can be manufactured at scale using CMOS-compatible processes. Compared with active components like spatial light modulators (SLMs) that have become widely used to generate trap arrays, metasurfaces provide an indefinitely scalable number of pixels, enabling large arrays of complex traps in a very small form factor, as well as reduced dynamic noise. To design metasurfaces that can generate three-dimensional bottle beams to serve as dark traps, we modified the Gerchberg-Saxton algorithm to enforce complex-amplitude profiles at the focal plane of the metasurface and to optimize the uniformity of the traps across the array. We fabricated and measured c-SOS metasurfaces that convert a Gaussian laser beam into arrays of bright traps, dark traps, and interleaved bright/dark traps.

Gerchberg-Saxton↗

Direct comparisons of X-ray scattering and atomistic molecular dynamics simulations for precise acid copolymers and ionomers

Designing acid- and ion-containing polymers for optimal proton, ion, or water transport would benefit profoundly from predictive models or theories that relate polymer structures with ionomer morphologies. Recently, atomistic molecular dynamics (MD) simulations were performed to study the morphologies of precise poly(ethylene-co-acrylic acid) copolymer and ionomer melts. Here, we present the first direct comparisons between scattering profiles, I(q), calculated from these atomistic MD simulations and experimental X-ray data for 11 materials. This set of precise polymers has spacers of exactly 9, 15, or 21 carbons between acid groups and has been partially neutralized with Li, Na, Cs, or Zn. In these polymers, the simulations at 120 °C reveal ionic aggregates with a range of morphologies, from compact, isolated aggregates (type 1) to branched, stringy aggregates (type 2) to branched, stringy aggregates that percolate through the simulation box (type 3). Excellent agreement is found between the simulated and experimental scattering peak positions across all polymer types and aggregate morphologies. The shape of the amorphous halo in the simulated I(q) profile is in excellent agreement with experimental I(q). We found that the modified hard-sphere scattering model fits both the simulation and experimental I(q) data for type 1 aggregate morphologies, and the aggregate sizes and separations are in agreement. Given the stringy structure in types 2 and 3, we develop a scattering model based on cylindrical aggregates. Both the spherical and cylindrical scattering models fit I(q) data from the polymers with type 2 and 3 aggregates equally well, and the extracted aggregate radii and inter- and intra-aggregate spacings are in agreement between simulation and experiment. Furthermore, these dimensions are consistent with real-space analyses of the atomistic MD simulations. By combining simulations and experiments, the ionomer scattering peak can be associated with the average distance between branches of type 2 or 3 aggregates. Furthermore, this direct comparison of X-ray scattering data to the atomistic MD simulations is a substantive step toward providing a comprehensive, predictive model for ionomer morphology, gives substantial support for this atomistic MD model, and provides new credibility to the presence of stringy, branched, and percolated ionic aggregates in precise ionomer melts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-focusing and stimulated Brillouin back-scattering of a long intense laser pulse in a finite temperature relativistic plasma

The nonlinear dynamics of electromagnetic waves propagating through a plasma considering the effects of relativistic mass and ponderomotive nonlinearities is investigated. The modified electron density distribution, the dispersion relation, and the spatial profiles of electromagnetic wave amplitude in the plasma are obtained. It is shown that the cut-off frequency decreases, and there is an intensity range in which the ponderomotive self-focusing takes place. In the upper limit of this range, the laser beam is defocused due to the relativistic ponderomotive force. In addition, the stability of electromagnetic waves to stimulated Brillouin scattering is studied, and the backscattered wave resulting from decay of high power electromagnetic beam is resolved in relativistic regime. The study of effects of electron density and temperature on the growth rate of backscattered wave has been shown that by increasing these effects, the growth rate of instability increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comprehensive comparisons of geodesic acoustic mode characteristics and dynamics between Tore Supra experiments and gyrokinetic simulations

In a dedicated collisionality scan in Tore Supra, the geodesic acoustic mode (GAM) is detected and identified with the Doppler backscattering technique. Observations are compared to the results of a simulation with the gyrokinetic code GYSELA. We found that the GAM frequency in experiments is lower than predicted by simulation and theory. Moreover, the disagreement is higher in the low collisionality scenario. Bursts of non harmonic GAM oscillations have been characterized with filtering techniques, such as the Hilbert-Huang transform. When comparing this dynamical behaviour between experiments and simulation, the probability density function of GAM amplitude and the burst autocorrelation time are found to be remarkably similar. In the simulation, where the radial profile of GAM frequency is continuous, we observed a phenomenon of radial phase mixing of the GAM oscillations, which could influence the burst autocorrelation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of microchannel plate based detectors for <25 keV x rays: Monte Carlo simulations and comparisons with experimental results

We present the results of Monte Carlo simulations of the microchannel plate (MCP) response to x rays in the 250 eV to 25 keV energy range as a function of both x-ray energy and impact angle and their comparisons with the experimental results from the X8A beamline at the National Synchrotron Light Source at Brookhaven National Laboratory. Incoming x rays interact with the lead glass of the microchannel plate, producing photoelectrons. Transport of the photoelectrons is neglected in this model, and it is assumed that photoelectrons deposit all their energy at the point they are created. This deposition leads to the generation of many secondary electrons, some fraction of which diffuse to the MCP pore surface where they can initiate secondary electron cascades in the pore under an external voltage bias. X-ray penetration through multiple MCP pore walls is increasingly important above 5 keV, and the effect of this penetration on MCP performance is studied. In agreement with past measurements, we find that the dependence of MCP sensitivity with angle relative to the pore bias changes from a cotangent dependence to angular independence and then proceeds to a secant dependence as the x-ray energy increases. We also find that with the increasing x-ray energy, the MCP gain sensitivity as a function of bias voltage decreases. The simulations also demonstrate that for x rays incident normal to the MCP surface, spatial resolution shows little dependence on the x-ray energy but degrades with the increasing x-ray energy as the angle of incidence relative to the surface normal increases. This agrees with experimental measurements. Simulation studies have also been completed for MCPs gated with a subnanosecond voltage pulse. In this work, we find that the optical gate profile width increases as the x-ray energy is increased above 5 keV, a consequence of increased x-ray penetration at energies >5 keV. Simulations of the pulsed dynamic range show that the dynamic range varies between ~100 and 1000 depending on x-ray energy and peak voltage.

47 OTHER INSTRUMENTATION↗

Effect of Nozzle Curvature on Supersonic Gas Jets Used in Laser-Plasma Acceleration

Supersonic gas jets produced by converging-diverging (C-D) nozzles are commonly used as targets for laser-plasma acceleration (LPA) experiments. A major point of interest for these targets is the gas density at the region of interaction where the laser ionizes the gas plume to create a plasma, providing the acceleration structure. Tuning the density profiles at this interaction region is crucial to LPA optimization. A "flat-top" density profile is desired at this line of interaction to control laser propagation and high energy electron acceleration, while a short high-density profile is often preferred for acceleration of lower-energy tightly-focused laser-plasma interactions. A particular design parameter of interest is the curvature of the nozzle's diverging section. We examine three nozzle designs with different curvatures: the concave "bell", straight conical and convex "trumpet" nozzles. We demonstrate that, at mm-scale distances from the nozzle exit, the trumpet and straight nozzles, if optimized, produce "flat-top" density profiles whereas the bell nozzle creates focused regions of gas with higher densities. An optimization procedure for the trumpet nozzle is derived and compared to the straight nozzle optimization process. We find that the trumpet nozzle, by providing an extra parameter of control through its curvature, is more versatile for creating flat-top profiles and its optimization procedure is more refined compared to the straight nozzle and the straight nozzle optimization process. Furthermore, we present results for different nozzle designs from computational fluid dynamics (CFD) simulations performed with the program ANSYS Fluent and verify them experimentally using neutral density interferometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental evidence of lower hybrid wave scattering in Alcator C-Mod due to scrape off layer density fluctuations

Novel experimental measurements of the lower hybrid (LH) wave electric field vector, E LH , obtained in Alcator C-Mod using dynamic Stark spectroscopy are presented. Three key results were obtained by comparison with 3D full-wave COMSOL simulations using the cold plasma dielectric tensor and reflectometry-measured density profiles. (1) As a function of space and LH power, the magnitude of E LH was found to agree with the full-wave simulation result. (2) The direction of E LH was found to have a substantial poloidal component and is in strong disagreement with the nearly radial full-wave simulation result. (3) The addition of scrape off layer density fluctuations to the density profile implemented in the full-wave simulations can be used to explain the E LH direction discrepancy. It is concluded that diffraction and scattering by turbulence-driven scrape off layer density fluctuations is the responsible mechanism behind the modification in the direction of E LH . Furthermore, these results further support the previously reported proposition that edge turbulence is an important mechanism behind the strong anomalous drop in LH current drive efficiency for diverted plasmas in the high-density regime.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗

Exploring the Frontiers of Energy Efficiency using Power Management at System Scale

In the face of surging power demands for exascale HPC systems, this work tackles the critical challenge of understanding the impact of software-driven power management techniques like Dynamic Voltage and Frequency Scaling (DVFS) and Power Capping. These techniques have been actively developed over the past few decades. By combining insights from GPU benchmarking to understand application power profiles, we present a telemetry data-driven approach for deriving energy savings projections. This approach has been demonstrably applied to the Frontier supercomputer at scale. Our findings based on three months of telemetry data indicate that, for certain resource-constrained jobs, significant energy savings (up to 8.5%) can be achieved without compromising performance. This translates to a substantial cost reduction, equivalent to 1438 MWh of energy saved. The key contribution of this work lies in the methodology for establishing an upper limit for these best-case scenarios and its successful application. This work enables HPC professionals to optimize the power-performance trade-off within constrained power budgets, not only for the exascale era but also beyond.

Karimi, Ahmad Maroof↗

Kraken2 Metagenomic Virus Database

The Database: Kraken2 [1] database built from a classification tree containing over 700k metagenomic viruses from JGI IMG/VR [2]. (1) Wood, D. E., Lu, J., and Langmead, B. (2019). Improved metagenomic analysis with Kraken 2. Genome Biol., 20(1), 1–13. doi: 10.1186/s13059-019-1891-0 (2) Paez-Espino D, Chen I-MA, Palaniappan K, Ratner A, Chu K, Szeto E, et al. IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses. Nucleic Acids Res. 2017;45:D457-65. For Paper: Title: A k-mer based approach for virus classification in metatranscriptomic and metagenomic samples identifies viral associations in the Populus phytobiome and autism brains Abstract Background Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an important role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which further limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples. Methods To identify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses. We then integrated the viral classification tree with the NCBI taxonomy for use with ParaKraken, a metagenomic/transcriptomic classifier. Results To illustrate the breadth of our utility for classifying viruses with ParaKraken, we analyzed data from a plant metagenome study identifying the differences between two Populus genotypes in three different compartments and on a human metatranscriptome study identifying the differences between Autism Spectrum Disorder patients and controls in post mortem brain biopsies. In the Populus study, we identified genotype and compartment specific viral signatures, while in the Autism study we identified a significant increased abundance of eight viral sequences in Autism brain biopsies. Conclusion Viruses represent an important aspect of the microbiome. The ability to classify viruses represents the first step in being able to better understand their role in the microbiome. The viral classification method presented here allows for more complete identification of viral sequences for use in identifying associations between viruses and the host and viruses and other microbiome members. Acknowledgements and Funding This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725. This research was also supported by the Plant-Microbe Interfaces Scientific Focus Area in the Genomic Science Program, the Office of Biological and Environmental Research (BER) in the U.S. Department of Energy Office of Science, and by the Department of Energy, Laboratory Directed Research and Development funding (ProjectID 8321), at the Oak Ridge National Laboratory. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the US DOE under contract DE-AC05-00OR22725. This research used resources of the Compute and Data Environment for Science (CADES).

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the hydrodynamic impact on the tool influence function during hemispherical subaperture optical polishing

To fabricate high-precision and accurate optics relative to the optical design surface, a high level of deterministic control of material removal (i.e., the tool influence function, TIF) during subaperture tool polishing is required. In this study, a detailed analysis of the pressure distribution, which is a key component of the TIF, has been performed using finite element analysis to couple together solid mechanics and fluid dynamics. Modeling experimental parameters of recently published work reveals that, when considering tool deformation, which in turn influences the fluid film thickness between the tool and workpiece, the effective pressure profile has a flat-top distribution. This flat-top pressure profile differs from the parabolic pressure distributions predicted by Hertzian mechanics. Furthermore, the shear contribution is shown here to be a key contributor to material removal, inducing the removal at the periphery of the contact edge and even outside the generally accepted contact area. Finally, the simulated fluid velocities provide evidence of mixed-mode contact polishing, supporting recent experimental findings that also suggest that onset of hydroplaning contributions lead to material removal drop-off.

36 MATERIALS SCIENCE↗