Surrogate-based partitioned methods for interface problems
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
As the penetration level of power electronics-interfaced renewables such as photovoltaics (PV) and wind has surged in modern electric grids, new operational risks caused by the dynamics of those inverter-based resources (IBRs) are emerging in parallel. Lessons learned from various grid events include that the impact of IBRs on system-level grid stability will become prominent along with the increase of renewables and that the short-timescale dynamic impacts of IBRs on grid stability are not fully captured by current commercial dynamic simulation tools [1] [2]. For example, IBRs can be controlled to mitigate those destabilizing interactions, but conventional phasor-domain tools (e.g. PSS/E, PSLF) often cannot capture that; likewise, the existing electromagnetic transient (EMT) simulation tools (e.g. PSCAD, EMTP) can simulate detailed IBR controls, but for large power systems with many IBRs, slow simulation speeds severely impede the ability to study dynamic events [3] [4]. Massively paralleling simulations using high-performance computing (HPC) can help address this, especially now that cloud-based HPC capability is widely available, but today s EMT tools are not HPC-compatible, and parallelization of dynamic simulation solvers is not trivial because each region can dynamically affect the others. Thus, dynamic simulation of grids with very large numbers of IBRs potentially poses a barrier to the ongoing energy transition.
Electrorefining is a controlled redox process used to regulate the behavior of ionic species. Through this process, metals can be deposited onto a cathode from an electrolyte solution in a controlled manner. The Mk-IV electrorefiner (Mk-IV ER) at Idaho National Laboratory is an engineering-scale, molten salt-based electrorefining cell that has been used for decades to recover metallic uranium from spent fuel. As a result, highly stable fission product chlorides have accumulated in the electrolyte. This accumulation results in changes to the salt’s properties, such as melting temperature, thermal conductivity, and density, as well as elevated product impurity and fissile materials criticality margin. These factors prompt the need for a salt regeneration process, such as melt-crystallization and species drawdown. This work focuses on providing a conceptual design to regenerate ER salt from used Mk-IV-ER salt in-situ, while minimizing salt waste volumes by concentrating the fission products in a final processed salt heal. We propose using a hot-finger crystallization apparatus design to fractionally crystallize salt in the Mk-IV-ER head space (or baffle space), allowing the collection of solid and liquid fractions. By using a cup-drain design, the used salt will be allowed to slowly solidify on the walls of a stainless-steel cup. The apparatus drain plug will then open to allow the liquid salt phase to drain to a lower cup, effectively separating the liquid phase from the solid phase. Under the hypothesis that the liquid phase salt concentrates the fission products, which is under examination in the accompanying work package, this separation allows the recovered solid salt to be reused while minimizing the high-level salt waste volume of used ER salt.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
We present an iterative generalisation of the quantum subspace expansion algorithm used with a Krylov basis. The iterative construction connects a sequence of subspaces via their lowest energy states. Diagonalising a Hamiltonian in a given Krylov subspace requires the same quantum resources in both the single step and sequential cases. We propose a variance-based criterion for determining a good iterative sequence and provide numerical evidence that these good sequences display improved numerical stability over a single step in the presence of finite sampling noise. Implementing the generalisation requires additional classical processing with a polynomial overhead in the subspace dimension. By exchanging quantum circuit depth for additional measurements the quantum subspace expansion algorithm appears to be an approach suited to near term or early error-corrected quantum hardware. Our work suggests that the numerical instability limiting the accuracy of this approach can be substantially alleviated in a parameter-free way.
Explore the source record for details and available documents.
This data set provides (1) the results of in-situ Sphagnum-peat hourly net ecosystem exchange (NEE) measured using a LICOR 8100 gas exchange system and (2) the component fluxes -- gross primary production (GPP) and ecosystem respiration (Re), derived using empirical regressions. NEE measurements were made from 6 June to 6 November 2014 and 20 March to 10 May 2015. Three 8100 chambers per dominant species (S. magellanicum or S. fallax) were placed in the S1 Bog in relatively open ground where there was no obvious hummock-hollow microtopography. The 8100 chambers were not located in the SPRUCE experimental enclosures. Measured site meteorological data used as model drivers are also included. All measurements were made at the 8.1-ha S1 bog forest site in northern Minnesota, 40 km north of Grand Rapids, in the USDA Forest Service Marcell Experimental Forest (MEF).
Plant growth is often constrained by the limited availability of resources in the microenvironment. Despite the continuous threat of attack from insect herbivores and pathogens, investment in defense represents a lost opportunity to expand photosynthetic capacity in leaves and absorption of nutrients and water by roots. To mitigate the metabolic expenditure on defense, plants have evolved inducible defense strategies. The plant hormone jasmonate (JA) is a key regulator of many inducible defenses. Synthesis of JA in response to perceived danger leads to the deployment of a variety of defensive structures and compounds, along with a potent inhibition of growth. Genetic studies have established an important role for JA in mediating tradeoffs between growth and defense. However, several gaps remain in understanding of how JA signaling inhibits growth, either through direct transcriptional control of JA-response genes or crosstalk with other signaling pathways. Here, we highlight recent progress in uncovering the role of JA in controlling growth-defense balance and its relationship to resource acquisition and allocation. We also discuss tradeoffs in the context of the ability of JA to promote increased leaf mass per area (LMA), which is a key indicator of leaf construction costs and leaf life span.
Recent studies of chondrites have found a tungsten isotopic anomaly between chondrules and matrix. Given the refractory nature of tungsten, this implies that W was carried into the solar nebula by at least two distinct families of pre-solar grains. The observed chondrule/matrix split requires that the distinct families were kept separate during the dust coagulation process, and that the two families of grain interacted with the chondrule formation mechanism differently. We take the co-existence of different families of solids in the same general orbital region at the chondrule-precursor size as given, and explore the requirements for them to have interacted with the chondrule formation process at significantly different rates. We show that this sorting of families of solids into chondrule- and matrix-destined dust had to have been at least as powerful a sorting mechanism as the relative settling of aerodynamically distinct grains at least two scale heights above the midplane. The requirement that the chondrule formation mechanism was correlated in some fashion with a dust-grain sorting mechanism argues strongly for spatially localized chondrule formation mechanisms such as turbulent dissipation in non-thermally ionized disk surface layers, and argues against volume-filling mechanisms such as planetesimal bow shocks.
M87, the active galaxy at the center of the Virgo cluster, is ideal for studying the interaction of a supermassive black hole (SMBH) with a hot, gas-rich environment. A deep Chandra observation of M87 exhibits an approximately circular shock front (13 kpc radius, in projection) driven by the expansion of the central cavity (filled by the SMBH with relativistic radio-emitting plasma) with projected radius ∼1.9 kpc. We combine constraints from X-ray and radio observations of M87 with a shock model to derive the properties of the outburst that created the 13 kpc shock. Principal constraints for the model are (1) the measured Mach number ( M ∼ 1.2), (2) the radius of the 13 kpc shock, and (3) the observed size of the central cavity/bubble (the radio-bright cocoon) that serves as the piston to drive the shock. We find that an outburst of ∼5 × 10{sup 57} erg that began about 12 Myr ago and lasted ∼2 Myr matches all the constraints. In this model, ∼22% of the energy is carried by the shock as it expands. The remaining ∼80% of the outburst energy is available to heat the core gas. More than half the total outburst energy initially goes into the enthalpy of the central bubble, the radio cocoon. As the buoyant bubble rises, much of its energy is transferred to the ambient thermal gas. For an outburst repetition rate of about 12 Myr (the age of the outburst), 80% of the outburst energy is sufficient to balance the radiative cooling.
An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.
Analyses of 15,314 electron velocity distribution functions (VDFs) within ±2 hr of 52 interplanetary (IP) shocks observed by the Wind spacecraft near 1 au are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar VDF, while both the halo and beam/strahl components were best fit to bi-kappa VDF. This is the first statistical study to show that the core electron distribution is better fit to a self-similar VDF than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The ranges of values defined by the lower and upper quartiles for the kappa exponents are κ ec ~ 5.40–10.2 for the core, κ eh ~ 3.58–5.34 for the halo, and κ eb ~ 3.40–5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents is s ec ~ 2.00–2.04, and those of asymmetric bi-self-similar core exponents are p ec ~ 2.20–4.00 for the parallel exponent and q ec ~ 2.00–2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.
Abstract A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low-energy core, tenuous, hot halo, and field-aligned beam/strahl are a statistically significant list of values illustrated with both histograms and tabular lists for reference and baselines in future work. Given the large statistics in this investigation, the beam/strahl fit results in the upstream are now the most comprehensive attempt to parameterize the beam/strahl electron velocity moments in the ambient solar wind. The median density, temperature, beta, and temperature anisotropy values for the core(halo)[beam/strahl] components, with subscripts ec(eh)[eb], of all fit results, respectively, are , , , and . This work will also serve as a 1 au baseline and reference for missions like Parker Solar Probe and Solar Orbiter .
This study investigates the interactions between cloud dynamics and aerosols in idealized large-eddy simulations (LES) of Arctic mixed-phase stratocumulus clouds (AMPS) observed at Oliktok Point, Alaska, in April 2015. This case was chosen because it allows the cloud to form in response to radiative cooling starting from a cloud-free state, rather than requiring the cloud ice an dliquid to adjust to an initial cloudy state. Sensitivity studies are used to identify whether there are buffering feedbacks that limit the impact of aerosol perturbations. The results of this study indicate that perturbations in ice nucleating particles (INPs) dominate over cloud condensation nuclei(CCN) perturbations; i.e., an equivalent fractional decrease in CCN and INPs results in an increase in the cloud-top long wave cooling rate, even though the droplet effective radius increases and the cloud emissivity decreases.The dominant effect of ice in the simulated mixed-phase cloud is a thinning rather than a glaciation, causing the mixed-phase clouds to radiate as a grey body and the radiative properties of the cloud to be more sensitive to aerosol perturbations. It is demonstrated that allowing prognostic CCN and INPs causes a layering of the aerosols, with increased concentrations of CCN above cloud top and increased concentrations of INPs at the base of the cloud-driven mixed layer. This layering contributes to the maintenance of the cloud liquid, which drives the dynamics of the cloud system.
Abstract. The onset of ice nucleation in mixed-phase clouds determines the lifetime and microphysical properties of ice clouds. In this work, we develop a novel method that differentiates between various phases of mixed-phase clouds, such as clouds dominated by pure liquid or pure ice segments, compared with those having ice crystals surrounded by supercooled liquid water droplets or vice versa. Using this method, we examine the relationship between the macrophysical and microphysical properties of Southern Ocean mixed-phase clouds at −40 to 0 °C (e.g. stratiform and cumuliform clouds) based on the in situ aircraft-based observations during the US National Science Foundation Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) flight campaign. The results show that the exchange between supercooled liquid water and ice crystals from a macrophysical perspective, represented by the increasing spatial ratio of regions containing ice crystals relative to the total in-cloud region (defined as ice spatial ratio), is positively correlated with the phase exchange from a microphysical perspective, represented by the increasing ice water content (IWC), decreasing liquid water content (LWC), increasing ice mass fraction, and increasing ice particle number fraction (IPNF). The mass exchange between liquid and ice becomes more significant during phase 3 when pure ice cloud regions (ICRs) start to appear. Occurrence frequencies of cloud thermodynamic phases show a significant phase change from liquid to ice at a similar temperature (i.e. −17.5 °C) among three types of definitions of mixed-phase clouds based on ice spatial ratio, ice mass fraction, or IPNF. Aerosol indirect effects are quantified for different phases using number concentrations of aerosols greater than 100 or 500 nm (N>100 and N>500, respectively). N>500 shows stronger positive correlations with ice spatial ratios compared with N>100. This result indicates that larger aerosols potentially contain ice-nucleating particles (INPs), which facilitate the formation of ice crystals in mixed-phase clouds. The impact of N>500 is also more significant in phase 2 when ice crystals just start to appear in the mixed phase compared with phase 3 when pure ICRs have formed, possibly due to the competing aerosol indirect effects on primary and secondary ice production in phase 3. The thermodynamic and dynamic conditions are quantified for each phase. The results show stronger in-cloud turbulence and higher updraughts in phases 2 and 3 when liquid and ice coexist compared with pure liquid or ice (phases 1 and 4, respectively). The highest updraughts and turbulence are seen in phase 3 when supercooled liquid droplets are surrounded by ice crystals. These results indicate both updraughts and turbulence support the maintenance of supercooled liquid water amongst ice crystals. Overall, these results illustrate the varying effects of aerosols, thermodynamics, and dynamics through various stages of mixed-phase cloud evolution based on this new method that categorizes cloud phases.
Explore the source record for details and available documents.