Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel processing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,171 records · Page 65

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

Tula: Optimizing Time, Cost, and Generalization in Distributed Large-Batch Training

Distributed training increases the number of batches processed per iteration either by scaling-out (adding more nodes) or scaling-up (increasing the batch-size). However, the largest configuration does not necessarily yield the best performance. Horizontal scaling introduces additional communication overhead, while vertical scaling is constrained by computation cost and device memory limits. Thus, simply increasing the batch-size leads to diminishing returns: training time and cost decrease initially but eventually plateaus, creating a knee-point in the time/cost vs. batch-size pareto curve. The optimal batch-size therefore depends on the underlying model, data and available compute resources. Large batches also suffer from worse model quality due to the well-known “generalization gap”. In this paper, we present Tula, an online service that automatically optimizes time, cost, and convergence quality for large-batch training of convolutional models. It combines parallel-systems modeling with statistical performance prediction to identify the optimal batchsize. Tula predicts training time and cost within 7.5−14% error across multiple models, and achieves up to 20× overall speedup and improves test accuracy by ≈9% on average over standard large-batch training on various vision tasks, thus successfully mitigating the generalization gap and accelerating training at the same time.

Tyagi, Sahil [ORNL] (ORCID:0009000783144745)↗

Scalable and accurate multi-GPU-based image reconstruction of large-scale ptychography data

Abstract While the advances in synchrotron light sources, together with the development of focusing optics and detectors, allow nanoscale ptychographic imaging of materials and biological specimens, the corresponding experiments can yield terabyte-scale volumes of data that can impose a heavy burden on the computing platform. Although graphics processing units (GPUs) provide high performance for such large-scale ptychography datasets, a single GPU is typically insufficient for analysis and reconstruction. Several works have considered leveraging multiple GPUs to accelerate the ptychographic reconstruction. However, most of these works utilize only the Message Passing Interface to handle the communications between GPUs. This approach poses inefficiency for a hardware configuration that has multiple GPUs in a single node, especially while reconstructing a single large projection, since it provides no optimizations to handle the heterogeneous GPU interconnections containing both low-speed (e.g., PCIe) and high-speed links (e.g., NVLink). In this paper, we provide an optimized intranode multi-GPU implementation that can efficiently solve large-scale ptychographic reconstruction problems. We focus on the maximum likelihood reconstruction problem using a conjugate gradient (CG) method for the solution and propose a novel hybrid parallelization model to address the performance bottlenecks in the CG solver. Accordingly, we have developed a tool, called PtyGer ( Pty chographic G PU(multipl e )-based r econstruction), implementing our hybrid parallelization model design. A comprehensive evaluation verifies that PtyGer can fully preserve the original algorithm’s accuracy while achieving outstanding intranode GPU scalability.

97 MATHEMATICS AND COMPUTING↗

Molecular underpinnings of ssDNA specificity by Rep HUH-endonucleases and implications for HUH-tag multiplexing and engineering

Abstract Replication initiator proteins (Reps) from the HUH-endonuclease superfamily process specific single-stranded DNA (ssDNA) sequences to initiate rolling circle/hairpin replication in viruses, such as crop ravaging geminiviruses and human disease causing parvoviruses. In biotechnology contexts, Reps are the basis for HUH-tag bioconjugation and a critical adeno-associated virus genome integration tool. We solved the first co-crystal structures of Reps complexed to ssDNA, revealing a key motif for conferring sequence specificity and for anchoring a bent DNA architecture. In combination, we developed a deep sequencing cleavage assay, termed HUH-seq, to interrogate subtleties in Rep specificity and demonstrate how differences can be exploited for multiplexed HUH-tagging. Together, our insights allowed engineering of only four amino acids in a Rep chimera to predictably alter sequence specificity. These results have important implications for modulating viral infections, developing Rep-based genomic integration tools, and enabling massively parallel HUH-tag barcoding and bioconjugation applications.

59 BASIC BIOLOGICAL SCIENCES↗

An innovative approach for atrazine electrochemical oxidation modelling: Process parameter effect, intermediate formation and kinetic constant assessment

Water reuse for irrigation activities is becoming a crucial worldwide challenge due to the depletion of water sources. Anyway, agricultural drainage can potentially contain dangerous contaminants such as metals, pesticides, and herbicides, including atrazine. To address the need for agriculture wastewater purification, we investigated atrazine removal from simulated wastewater by electro-oxidation using platinum-coated titanium electrodes on a lab-scale experimental apparatus. The effects of electrolyte composition and concentration, i.e. ionic strength and applied current density on atrazine removal, were investigated. The results demonstrated that the electrochemical oxidation of the herbicide occurred through two routes, depending on the presence or absence of oxidizing chlorine species. The generation of intermediates during the treatment was monitored and quantified by evaluating the effect of an inert electrolyte (NaClO 4 ) versus an oxidizable chlorine species (NaCl). In both experimental conditions, five intermediates were identified, including desethyl-atrazine (DEA), hydroxyatrazine (ATZ-OH), desisopropyl-atrazine (DIA) and desethyl-desisopropyl-atrazine (DEDIA). A degradation mechanism and a model for describing hydroxyl radicals and active chlorine species contributions at ATZ oxidation were also proposed. Intermediate evolution profiles suggest that ATZ degradation can be considered as a series–parallel reaction system. Finally, the energy requirement assessment for ATZ removal was carried out. The highest ATZ removal (≅98%) was achieved with NaCl = 0.08 M, J = 60 A/m −2 , and E C = 5.83 kWh m −3 . Results highlight that atrazine removal was improved when an active chlorine species (NaCl) was present in the water solution. Moreover, the addition of chlorine species during electro-oxidation is an energy-saving strategy. Collectively, electro-oxidation technique can be efficiently applied to treat polluted water in order to meet the needs of recycling water quality and reduce resource consumption.

Electro-chemical oxidation↗

Mn(II)-induced phase transformation of Mn(IV) oxide in seawater

Manganese (Mn) oxides are key components of oceanic and lacustrine Mn nodules and influence metal cycling through oxidation and adsorption processes. Layered Mn oxides (LMOs) are the most common minerals in these nodules and the immediate products of microbially mediated Mn(II) oxidation by O 2 . LMOs can transform into tunneled Mn oxides (TMOs), Mn oxyhydroxides (MnOOH), or Mn(II,III) phase (Mn 3 O 4 ). LMOs often concur with Mn(II) in the environment and the adsorption and oxidation of Mn(II) by LMOs can greatly promote the transformation of LMOs to those phases. However, the Mn(II)-promoted transformation of LMOs in seawater—rich in various cations (300 mM Na + , 10 mM K + , 50 mM Ca 2+ , and 10 mM Mg 2+ ) remains poorly understood. We examined the transformation of δ-MnO 2 in artificial seawater (pH 8.2) under anoxic conditions with the Mn(II)/MnO 2 ratio (r) ranging from 0.08 to 3.83. To assess the effect of ionic strength (IS), parallel experiments were conducted in a mixed 530 mM NaCl and 10 mM KCl solution (having seawater ionic strength but without Ca 2+ and Mg 2+ ) and in 100 mM NaCl solution as a control. At low r (0.08), δ-MnO 2 transformed into triclinic birnessite and a 4 × 4 TMO in 100 mM NaCl solution, which, however, was suppressed in seawater due to strong interactions of Ca 2+ /Mg 2+ with δ-MnO 2 . In the mixed 530 mM NaCl and 10 mM KCl solution (the same ionic strength as of seawater), the transformation occurred extensively but the products had lower crystallinity compared to in 100 mM NaCl solution. At the high Mn(II)/MnO 2 ratios (0.5 ≤ r ≤ 3.83), δ-MnO 2 transformed extensively into MnOOH phases and hausmannite (Mn 3 O 4 ) in 100 mM NaCl solution. The seawater suppressed the transformation, but the suppression became weaker with increasing Mn(II)/MnO 2 ratio. For example, the transformation was completely suppressed at r = 0.5 but essentially negligible at r = 3.83. The suppression at these high Mn(II)/MnO 2 ratios was mainly ascribed to the influence of Ca 2+ and Mg 2+ rather than of the high IS, and the weaker suppression at the higher Mn(II)/MnO 2 ratio suggests stronger competition of Mn(II) with Ca 2+ /Mg 2+ for interacting with δ-MnO 2 . Moreover, the composition and crystallinity of the transformation products (i.e., the relative abundance of MnOOH (α, β, and γ) and Mn 3 O 4 ) were influenced by both the high ionic strength and the presence of Ca 2+ and Mg 2+ . Therefore, even though Ca 2+ and Mg 2+ concentrations are much lower than Na + in seawater, their impacts on the transformation are dominant. Our study explains why MnOOH phases, hausmannite, and TMOs are less common than LMOs in oceanic environments, partially because seawater chemistry suppresses their formation. LMOs are the most reactive for metal adsorption and oxidation among all Mn oxides. Thus, the high stability of LMOs in an oceanic environment confers the high impacts of Mn oxides on metal cycling in the ocean.

Divalent manganese↗

MECHANISMS OF RADIATION CARCINOGENESIS: WHAT IS REALLY INDUCED?

Abstract It has been difficult to understand why the relative risk for cancer decreases with an increase in time since an exposure to radiation. It was recently recognized that this decline can be explained by a parallel shift of the age-related cancer mortality curve toward younger ages. In fact, it has been known for many years that mouse survival curves exhibit a parallel shift toward younger ages following an exposure to radiation, but it was not recognized that the mutation induction theory is incompatible with this parallel shift. This is because a parallel shift in the survival curve implies that all the irradiated individuals are affected, but the mutation induction theory assumes that only a fraction of the irradiated individuals is affected following an exposure to radiation. Thus, it seems likely that the target of radiation action, which leads to carcinogenesis, is not restricted to epithelial cells but includes all of the surrounding cells leading to an altered microenvironment. Since it is repeatedly observed that radiation-exposed normal tissues can stimulate transplanted or spontaneously arising tumor cells to grow faster, worsen the malignant phenotypes and finally kill the host earlier than usual, an exposure to radiation seems most likely to cause tissue inflammation, which creates conditions favorable for the growth of spontaneously arising tumor cells. This new concept suggests that it might be possible to attenuate the extent of radiation carcinogenesis by intervening in tissue inflammatory processes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Status of the laser manipulations of H- beam at J-PARC

The laser manipulations of H- ion beam by single or double neutralization is a very promising technique and highly essential to utilized in accelerator processes such as stripping, pulse chopping, collimation, extraction, and beam diagnostics for the present and future high-intensity proton accelerators. At J-PARC, we are preparing for a POP (Proof-of-Principle) demonstration 400 MeV H- stripping by using only lasers. A prototype YAG laser system and a laser cavity system to reduce the laser power are being developed through 3 MeV H- neutralization studies. Fermilab utilizes H- neutralization at 0.75 MeV by establishing a laser Notcher system for a gap in the H- pulse needed for a clean beam extraction from the ring. To minimize the laser power and maximize the interaction efficiency are common issue at both laboratories. Under the US-Japan collaboration for high-intensity neutrino beam, we are closely working to establish recycling/reusing of seed laser pulses to reduce a size of the laser system. In parallel, we have also developed non-destructive beam diagnostic systems at lower H- energy of 3 MeV, which will be implemented to the 400 MeV as well as easily applicable to the Fermilab linac including PIP-II linac. Installation of the laser system at J-PARC linac for 400 MeV H- stripping is ongoing to start the POP experimental study in 2023.

43 PARTICLE ACCELERATORS↗

HPDR: High-Performance Portable Scientific Data Reduction Framework

The rapid growth in scientific data generation is outpacing advancements in computing systems necessary for efficient storage, transfer, and analysis, particularly in the context of exascale computing. With the deployment of first-generation exascale computing systems and next-generation experimental facilities, this gap is widening and necessitates effective data reduction techniques to manage enormous data volumes. Over the past decade, various data reduction methods, including lossless compression, error-controlled lossy compression, and data refactoring, have been developed to accelerate I/O in scientific workflows. Despite significant reductions in data volume, these methods introduce considerable computational overhead, which can become the new bottleneck in data processing. To mitigate this, GPU-accelerated data reduction algorithms have been introduced. However, challenges remain in their integration into exascale workflows, including limited portability across different GPU architectures, substantial memory transfer overhead, and reduced scalability on dense multi-GPU systems. To address these challenges, we propose HPDR, a high-performance and portable data reduction framework. HPDR is designed to enable the execution of state-of-the-art reduction algorithms across diverse processor architectures while reducing memory transfer overhead to 2.3 % of the original, resulting in up to 3.5× faster throughput compared to existing solutions. It also achieves up to 96% of the theoretical speedup in multi-GPU settings. In addition, evaluations on accelerating I/O operations at scale up to 1,024 nodes of the Frontier supercomputer demonstrate that HPDR can achieve up to 103 TB/s reduction throughput, providing up to 4× acceleration in parallel I/O performance compared to existing data reduction routines. This work highlights the potential of HPDR to significantly enhance data reduction efficiency in exascale computing environments.

Chen, Jieyang [University of Oregon]↗

Parallel Grid-Forming Inverter-Driven Black Start: Power-Hardware-in-the-Loop Validation: Preprint

With power systems encountering increasing deployment levels of inverter-based resources (IBRs), system restoration using grid-forming (GFM) IBRs has gained attention. Engineered to establish grid voltages in the absence of a stiff grid, black-start-capable GFM IBRs are expected to enhance power system resilience by playing a critical role in bottom-up system restoration. This paper experimentally studies the feasibility of the novel approach through power-hardware-in-the-loop (PHIL) testing and validation with a commercial GFM inverter. The PHIL test setup demonstrates an inverter-driven black start of a 5-MW unbalanced distribution feeder where two GFM inverters collectively black start the feeder: one is a commercial hardware GFM inverter interacting through the PHIL interface, and the other is a software GFM emulated by a real-time simulation with full electromagnetic transient (EMT) models. Both GFM inverters are equipped with negative-sequence voltage control to suppress the voltage imbalance resulting from the unbalanced loading, allowing us to study the dynamic interactions between the GFMs with their control parameters unknown. To evaluate the dynamic behavior of the GFM inverters under the entire black-start process, the EMT model of the distribution system details the transformer and motor dynamics to emulate their inrush and startup behaviors. It abstracts conventional grid-connected inverter dynamics, i.e., grid-following inverters. Oscilloscope measurements of the hardware GFM inverter are also presented for critical steps. Takeaways for further study and field deployment are provided.

grid-forming inverter↗

A fast and large bandwidth superconducting variable coupler

Variable microwave-frequency couplers are highly useful components in classical communication systems and likely will play an important role in quantum communication applications. Conventional semiconductor-based microwave couplers have been used with superconducting quantum circuits, enabling, for example, the in situ measurements of multiple devices via a common readout chain. However, the semiconducting elements are lossy and furthermore dissipate energy when switched, making them unsuitable for cryogenic applications requiring rapid, repeated switching. Superconducting Josephson junction-based couplers can be designed for dissipation-free operation with fast switching and are easily integrated with superconducting quantum circuits. These enable on-chip, quantum-coherent routing of microwave photons, providing an appealing alternative to semiconductor switches. Here, we present and characterize a chip-based broadband microwave variable coupler, tunable over 4-8GHz with over 1.5GHz instantaneous bandwidth, based on the superconducting quantum interference device with two parallel Josephson junctions. The coupler is dissipation-free and features large on-off ratios in excess of 40dB, and the coupling can be changed in about 10ns. The simple design presented here can be readily integrated with superconducting qubit circuits and can be easily generalized to realize a four- or more port device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Stimulation of soil respiration by elevated CO 2 is enhanced under nitrogen limitation in a decade-long grassland study

Whether and how CO 2 and nitrogen (N) availability interact to influence carbon (C) cycling processes such as soil respiration remains a question of considerable uncertainty in projecting future C–climate feedbacks, which are strongly influenced by multiple global change drivers, including elevated atmospheric CO 2 concentrations (eCO 2 ) and increased N deposition. However, because decades of research on the responses of ecosystems to eCO 2 and N enrichment have been done largely independently, their interactive effects on soil respiratory CO 2 efflux remain unresolved. Here, we show that in a multifactor free-air CO 2 enrichment experiment, BioCON (Biodiversity, CO 2 , and N deposition) in Minnesota, the positive response of soil respiration to eCO 2 gradually strengthened at ambient (low) N supply but not enriched (high) N supply for the 12-y experimental period from 1998 to 2009. In contrast to earlier years, eCO 2 stimulated soil respiration twice as much at low than at high N supply from 2006 to 2009. In parallel, microbial C degradation genes were significantly boosted by eCO 2 at low but not high N supply. Further, incorporating those functional genes into a coupled C–N ecosystem model reduced model parameter uncertainty and improved the projections of the effects of different CO 2 and N levels on soil respiration. If our observed results generalize to other ecosystems, they imply widely positive effects of eCO 2 on soil respiration even in infertile systems.

54 ENVIRONMENTAL SCIENCES↗

Pretraining Billion-Scale Geospatial Foundational Models on Frontier

As AI workloads increase in scope, generalization capability becomes challenging for small task-specific models and their demand for large amounts of labeled training samples increases. On the contrary, Foundation Models (FMs) are trained with internet-scale unlabeled data via self-supervised learning and have been shown to adapt to various tasks with minimal fine-tuning. Although large FMs have demonstrated significant impact in natural language processing and computer vision, efforts toward FMs for geospatial applications have been restricted to smaller size models, as pretraining larger models requires very large computing resources equipped with state-of-the-art hardware accelerators. Current satellite constellations collect 100+TBs of data a day, resulting in images that are billions of pixels and multimodal in nature. Such geospatial data poses unique challenges opening up new opportunities to develop FMs. We investigate billion scale FMs and HPC training profiles for geospatial applications by pretraining on publicly available data. We studied from end-to-end the performance and impact in the solution by scaling the model size. Our larger 3B parameter size model achieves up to 30% improvement in top1 scene classification accuracy when comparing a 100M parameter model. Moreover, we detail performance experiments on the Frontier supercomputer, America's first exascale system, where we study different model and data parallel approaches using PyTorch's Fully Sharded Data Parallel library. Specifically, we study variants of the Vision Transformer architecture (ViT), conducting performance analysis for ViT models with size up to 15B parameters. By discussing throughput and performance bottlenecks under different parallelism configurations, we offer insights on how to leverage such leadership-class HPC resources when developing large models for geospatial imagery applications.

Tsaris, Aristeidis (aris)↗

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct Functionalization of Established 3D-Printed Aza-Michael Liquid Crystal Elastomers with Donor–Acceptor Stenhouse Adducts

Extrusion 3D printing has advanced the manufacturing of complex liquid crystal elastomer (LCE) architectures. In parallel, donor–acceptor Stenhouse adducts (DASAs), a class of white-light-responsive photoswitches, have enabled both photochemical and photothermal LCE actuation. However, DASA–LCEs have yet to be extruded and 3D-printed. Two key challenges exist: DASA’s inherent sensitivity to heat and radicals can lead to degradation during ink preparation and printing, and small changes in the concentration of the added DASA component impact the properties of the extrudable ink, requiring reoptimization of well-established 3D-printing protocols. To overcome these challenges, we present a post-printing functionalization strategy that circumvents these limitations. Residual secondary amines, inherent to inks synthesized via standard aza-Michael addition, serve as active sites for covalent attachment of DASA photoresponsive groups following printing and cross-linking. Our method means that DASAs can be directly grafted onto 3D-printed aza-Michael LCEs without modifying the ink formulation or processing. The resulting DASA–LCEs exhibit wavelength tunability within the visible range and a variety of photothermal and photochemical responses. The post-functionalization can occur within 2 min and enables spatial control of the DASA concentration, producing films with tunable color gradients and locally varied photothermal and photochemical responses under visible light. In conclusion, this approach enables the rapid fabrication of DASA-based light-responsive LCEs using established ink formulations with the potential for the design of complex 3D architectures.

3D printing↗

Stochastic evaluation of fourth-order many-body perturbation energies

A scalable, stochastic algorithm evaluating the fourth-order many-body perturbation (MP4) correction to energy is proposed. Three hundred Goldstone diagrams representing the MP4 correction are computer generated and then converted into algebraic formulas expressed in terms of Green’s functions in real space and imaginary time. They are evaluated by the direct (i.e., non-Markov, non-Metropolis) Monte Carlo (MC) integration accelerated by the redundant-walker and control-variate algorithms. The resulting MC-MP4 method is efficiently parallelized and is shown to display O(n 5.3 ) size-dependence of cost, which is nearly two ranks lower than the O(n 7 ) dependence of the deterministic MP4 algorithm. Furthermore, it evaluates the MP4/aug-cc-pVDZ energy for benzene, naphthalene, phenanthrene, and corannulene with the statistical uncertainty of 10 mE h (1.1% of the total basis-set correlation energy), 38 mE h (2.6%), 110 mE h (5.5%), and 280 mE h (9.0%), respectively, after about 10 9 MC steps.

74 ATOMIC AND MOLECULAR PHYSICS↗

Project Title: Adaptive Vertical Grid Enhancement for E3SM (Final Technical Report)

During the first phase of the project (2018/05-2021/05), we proved the hypothesis that, without tuning, high vertical resolution improves the representation of unresolved, parameterized low altitude clouds such as stratocumulus clouds in global climate models. We evaluated the effectiveness of the newly developed computational framework, Framework for Improvement by Vertical Enhancement (FIVE) by implementing it into E3SM. FIVE aims to improve representation of the target cloud, low level cloud for this project, with use of a locally high vertical resolution for computation of selected processes. FIVE was assessed via direct comparison among E3SM, E3SM-FIVE, and satellite data as well as a newly developed ensemble artificial neural network methodology. We demonstrated that E3SM-FIVE remarkably improves marine stratocumulus representation. To advance the FIVE methodology, we took a first step toward the adaptive vertical grid in the FIVE method by modifying the E3SM single column model. For future application of FIVE to the global storm resolving model with a few km horizontal grid resolution, we developed the vertical-meridional two-dimensional (2D) Hadley circulation model, which allows one to run simulations with much higher vertical resolution and significantly longer duration than current three-dimensional global storm resolving model. During the second phase of the project (2021/05-2023/11; period for this technical report), we sought a way to appropriately transform FIVE to a regionalized version of FIVE, which allows one to select locations with/without FIVE. In parallel, we demonstrated that E3SM-FIVE with a 25-km mesh significantly reduced the bias associated with coastal stratocumulus. Furthermore, we showed that the same improvement was able to be obtained with the use of both the Regional Refined Model (RRM) targeted to one of the stratocumulus regions with 25-km resolution and FIVE. We explored a new implementation strategy of FIVE with our 2D Hadley circulation model. In parallel we applied the 2D Hadley circulation model to study the double intertropical convergence zone bias with storm resolving model resolution and aerosol-cloud interactions within the Hadley circulation with near large eddy simulation resolution. Last, we have been investigating Arctic mixed-phase stratocumulus clouds with large eddy simulations in order to improve their representation in global models. During the course of the entire project, 6 journal articles were published, and 3 manuscripts are currently in preparation. We gave a total of 30 presentations at various domestic and international conferences. In addition we gave 6 seminar talks at universities and a national laboratory. E3SM- FIVE is available as a branch of E3SM.

54 ENVIRONMENTAL SCIENCES↗

Optimization and Experimental Validation of Annular Finned PCM-HX for a Domestic Hot Water Heater Application

The load profile for domestic water heating is time-dependent and can result in high energy demand during peak operating times. Shifting this peak load can have significant environmental and economic impacts. Phase change material (PCM)-based thermal energy storage (TES) is a potentially useful technology for peak load shifting in domestic hot water (DHW) applications thanks to its high latent heat and energy density. In this study, an annular finned-tube PCM-HX design concept was optimized for a load-shifting TES unit to meet the Department of Energy standard for a medium-usage DHW heater using a resistance-capacitance model (RCM) integrated with a Multi-Objective Genetic Algorithm. The optimized design comprised 70 identical annular finned-tube PCM-HX units connected in parallel and utilizing RT62HC as the PCM. A single PCM-HX unit was prototyped and tested in a vertically oriented setup with upward heat transfer fluid (HTF) flow. The hot water supply time was defined based on a cutoff temperature of 51.7°C. The as-designed mass flow rate (1.5 g/s) was tested to assess the performance of the prototyped PCM-HX unit for RCM validation. For the experimental investigation, RTD sensor bundles measured HTF temperature at the PCM-HX inlet and outlet, and a Coriolis flow meter accurately measured the HTF mass flow rate. The simulated discharging power underpredicted the experimental result by about 12%, and the simulated hot water supply time underpredicted the experimental result by approximately 13% for the as-designed mass flow rate (1.5 g/s). The average deviation of the hot water supply temperature between the experimental and RCM results during the complete PCM solidification process was 1.3 K for the as-designed mass flow rate. The overall good agreement between the experimental and RCM results provides confidence that computationally efficient models such as RCM can be utilized for design optimization of PCM-HXs.

42 ENGINEERING↗