SEARCH · Search NASA

Results for “constraint reduction”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A Scalable Space-Time Domain Decomposition Approach for Solving Large Scale Nonlinear Regularized Inverse Ill Posed Problems in 4D Variational Data Assimilation

We address the development of innovative algorithms designed to solve the strong-constraint Four Dimensional Variational Data Assimilation (4DVar DA) problems in large scale applications. We present a space-time decomposition approach which employs the whole domain decomposition, i.e. both along the spacial and temporal direction in the overlapping case, and the partitioning of both the solution and the operator. Starting from the global functional defined on the entire domain, we get to a sort of regularized local functionals on the set of sub domains providing the order reduction of both the predictive and the Data Assimilation models. The algorithm convergence is developed. Performance in terms of reduction of time complexity and algorithmic scalability is discussed on the Shallow Water Equations on the sphere. The number of state variables in the model, the number of observations in an assimilation cycle, as well as numerical parameters as the discretization step in time and in space domain are defined on the basis of discretization grid used by data available at repository Ocean Synthesis/Reanalysis Directory of Hamburg University.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Safety Basis Compliant Flexible Operations [Capstone Project]

Compliant nuclear operations must be weighed against the radiological risk to the worker, public, and the environment. This project set out to look at one facility at Los Alamos National Laboratory (LANL or the Laboratory) and compare the risk with the operational restrictions associated with the Technical Safety Requirements (TSRs). The hope was to identify controls that could be reduced or simplified while maintaining compliance with 10 CFR 830, Nuclear Safety Management. This would result in increasing facility operational flexibility and efficiency. Controls were identified to be simplified or eliminated based on risk. This identification was initially envisioned from the perspective of reanalysis of the accident and hazard analysis consequences. Additionally, it was envisioned that fire testing could show that the consequences of a fire are not as dramatic as analyzed. Specific facility controls were identified that cause operational inefficiencies. Based on data collection, two controls were identified for further evaluation. The first control was combustible loading within the waste storage buildings, and the second control is related to the seismic power cutoff switch. In evaluating the combustible loading control, the thought is that the fire suppression system (FSS) provides adequate protection and the combustible loading controls do not need elevation to a TSR control. The FSS, along with the Fire Protection Program, would provide protection to prevent a fire from impacting large amounts of stored waste material. The seismic switch was more of a challenge. When reviewing the facility operation, it was determined that the facility does not operate anywhere near the established material-at-risk (MAR – radioactive material) limits. By adjusting these limits, a reduced consequence could be achieved and hence, a significant risk reduction achieved. A new operating mode would be used to achieve this reduction. A reduced MAR limit was successfully used as a test case for snow removal during this project. This project expands on that trial case as related to combustible loading and the seismic power cutoff switch TSRs. This will increase the efficiency of operations with reduced TSR controls while reducing overall facility risk. This is a win-win situation for Triad National Security, LLC (Triad) and the U.S. Department of Energy (DOE). Due to schedule and budgetary constraints, not all envisioned work scope was accomplished. However, it is maintained in this report for future efforts.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2026

Memory-efficient nonsmooth dynamic optimization using adaptive randomized compression

Dynamic optimization problems arise in many applications including flow control, full waveform inversion, and medical imaging. These problems are plagued by significant computational challenges. One such challenge — and the focus of this work — is the memory limitation induced by the size of the underlying dynamical system. In particular, the entire dynamic trajectory is required for derivative computation and therefore must be stored or recomputed using, e.g., checkpointing. Although recent work demonstrated the use of adaptive randomized sketching to overcome the memory challenge, that work only applies to smooth unconstrained problems, prohibiting its use for nonsmooth regularized and constrained problems. The inclusion of nonsmooth regularizers and constraints is critical as they often arise in an attempt to preserve certain physical properties or to promote sparsity. To solve these problems, we introduce a trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function that leverages randomized sketching to compress the dynamical system trajectories and adaptively adjust the sketch rank to satisfy a gradient inexactness condition. We prove convergence of this algorithm and demonstrate that it achieves substantial memory reduction on three discretized PDE-constrained optimization applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Power generation-cooling water Nexus: Impacts of cooling water shortage on power system operation - a simulation case study in Illinois, U.S

Cooling water shortage, frequently attributed to drought and heat waves, poses a significant threat to the operations of thermoelectric power plants and further poses a challenge for the entire power system and environmental stakeholders. Recognizing the critical nexus between power generation and cooling water availability and the potential ability of power generations to adjust generation schedules during cooling water shortages, this paper introduces a security-constrained unit commitment and economic dispatch model considering water-energy nexus. In specific, the model is augmented with a unit-level cooling water requirement (CWR) model and multi-level cooling water availability (CWA) constraints. The unit-level CWR model quantifies the cooling water withdrawal per MWh of power generation, taking into account factors such as thermoelectric generation technologies, cooling system technologies, and environmental parameters. The multi-level CWA constraints incorporate pump-level, plant-level, watershed-level, and forced minimum power constraints, utilizing data derived from actual-based cooling water shortage scenarios. Using a simulation case study in Illinois, United States, this research examines the reliability, economic, and environmental implications of cooling water shortages on power system operations. The results show that Illinois may experience 10-15% daily load curtailment and severe congestion between certain regions from the east to central during cooling water shortages, while once-through and wet-tower units experience a 52% and 17% reduction in power generation. In conclusion, overall cooling water withdrawal decreases by 24-38% as severity intensifies.

Cooling water shortage↗

DOE OSTI2020

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

REIMAGINING HEAT EXCHANGERS FOR NEXT GENERATION ENVIRONMENTAL SYSTEMS

Air-to-refrigerant heat exchangers (HXs) are essential components in space conditioning, refrigeration, and power systems, and recent efforts have focused on making these devices more compact, reducing refrigerant charge and lowering manufacturing costs. Historically, HX innovation has been limited by available computational resources, design tools, and manufacturing constraints. The best available technologies utilize tube-fin and micro- or macro-channel tubes with fins, which are not necessarily the optimal designs achievable with current technology. In this paper, we highlight the latest advancements in air-to-refrigerant HXs, specifically emphasizing innovations achieved through shape and topology optimization. A multi-scale design optimization approach is introduced, alongside similar methods in literature, which enable highly sophisticated shape-optimized tube designs with more than 50% reduction in size and 25% reduction in refrigerant charge, essential for A3 refrigerant charge limit compliance. The frameworks integrate traditional heat and mass transfer science with state-of-the-art machine learning, genetic algorithms, and adjoint algorithms to create novel designs. While many of these innovative designs may not be manufacturable using conventional methods, they allow us explore the boundaries of what is possible. These novel air-to-refrigerant HXs are key enablers for ultra-low-refrigerant charge heat pump and refrigeration systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Energy efficiency in industrial drying: A hybrid ultrasonic system with a novel dynamic optimization framework

Drying processes are among the most energy-consuming operations in industrial and manufacturing settings, demanding strategic selection, design, and control for enhanced efficiency. Advancing drying technologies is critical for improving sustainability, lowering energy use, reducing carbon emissions, and minimizing waste. This study explores two innovative strategies aimed at transforming drying processes into sustainable, low-carbon systems by reducing energy consumption, minimizing waste, and maintaining a strong emphasis on preserving product quality. The first strategy showcases a sub-pilot scale hybrid ultrasonic-convective dryer for agrifood products. This technology, powered by electricity (process electrification), integrates non-thermal ultrasonic dehydration with convective heating and is presented as a sustainable and energy-efficient solution that enhances eco-friendly practices. The second strategy involves introducing and implementing a novel, multiobjective, mixed integer dynamic optimization technique to determine the optimal time-dependent process parameter values for the drying operation. This optimization technique yields operating conditions that are piecewise constant in time aiming to maximize the energy efficiency of the hybrid ultrasonic-convective dryer while ensuring strict adherence to product quality constraints. By adopting the hybrid ultrasonic-convective dryer, a notable 35% improvement in energy efficiency was achieved compared to conventional hot-air drying systems for drying apple slices. The proposed optimization framework further enhanced energy efficiency by nearly 14% over the most efficient process on the identical testbed, under static operating conditions. The reported enhancements have been experimentally validated. Regarding drying time (thereby improving production yield), the developed hybrid ultrasonic-convective dryer demonstrates as much as a 41% reduction in total processing time, which is further optimized by an additional 10% using our proposed optimization framework. The research outcomes have profound implications for the design and operation of drying systems, encompassing crucial aspects such as process electrification, cost-effectiveness, energy savings, time efficiency, product yield, product quality, and process automation.

Dynamic optimization↗

DOE OSTI2025

Moltensaltpropnet

MoltenSaltPropnet is a physics-informed machine learning framework that aims to predict the thermophysical properties of molten fluoride and chloride salt mixtures, which are crucial for the design and safety of Generation IV molten salt reactors. The code processes data from the Molten-Salt Thermal Properties Database (MSTDB-TP) and the Janz compendium, converting critically evaluated correlations into fast, differentiable surrogate models for density, viscosity, thermal conductivity, and heat capacity across 448 distinct salt systems. The implementation consists of several key components: 1. Data Curation: The code parses and cleans the raw data, normalizing elemental mole fractions and extracting relevant regression coefficients for various thermophysical properties. 2. Feature Engineering: It generates fixed-length numerical descriptors that encapsulate the composition and temperature, incorporating polynomial interaction terms and dimensionality-reduction techniques to optimize model performance. 3. Coefficient Learning: Four different machine learning architectures are employed: a deep residual network (ResNet), a Kolmogorov–Arnold network (KAN), a sparsity-inducing neural network (SNN), and classical regression models. Each model learns to predict coefficients that define the temperature-dependent correlations for the thermophysical properties. 4. Property Reconstruction: The predicted coefficients are used to compute temperature-dependent property values, ensuring positivity and monotonic trends through a composite loss function that enforces physical constraints. 5. User Interface: An open-source web application enables users to filter the database, train task-specific models, and visualize the results, allowing for rapid exploration of candidate salt mixtures. MoltenSaltPropnet bridges the gap between limited experimental data and high-fidelity reactor simulations, providing a powerful tool for researchers in the field of molten salt reactors and advanced nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗

DOE OSTI2023

Optimization of Operating Profit for the FEED of the PZAS Process for CO 2 Capture on a 460 MW NGCC

A rating model of PZAS (Piperazine with the Advanced Stripper) for a 460 MW natural gas combined cycle was developed from vendor-provided equipment performance data. Surrogate rating models were built and used to optimize variable profit of the system subject to natural gas and electricity price at a fixed value for CO 2 of $\$80$/tonne. Several process constraints were considered between ambient temperatures of 4.4 and 40.5 °C. Lower lean loading resulted in lower reboiler heat duty of about 2.2 GJ/tonne. Lean loading at 0.2 mol/equiv N was sufficient at an ambient temperature of 4.4 °C to maximize variable profit. Solvent precipitation was avoided even at low lean loading (0.18 mol/mol) in the absence of a trim cooler for the lean solvent between the absorber and the cold cross exchanger. Optimal CO 2 removal was always higher than the design case removal of 90% and achieved a peak of 93.6% for 4.4 ºC ambient temperature. Lean solvent rate was a significant driver of profitability. The maximum capacity of the compressor was a significant constraint on profitability. Increasing compressor capacity by 20–30% could increase variable profit by about 21% in the case of 0.18 lean loading at 4.4 °C ambient temperature. Furthermore, this increase in profit was lower at higher ambient temperature and higher lean loading. The breakeven gas price was $\$4.6$/MMBtu at high lean loading of 0.22 mol/mol. Below this value, profit was maximized by increasing the boiler heat rate. Above this value, profit could not be increased by burning additional natural gas, requiring a reduction in boiler heat rate and consequently, specific heat duty. At lower lean loading (0.2, 0.18 mol/mol), this transition in heat duty was not observed due to the lower heat duty. As the ambient temperature increased, a similar trend in heat duty was observed at 0.22 lean loading, but the transition to lower heat duty occurred earlier at $\$4$/MMBtu. This was because of the lower CO 2 production at higher ambient temperature, making the process unprofitable at a gas price higher than $\$3$/MMBtu unless the heat duty of the stripper and natural gas boiler is steeply decreased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Glacial‐Interglacial and Millennial‐Scale Changes in Nitrous Oxide Emissions Pathways and Source Regions

During the transition from the Last Glacial Maximum (LGM) to the Holocene, the atmospheric N 2 O mole fraction increased by 80 nmol mol −1 . Using ice core measurements of N 2 O isotopomer ratios, we show that this increase was driven by increases in both nitrification and denitrification, with the relative partitioning between both production pathways depending on the assumed isotopic end-member source signatures. Similarly, we also attribute a 35 nmol mol −1 N 2 O mole fraction increase during the Heinrich Stadial 4/Dansgaard Oeschger 8 (HS4/DO8) millennial-scale event to increases in both N 2 O production pathways. In contrast, the 25 nmol mol −1 N 2 O mole fraction decrease during the Younger Dryas was driven almost exclusively by a decrease in nitrification. The deglacial and HS4/DO8 increases in N 2 O production occurred in both marine and terrestrial environments, with the terrestrial source responding faster to warming by about two centuries. Constraints on changes in nitrification and denitrification emissions are robust and consistent with previous studies showing the sensitivity of N 2 O emissions to abrupt Northern Hemisphere warming. This study demonstrates for the first time the importance of both denitrification and nitrification pathways in driving source changes. Absolute emissions are more uncertain due to uncertainty about source isotopomer signatures. For instance, the contribution of denitrification to emissions at the LGM shifts from (65 ± 10) % to (91 ± 6) % when factoring in isotope enrichment due to partial reduction of N 2 O to N 2 during denitrification. Reducing uncertainty in source signatures will increase the power of ice core N 2 O isotope records in deducing environmental change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Immunogenicity and protective efficacy of a rhesus adenoviral vaccine targeting conserved COVID-19 replication transcription complex

The COVID-19 pandemic marks the third coronavirus pandemic this century (SARS-CoV-1, MERS, SARS-CoV-2), emphasizing the need to identify and evaluate conserved immunogens for a pan-sarbecovirus vaccine. Here we investigate the potential utility of a T-cell vaccine strategy targeting conserved regions of the sarbecovirus proteome. We identified the most conserved regions of the sarbecovirus proteome as portions of the RNA-dependent RNA polymerase (RdRp) and Helicase proteins, both of which are part of the coronavirus replication transcription complex (RTC). Fitness constraints suggest that as SARS-CoV-2 continues to evolve these regions may better preserve cross-reactive potential of T-cell responses than Spike, Nucleocapsid, or Membrane proteins. We sought to determine if vaccine-elicited T-cell responses to the highly conserved regions of the RTC would reduce viral loads following challenge with SARS-CoV-2 in mice using a rhesus adenovirus serotype 52 (RhAd52) vector. The RhAd52.CoV.Consv vaccine generated robust cellular immunity in mice and led to significant reductions in viral loads in the nasal turbinates following challenge with a mouse-adapted SARS-CoV-2. These data suggest the potential utility of T-cell targeting of conserved regions for a pan-sarbecovirus vaccine.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Information content of higher order galaxy correlation functions

The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance–redshift relationship. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence, all the cosmological information can be extracted from their two-point correlation function. Subsequent non-linear evolution under gravity, as well as halo and then galaxy formation, generates higher order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field’s power spectrum is equal to that of all the higher order functions of the final, non-linear field. This claim implies that reconstruction of the initial density field from the non-linear field renders analysis of higher order correlation functions of the latter superfluous. We show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the two and three-point correlation functions can, in some cases, exceed those offered by the initial power spectrum. We provide a mathematical justification for this claim and demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot-noise, taking modes up to k max = 0.2 h Mpc −1 , using the bispectrum alone offers a factor of 2 reduction in the variance on the cosmic distance scale relative to that available from the linear power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

A Comfort Model Simplification for Tight Integration with Grid Service Optimizations

Localized generation, storage, and intelligent power electronics are increasingly being integrated into building. The excess storage, thermal, and generation capacity can be used to develop new markets for ancillary services such as peak reduction, voltage and frequency support, or load flattening. Widespread participation in these markets will make the power grid more resilient. However, widespread participation in these markets will require a simple and intuitive control over the excess capacity by the building occupants. Most importantly for adopting these technologies, participation in these markets should not require sacrificing occupant comfort. Many of the proposed control and optimization technologies for supplying the ancillary services do not explicitly incorporate occupant comfort into their models and when they do, generally utilize a rigid constraint on the indoor air temperature set a priori. This reduces flexibility in using excess thermal capacity. It also cannot take into account uncertainties in local factors that affect thermal comfort such as clothing insulation or relative humidity. In this paper we present and justify some assumptions that simplify the standard predictors of thermal comfort. We also develop a regression model of the thermal comfort that is linear with respect to the indoor air temperature. Next we develop the concept of thermal comfort variation that eliminates the dependence of the thermal comfort on uncertain thermal comfort factors. Then, we present a quadratic program utilizing the thermal comfort variation in both the objective function and as a constraint. Finally, we show the results of a time-of-use cost optimization that utilizes the simplified thermal comfort model.

Melin, Alexander↗

DOE OSTI2020

Observed Sensitivity of Low-Cloud Radiative Effects to Meteorological Perturbations over the Global Oceans

Understanding how marine low clouds and their radiative effects respond to changing meteorological conditions is crucial to constrain low-cloud feedbacks to greenhouse warming and internal climate variability. In this paper, we use observations to quantify the low-cloud radiative response to meteorological perturbations over the global oceans to shed light on physical processes governing low-cloud and planetary radiation budget variability in different climate regimes. We assess the independent effect of perturbations in sea surface temperature, estimated inversion strength, horizontal surface temperature advection, 700-hPa relative humidity, 700-hPa vertical velocity, and near-surface wind speed. Stronger inversions and stronger cold advection greatly enhance low-level cloudiness and planetary albedo in eastern ocean stratocumulus and midlatitude regimes. Warming of the sea surface drives pronounced reductions of eastern ocean stratocumulus cloud amount and optical depth, and hence reflectivity, but has a weaker and more variable impact on low clouds in the tropics and middle latitudes. By reducing entrainment drying, higher free-tropospheric relative humidity enhances low-level cloudiness. At low latitudes, where cold advection destabilizes the boundary layer, stronger winds enhance low-level cloudiness; by contrast, wind speed variations have weak influence at midlatitudes where warm advection frequently stabilizes the marine boundary layer, thus inhibiting vertical mixing. These observational constraints provide a framework for understanding and evaluating marine low-cloud feedbacks and their simulation by models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Mooring System Cost Estimates for Wave Energy Farms in Shared Mooring Arrays

As wave energy converters (WECs) become more advanced and cost-efficient, so too must their mooring systems. A key question in the development plans for WECs is the cost of the mooring system, particularly for large wave farms. WEC devices deployed in a WEC farm array, where each device can be connected by shared mooring lines in various array layouts, have potential to reduce mooring system costs significantly. This paper presents the modeling and designing of mooring systems for large WEC arrays and the calculation of the cost of each mooring system to determine the change in cost as the number of WECs in a farm increases. A baseline mooring system for a single floating oscillating water column (OWC) WEC was developed for this analysis. The mooring system utilizes four anchored mooring lines connected to a square assembly of wire rope mooring lines, supported by four floating buoys, and attached to the floating WEC by four polyester rope mooring lines. This assembly, referred to as a floating cell, can be tiled to form various rectangular WEC arrays. The objective of this analysis is to determine how the mooring system cost changes as more WECs are added to an array layout, each with their own interconnected floating cell. To do this, complete mooring systems need to be designed for each WEC array layout. To narrow down the design space of a WEC array mooring system, a couple assumptions were made. It was assumed that the floating cell parameters of the baseline design were to stay constant across all floating cells in the WEC array. It was also assumed that the anchored mooring lines would be of the baseline configuration, a predominantly chain mooring line with a short section of polyester rope near the fairlead, and a drag-embedment anchor. These anchored lines were assumed to extend from the outer edges of the WEC array, inline with the headings of the wire rope mooring lines of the floating cells, or diagonal if extending from a corner of the WEC array. Full mooring systems were designed for 2xN and NxN WEC array layouts and efficiently simulated in the mooring dynamics simulation tool, MoorDyn, to ensure all dynamic constraints were met. The system costs were calculated and then refined by shortening unnecessary chain line lengths and reducing the chain diameters of the downstream anchored mooring lines. It was found that, in general, mooring system costs per WEC decrease when WECs are installed in an array. Compared to the baseline mooring system for a single WEC, the 2x3 mooring system array had the lowest mooring system cost per WEC, reducing the cost per WEC by 59%. The 3x3 and 4x4 array mooring systems also saw significant reductions in cost per WEC but had negligible cost savings between the two designs, primarily because the larger 4x4 mooring system requires larger chain diameters, which increases cost. These results provide an interesting glimpse into modeling, designing, and calculating the cost of mooring systems for large WEC arrays.

cost↗

DOE OSTI2022

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗