Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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,135 records · Page 63

Higher-order Zeno sequences

The quantum Zeno effect typically refers to freezing the dynamics of a quantum system through frequent observations. In general, quantum Zeno dynamics is obtained with an error of order 𝒪⁢(1/𝑁), where 𝑁 is the number of projective measurements performed within a fixed evolution time. In this work, we develop higher-order Zeno sequences that achieve faster convergence to Zeno dynamics, yielding an improved error scaling of 𝒪⁢(1/𝑁 2⁢𝑘 ), where 𝑘 describes the order of the Zeno sequence. This is achieved by relating higher-order Zeno sequences to higher-order Trotter formulas that achieve similar convergence behavior. We leverage this relation to develop higher-order Zeno sequences for different manifestations of the quantum Zeno effect, including frequent projective measurements and unitary kicks. We go on to discuss achieving quantum Zeno dynamics through periodic control fields of high frequency. We explicitly develop control fields that yield a second-order type improvement in the Zeno error scaling and present shorter Zeno sequences. Finally, we discuss the connection to randomized and Uhrig dynamical decoupling to develop more efficient implementations in the weak-coupling regime.

Quantum Zeno dynamics↗

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Grand-scale Atmospheric Imaging Apparatus (GAIA) and Wind LiDAR Multi-scale Measurements in the Atmospheric Surface Layer

Understanding the organization and dynamics of turbulence structures in the atmospheric surface layer (ASL) is important for fundamental and applied research in different fields, including weather prediction, snow settling, particle and pollutant transport, and wind energy. The main challenges associated with probing and modeling turbulence in the ASL are: i) the broad range of turbulent scales associated with the different eddies present in high Reynolds-number boundary layers ranging from the viscous scale ($\mathcal{O}$(mm)) up to large energy-containing structures ($\mathcal{O}$(km)); ii) the non-stationarity of the wind conditions and the variability associated with the daily cycle of the atmospheric stability; iii) the interactions among eddies of different sizes populating different layers of the ASL, which contribute to momentum, energy, and scalar turbulent fluxes. Creative and innovative measurement techniques are required to probe near-surface turbulence by generating spatio-temporally-resolved data in the proximity of the ground and, at the same time, covering the entire ASL height with large enough streamwise extent to characterize the dynamics of larger eddies evolving aloft. To this aim, the U.S. National Science Foundation sponsored the development of the Grand-scale Atmospheric Imaging Apparatus (GAIA) enabling super-large snow particle image velocimetry (SLPIV) in the near-surface region of the ASL. This inaugural version of GAIA provides a comprehensive measuring system by coupling SLPIV and two scanning Doppler LiDARs to probe the ASL at an unprecedented resolution. A field campaign performed in 2021–2022 and its preliminary results are presented herein elucidating new research opportunities enabled by the GAIA measuring system.

17 WIND ENERGY↗

Testing of Helium-Cooled Metal Molybdenum Disk Target

NorthStar Medical Radioisotopes LLC is planning to produce an important medical radioisotope, molybdenum-99 (Mo-99), through photonuclear reaction on molybdenum-100 (Mo-100). In this approach, molybdenum metal will be bombarded with a 40-MeV electron beam. Because enriched Mo-100 is expensive, it is desired to use as much beam power as possible to achieve maximum production yield and minimize target mass. This objective leads to very high beam power density (heat deposition in the target), which sets challenging requirements for cooling. Together with scientists at Los Alamos National Laboratory, a team at Argonne National Laboratory has developed and demonstrated a cooling approach using pressurized helium, which allows for efficient heat removal. One of the main challenges in this approach is the management of the heat load on the target window. The target window separates the high-pressure helium inside the target from the vacuum in the beamline, so it is constantly under stress from differential pressure. Also, the window is cooled only by the helium gas flowing on one side, making the window cooling even more challenging. High heat deposition in the target disks also imposes a strict requirement on performance of the helium cooling system and thickness of the target disks. The target disks are produced from metal powder via a press-and-sinter process. The resulting disks do not possess tensile strength as high as solid molybdenum and might not survive the vibration from the high-velocity helium coolant and high thermal stress from beam heating. An Argonne team of scientists performed a series of tests at Argonne’s Low Energy Accelerator Facility (LEAF) [3-7]. This report describes two series of tests for scale down production target designs that utilize full-scale 29 millimeters diameter, 0.75 mm thick press-and-sintered disks. We performed two thermal tests with different beam parameters and configurations of the disk laminations. We compared the results of the window temperature measurements and cooling system parameters obtained in the experiments with those predicted by analytical calculations and Computation Flow Dynamic (CFD) simulations. Results of the experiments and calculations are presented below.

36 MATERIALS SCIENCE↗

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The GREEN ‘omics of Nutrient Feedbacks to Soil Warming

The GREEN ‘omics of Nutrient Feedbacks in Soil project advanced the DOE Biological and Environmental Research (BER) mission by developing and applying isotope-enabled ’omics tools to understand how soil microbes regulate carbon and nutrient cycling. Guided by the Growth Rate, growth Efficiency, and stoichiometry of Essential Nutrients (GREEN ’omics) framework, the project aimed to build a predictive, systems-level understanding of microbial traits that control ecosystem biogeochemistry. In a collaboration among Northern Arizona University (lead), West Virginia University, Lawrence Livermore National Laboratory, and Pacific Northwest National Laboratory, we combined quantitative stable isotope probing (qSIP), Chip-SIP, NanoSIMS, and genome-resolved metagenomics across long-term experiments in Arctic, boreal, temperate, and tropical ecosystems. The project produced three key outcomes: 1) We showed that community-weighted temperature sensitivities of bacterial growth (Q10) can predict ecosystem-scale soil respiration responses across diverse soils. 2) We provided the first in situ evidence for density-dependent population dynamics in soil bacteria and demonstrated that nutrient additions intensify competition, concentrating carbon use into fewer taxa. 3) We improved and extended isotope-enabled ’omics methods by quantifying qSIP measurement error to guide experimental design and coupling SIP with genome-resolved metagenomics to reveal cross-kingdom interactions among bacteria, fungi, and viruses. Together, these results show that a small number of microbial traits and taxa exert disproportionate control over soil carbon and nutrient cycling, providing critical data and methods to improve representation of microbial processes in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Non-growing Season Plant Nutrient Uptake Controls Arctic Tundra Vegetation Composition: Modeling Archive

This Modeling Archive is in support of a NGEE-Arctic publication: Riley et al. (2021) “Non-growing season plant nutrient uptake controls Arctic tundra vegetation composition under future climate” (Riley et al., 2021). The dataset contains “ecosys” model outputs reported in Riley et al. (2021) that examines the role of non-growing season (NGS) nutrient dynamics on 21st century vegetation composition. The study highlighted the importance of these nutrient dynamics on the emergent plant functional type distributions, focusing on their role in shrub expansion. The study showed that ignoring NGS nutrient dynamics led to a dramatic under-prediction of shrub expansion by year 2100, motivating the need to include these processes in Earth System Models. Included are modeled net primary production (NPP) and nitrogen (N) uptake for two scenarios: (1) baseline and (2) no NGS N uptake. The output is provided (1) at 25 km resolution across the North American tundra; (2) for years 2012 through 2100; and (3) PFT-specific. Dataset contains 20 NetCDF files in a zip file and one pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Plant physiology, shrub size, thaw depth and soil water content, Seward Peninsula, Alaska, 2023

Stomatal response to changing irradiance measurements (m curves), CO2 response (ACi) and light response (AQ) curves on common tussock tundra shrub, graminoid and forb species. Data were collected in July 2023 from three sites on the Seward Peninsula, Kougarok Mile 64, Teller Mile 27 and Council. Measurements were made with LI-COR LI-6800 gas exchange systems on 18 species. For samples from Teller and Kougarok, thaw depth, soil moisture content and dGPS locations are also provided for most samples. Canopy height and canopy diameter are provided for shrubs measured at the Teller site. The data package files include data files, metadata files and the complete instrument output for all gas exchange measurements in .csv format, and a pdf describing the experimental protocols. See the related data package NGA508 for foliar trait data (leaf mass per area, leaf nitrogen concentration).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Sun-induced chlorophyll fluorescence/eddy covariance measurements at NGEE Arctic Council site, Seward Peninsula, Alaska, 2022

Sun-induced chlorophyll fluorescence (SIF) is emitted from the chlorophyll antenna complexes during the light harvesting and photophysical reactions of photosynthesis. It is directly coupled with and therefore can be used to estimate photosynthetic electron transport and CO2 assimilation. SIF, when observed at the ecosystem level, provides information about ecosystem functions that complements the net ecosystem exchange (NEE) of CO2, CH4, water vapor, and sensible heat. NEE observations have been conducted with the eddy covariance (EC) approach at the NGEE Arctic Council, Alaska site for years. However, concurrent SIF observations have not been available. In June 2022, we installed the fluorescence automated measurement equipment (FAME) at the NGEE Arctic Council site. The zip file contains 3,187 *.dat files that can be opened with a text editor or in *.csv.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Feasibility of Hybrid Electro- and Acousto-Dynamical Systems for Quantum Optical Networks

We present results on an exploratory research program that aims to study the feasibility of novel, sustained, dense, and efficient quantum information generation, storage, retrieval, relay and distribution devices and systems towards future advanced quantum networking architectures. The pathfinder program is anchored on original ideas and device conceptual designs that turn the phononic-sourced challenges into an opportunity towards optimized transduction chains over a broad range of length and energy/temperature scales. The proposal involves conceptual design and feasibility towards integration, commissioning and benchmarking of challenging transduction devices in a scaled quantum network setup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Plant Physiology, Alder and Willow Species, Seward Peninsula, Alaska, 2019

Photosynthetic parameters (Vcmax and Jmax), CO2 response (ACi) curves and dark respiration (Rdark) measurements of tundra shrub species. Data were collected at the NGEE-Arctic Kougarok and Teller study sites on the Seward Peninsula, AK, USA in July 2019. Four species were measured using LI-COR LI-6400 portable gas exchange systems: Alnus viridis ssp. fruticosa, Salix richardsonii and Salix pulchra. The data package files include data files, metadata files, dGPS locations and the complete instrument output for all gas exchange measurements in .csv format. See the related data packages for foliar trait data (leaf mass per area, leaf nitrogen concentration), and leaf spectral reflectance data collected simultaneously on the same shrubs (NGA210, NGA212).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

FOCAL Campaign I: Advanced Wind Turbine Control Strategies

Campaign I of the Floating Offshore-wind Controls Advanced Laboratory Experimental Program (FOCAL) aims to generate a dataset enabling the validation of aerodynamic performance of a scaled turbine mounted on a rigid tower in a fixed condition. The turbine considered in the FOCAL testing campaigns is the IEA-Wind 15MW Reference Wind Turbine. This scaled model is capable of simulating advanced blade-pitch control strategies in a high-quality wind field. The turbine is fully instrumented to record a variety of parameters in real time such as structural loads and dynamics. The test data considered was generated at the University of Maine's Harold Alfond Wind and Wave (W2) testing facility. The Load Cases (LC) considered in this testing campaign are as follows: LC 1.X - Constant wind and blade-pitch with varying rotor speeds LC 2.X - Constant wind and rotor speed with varying blade-pitch LC 3.X - Varying wind with active closed-loop control Detailed properties on the modeled system are found in the following reference: Lenfest E., Floating Offshore-wind Controls Advanced Laboratory (FOCAL) Experimental Program - Campaign I: 1:70 Model-scale Testing of the IEA-Wind 15MW Reference Turbine. UMaine ASCC Report Number 23-40-1183. Details on the results of the verification and validation are found in the following reference: Mendoza, Nicole et al., "Verification and Validation of Model-Scale Turbine Performance and Control for the IEA Wind 15 MW Reference Wind Turbine," Energies, vol. 15, no. 20, 2022, https://doi.org/10.3390/en15207649.

17 WIND ENERGY↗

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗

Validation of finite element analysis strategy to investigate acoustic levitation in a two-axis acoustic levitator

A two-axis acoustic levitator can be used to generate a standing pressure wave capable of levitating solid and liquid particles at appropriate input conditions. This work proposes a simulation framework to investigate the two-axis levitation particle stability using a commercial, computational fluid dynamics software based on the harmonic solution to the acoustic wave equation. The simulation produced predictions of the standing wave that include a strong "+" shaped pattern of nodes and anti-nodes that are aligned with the levitator axes. To verify the simulation, a levitator was built and used to generate the standing wave. The field was probed with a microphone and a motorized-scanning system. After scaling the simulated pressure to the measured pressure, the magnitudes of the sound pressure level at corresponding high-pressure locations were different by no more than 5%. This is the first time a measurement of a two-axis levitator standing pressure wave has been presented and shown to verify simulations. Furthermore, as an additional verification, the authors consulted high speed camera measurements of a reference-levitator transducer, which was found to have a maximum peak-to-peak displacement of 50 +/- 5 mu m. The reference-levitator is known to levitate water at 160 dB. The system for this work was simulated to match the operation of the reference-levitator so that it produced sound pressure levels of 160 dB. This pressure was achieved when the transducer maximum peak-to-peak displacement was 50.8 mu m. The agreement between the two levitators' displacements provides good justification that the modeling approach presented here produces reliable results.

2 axis levitator↗

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗

Characterizing Output Bottlenecks of a Production Supercomputer: Analysis and Implications

This article studies the I/O write behaviors of the Titan supercomputer and its Lustre parallel file stores under production load. The results can inform the design, deployment, and configuration of file systems along with the design of I/O software in the application, operating system, and adaptive I/O libraries.We propose a statistical benchmarking methodology to measure write performance across I/O configurations, hardware settings, and system conditions. Moreover, we introduce two relative measures to quantify the write-performance behaviors of hardware components under production load. In addition to designing experiments and benchmarking on Titan, we verify the experimental results on one real application and one real application I/O kernel, XGC and HACC IO, respectively. These two are representative and widely used to address the typical I/O behaviors of applications.In summary, we find that Titan’s I/O system is variable across the machine at fine time scales. This variability has two major implications. First, stragglers lessen the benefit of coupled I/O parallelism (striping). Peak median output bandwidths are obtained with parallel writes to many independent files, with no striping or write sharing of files across clients (compute nodes). I/O parallelism is most effective when the application—or its I/O libraries—distributes the I/O load so that each target stores files for multiple clients and each client writes files on multiple targets in a balanced way with minimal contention. Second, our results suggest that the potential benefit of dynamic adaptation is limited. In particular, it is not fruitful to attempt to identify “good locations” in the machine or in the file system: component performance is driven by transient load conditions and past performance is not a useful predictor of future performance. For example, we do not observe diurnal load patterns that are predictable.

97 MATHEMATICS AND COMPUTING↗

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗

Extreme-scale stochastic optimization and simulation via learning-enhanced decomposition and parallelization (Final Technical Report)

Stochastic optimization and simulation models ubiquitously arise in designing and operating complex service/engineering systems. They can be extreme in scale due to high-dimensional data and decisions, and can also involve decisions made sequentially in response to newly revealed data, both causing significant computational challenge. The objective of this research is to explore a unified framework that integrates machine learning with discrete optimization and risk-averse modeling, to improve the efficiency of decomposition paradigms for stochastic optimization and simulations at extreme scale. The models we consider represent a broad class of complex decision-making problems, where 0-1 or continuous decisions are made before and/or after knowing multiple sources of uncertainties that could be correlated. We will employ machine learning methods to dynamically decide and prioritize computational procedures, including cut generation, branching, and bounding of the optimal objective. Furthermore, the research will shed new lights on the traditional decomposition algorithms for extreme-scale computing. Deliverables of the research include new modeling and computational methods for advancing the state-of-the-art research in optimization and simulation, bringing many relevant risk-averse, data-driven optimization problems in practice within the range of tractability. Examples include distributed computing server scheduling and sensor deployment for monitoring critical infrastructures. Success in this effort will enable progress in solving multiple extreme-scale problems in the complex system design and operations arising from DoE missions in energy, environment, and national security.

24 POWER TRANSMISSION AND DISTRIBUTION↗