Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic processes”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2021

The Urban Observatory: A Multi-Modal Imaging Platform for the Study of Dynamics in Complex Urban Systems

We describe an “Urban Observatory” facility designed for the study of complex urban systems via persistent, synoptic, and granular imaging of dynamical processes in cities. An initial deployment of the facility has been demonstrated in New York City and consists of a suite of imaging systems—both broadband and hyperspectral—sensitive to wavelengths from the visible (∼400 nm) to the infrared (∼13 micron) operating at cadences of ∼0.01–30 Hz (characteristically ∼0.1 Hz). Much like an astronomical survey, the facility generates a large imaging catalog from which we have extracted observables (e.g., time-dependent brightnesses, spectra, temperatures, chemical species, etc.), collecting them in a parallel source catalog. We have demonstrated that, in addition to the urban science of cities as systems, these data are applicable to a myriad of domain-specific scientific inquiries related to urban functioning including energy consumption and end use, environmental impacts of cities, and patterns of life and public health. We show that an Urban Observatory facility of this type has the potential to improve both a city’s operations and the quality of life of its inhabitants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

DOE NETL Gas Hydrate Project FWP 65213: Kinetic Parameters for the Exchange of Hydrate Formers (Final Report)

This report documents the research activities sponsored by the U.S. Department of Energy, National Energy Technology Laboratory (NETL), by the Pacific Northwest National Laboratory (PNNL) for the gas hydrate project entitled, “Kinetic Parameters for the Exchange of Hydrate Formers.” This project started June 1, 2013 and concluded December 31, 2018, and comprised five budget periods, under FWP-65213. The overarching goal of this project was to gain an improved understanding of the dynamic processes of gas hydrate accumulations in geologic media by combining laboratory studies, numerical simulation, and analysis of shipboard infrared imaging of hydrate core samples. This project comprised four principal components: 1) fundamental laboratory investigations, 2) numerical simulator development and verification, 3) hydrate core characterization and analysis, and 4) applied laboratory and numerical investigations. Laboratory and field experiments have demonstrated the four core methodologies for producing natural gas hydrates from geologic formations: 1) depressurization, 2) thermal stimulation, 3) inhibitor injection, and 4) guest molecule exchange. Converting these methodologies into viable commercial technologies remains challenging, either from an economic or technical perspective. This project was largely focused with the technical challenges associated with the guest molecule exchange method, which is often referred to as gas swapping. Previous research by PNNL had shown that gas swapping with pure CO 2 was not feasible, as first originally thought, principally due to the formation of secondary hydrates, which clogged or reduced the effective permeability of the hydrate-bearing formation. This research was directed at discovering the potential for using mixtures of gases to affect the gas swapping methodology at commercial scales. The Ignik Sikumi Gas Hydrate Field Trial demonstrated the injection of a N 2 -CO 2 gas mixture into a natural gas hydrate bearing layer, and exchange behavior with the CH 4 hydrate former or guest molecules. The original concept behind the use of a N 2 -CO 2 was to maintain gas conditions at the reservoir state. Pure CO 2 would have existed under supercritical conditions, with an associated density closer to a liquid than gas. One of the key outcomes from this research project has been the discovery, through laboratory experiments and numerical simulation, that gas mixtures, especially those with components having different gas hydrate equilibria conditions provide a control option for the guest molecule exchange methodology. For example, the dissociation or formation of hydrate could be controlled by adjusting the composition of the gas mixture. This research project was dynamic in that the discoveries during each budget period would set the course of research during successive budget periods. Therefore, this report is arranged in chronological order progressing through the five budget periods of the project.

03 NATURAL GAS↗
DOE OSTI2018

Insights on the Alumina–Water Interface Structure by Direct Comparison of Density Functional Simulations with X-ray Reflectivity

Density functional theory molecular dynamics (DFT-MD) simulations are frequently used to predict the interfacial structures and dynamical processes at solid-water interfaces in efforts to gain a deeper understanding of these systems. However, the accuracy of these predictions of interfacial structure has not been rigorously quantified. Here, direct comparisons between large-scale DFT-MD simulations and high-resolution X-ray reflectivity (XR) measurements of the well-defined Al2O3(001)-water interface reveal the relative accuracy of these two methods to describe interfacial structure, a comparison that is enabled by XR’s high sensitivity to atomic-scale displacements. The DFT-MD simulated and XR model-fit structures are qualitatively similar, but XR signals calculated directly from the DFT-MD predictions deviate significantly from the experimental data, revealing discrepancies in these two approaches. Differences in the derived interfacial Al2O3 relaxation profiles of ~0.02 A within the top 5 layers are significant to XR but at the limit of the accuracy of DFT. Further differences are found in the surface hydration layer with a simulated average water layer height ~0.2 A higher than that observed experimentally. This is outside the accuracy of both XR and DFT and is not improved by the inclusion of a phenomenological correction for hydrogen bonding (e.g., Grimme).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Modeling Mesoscale Processes of Scalable Synthesis (Final Report)

The main objective of the center is to develop tools that can help for the understanding of equilibrium and non-equilibrium properties of nanoparticles in aqueous electrolytes. To achieve this goal, we need to develop new computational methodologies that can integrate molecular level description such as the electronic structure analysis to mesoscale description such as generalized hydrodynamics. The challenge is to develop systematic coarse graining techniques that can adaptively choose the level of description needed for the particular task. For example, if bonding breaking is involved, then one should use appropriate models from quantum mechanics to accurately model the time scale and molecular mechanism for the breaking process. However, to model transport process or the process of self-assembly, much less detailed models can be used. There are two main components to this task. One is to develop accurate models that bridge different levels of description. The classical example is the QM/MM model that bridges quantum mechanics and molecular mechanics for studying chemical reactions of large molecules. This work was awarded the Nobel prize this year. For our purpose, we also need to develop methods that bridge molecular and continuum levels of description. Another main task is to develop systematic coarse-graining methods that can give rise to accurate mesoscale and macroscale models, starting from more detailed models at the atomic or molecular scale. None of these ideas are new. By now they have been fashionable topics for quite some time. But neither problem is satisfactorily resolved. The main problem is that uncontrolled approximations are used. As a result, we have very little idea about the validity and accuracy of these models. The main purpose of the center is to exploit the Mori-Zwanzig formalism to achieve the goals stated above. The Mori-Zwanzig formalism is a general tool for eliminating degrees of freedom from existing models. What is attractive about this formalism is its generality. The flip side is that it is also somewhat of a tautology. In order to make it a practical tool, one has to make approximations and this is where the trouble comes. Instead of focusing on Mori-Zwanzig formalism, the PI proposes to study the closely related renormalization group formalism (RG). Like Mori-Zwanzig, RG is a general strategy for eliminating small scale degrees of freedom in a model. Unlike Mori-Zwanzig, RG has the advantage that at least for some problems, one can write a PDE (partial differential equation) type of model for the exact RG flow, without any approximation. This gives us a solid starting point for discussing and analyzing approximations. The PI studied approximate RG flows for turbulent diffusion and related problems. However, extending it to more practical problems proved to be hard. The second line that the PI has pursued is to develop special numerical algorithms for studying particularly important events for dynamics processes in meso-scale systems. The PI succeeded in doing so for the problem of understanding phase transition in metallic systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Time-dependent oral absorption models

The plasma concentration-time profiles following oral administration of drugs are often irregular and cannot be interpreted easily with conventional models based on first- or zero-order absorption kinetics and lag time. Six new models were developed using a time-dependent absorption rate coefficient, ka(t), wherein the time dependency was varied to account for the dynamic processes such as changes in fluid absorption or secretion, in absorption surface area, and in motility with time, in the gastrointestinal tract. In the present study, the plasma concentration profiles of propranolol obtained in human subjects following oral dosing were analyzed using the newly derived models based on mass balance and compared with the conventional models. Nonlinear regression analysis indicated that the conventional compartment model including lag time (CLAG model) could not predict the rapid initial increase in plasma concentration after dosing and the predicted Cmax values were much lower than that observed. On the other hand, all models with the time-dependent absorption rate coefficient, ka(t), were superior to the CLAG model in predicting plasma concentration profiles. Based on Akaike's Information Criterion (AIC), the fluid absorption model without lag time (FA model) exhibited the best overall fit to the data. The two-phase model including lag time, TPLAG model was also found to be a good model judging from the values of sum of squares. This model also described the irregular profiles of plasma concentration with time and frequently predicted Cmax values satisfactorily. A comparison of the absorption rate profiles also suggested that the TPLAG model is better at prediction of irregular absorption kinetics than the FA model. In conclusion, the incorporation of a time-dependent absorption rate coefficient ka(t) allows the prediction of nonlinear absorption characteristics in a more reliable manner.

Non-NASA Center↗
DOE OSTI2017

Sea-ice thickness from field measurements in the northwestern Barents Sea: SEA-ICE THICKNESS, NW BARENTS SEA

The Barents Sea is one of the fastest changing regions of the Arctic, and has experienced the strongest decline in winter-time sea-ice area in the Arctic, at 22364% decade21. Sea-ice thickness in the Barents Sea is not well studied. We present two previously unpublished helicopter-borne electromagnetic(HEM) ice thickness measurements from the northwestern Barents Sea acquired in March 2003 and 2014.The HEM data are compared to ice thickness calculated from ice draft measured by ULS deployed between1994 and 1996. These data show that ice thickness varies greatly from year to year; influenced by the thermodynamic and dynamic processes that govern local formation vs long-range advection. In a year with a large inflow of sea-ice from the Arctic Basin, the Barents Sea ice cover is dominated by thick multiyear ice;as was the case in 2003 and 1995. In a year with an ice cover that was mainly grown in situ, the ice will be thin and mechanically unstable; as was the case in 2014. The HEM data allow us to explore the spatial and temporal variability in ice thickness. In 2003 the dominant ice class was more than 2 years old; and modal sea-ice thickness varied regionally from 0.6 to 1.4 m, with the thinner ice being either irst-year ice, or multiyear ice which had come into contact with warm Atlantic water. In 2014 the ice cover was predominantly locally grown ice less than 1 month old (regional modes of 0.5–0.8 m). These two situations represent two extremes of a range of possible ice thickness distributions that can present very different conditions for shipping traffic; or have a different impact on heat transport from ocean to atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Radiation belt seed population and its association with the relativistic electron dynamics: A statistical study: Radiation Belt Seed Population

Using the Van Allen Probes data, we study the radiation belt seed population and it associated with the relativistic electron dynamics during 74 geomagnetic storm events. Based on the flux changes of 1 MeV electrons before and after the storm peak, these storm events are divided into two groups of “non-preconditioned” and “preconditioned”. The statistical study shows that the storm intensity is of significant importance for the distribution of the seed population (336 keV electrons) in the outer radiation belt. However, substorm intensity can also be important to the evolution of the seed population for some geomagnetic storm events. For non-preconditioned storm events, the correlation between the peak fluxes and their L-shell locations of the seed population and relativistic electrons (592 keV, 1.0 MeV, 1.8 MeV, and 2.1 MeV) is consistent with the energy-dependent dynamic processes in the outer radiation belt. For preconditioned storm events, the correlation between the features of the seed population and relativistic electrons is not fully consistent with the energy-dependent processes. It is suggested that the good correlation between the radiation belt seed population and ≤1.0 MeV electrons contributes to the prediction of the evolution of ≤1.0 MeV electrons in the Earth’s outer radiation belt during periods of geomagnetic storms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Low-energy (< 200 eV) electron acceleration by ULF waves in the plasmaspheric boundary layer: Van Allen Probes observation

Here, we report observational evidence of cold plamsmaspheric electron (< 200 eV) acceleration by ultra-low-frequency (ULF) waves in the plasmaspheric boundary layer on 10 September 2015. Strongly enhanced cold electron fluxes in the energy spectrogram were observed along with second harmonic mode waves with a period of about 1 minute which lasted several hours during two consecutive Van Allen Probe B orbits. Cold electron (<200 eV) and energetic proton (10-20 keV) bi-directional pitch angle signatures observed during the event are suggestive of the drift-bounce resonance mechanism. The correlation between enhanced energy fluxes and ULF waves leads to the conclusions that plasmaspheric dynamics is strongly affected by ULF waves. Van Allen Probe A and B, GOES 13, GOES 15 and MMS 1 observations suggest ULF waves in the event were strongest on the dusk-side magnetosphere. Measurements from MMS 1 contain no evidence of an external wave source during the period when ULF waves and injected energetic protons with a bump-on-tail distribution were detected by Van Allen Probe B. This suggests that the observed ULF waves were probably excited by a localized drift-bounce resonant instability, with the free energy supplied by substorm-injected energetic protons. The observations by Van Allen Probe B suggest that energy transfer between particle species in different energy ranges can take place through the action of ULF waves, demonstrating the important role of these waves in the dynamical processes of the inner magnetosphere.

58 GEOSCIENCES↗
DOE OSTI2018

The Challenge of Identifying Controls on Cloud Properties and Precipitation Onset for Cumulus Congestus Sampled During MC3E

Here, the controls on precipitation onset and the transition from shallow cumulus to congestus are explored using a suite of 16 large–eddy simulations based on the 25 May 2011 event from the Midlatitude Continental Convective Clouds Experiment (MC3E). The thermodynamic variables in the model are relaxed at various timescales to observationally constrained temperature and moisture profiles in order to better reproduce the observed behavior of precipitation onset and total precipitation. Three of the simulations stand out as best matching the precipitation observations and also perform well for independent comparisons of cloud fraction, precipitation area fraction, and evolution of cloud top occurrence. All three simulations exhibit a destabilization over time, which leads to a transition to deeper clouds, but the evolution of traditional stability metrics by themselves is not able to explain differences in the simulations. Conditionally sampled cloud properties (in particular, mean cloud buoyancy), however, do elicit differences among the simulations. The inability of environmental profiles alone to discern subtle differences among the simulations and the usefulness of conditionally sampled model quantities argue for hybrid observational/modeling approaches. These combined approaches enable a more complete physical understanding of cloud systems by combining observational sampling of time–varying three–dimensional meteorological quantities and cloud properties, along with detailed representation of cloud microphysical and dynamical processes from numerical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using Diel Solute Signals to Assess Ecohydrological Processing in Lotic Systems

Lotic systems are a prominent ecohydrological interface within which complex, coupled interactions occur between ecological, hydrological, and geochemical processes. This chapter reviews the ecological processes driving diel signals in oxygen as well as major nutrients and other trace elements. It discusses how diel signals can be used to evaluate functioning and interactions within stream ecosystems. The chapter considers potential future research directions and outstanding questions still to be addressed by or about diel solute signals and associated processes. High-resolution sensor-based time series and associated tools and conceptual advances are being integrated into environmental management. Quantitative tools to disentangle diel biogeochemically-driven signals from other sources of variability in continuous time series are growing. Theoretical models of diel process dynamics are helping to spur hypothesis generation and new ways of approaching diel signals in lotic systems.

Kurz, Marie↗
DOE OSTI2015

On the Localized Nature of the Structural Transformations of Li2MnO 3 Following Electrochemical Cycling

Although the Li-excess layered-oxide Li2MnO3 has a high theoretical capacity, structural transformations within the oxide during electrochemical cycling lead to relatively low experimental capacities, hindering its use in practical applications. Here, aberration-corrected scanning transmission electron microscopy/electron energy loss spectroscopy and high-resolution X-ray diffraction are used to characterize the oxide following electrochemical cycling. Microscopy reveals the coexistence of regions with local monoclinic, spinel, and rock-salt symmetries, indicating localized and inhomogeneous structural evolutions. Crystal structure transformations are observed both at the particle surface and in the bulk. At the surface, these transformed regions resemble spinel Mn3O4 or rock-salt MnO, consistent with oxygen loss. In the bulk, the regions resemble defect spinels, such as the layered-spinel LixMn4/3O4, which suggest a partial phase transformation consistent with oxygen retention. Both microscopy and diffraction data of the cycled sample indicate areas of pristine Li2MnO3; the presence of such areas, in close proximity to LixMn4/3O4 areas, suggests that the layered to spinel structure transformation is partially reversible. Spinel, disordered rock salt, and pristine areas are also observed in Li2MnO3 samples intentionally damaged by electron beam irradiation. This observation indicates that the dynamic processes resulting in phase transformations can be studied for a variety of oxide systems by a judicious selection of irradiation conditions.

aberration-corrected scanning transmission electro↗
DOE OSTI2017

Mechanistic Insight into Nanoparticle Surface Adsorption by Solution NMR Spectroscopy in an Aqueous Gel

Abstract Engineering nanoparticle (NP) functions at the molecular level requires a detailed understanding of the dynamic processes occurring at the NP surface. Herein we show that a combination of dark‐state exchange saturation transfer (DEST) and relaxation dispersion (RD) NMR experiments on gel‐stabilized NP samples enables the accurate determination of the kinetics and thermodynamics of adsorption. We used the former approach to describe the interaction of cholic acid (CA) and phenol (PhOH) with ceria NPs with a diameter of approximately 200 nm. Whereas CA formed weak interactions with the NPs, PhOH was tightly bound to the NP surface. Interestingly, we found that the adsorption of PhOH proceeds via an intermediate, weakly bound state in which the small molecule has residual degrees of rotational diffusion. We believe the use of aqueous gels for stabilizing NP samples will increase the applicability of solution NMR methods to the characterization of nanomaterials.

Egner, Timothy K.↗
DOE OSTI2017

Mechanistic Insight into Nanoparticle Surface Adsorption by Solution NMR Spectroscopy in an Aqueous Gel

Engineering nanoparticle (NP) functions at the molecular level requires a detailed understanding of the dynamic processes occurring at the NP surface. Herein we show that a combination of dark-state exchange saturation transfer (DEST) and relaxation dispersion (RD) NMR experiments on gel-stabilized NP samples enables the accurate determination of the kinetics and thermodynamics of adsorption. We used the former approach to describe the interaction of cholic acid (CA) and phenol (PhOH) with ceria NPs with a diameter of approximately 200 nm. Whereas CA formed weak interactions with the NPs, PhOH was tightly bound to the NP surface. Interestingly, we found that the adsorption of PhOH proceeds via an intermediate, weakly bound state in which the small molecule has residual degrees of rotational diffusion. Here we believe the use of aqueous gels for stabilizing NP samples will increase the applicability of solution NMR methods to the characterization of nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Visualizing Assembly Dynamics of All‐Liquid 3D Architectures

Abstract To better exploit all‐liquid 3D architectures, it is essential to understand dynamic processes that occur during printing one liquid in a second immiscible liquid. Here, the interfacial assembly and transition of 5,10,15,20‐tetrakis(4‐sulfonatophenyl) porphyrin (H 6 TPPS) over time provides an opportunity to monitor the interfacial behavior of nanoparticle surfactants (NPSs) during all‐liquid printing. The formation of J‐aggregates of H 4 TPPS 2− at the interface and the interfacial conversion of the J‐aggregates of H 4 TPPS 2− to H‐aggregates of H 2 TPPS 4− is demonstrated by interfacial rheology and in situ atomic force microscopy. Equally important are the chromogenic changes that are characteristic of the state of aggregation, where J‐aggregates are green in color and H‐aggregates are red in color. In all‐liquid 3D printed structures, the conversion in the aggregate state with time is reflected in a spatially varying change in the color, providing a simple, direct means of assessing the aggregation state of the molecules and the mechanical properties of the assemblies, linking a macroscopic observable (color) to mechanical properties.

Gu, Pei‐Yang↗
DOE OSTI2022

Volcanic electrification: recent advances and future perspectives

Abstract The electrification of volcanic plumes has been described intermittently since at least the time of Pliny the Younger and the 79 AD eruption of Vesuvius. Although sometimes disregarded in the past as secondary effects, recent work suggests that the electrical properties of volcanic plumes reveal intrinsic and otherwise inaccessible parameters of explosive eruptions. An increasing number of volcanic lightning studies across the last decade have shown that electrification is ubiquitous in volcanic plumes. Technological advances in engineering and numerical modelling, paired with close observation of recent eruptions and dedicated laboratory studies (shock-tube and current impulse experiments), show that charge generation and electrical activity are related to the physical, chemical, and dynamic processes underpinning the eruption itself. Refining our understanding of volcanic plume electrification will continue advancing the fundamental understanding of eruptive processes to improve volcano monitoring. Realizing this goal, however, requires an interdisciplinary approach at the intersection of volcanology, atmospheric science, atmospheric electricity, and engineering. Our paper summarizes the rapid and steady progress achieved in recent volcanic lightning research and provides a vision for future developments in this growing field.

58 GEOSCIENCES↗
DOE OSTI2018

A peridynamic failure analysis of fiber-reinforced composite laminates using finite element discontinuous Galerkin approximations

Here, this paper presents a discontinuous Galerkin weak form for bond-based peridynamic models to predict the damage of fiber-reinforced composite laminates. To represent the anisotropy of a laminate in a peridynamic model, a lamina is simplified as a transversely isotropic medium under a plane stress condition. The laminated structure is modeled by stacking the surface mesh layers along the thickness direction according to the laminate sequence. To avoid a mesh dependence on either the fiber orientation or the discretization, the spherical harmonic expansion theory is employed to construct a function for the micro-elastic modulus in terms of the bond-fiber angle. The laminate material is decomposed into an isotropic matrix material part and a transversely isotropic material part. The bond stiffness can be evaluated using the engineering material constants, based on the equivalence between the elastic energy density in the peridynamic theory and the elastic energy density in the classic continuum mechanics theory. Benchmark tests are conducted to verify the proposed model. Numerical results illustrate that the convergence of simulations with different horizon sizes and meshes can be achieved. In terms of damage analysis, the proposed model can capture the dynamic process of the complex coupling of the inner-layer and delamination damage modes.

42 ENGINEERING↗