Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local dynamic operator”

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 847 records · Page 47

Passive Islanding Detection of Inverter-Based Resources in a Noisy Environment

Islanding occurs when a load is energized solely by local generators and can result in frequency and voltage instability, changes in current, and poor power quality. Poor power quality can interrupt industrial operations, damage sensitive electrical equipment, and induce outages upon the resynchronization of the island with the grid. This study proposes an islanding detection method employing a Duffing oscillator to analyze voltage fluctuations at the point of common coupling (PCC) under a high-noise environment. Unlike existing methods, which overlook the noise effect, this paper mitigates noise impact on islanding detection. Power system noise in PCC measurements arises from switching transients, harmonics, grounding issues, voltage sags and swells, electromagnetic interference, and power quality issues that affect islanding detection. Transient events like lightning-induced traveling waves to the PCC can also introduce noise levels exceeding the voltage amplitude by more than seven times, thus disturbing conventional detection techniques. The noise interferes with measurements and increases the nondetection zone (NDZ), causing failed or delayed islanding detection. The Duffing oscillator nonlinear dynamics enable detection capabilities at a high noise level. The proposed method is designed to detect the PCC voltage fluctuations based on the IEEE standard 1547 through the Duffing oscillator. For the voltages beyond the threshold, the Duffing oscillator phase trajectory changes from periodic to chaotic mode and sends an islanded operation command to the inverter. The proposed islanding detection method distinguishes switching transients and faults from an islanded operation. Experimental validation of the method is conducted using a 3.6 kW PV setup.

Energy & Fuels↗

Synthesis of historical reservoir operations from 1980 to 2020 for the evaluation of reservoir representation in large-scale hydrologic models

Abstract. All the major river systems in the contiguous United States (CONUS) (and many in the world) are impacted by dams, yet reservoir operations remain difficult to quantify and model due to a lack of data. Reservoir operation data are often inaccessible or distributed across many local operating agencies, making the acquisition and processing of data records quite time-consuming. As a result, large-scale models often rely on simple parameterizations for assumed reservoir operations and have a very limited ability to evaluate how well these approaches match actual historical operations. Here, we use the first national dataset of historical reservoir operations in the CONUS domain, ResOpsUS, to analyze reservoir storage trends and operations in more than 600 major reservoirs across the US. Our results show clear regional differences in reservoir operations. In the eastern US, which is dominated by flood control storage, we see storage peaks in the winter months with sharper decreases in the operational range (i.e., the difference between monthly maximum and minimum storage) in the summer. In the more arid western US where storage is predominantly for irrigation, we find that storage peaks during the spring and summer with increases in the operational range during the summer months. The Lower Colorado region is an outlier because its seasonal storage dynamics more closely mirrored those of flood control basins, yet the region is classified as arid, and most reservoirs have irrigation uses. Consistent with previous studies, we show that average annual reservoir storage has decreased over the past 40 years, although our analyses show a much smaller decrease than previous work. The reservoir operation characterizations presented here can be used directly for development or evaluation of reservoir operations and their derived parameters in large-scale models. We also evaluate how well historical operations match common assumptions that are often applied in large-scale reservoir parameterizations. For example, we find that 100 dams have maximum storage values greater than the reported reservoir capacity from the Global Reservoirs and Dams database (GRanD). Finally, we show that operational ranges have been increasing over time in more arid regions and decreasing in more humid regions, pointing to the need for operating policies which are not solely based on static values.

Geology↗

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computer network environment planning and analysis

The GSFC Computer Network Environment provides a broadband RF cable between campus buildings and ethernet spines in buildings for the interlinking of Local Area Networks (LANs). This system provides terminal and computer linkage among host and user systems thereby providing E-mail services, file exchange capability, and certain distributed computing opportunities. The Environment is designed to be transparent and supports multiple protocols. Networking at Goddard has a short history and has been under coordinated control of a Network Steering Committee for slightly more than two years; network growth has been rapid with more than 1500 nodes currently addressed and greater expansion expected. A new RF cable system with a different topology is being installed during summer 1989; consideration of a fiber optics system for the future will begin soon. Summmer study was directed toward Network Steering Committee operation and planning plus consideration of Center Network Environment analysis and modeling. Biweekly Steering Committee meetings were attended to learn the background of the network and the concerns of those managing it. Suggestions for historical data gathering have been made to support future planning and modeling. Data Systems Dynamic Simulator, a simulation package developed at NASA and maintained at GSFC was studied as a possible modeling tool for the network environment. A modeling concept based on a hierarchical model was hypothesized for further development. Such a model would allow input of newly updated parameters and would provide an estimation of the behavior of the network.

Dalphin, John F.↗

Applications of Data Assimilation to Analysis of the Ocean on Large Scales

It is commonplace to begin talks on this topic by noting that oceanographic data are too scarce and sparse to provide complete initial and boundary conditions for large-scale ocean models. Even considering the availability of remotely-sensed data such as radar altimetry from the TOPEX and ERS-1 satellites, a glance at a map of available subsurface data should convince most observers that this is still the case. Data are still too sparse for comprehensive treatment of interannual to interdecadal climate change through the use of models, since the new data sets have not been around for very long. In view of the dearth of data, we must note that the overall picture is changing rapidly. Recently, there have been a number of large scale ocean analysis and prediction efforts, some of which now run on an operational or at least quasi-operational basis, most notably the model based analyses of the tropical oceans. These programs are modeled on numerical weather prediction. Aside from the success of the global tide models, assimilation of data in the tropics, in support of prediction and analysis of seasonal to interannual climate change, is probably the area of large scale ocean modeling and data assimilation in which the most progress has been made. Climate change is a problem which is particularly suited to advanced data assimilation methods. Linear models are useful, and the linear theory can be exploited. For the most part, the data are sufficiently sparse that implementation of advanced methods is worthwhile. As an example of a large scale data assimilation experiment with a recent extensive data set, we present results of a tropical ocean experiment in which the Kalman filter was used to assimilate three years of altimetric data from Geosat into a coarsely resolved linearized long wave shallow water model. Since nonlinear processes dominate the local dynamic signal outside the tropics, subsurface dynamical quantities cannot be reliably inferred from surface height anomalies. Because of its potential for large scale synoptic coverage of the deep ocean, acoustic travel time data should be a natural complement to satellite altimetry. Satellite data give us vertical integrals associated with thermodynamic and dynamic processes.

Miller, Robert N.↗

Feeding Ten Billion People Is Possible Within Four Terrestrial Planetary Boundaries

Global agriculture puts heavy pressure on planetary boundaries, posing the challenge to achieve future food security without compromising Earth system resilience. On the basis of process-detailed, spatially explicit representation of four interlinked planetary boundaries (biosphere integrity, land-system change, freshwater use, nitrogen flows) and agricultural systems in an internally consistent model framework, we here show that almost half of current global food production depends on planetary boundary transgressions. Hotspot regions, mainly in Asia, even face simultaneous transgression of multiple underlying local boundaries. If these boundaries were strictly respected, the present food system could provide a balanced diet (2,355 kcal per capita per day) for 3.4 billion people only. However, as we also demonstrate, transformation towards more sustainable production and consumption patterns could support 10.2 billion people within the planetary boundaries analysed. Key prerequisites are spatially redistributed cropland, improved water–nutrient management, food waste reduction and dietary changes. Adoption of the Sustainable Development Goals by all nations in 2015 is the first ever commitment to a world development path that safeguards the stability of the Earth system as a prerequisite for meeting universal human standards1. The longstanding challenge of achieving food security through sustainable agriculture is particularly acute in this context as world agriculture is a leading cause for the current transgressions of multiple planetary boundaries (PBs) globally and regionally2–5. The PB framework is a comprehensive scientific attempt to synoptically define our planet’s biogeophysical limits to anthropogenic interference. It suggests bounds to nine interacting processes that together delineate a Holocene-like Earth system state. The Holocene is chosen as the reference state as it is the only period known to provide a safe operating space for a world population of several billion people, and according to a precautionary principle, the PBs are set in sufficient distance from processes that may critically undermine Earth system resilience and global sustainability. A challenging question, thus, is whether human development goals such as food security can be met while maintaining multiple PBs along with their subglobal manifestations. Further PB transgressions could jeopardize the chances of providing sufficient food for a world population projected to be wealthier and reach >9 billion by 2050. This conundrum portrays a tradeoff between Earth’s biophysical carrying capacity and humankind’s rising food demand, calling in response for radical rethinking of food production and consumption patterns6–9. Yield gap closures, avoidance of excessive input use, shifts towards less resource-demanding diets, food waste reductions and efficient international trade are crucial options for sustainably increasing the food supply10–15. For example, enhancing water-use efficiency on irrigated and rain-fed farms can triple or quadruple crop yields in low-performing systems, suggesting possible global gains of >20% (ref. 16). Even higher gains appear feasible through globally optimized configurations of the land-use pattern17, and cutting food losses by half could generate food for another billion people18. Thus, collective large-scale implementation of such options could sustain food for a further growing world population19. Yet achieving this within a safe operating space as defined by PBs requires not only a halt to but actually a reversal of existing PB transgressions. Previous studies suggest that such a reconciliation might be possible, but these were based on aggregate representations of PBs (not accounting for the spatial patterns of limits, transgressions and interactions) or considered only one boundary in isolation17,20–23. Here, we systematically quantify to what extent current food production depends on local to global transgressions of the PBs for biosphere integrity, land-system change, freshwater use and nitrogen (N) flows, along with the potential of a range of solutions to avoid these transgressions and still increase food supply (Table 1). To this end, we configured an internally consistent process-based model of the terrestrial biosphere including agriculture (LPJmL) with multiple spatially distributed PBs and their interactions. LPJmL is among the longest-established and best-evaluated biosphere models, showing robust performance regarding simulation of, for example, carbon, water and crop yield dynamics (Supplementary Figs. 1 and 2 and Supplementary Table 1; see ref. 24 for a comprehensive benchmarking and Supplementary Methods for more detail on model evaluations). In principle following established definitions4, we refine the computation of some PBs with respect to their regional patterns and interactions (Methods), providing globally gridded precautionary limits to human interference with the Earth system at a level of great detail. In particular, we account for the evidence that many PBs need to be represented spatially explicitly4 to cover their

Gerten, Dieter↗

Non-Markovian reduced dynamics based upon a hierarchical effective-mode representation

A reduced dynamics representation is introduced which is tailored to a hierarchical, Mori-chain type representation of a bath of harmonic oscillators which are linearly coupled to a subsystem. We consider a spin-boson system where a single effective mode is constructed so as to absorb all system-environment interactions, while the residual bath modes are coupled bilinearly to the primary mode and among each other. Using a cumulant expansion of the memory kernel, correlation functions for the primary mode are obtained, which can be suitably approximated by truncated chains representing the primary-residual mode interactions. A series of reduced-dimensional bath correlation functions is thus obtained, which can be expressed as Fourier-Laplace transforms of spectral densities that are given in truncated continued-fraction form. For a master equation which is second order in the system-bath coupling, the memory kernel is re-expressed in terms of local-in-time equations involving auxiliary densities and auxiliary operators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

Mesoscale weather conditions can significantly affect the space launch and landing operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). During the summer months, land-sea interactions that occur across KSC and CCAFS lead to the formation of a sea breeze, which can then spawn deep convection. These convective processes often last 60 minutes or less and pose a significant challenge to the forecasters at the National Weather Service (NWS) Spaceflight Meteorology Group (SMG). The main challenge is that a "GO" forecast for thunderstorms and precipitation at the Shuttle Landing Facility is required at the 90 minute deorbit decision for End Of Mission (EOM) and at the 30 minute Return To Launch Site (RTLS) decision. Convective initiation, timing, and mode also present a forecast challenge for the NWS in Melbourne, FL (MLB). The NWS MLB issues such tactical forecast information as Terminal Aerodrome Forecasts (TAF5), Spot Forecasts for fire weather and hazardous materials incident support, and severe/hazardous weather Watches, Warnings, and Advisories. Lastly, these forecasting challenges can also affect the 45th Weather Squadron (45 WS), which provides comprehensive weather forecasts for shuttle launch, as well as ground operations, at KSC and CCAFS. The need for accurate mesoscale model forecasts to aid in their decision making is crucial. This study specifically addresses the skill of different model configurations in forecasting warm season convective initiation. Numerous factors influence the development of convection over the Florida peninsula. These factors include sea breezes, river and lake breezes, the prevailing low-level flow, and convergent flow due to convex coastlines that enhance the sea breeze. The interaction of these processes produces the warm season convective patterns seen over the Florida peninsula. However, warm season convection remains one of the most poorly forecast meteorological parameters. To determine which configuration options are best to address this specific forecast concern, the Weather Research and Forecasting (WRF) model, which has two dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM) was employed. In addition to the two dynamical cores, there are also two options for a "hot-start" initialization of the WRF model - the Local Analysis and Prediction System (LAPS; McGinley 1995) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS; Brewster 1996). Both LAPS and ADAS are 3- dimensional weather analysis systems that integrate multiple meteorological data sources into one consistent analysis over the user's domain of interest. This allows mesoscale models to benefit from the addition of highresolution data sources. Having a series of initialization options and WRF cores, as well as many options within each core, provides SMG and MLB with considerable flexibility as well as challenges. It is the goal of this study to assess the different configurations available and to determine which configuration will best predict warm season convective initiation.

Watson, Leela R.↗

Orbital ice: An exact Coulomb phase on the diamond lattice

We demonstrate the existence of an orbital Coulomb phase as the exact ground state of a p-orbital exchange Hamiltonian on the diamond lattice. The Coulomb phase is an emergent state characterized by algebraic dipolar correlations and a gauge structure resulting from local constraints (ice rules) of the underlying lattice models. For most ice models on the pyrochlore lattice, these local constraints are a direct consequence of minimizing the energy of each individual tetrahedron. On the contrary, the orbital ice rules are emergent phenomena resulting from the quantum orbital dynamics. We show that the orbital ice model exhibits an emergent geometrical frustration by mapping the degenerate quantum orbital ground states to the spin-ice states obeying the 2-in-2-out constraints on the pyrochlore lattice. We also discuss possible realization of the orbital ice model in optical lattices with p-band fermionic cold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison of Energetic Electron Fluxes Measured by GPS and THEMIS Spacecraft in the Inner Magnetosphere

Key elements of space weather models are energetic electron fluxes in the inner magnetosphere and the outer radiation belt. Flux depletion is driven by various loss processes: scattering into atmosphere, magnetopause shadowing. Flux enhancement is driven by various acceleration processes: local wave-particle interactions, radial transport, plasma sheet injections. Many of these processes operate on ~ hour timescales. Such mesoscale flux variations are not well traced by equatorial spacecraft with much longer orbits. Energetic electron detectors onboard the Global Positioning System (GPS) constellation provide a unique opportunity for probing such ~ hour-scale flux variations. Measurements from up to 23 identically instrumented GPS satellites cover a wide energy and L-shell range with a subhour time resolution. However, their orbits are inclined and thus all measurements at L-shell >4.3 are off-equatorial. Here, in this report, we present a comparison of equatorial THEMIS and nonequatorial GPS measurements of omnidirectional ≤600 keV electron fluxes. Such a comparison allows us to derive coefficients for using off-equatorial GPS fluxes to infer the equatorial values. These coefficients depend on particle energy and L-shell. We demonstrate a new data set derived from GPS measurements and discuss how it can be used to investigate mesoscale dynamics of energetic electron fluxes in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine Learning-Driven Quantification of CO2 Plume Dynamics at Illinois Basin Decatur Project Sites Using Microseismic Data

This study utilizes machine learning to quantify CO2 plume extents by analyzing microseismic data from the Illinois Basin Decatur Project (IBDP). Leveraging a unique dataset of well logs, microseismic records, and CO2 injection metrics, this work aims to predict the temporal evolution of subsurface CO2 saturation plumes. The findings illustrate that machine learning can predict plume dynamics, revealing vertical clustering of microseismic events over distinct time periods within certain proximities to the injection well, consistent with an invasion percolation model. The buoyant CO2 plume partially trapped within sandstone intervals periodically breaches localized barriers or baffles, which act as leaky seals and impede vertical migration until buoyancy overcomes gravity and capillary forces, leading to breakthroughs along vertical zones of weakness. Between different unsupervised clustering techniques, K-Means and DBSCAN were applied and analyzed in detail, where K-means outperformed DBSCAN in this specific study by indicating the combination of the highest Silhouette Score and the lowest Davies–Bouldin Index. The predictive capability of machine learning models in quantifying CO2 saturation plume extension is significant for real-time monitoring and management of CO2 sequestration sites. The models exhibit high accuracy, validated against physical models and injection data from the IBDP, reinforcing the viability of CO2 geological sequestration as a climate change mitigation strategy and enhancing advanced tools for safe management of these operations.

Iyegbekedo, Ikponmwosa↗

Wind direction estimation using SCADA data with consensus-based optimization

Abstract. Wind turbines in a wind farm typically operate individually to maximize their own performance and do not take into account information from nearby turbines. To enable cooperation to achieve farm-level objectives, turbines will need to use information from nearby turbines to optimize performance, ensure resiliency when other sensors fail, and adapt to changing local conditions. A key element of achieving a more efficient wind farm is to develop algorithms that ensure reliable, robust, real-time, and efficient operation of wind turbines in a wind farm using local sensor information that is already being collected, such as supervisory control and data acquisition (SCADA) data, local meteorological stations, and nearby radars/sodars/lidars. This article presents a framework for developing a cooperative wind farm that incorporates information from nearby turbines in real time to better align turbines in a wind farm. SCADA data from multiple turbines can be used to make better estimates of the local inflow conditions at each individual turbine. By incorporating measurements from multiple nearby turbines, a more reliable estimate of the wind direction can be obtained at an individual turbine. The consensus-based approach presented in this paper uses information from nearby turbines to estimate wind direction in an iterative way rather than aggregating all the data in a wind farm at once. Results indicate that this estimate of the wind direction can be used to improve the turbine's knowledge of the wind direction. This estimated wind direction signal has implications for potentially decreasing dynamic yaw misalignment, decreasing the amount of time a turbine spends yawing due to a more reliable input to the yaw controller, increasing resiliency to faulty wind-vane measurements, and increasing the potential for wind farm control strategies such as wake steering.

17 WIND ENERGY↗

PheniX: A New Vision for the Hard X-ray Sky

We are proposing a mission devoted to high energy X-ray astronomy that is based on a focusing telescope operating in the 1-200 keV energy range but optimized for the hard X-ray range. The main scientific topics concern: Physics of compact objects: The proximity of compact objects provides a unique laboratory to study matter and radiation in extreme conditions of temperature and density in strong gravitational environment. The emission of high energy photons from these objects is far from being understood. The unprecedented sensitivity in the high energy domain will allow a precise determination of the non-thermal processes at work in the vicinity of compact objects. The full 1-200 keV energy coverage will be ideal to disentangle the emission processes produced in the spacetime regions most affected by strong-gravity, as well as the physical links: disk-thermal emission-iron line-comptonisation-reflection-non-thermal emission-jets. Neutron stars-magnetic field-cyclotron lines: Time resolved spectroscopy (and polarimetry) at ultra-high sensitivity of AXP, milliseconds pulsars and magnetars will give new tools to study the role of the synchrotron processes at work in these objects. Cyclotron lines-direct measurement of magnetic filed-equation of state constraints-short bursts-giant flares could all be studied with great details. AGN: The large sensitivity improvement will provide detailed spectral properties of the high energy emission of AGN's. This will give a fresh look to the connection between accretion and jet emission and will provide a new understanding of the physical processes at work. Detection of high-redshift active nuclei in this energy range will allow to introduce an evolutionary aspect to high-energy studies of AGN, probing directly the origin of the Cosmic X-ray Background also in the non-thermal range (> 20 keV). Element formation-Supernovae: The energy resolution achievable for this mission (<0.5 keV) and a large high energy effective area are ideally suited for the 44Ti line study (68 and 78 keV). This radioactive nuclei emission will give an estimate of their quantities and speed in their environment. In addition the study of the spatial structure and spectral emission of SNR will advance our knowledge of the dynamics of supernovae explosions, of particles acceleration mechanisms and how the elements are released in the interstellar medium. Instrumental design: The progress of X-ray focusing optics techniques allows a major step in the instrumental design: the collecting area becomes independent of the detection area. This drastically reduces the instrumental background and will open a new era. The optics will be based on depth-graded multi-layer mirrors in a Wolter I configuration. To obtain a significant effective area in the hundred of keV range a focal length in the 40-50 meters range (attainable with a deployable mast) is needed. In addition such a mission could benefit from recent progress made on mirror coating. We propose to cover the 1-200 keV energy range with a single detector, a double-sided Germanium strip detector operating at 80 K. The main features will be: (a) good energy resolution (.150 keV at 5 keV and <.5 keV at 100 keV), (b) 3 dimensional event localization with a low number of electronic chains, (c) background rejection by the 3D localization, (d) polarisation capabilities in the Compton regime.

PHENIX↗

A RADIATION TRANSFER SOLVER FOR ATHENA USING SHORT CHARACTERISTICS

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code that solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-local thermodynamic equilibrium (LTE) effects. The module is based on well known and well tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator splitting: we describe this approach in detail, including a new constraint on the time step for stability due to radiation diffusion modes. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g., LTE gray) problems can be comparable to or less than a single time step of Athena's MHD integrators, and only few times more expensive than that for more general (non-LTE) problems.

79 ASTRONOMY AND ASTROPHYSICS↗

Disorder Dynamics in Battery Nanoparticles During Phase Transitions Revealed by Operando Single-Particle Diffraction

Structural and ion-ordering phase transitions limit the viability of sodium-ion intercalation materials in grid scale battery storage by reducing their lifetime. However, the combination of phenomena in nanoparticulate electrodes creates complex behavior that is difficult to investigate, especially on the single-nanoparticle scale under operating conditions. In this work, operando single-particle X-ray diffraction (oSP-XRD) is used to observe single-particle rotation, interlayer spacing, and layer misorientation in a functional sodium-ion battery. oSP-XRD is applied to Na 2/3 [Ni 1/3 Mn 2/3 ]O 2 , an archetypal P2-type sodium-ion-positive electrode material with the notorious P2-O2 phase transition induced by sodium (de)intercalation. Here, it is found that during sodium extraction, the misorientation of crystalline layers inside individual particles increases before the layers suddenly align just prior to the P2-O2 transition. The increase in the long-range order coincides with an additional voltage plateau signifying a phase transition prior to the P2-O2 transition. To explain the layer alignment, a model for the phase evolution is proposed that includes a transition from localized to correlated Jahn-Teller distortions. The model is anticipated to guide further characterization and engineering of sodium-ion intercalation materials with P2-O2 type transitions. oSP-XRD, therefore, opens a powerful avenue for revealing complex phase behavior in heterogeneous nanoparticulate systems.

25 ENERGY STORAGE↗

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Markov state models from short non-equilibrium simulations—Analysis and correction of estimation bias

To begin, many state-of-the-art methods for the thermodynamic and kinetic characterization of large and complex biomolecular systems by simulation rely on ensemble approaches, where data from large numbers of relatively short trajectories are integrated. In this context, Markov state models (MSMs) are extremely popular because they can be used to compute stationary quantities and long-time kinetics from ensembles of short simulations, provided that these short simulations are in “local equilibrium” within the MSM states. However, over the last 15 years since the inception of MSMs, it has been controversially discussed and not yet been answered how deviations from local equilibrium can be detected, whether these deviations induce a practical bias in MSM estimation, and how to correct for them. In this paper, we address these issues: We systematically analyze the estimation of MSMs from short non-equilibrium simulations, and we provide an expression for the error between unbiased transition probabilities and the expected estimate from many short simulations. We show that the unbiased MSM estimate can be obtained even from relatively short non-equilibrium simulations in the limit of long lag times and good discretization. Further, we exploit observable operator model (OOM) theory to derive an unbiased estimator for the MSM transition matrix that corrects for the effect of starting out of equilibrium, even when short lag times are used. Finally, we show how the OOM framework can be used to estimate the exact eigenvalues or relaxation time scales of the system without estimating an MSM transition matrix, which allows us to practically assess the discretization quality of the MSM. Applications to model systems and molecular dynamics simulation data of alanine dipeptide are included for illustration. The improved MSM estimator is implemented in PyEMMA of version 2.3.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum and classical behavior in interacting bosonic systems

It is understood that in free bosonic theories, the classical field theory accurately describes the full quantum theory when the occupancy numbers of systems are very large. However, the situation is less understood in interacting theories, especially on time scales longer than the dynamical relaxation time. Recently there have been claims that the quantum theory deviates spectacularly from the classical theory on this time scale, even if the occupancy numbers are extremely large. Furthermore, it is claimed that the quantum theory quickly thermalizes while the classical theory does not. The evidence for these claims comes from noticing a spectacular difference in the time evolution of expectation values of quantum operators compared to the classical micro-state evolution. If true, this would have dramatic consequences for many important phenomena, including laboratory studies of interacting BECs, dark matter axions, preheating after inflation, etc. In this work we critically examine these claims. We show that in fact the classical theory can describe the quantum behavior in the high occupancy regime, even when interactions are large. The connection is that the expectation values of quantum operators in a single quantum micro-state are approximated by a corresponding classical ensemble average over many classical micro-states. Furthermore, by the ergodic theorem, a classical ensemble average of local fields with statistical translation invariance is the spatial average of a single micro-state. So the correlation functions of the quantum and classical field theories of a single micro-state approximately agree at high occupancy, even in interacting systems. Furthermore, both quantum and classical field theories can thermalize, when appropriate coarse graining is introduced, with the classical case requiring a cutoff on low occupancy UV modes. We discuss applications of our results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗