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,153 records · Page 64

Prediction of divertor heat flux width for ITER using BOUT++ transport and turbulence module

Investigation of turbulent transport dynamics in scrape-off-layer (SOL) and divertor heat flux width prediction is performed for ITER. Both BOUT++ transport and BOUT++ turbulence codes are applied to capture the physics on different temporal scales. Simulations start with an ITER 15MA baseline scenario profile generated by CORSICA (Kim et al 2015 Paper ITER_D_R9T8J9 v1.1). In BOUT++ transport code, the plasma parameters ( n i , T i , T e ) and radial electric () profiles are evolved to steady state. Here, the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients are calculated by inverting the plasma profiles inside the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of zonal shear collapse with hydrodynamic electrons

This paper introduces a theory for the collapse of the edge zonal shear layer, as observed at the density limit at low β. We investigate the scaling of the transport and mean profiles with the adiabaticity parameter α, with special emphasizes on fluxes relevant to zonal flow (ZF) generation. We show that the adiabaticity parameter characterizes the strength of production of zonal flows and so determines the state of turbulence. A 1D reduced model that self-consistently describes the spatiotemporal evolution of the mean density $\overline{n}$, the azimuthal flow $\overline{v}_{y}$, and the turbulent potential enstrophy ε=$\langle$$\tilde{n}$–∇ 2 $\tilde{\phi}$) 2 /2$\rangle$—related to fluctuation intensity—is presented. Quasi-linear analysis determines how the particle flux Γ n and vorticity flux Π=–χ y ∇ 2 v y +Π res scale with α, in both hydrodynamic and adiabatic regimes. As the plasma response passes from adiabatic (α>1) to hydrodynamic (α<1), the particle flux Γ n is enhanced and the turbulent viscosity χ y increases. However, the residual flux Π res —which drives the flow—drops with α. As a result, the mean vorticity gradient ∇ 2 $\overline{v}$ y =Π res /χ y —representative of the strength of the shear—also drops. The shear layer then collapses and turbulence is enhanced. The collapse is due to a decrease in ZF production, not an increase in damping. A physical picture for the onset of collapse is presented. The findings of this paper are used to motivate an explanation of the phenomenology of low β density limit evolution. A change from adiabatic (α=k$^{2}_{z}$v$^{2}_{th}$/(||ω||ν ei )>1) to hydrodynamic (α<1) electron dynamics is associated with the density limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NOVEL METHODS FOR EXPERIMENTAL CHARACTERIZATION OF 3D SUPERCONDUCTING LINAC BEAM DYNAMICS

This paper describes an approach to measure initial Twiss parameters in transverse and longitudinal directions at the entrance of a linac with independent short accelerating cavities. For the transverse plane the usual technique of transverse profiles is used, and for the longitudinal direction a recently developed non-intercepting method is applied. The new method is based on a beam position monitor amplitudes analysis. The applicability of the methods are discussed and demonstrated on an example of the Spallation Neutron Source superconducting linac.

Shishlo, Andrei P↗

New insights into the nonadiabatic state population dynamics of model proton-coupled electron transfer reactions from the mixed quantum-classical Liouville approach

In a previous study [F. A. Shakib and G. Hanna, J. Chem. Phys. 141, 044122 (2014)], we investigated a model proton-coupled electron transfer (PCET) reaction via the mixed quantum-classical Liouville (MQCL) approach and found that the trajectories spend the majority of their time on the mean of two coherently coupled adiabatic potential energy surfaces. This suggested a need for mean surface evolution to accurately simulate observables related to ultrafast PCET processes. In this study, we simulate the time-dependent populations of the three lowest adiabatic states in the ET-PT (i.e., electron transfer preceding proton transfer) version of the same PCET model via the MQCL approach and compare them to the exact quantum results and those obtained via the fewest switches surface hopping (FSSH) approach. We find that the MQCL population profiles are in good agreement with the exact quantum results and show a significant improvement over the FSSH results. All of the mean surfaces are shown to play a direct role in the dynamics of the state populations. Interestingly, our results indicate that the population transfer to the second-excited state can be mediated by dynamics on the mean of the ground and second-excited state surfaces, as part of a sequence of nonadiabatic transitions that bypasses the first-excited state surface altogether. This is made possible through nonadiabatic transitions between different mean surfaces, which is the manifestation of coherence transfer in MQCL dynamics. We also investigate the effect of the strength of the coupling between the proton/electron and the solvent coordinate on the state population dynamics. Drastic changes in the population dynamics are observed, which can be understood in terms of the changes in the potential energy surfaces and the nonadiabatic couplings. Finally, we investigate the state population dynamics in the PT-ET (i.e., proton transfer preceding electron transfer) and concerted versions of the model. The PT-ET results confirm the participation of all of the mean surfaces, albeit in different proportions compared to the ET-PT case, while the concerted results indicate that the mean of the ground- and first-excited state surfaces only plays a role, due to the large energy gaps between the ground- and second-excited state surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗

Controls and relationships of soil organic carbon abundance and persistence vary across pedo–climatic regions

One of the largest uncertainties in the terrestrial carbon cycle is the timing and magnitude of soil organic carbon (SOC) response to climate and vegetation change. This uncertainty prevents models from adequately capturing SOC dynamics and challenges the assessment of management and climate change effects on soils. Reducing these uncertainties requires simultaneous investigation of factors controlling the amount (SOC abundance) and duration (SOC persistence) of stored C. We present a global synthesis of SOC and radiocarbon profiles (n Profile = 597) to assess the timescales of SOC storage. We use a combination of statistical and depth–resolved compartment models to explore key factors controlling the relationships between SOC abundance and persistence across pedo–climatic regions and with soil depth. This allows us to better understand (i) how SOC abundance and persistence covary across pedo–climatic regions and (ii) how the depth dependence of SOC dynamics relates to climatic and mineralogical controls on SOC abundance and persistence. We show that SOC abundance and persistence are differently related; the controls on these relationships differ substantially between major pedo–climatic regions and soil depth. For example, large amounts of persistent SOC can reflect climatic constraints on soils (e.g., in tundra/polar regions) or mineral absorption, reflected in slower decomposition and vertical transport rates. In contrast, lower SOC abundance can be found with lower SOC persistence (e.g., in highly weathered tropical soils) or higher SOC persistence (e.g., in drier and less productive regions). We relate variable patterns of SOC abundance and persistence to differences in the processes constraining plant C input, microbial decomposition, vertical C transport and mineral SOC stabilization potential. This process–oriented grouping of SOC abundance and persistence provides a valuable benchmark for global C models, highlighting that pedo–climatic boundary conditions are crucial for predicting the effects of climate change and soil management on future C abundance and persistence.

54 ENVIRONMENTAL SCIENCES↗

Shallow Cumulus Entrainment Dynamics in a Sheared Environment

Abstract Vertical wind shear has long been known to tilt convective towers and reduce thermal ascent rates. The purpose of this study is to better understand the physical mechanisms responsible for reduced ascent rates in shallow convection. In particular, the study focuses on cloud-edge mass flux to assess how shear impacts mass-flux profiles of both the ensemble and individual clouds of various depths. A compositing algorithm is used to distill large-eddy simulation (LES) output to focus on up- and down-shear cloud edges that are not influenced by complex cloud geometry or nearby clouds. A direct entrainment algorithm is used to estimate the mass flux through the cloud surface. We find that the dynamics on the up- and down-shear sides are fundamentally different, with the entrainment of environmental momentum and dilution of buoyancy being primarily responsible for the reduced down-shear ascent rates. Direct estimates of fluid flow through the cloud interface indicate a counter-shear organized flow pattern that entrains on the down-shear side and detrains on the up-shear side, resulting from the subcloud shear being lifted into the cloud layer by the updraft. In spite of organized regions of entrainment and detrainment, the overall net lateral mass flux remains unchanged with respect to the no shear run, with weak detrainment present throughout cloud depth.

Meteorology & Atmospheric Sciences↗

A Data-Driven Democratized Control Architecture for Regional Transmission Operators

As probably the most complicated and critical infrastructure system, U.S. power grids become increasingly vulnerable to extreme events such as cyber-attacks and severe weather, as well as higher DER penetrations and growing information mismatch among system operators, utilities (transmission or generation owners), and end-users. This paper proposes a data-driven democratized control architecture considering two democratization pathways to assist transmission system operators, with a targeted use case of developing online proactive islanding strategies. Detailed discussions on load capability profiling at transmission buses and disaggregation of DER generations are provided and illustrated with real-world utility data. By Combining network and operational constraints, transmission system operators can be equipped with new tools built on top of this architecture, to derive accurate, proactive, and strategic islanding decisions to incorporate the wide range of dynamic portfolios and needs when facing extreme events or unseen grid contingencies.

Power System Reliability, Islanding, Democratized ↗

Dynamic PET simulator via tomographic emission projection for kinetic modeling and parametric image studies

Purpose: To develop and evaluate a fast and simple tool called dPETSTEP (Dynamic PET Simulator of Tracers via Emission Projection), for dynamic PET simulations as an alternative to Monte Carlo (MC), useful for educational purposes and evaluation of the effects of the clinical environment, postprocessing choices, etc., on dynamic and parametric images. Methods: The tool was developed in MATLAB using both new and previously reported modules of PETSTEP (PET Simulator of Tracers via Emission Projection). Time activity curves are generated for each voxel of the input parametric image, whereby effects of imaging system blurring, counting noise, scatters, randoms, and attenuation are simulated for each frame. Each frame is then reconstructed into images according to the user specified method, settings, and corrections. Reconstructed images were compared to MC data, and simple Gaussian noised time activity curves (GAUSS). Results: dPETSTEP was 8000 times faster than MC. Dynamic images from dPETSTEP had a root mean square error that was within 4% on average of that of MC images, whereas the GAUSS images were within 11%. The average bias in dPETSTEP and MC images was the same, while GAUSS differed by 3% points. Noise profiles in dPETSTEP images conformed well to MC images, confirmed visually by scatter plot histograms, and statistically by tumor region of interest histogram comparisons that showed no significant differences (p < 0.01). Compared to GAUSS, dPETSTEP images and noise properties agreed better with MC. Conclusions: The authors have developed a fast and easy one-stop solution for simulations of dynamic PET and parametric images, and demonstrated that it generates both images and subsequent parametric images with very similar noise properties to those of MC images, in a fraction of the time. They believe dPETSTEP to be very useful for generating fast, simple, and realistic results, however since it uses simple scatter and random models it may not be suitable for studies investigating these phenomena. dPETSTEP can be downloaded free of cost from https://github.com/CRossSchmidtlein/dPETSTEP.

60 APPLIED LIFE SCIENCES↗

Revealing the dynamics of ultrarelativistic non-equilibrium many-electron systems with phase space tomography

The description of physical processes with many-particle systems is a key approach to the modeling of numerous physical systems. For example in storage rings, where ultrarelativistic particles are agglomerated in dense bunches, the modeling and measurement of their phase-space distribution is of paramount importance: at any time the phase-space distribution not only determines the complete space-time evolution but also provides fundamental performance characteristics for storage ring operation. Here, we demonstrate a non-destructive tomographic imaging technique for the 2D longitudinal phase-space distribution of ultrarelativistic electron bunches. For this purpose, we utilize a unique setup, which streams turn-by-turn near-field measurements of bunch profiles at MHz repetition rates. To demonstrate the feasibility of our method, we induce a non-equilibrium state and show that the phase-space distribution microstructuring as well as the phase-space distribution dynamics can be observed in great detail. Our approach offers a pathway to control ultrashort bunches and supports, as one example, the development of compact accelerators with low energy footprints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tully-Fisher relation, galactic rotation curves and dissipative mirror dark matter

If dark matter is dissipative then the distribution of dark matter within galactic halos can be governed by dissipation, heating and hydrostatic equilibrium. Previous work has shown that a specific model, in the framework of mirror dark matter, can explain several empirical galactic scaling relations. It is shown here that this dynamical halo model implies a quasi-isothermal dark matter density, ρ(r) ≅ ρ{sub 0}r{sub 0}{sup 2}/(r{sup 2}+r{sub 0}{sup 2}), where the core radius, r{sub 0}, scales with disk scale length, r{sub D}, via r{sub 0}/kpc ≈ 1.4(r{sub D}/kpc). Additionally, the product ρ{sub 0}r{sub 0} is roughly constant, i.e. independent of galaxy size (the constant is set by the parameters of the model). The derived dark matter density profile implies that the galactic rotation velocity satisfies the Tully-Fisher relation, L{sub B}∝v{sup 3}{sub max}, where v{sub max} is the maximal rotational velocity. Examples of rotation curves resulting from this dynamics are given.

79 ASTRONOMY AND ASTROPHYSICS↗

Proteomic profile changes associated with diminished expression of T-cell intracellular antigens reveal a hormesis response

T-cell intracellular antigen (TIA) proteins function as regulators of cell homeostasis by controlling global gene expression in response to dynamic regulatory changes and environmental stress. Here, we used two-dimensional differential in-gel electrophoresis (2D-DIGE) and mass spectrometry (MALDI-TOF/TOF) to identify protein changes associated with the down-regulated expression of TIA proteins. We detected 30 differentially expressed proteins (DEPs), 24 of which were identified, and some of these DEPs were validated by western blotting. In silico analysis showed that DEPs were associated with metabolic processes, detoxification and proteostasis. We mapped the DEPs to the available biological pathways and networks, which included the metabolism of small molecules such as sugars, lipids, amino acids, and nucleotides. Our findings support previous studies and suggest that low expression of TIA proteins might act as a potential adaptive switch to link gene expression reprogramming to a proliferative phenotype mediated by a hormesis phenomenon.

60 APPLIED LIFE SCIENCES↗

Temporal dynamics of the multi-omic response to endurance exercise training

Regular exercise promotes whole-body health and prevents disease, but the underlying molecular mechanisms are incompletely understood. Here, the Molecular Transducers of Physical Activity Consortium profiled the temporal transcriptome, proteome, metabolome, lipidome, phosphoproteome, acetylproteome, ubiquitylproteome, epigenome and immunome in whole blood, plasma and 18 solid tissues in male and female Rattus norvegicus over eight weeks of endurance exercise training. The resulting data compendium encompasses 9,466 assays across 19 tissues, 25 molecular platforms and 4 training time points. Thousands of shared and tissue-specific molecular alterations were identified, with sex differences found in multiple tissues. Temporal multi-omic and multi-tissue analyses revealed expansive biological insights into the adaptive responses to endurance training, including widespread regulation of immune, metabolic, stress response and mitochondrial pathways. Many changes were relevant to human health, including non-alcoholic fatty liver disease, inflammatory bowel disease, cardiovascular health and tissue injury and recovery. The data and analyses presented in this study will serve as valuable resources for understanding and exploring the multi-tissue molecular effects of endurance training and are provided in a public repository (https://motrpac-data.org/).

59 BASIC BIOLOGICAL SCIENCES↗

Poly(N-isopropylacrylamide)-grafted gold nanoparticles at the vapor/water interface

Hypothesis: Grafting nanoparticles surfaces with water-soluble polymers modify interparticle interactions that are pivotal for assembling them into ordered phases. By manipulating salt concentrations of gold nanoparticles (AuNPs) that are grafted with poly(N-isopropylacrylamide) (PNIPAM-AuNPs), we hypothesize that various aggregated phases form at the suspension/vapor interface or in the bulk that depend on the molecular weight (MW) of PNIPAM and on salt concentrations. Experiments: AuNPs are grafted with thiolated PNIPAM of molecular weights of 3 or 6 kDa, and grafting is confirmed by dynamic light scattering. Liquid-surfaces X-ray reflectivity and grazing incidence small-angle X-ray scattering are used to determine the density profiles of the suspension/vapor interface and their inplane structure as salt is added to the suspensions. Findings: We observe that surface enrichment is induced by adding NaCl to the suspensions, and that at low salt concentrations, the monoparticle layer formed is dispersed, and above a threshold salt concentration, depending on MW of PNIPAM, the PNIPAM-AuNPs order in a hexagonal structure. We show that the lattice constant of the two-dimensional hexagonal structure varies with salt concentration, and more significantly with MW of PNIPAM.

36 MATERIALS SCIENCE↗

Digital force prediction for milling

This paper presents a fully digital, integrated approach for milling force prediction with arbitrary end mill-work material combinations. The approach includes: 1) structured light scanning to identify the end mill’s cutting edge macro-geometry along the tool axis; 2) structured light scanning to measure the cutting edge cross-sectional rake and relief profiles; 3) finite element analysis of orthogonal cutting to determine the force model coefficients that relate the force to chip area using the work material’s constitutive model and tool’s rake and relief profiles; and 4) time domain simulation with inputs that include the measured cutting edge macro-geometry, finite element-based force model, and measured structural dynamics. Milling force predictions are compared to in-process measurements to validate the method.

42 ENGINEERING↗