Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

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,567 records · Page 87

Hybrid-Electric Aero-Propulsion Controls Laboratory: Overview and Capability

A hardware-in-the-loop (HIL) laboratory is developed to investigate control technologies for electrified aircraft propulsion (EAP). The laboratory emulates a propulsion system by reproducing the mechanical shaft interface to the electrical power system in hardware. The experimental electric power system includes supercapacitor energy storage and a dynamically variable electrical load. A novel method of scaling power and inertia is provided in software to accurately reproduce the transient, off-design turbomachinery performance dynamics without including actual turbomachinery. An overview of how the HIL system can accommodate a broad range of EAP architectures, including power extraction and insertion, with capability for transient energy management is described. The real-time system operates in the 100kW power class and is instrumented to emulate turbomachinery – power system interactions. The platform is an agile, flexible laboratory for low-cost, risk-reduction development and testing of propulsion control, operability, and energy management technologies.

controls↗

High Speed Data Investigation of an LOX/LNG Combustor Assisted by Non-linear Acoustic Simulations

During testing of a regenerative cooled liquid oxygen/liquid natural gas combustor large amplitude combustion instabilities were observed in half of the series conducted. These instabilities produced an unsteady pressure Root Mean Square (RMS) value greater than the chamber steady state pressure. Quick look investigations of the test data presented some standard indications of non-linear tangential instability such as a non-sinusoidal time series, and larger peaks than trough. Also observed in the data quick look were some atypical findings including a first harmonic amplitude larger than the fundamental, and numerous side bands on the fundamental and subsequent harmonics. A simplified model of the combustor was created using a commercially available finite element software. The first tangential acoustic mode was then simulated using a non-linear wave equation in conjunction with a time explicit discontinuous Gelarkin flow solver. Numerous simulations were evaluated providing valuable insight and reproducing the atypical findings in the data. Side bands were created through non-linear signal modulation like an AM radio. This happens due to the non-linear steepening and relaxing that occurs in large amplitude oscillations. The simulation was also able to reproduce a signal with a first harmonic amplitude larger than the fundamental. Investigating the pressure signal at multiple locations across the chamber head end revealed radial dependent frequency content. A literature review of the topic revealed corroboration of these effect in idealized geometries.

Combustion Stability↗

Improving “Domain-Relevant Metadata Requirements” for Supporting Open-Source Science Initiative

Implementation of the NASA Open-Source Science Initiative (OSSI) requires sharing of all relevant information to ensure “open reproducible science” [1]. However, there are several challenges in applying the OSSI to airborne field campaigns focused on atmospheric composition, which often involve a wide variety of in-situ measurements for trace gases, aerosol and cloud properties, meteorological parameters, and radiation fields. To ensure open reproducibility from airborne field campaigns, it is essential to obtain detailed measurement descriptions, which include the detection principle, sample procedure and treatment, and data processing and correction method. The challenge is that some information, e.g., sampling procedure and treatment, may be instrument-specific and campaign or platform-dependent. The data processing may also involve empirical corrections which may evolve over time. In addition, these details (especially operation- or campaign-specific ones) are often not given in journal publications. Given these issues, there is a need to leverage and improve the current “domain-relevant metadata requirements” to represent the measurement description in standardized metadata. These requirements can then facilitate systematic collection of measurement specific metadata and serve as a foundation to develop tools for making the information accessible and data more interoperable and usable or reusable. Here we show a review of existing metadata collections, use cases, and needs for new standards.

Sean Leavor↗

Developing a Vision for Maturing the Heliophysics Infrastructure towards Open Science

In the dawn of open science and the upcoming requirements, we speak about the existing state of Heliophysics infrastructure and detail the evolution required to address capability or interconnection shortcomings. Such a daunting barrier calls for an analysis ecosystem with multi-faceted capability. We propose such an ecosystem, called DIARieS, to be built upon five conceptual pillars: Discovery, Implementation, Analysis, Reproducibility, and Sharing of results. The combination of these concepts in a single platform will enable users to more intuitively combine recent advances in technology to create ‘DIARieS’ of their workflows, which can be easily made open to others in the community. The DIARieS ecosystem will also increase our efficiency by streamlining our various workflow processes, including automatic incorporation of the impending requirements of open science. The various components of the ecosystem will simplify software installation and data implementation, including automatically generated citation lists based on the components included. Automatic containerization and version control of the ecosystem will make the custom workflows easily reproducible. Employing widget technology will ease the difficulty of producing publication and commercial quality visualizations and applying common analyses techniques. Incorporating multiple technologies will streamline the various sharing methods common in our work environments today. Overall, the totality of capabilities to be offered by this analysis ecosystem will drastically simplify the application of open science principles to our work in addition to improving our efficiency and ease of collaboration.

Infrastructure↗

Interpreting Solar Wind Turbulent Spectra beyond Taylor’s Hypothesis

In this Letter we apply a methodology, recently proposed by Bourouaine & Perez (BP19), to interpret solar wind turbulent power spectra beyond the Taylor approximation (TA). The turbulent power spectra were measured using Helios spacecraft data near 0.6 au. We use the model proposed in BP19 to reproduce the field-perpendicular power spectrum E(k(sub ⊥)) of antisunward Alfvénic fluctuations in the plasma frame (where k(sub ⊥) is the field-perpendicular wavenumber) from the corresponding measured frequency power spectrum P(sub sc)(ω,θ(sub b)) along the sampling angle ,θ(sub b), which is the angle between the local magnetic field and the sampling direction. Here ω = 2πf and f is the frequency of the time signal. Interestingly enough, we found that for all corresponding measured frequency power spectrum P(sub sc)(ω,θ(sub b)) the reproduced field-perpendicular power spectrum E(k(sub ⊥)) is the same and independent of the considered sampling angle θ(sub b). This finding is consistent with the fact that the analyzed turbulence is strong and highly anisotropic with k(sub ∥) ≪ k(sub ⊥) (where k(sub ∥) is the field-parallel wavenumber). Furthermore, for this specific time signal we found that the commonly used TA is still approximately valid with the important difference that a broadening in k(sub ⊥) for each angular frequency ω is present. This broadening can be described in the context of the methodology proposed in BP19.

Solar magnetic fields↗

Ensemble Simulations of the 2012 July 12 Coronal Mass Ejection with the Constant-turn Flux Rope Model

Flux-rope-based magnetohydrodynamic modeling of coronal mass ejections (CMEs) is a promising tool for prediction of the CME arrival time and magnetic field at Earth. In this work, we introduce a constant-turn flux rope model and use it to simulate the 2012 July 12 16:48 CME in the inner heliosphere. We constrain the initial parameters of this CME using the graduated cylindrical shell (GCS) model and the reconnected flux in posteruption arcades. We correctly reproduce all the magnetic field components of the CME at Earth, with an arrival time error of approximately 1 hr. We further estimate the average subjective uncertainties in the GCS fittings by comparing the GCS parameters of 56 CMEs reported in multiple studies and catalogs. We determined that the GCS estimates of the CME latitude, longitude, tilt, and speed have average uncertainties of 5°.74, 11°.23, 24°.71, and 11.4%, respectively. Using these, we have created 77 ensemble members for the 2012 July 12 CME. We found that 55% of our ensemble members correctly reproduce the sign of the magnetic field components at Earth. We also determined that the uncertainties in GCS fitting can widen the CME arrival time prediction window to about 12 hr for the 2012 July 12 CME. On investigating the forecast accuracy introduced by the uncertainties in individual GCS parameters, we conclude that the half-angle and aspect ratio have little impact on the predicted magnetic field of the 2012 July 12 CME, whereas the uncertainties in longitude and tilt can introduce relatively large spread in the magnetic field predicted at Earth.

Talwinder Singh↗

Earth-System-Model Evaluation of Cloud and Precipitation Occurrence for Supercooled and Warm Clouds Over the Southern Ocean's Macquarie Island

Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

cloud feedbacks↗

Cold Traps of Hypervolatiles in the Protosolar Nebula at the Origin of the Peculiar Composition of Comet C/2016 R2 (PanSTARRS)

Recent observations of the long-period comet C/2016 R2 (PanSTARRS; hereafter R2) indicate an unusually high N2/CO abundance ratio, typically larger than ∼0.05, and at least 2–3 times higher than the one measured in 67P/Churyumov–Gerasimenko. Another striking compositional feature of this comet is its heavy depletion in H2O (H2O/CO ∼ 0.32%), compared to other comets. Here we investigate the formation circumstances of a generic comet whose composition reproduces these two key features. We first envisage the possibility that this comet agglomerated from clathrates, but we find that such a scenario does not explain the observed low water abundance. We then alternatively investigate the possibility that the building blocks of R2 agglomerated from grains and pebbles made of pure condensates via the use of a disk model describing the radial transport of volatiles. We show that N2/CO ratios reproducing the value estimated in this comet can be found in grains condensed in the vicinity of the CO and N2 ice lines. Moreover, high CO/H2O ratios (>100 times the initial gas-phase value) can be found in grains condensed in the vicinity of the CO ice line. If the building blocks of a comet assembled from such grains, they should present N2/CO and CO/H2O ratios consistent with the measurements made in R2's coma. Our scenario indicates that R2 formed in a colder environment than the other comets that share more usual compositions. Our model also explains the unusual composition of the interstellar comet 2l/Borisov.

Olivier Mousis↗

The Hydrological Cycle and Ocean Circulation of the Maritime Continent in the Mid-Pliocene: Results From PlioMIP2

The Maritime Continent (MC) forms the western boundary of the tropical Pacific Ocean, and relatively small changes in this region can impact the climate locally and remotely. In the mid-Piacenzian warm period of the Pliocene (mPWP; 3.264 to 3.025 Ma) atmospheric CO2 concentrations were ∼ 400 ppm, and the subaerial Sunda and Sahul shelves made the land–sea distribution of the MC different to today. Topographic changes and elevated levels of CO2, combined with other forcings, are therefore expected to have driven a substantial climate signal in the MC region at this time. By using the results from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), we study the mean climatic features of the MC in the mPWP and changes in Indonesian Throughflow (ITF) with respect to the preindustrial. Results show a warmer and wetter mPWP climate of the MC and lower sea surface salinity in the surrounding ocean compared with the preindustrial. Furthermore, we quantify the volume transfer through the ITF; although the ITF may be expected to be hindered by the subaerial shelves, 10 out of 15 models show an increased volume transport compared with the preindustrial. In order to avoid undue influence from closely related models that are present in the PlioMIP2 ensemble, we introduce a new metric, the multi-cluster mean (MCM), which is based on cluster analysis of the individual models. We study the effect that the choice of MCM versus the more traditional analysis of multi-model mean (MMM) and individual models has on the discrepancy between model results and data. We find that models, which reproduce modern MC climate well, are not always good at simulating the mPWP climate anomaly of the MC. By comparing with individual models, the MMM and MCM reproduce the preindustrial sea surface temperature (SST) of the reanalysis better than most individual models and produce less discrepancy with reconstructed sea surface temperature anomalies (SSTA) than most individual models in the MC. In addition, the clusters reveal spatial signals that are not captured by the MMM, so that the MCM provides us with a new way to explore the results from model ensembles that include similar models.

Hydrological cycle↗

The Statistical Emulators of GGCMI Phase 2: Responses of Year-to-Year Variation of Crop Yield to CO 2 , Temperature, Water and Nitrogen Perturbations

Understanding the impact of climate change on year-to-year variation of crop yield is critical to global food stability and security. While crop model emulators are believed to be lightweight tools to replace the models, few emulators have been developed to capture such interannual variation of crop yield in response to climate variability. In this study, we developed a statistical emulator with a machine learning algorithm to reproduce the response of year-to-year variation of four crop yields to CO 2 (C), temperature (T), water (W), and nitrogen (N) perturbations defined in the Global Gridded Crop Model Intercomparison Project (GGCMI) phase 2. The emulators were able to explain more than 52 % of the variance of simulated yield and performed well in capturing the year-to-year variation of global average and gridded crop yield over current croplands in the baseline. With the changes in CO 2 –temperature–water–nitrogen (CTWN) perturbations, the emulators could reproduce the year-to-year variation of crop yield well over most current cropland. The variation of R and the mean absolute error was small under the single CTWN perturbations and dual-factor perturbations. These emulators thus provide statistical response surfaces of yield, including both its mean and interannual variability, to climate factors. They could facilitate spatiotemporal downscaling of crop model simulation, projecting the changes in crop yield variability in the future and serving as a lightweight tool for multi-model ensemble simulation. The emulators enhanced the flexibility of crop yield estimates and expanded the application of large-ensemble simulations of crop yield under climate change.

Climate change↗

Bridging the Gap: Enhancing Prominence and Provenance of NASA Datasets in Research Publications

Attribution of datasets that were used to generate research results described in peer-reviewed publications to the original source of these datasets (which are often archived at NASA Earth Science data centers) has been very challenging. Even though the data citation standard of citing datasets as research artifacts and citing them with Digital Object Identifiers (DOIs) was introduced over a decade ago, most authors do not properly reference the data used in their studies and merely mention them in the text. The lack of proper citations of datasets makes the peer-reviewed publication less transparent, imperils reproducibility, and impedes open science. We offer an open-source publication management methodology and a tool that can help to enhance usage-based data discovery, prominence, and provenance of the data; reproducibility of the research results; and potentially increase the return on investment on NASA-funded research.

open-source↗

Boundary-Layer Cloud Modeling Challenges on the North Slope of Alaska

The accurate modeling and prediction of cloud base heights is critical for energy balance calculations and aviation operations, alike. Low-level (i.e., boundary-layer) Arctic clouds can be difficult to model, making prediction of formation and dissipation challenging. Primarily mixed-phase, these clouds typically contain low quantities of supercooled liquid water and often slowly precipitate relatively small amounts of moderately and heavily rimed snow particles. While this appears to be the predominant cloudy state on the North Slope of Alaska (NSA), the delicate balance of microphysical, dynamical, radiative, surface coupling, and advective processes can rapidly shift to heavy snow (with various degrees of riming) or to a complete dissipation of the cloud layer without any precipitation, depending on the dominant processes. Here we strive to disentangle these various processes. First, we compare the predictive performances of four different numerical weather models in forecasting the presence and base-heights of low-level clouds: the High-Resolution Rapid Refresh - Alaska (HRRR-AK) model, the Polar Weather Research and Forecasting (Polar WRF) model, the Unified Model (UM), and the European Centre for Medium-range Weather Forecasting (ECMWF) model. Initial results comparing model output at two U.S. Department of Energy Atmospheric Radiation Measurement (AMT) NSA sites, during the fall season in 2019 and 2022, show that the UM slightly outperforms the HRRR-AK in terms of accurately forecasting the presence of a low-level cloud layer (89% of the time). All models have a significant bias of 300 to 800 meters in forecasting cloud base height (lower than is observed); however, the UM and ECMWF models have the lowest biases. Finally, a case study for a particularly challenging April 2017 thin-cloud event is presented, wherein we compare the performance of four different bulk microphysical parameterization schemes using a higher-resolution large eddy simulation (LES) model, the WRF-LES. Initial results show that the Thompson scheme was the only one able to reproduce and sustain a substantial supercooled liquid layer, but it was unable to reproduce the transition from a deep, liquid-rich cloud to a thin layer with moderately and heavily rimed precipitation. This is the first step in linking simulated LES-scale riming processes with those parameterized at a coarser mesoscale model scale. This has important implications for forecasting low-level clouds in an operational environment, given the efficiency of the riming process.

cloud base heights↗

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

Solar energetic particles↗

Azimuthal Size Scales of Solar Wind Periodic Density Structures

Periodic Density Structures (PDSs) are quasi-periodic variations of solar wind density ranging from a few minutes to a few hours. PDSs advect with the solar wind and have radial length scales (Lx) of tens to several thousand Mm, thus belonging to the class of “mesoscale structures”. Current interplanetary multi-spacecraft observations are not at spatial separations capable of directly measuring the 3D size scale of PDSs or other mesoscale structures. Instead, previous investigations estimated characteristic spatial scales in solar wind parameters using cross-correlation and/or coherence analysis applied to multi-spacecraft observations. For solar wind density and interplanetary magnetic field (IMF) intensity, the reported size scales perpendicular to the Sun-Earth line (Ly) ranged between ≈30 and ≈200 Earth Radii (RE). Here, we implemented a similar approach on the same parameters but focused on high density slow solar wind intervals with PDSs observed by the Wind and ARTEMIS-P1 spacecraft. Additionally, this is the first statistical study of the IMF intensity periodicities in relation to PDSs. We identified intervals in which the two spacecraft observed the same periodicity obtaining two PDS groups based on their radial length scale: Lx1 ≈86RE and Lx2 ≈35RE. Then, we classified the events based on the periodic variations’ coherence level. Reproducing the results with simulations of the PDSs’ transit, we inferred the Ly order of magnitudes for the two PDS groups: Ly1 ≈340RE and Ly2 ≈187RE. Knowing the PDSs’ size scales is fundamental to constrain models aimed at reproducing these structures and is critical to better understand the PDS-magnetosphere coupling.

Interplanetary physics↗

Identification of Sea Surface Temperature and Sea Surface Salinity Fronts along the California Coast: Application Using Saildrone and Satellite Derived Products

Coastal upwelling regions are one of the most dynamic areas of the world’s oceans. The California and Baja California Coasts are impacted by both coastal upwelling and the California Current, leading to frontal activity that is captured by gradients in both Sea Surface Temperature (SST) and Sea Surface Salinity (SSS). Satellite data are a great source of spatial data to study fronts. However, biases near coastal areas and coarse resolutions can impair its usefulness in upwelling areas. In this work gradients in SST from NASA Multi-Scale Ultra-High Resolution (MUR) and in two SSS products derived from the Soil Moisture Active Passive (SMAP) NASA mission are compared directly with gradients derived from the Saildrone uncrewed vehicles to validate the gradients as well as to assess their ability to detect known frontal features. The three remotely sensed data sets (MURSST/JPL, SMAP SSS/RSS, SMAP SSS) were co-located with the Saildrone data prior to the calculation of the gradients. Wavelet analysis is used to determine how well the satellite derived SST and SSS products are reproducing the Saildrone derived gradients. Overall results indicate the remote sensing products are reproducing features of known areas of coastal upwelling. Differences between the SST and SSS gradients are mainly associated with the limitations of the microwave derived SSS coverage near land and its reduced spatial resolution. The results are promising for using remote sensing data sets to monitor frontal structure along the California Coast and the application to long term changes in coastal upwelling and dynamics.

Sea Surface Salinity↗

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

Systems and methods for synchronizing converter modules

Various examples are directed to a converter system comprising first and second series-connected converter modules and a synchronization circuit. The synchronization circuit may modulate a reference signal onto a carrier signal to generate a synchronization current signal and the synchronization current signal to an output current of the converter system to generate an aggregated output current. A first converter module may receive the aggregated output current from a first current sensor and generate a first reproduced synchronization signal at least in part from the aggregated output current. A first switch control signal for switching at least one switch at the first converter may be generated based at least in part on the first reproduced synchronization signal.

Cox, Robert Williams↗

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗