Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial transferability”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Efficient generation and amplification of intense vortex and vector laser pulses via strongly-coupled stimulated Brillouin scattering in plasmas

The past decade has seen tremendous progress in the production and utilization of vortex and vector laser pulses. Although both are considered as structured light beams, the vortex lasers have helical phase fronts and phase singularities, while the vector lasers have spatially variable polarization states and polarization singularities. In contrast to the vortex pulses that carry orbital angular momentum (OAM), the vector laser pulses have a complex spin angular momentum (SAM) and OAM coupling. Despite many potential applications enabled by such pulses, the generation of high-power/-intensity vortex and vector beams remains challenging. Here, we demonstrate using theory and three-dimensional simulations that the strongly-coupled stimulated Brillouin scattering (SC-SBS) process in plasmas can be used as a promising amplification technique with up to 65% energy transfer efficiency from the pump beam to the seed beam for both vortex and vector pulses. We also show that SC-SBS is strongly polarization-dependent in plasmas, enabling an all-optical polarization control of the amplified seed beam. Additionally, the interaction of such structured lasers with plasmas leads to various angular momentum couplings and decouplings that produce intense new light structures with controllable OAM and SAM. This scheme paves the way for novel optical devices such as plasma-based amplifiers and light field manipulators.

Physics↗
DOE OSTI2021

Design and engineering challenges of a multi-energy hard x-ray camera for long-pulse profile measurements at WEST tokamak

The WEST tokamak has recently been prepared for long pulse operation with a water-cooled full first wall. Heating is provided by radiofrequency systems, including Lower Hybrid Current Drive (LHCD). The Princeton Plasma Physics Laboratory has developed a multi-energy hard x-ray camera for profile measurements at WEST. The camera is based on a new generation of 2D pixel array detectors that allow the energy threshold to be independently set for each pixel. The diagnostic will provide spatial, temporal and energy resolved measurements of the hard x-ray emission from the full plasma cross-section, investigating several physical quantities such as the electron temperature from continuum emission, the fast electron tail density produced by radiofrequency current drive and runaway electrons, as well as characteristic tungsten x-rays due to beam-target emission at the edge. Overall, this work describes the engineering challenges that the WEST long pulse scenario poses for this diagnostic and how the design addresses and solves them. Vacuum, thermal stress and heat transfer calculations are presented and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Underlying Mechanisms of Hot Carrier-Driven Reactivity on Bimetallic Nanostructures

Bimetallic nanostructures exhibit unique catalytic activity and selectivity that is not evident for their monometallic analogs. Such nanostructures contain plasmonic metals, such as gold or silver, that afford highly efficient harvesting of electromagnetic radiation and its conversion into hot carriers. These highly energetic species are transferred to the catalytic metal sub-component of the bimetallic nanostructure, where a large spectrum of chemical reactions may be catalyzed. The strength of the electric field and the interplay between catalytic and plasmonic metals at the nanoscale are thus critically important for the catalytic activity of bimetallic nanostructures. In this study, we investigate the relationship between the catalytic activity and local electric fields sustained on the surface of gold-palladium (Au@PdNPs) and gold-platinum (Au@PtNPs) nanoplates using tip-enhanced Raman spectroscopy (TERS). We image the spatially varying magnitudes of rectified (DC) local electric fields on the surface of these nanostructures and compare them to fields sustained on the surface of monometallic nanoplates. We find substantially larger electric field magnitudes on Au@PdNPs and Au@PtNPs as compared to their monometallic analogs. These findings suggest that catalytic efficiency of bimetallic nanostructures may be mediated and potentially tuned through precise control of electric fields sustained on their surfaces.

Li, Zhandong↗
DOE OSTI2025

Prolonging All-Optical Molecular Electron Spin Coherence in the Tissue Transparency Window

Coherent electron spin states within paramagnetic molecules hold significant potential for microscopic quantum sensing. However, all-optical coherence measurements amenable to high spatial and temporal resolution under ambient conditions remain a significant challenge. Here we conduct room-temperature, picosecond time-resolved Faraday ellipticity/rotation (TRFE/R) measurements of the electron spin decoherence time T 2 * in [IrBr 6 ] 2- . Decoherence is strongly sensitive to solution phase viscosity, pointing to molecular tumbling as an important decoherence mechanism. Accordingly, immobilization of [IrBr 6 ] 2- molecules in thin polymer films results in an order-of-magnitude increase in coherence lifetime and significantly greater magnetic field sensitivity. Here, by tuning energies of ligand-to-metal charge transfer (LMCT) states, TRFE/R enables spin initialization and readout in the tissue transparency window, paving the way toward all-optical, ultrafast molecular electron spin coherence imaging in biological systems.

Electron paramagnetic resonance spectroscopy↗
DOE OSTI2021

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗
DOE OSTI2023

Measurement of 2D density profiles using a second-harmonic, dispersion interferometer

A second-harmonic, dispersion interferometer is used to image large-area (≃5 cm 2 ) plasma-jet and gas-jet density profiles. Achromatic telescopes magnify the diameters of the primary-laser beam (1064 nm) and its second-harmonic (532 nm) before probing the sample and de-magnify the beam diameters after the sample, where the primary beam transfers its phase change to a second, second-harmonic beam, allowing the sample’s dispersive-phase change to be measured between two, orthogonally polarized second harmonic beams. The telescopes produce an azimuthally symmetric, dispersive-phase shift in the sample + background phase-change image and in the background phase-change image, which is removed by digital subtraction. The interferometer’s performance was verified using standard-optical components as dispersive elements (BK7 lenses and wedge plates), resolving a minimum, phase-change sensitivity of ΔΦ min ≳ 15 mrad and spatial resolution of Δx res ≃ 100 μm. The phase change produced by unknown-density objects (a pulsed-plasma-jet and a pulsed-gas-jet) was measured, and their data were used to recover the original, 2D density profiles using an inverse Abel transform: peak-number density, N gas ≃ 6 × 10 20 cm –3 and N e ≃ 5 × 10 16 cm –3 ; line-integrated density, ∫N gas dl ~ 2 × 10 19 cm –2 and ∫N e dl ~ 1 × 10 16 cm –2 . The techniques and methods developed here are scalable to even larger probe-beam diameters and frame-capture rates, leading to a diagnostic capability that is well-suited for applications involving the real-time measurement of density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Permeability-controlled migration of induced seismicity to deeper depths near Venus in North Texas

Abstract Migration of seismic events to deeper depths along basement faults over time has been observed in the wastewater injection sites, which can be correlated spatially and temporally to the propagation or retardation of pressure fronts and corresponding poroelastic response to given operation history. The seismicity rate model has been suggested as a physical indicator for the potential of earthquake nucleation along faults by quantifying poroelastic response to multiple well operations. Our field-scale model indicates that migrating patterns of 2015–2018 seismicity observed near Venus, TX are likely attributed to spatio-temporal evolution of Coulomb stressing rate constrained by the fault permeability. Even after reducing injection volumes since 2015, pore pressure continues to diffuse and steady transfer of elastic energy to the deep fault zone increases stressing rate consistently that can induce more frequent earthquakes at large distance scales. Sensitivity tests with variation in fault permeability show that (1) slow diffusion along a low-permeability fault limits earthquake nucleation near the injection interval or (2) rapid relaxation of pressure buildup within a high-permeability fault, caused by reducing injection volumes, may mitigate the seismic potential promptly.

58 GEOSCIENCES↗
DOE OSTI2024

Time reflection of light from a quantum perspective and vacuum entanglement

If a boundary between two static media is moving with a constant superluminal velocity, or there is a sudden change of the refractive index with time, this yields generation of entangled pairs of photons out of vacuum propagating in the opposite directions. Here we show that during this process, entanglement of Minkowski vacuum is transferred to the entanglement of the generated photon pairs. If initially an electromagnetic pulse is present in the medium the photon generation is stimulated into the pulse mode, and since photons are created as entangled pairs the counter-propagating photon partners produce a pulse moving in the opposite direction, which is known as time reflection. Thus, time reflection occurs due to stimulated generation of the entangled photon pairs out of entangled vacuum and no photons in the original pulse are in fact being reflected. This is different from the mechanism of light reflection from spatial inhomogeneities for which no photons are generated.

Optics↗
DOE OSTI2022

Coastal Acoustic Buoy for Offshore Wind: Project Synthesis

The population of North Atlantic right whales is critically endangered and their habitat overlaps with offshore windfarm leases. It is therefore imperative that effective mitigation strategies be used to avoid impacts on right whales during the construction of offshore windfarms. The Department of Energy issued FOA Number DE-FOA-0001924 to encourage the development of technology that could monitor large exclusion zones for right whales in order to mitigate potential impact of construction noise on right whales. This report summarizes the past two years of the development and evaluation of the Coastal Acoustic Buoy for Offshore Wind (CABOW) project which aimed to develop technology to monitor large exclusion zones for North Atlantic right whales. Over the course of the project SMRU Consulting have implemented a rigorous design process including comparison of different approaches (e.g., single sensor vs multiple sensors), as well as consideration of placement and timing of acoustic monitoring. We have evaluated critical components of the CABOW system including, reliability, detection range, and bearing accuracy in areas adjacent to offshore windfarm leases in Maryland by conducting 3,536 playbacks of simulated right whale upcalls. The maximum call detection range was 7.5 km when noise was 99 dB re 1µPa rms (50-225 Hz), but this reduced to < 1 km when ambient noise levels were high. Our detection probability in the field was measured as a function of range as well as the source-to-noise level ratio allowing us to build a model to predict the probability of detection under various scenarios (sample size: 3,536 calls x 5 buoys = 17,680). The median bearing error was -0.25° but this is likely an underestimate of error due to experimental design. Using the published recall and precision of the two detectors we implemented in the CABOW system, we estimate that at a recall of 80%, our precision was > 80%, within the range of what we were aiming for in this project. To estimate our exclusion zone false negative and false positive rates, we built a simulation model using the empirical data from our field trial. We modelled three to nine CABOW units placed on the 10 km exclusion zone and estimated our false negative rate to be 1% or less (which was our project goal) and our false positive rate to be between seven and nine percent, slightly above our goal of 5%. However, we also modelled an equivalent PAM system that does not have bearing capabilities and found the false positive rates for that system to be six to eight times higher than the CABOW rate. This higher false positive rate of PAM systems with low spatial information could have significant cost repercussions for offshore wind developers by adding work shutdowns or delays without providing additional protection for right whales. The model we built allows us to explore the placement of PAM systems under various scenarios and will thus help facilitate planning of PAM mitigation systems to meet NOAA Incidental Harassment Authorizations for specific windfarms. We achieved an average system uptime of 98.3%, just below our goal of 99%. The issues that caused these short losses of data have been identified and fixed. Right whale detections and audio clips were typically transferred via radio from the buoys to the base station in two to four seconds. We therefore believe we have developed a highly robust real-time PAM system. Based on the above, we feel we have achieved the stated funding goal of developing a cost-effective and robust real-time PAM system that enables the monitoring of large exclusion zones for right whales during the constructions of offshore windfarms. This should lead to decreased costs and risks for the offshore wind sector while providing robust mitigation for right whales. It is important to state that PAM mitigation will need to be implemented with other mitigation strategies (e.g., visual observers) to provide a complete mitigation strategy to ensure that any effects on right whales from offshore windfarm construction is minimized.

17 WIND ENERGY↗
DOE OSTI2022

Representative-cell-based Modeling of HTTF in SAM

Although the current nuclear power market is primarily occupied by the light water reactors (LWRs), the TRISO-fueled helium-cooled graphite-moderated high temperature gas reactors (HTGRs) are drawing growing attention as the nuclear power industry marches towards more advanced systems. Compared to that of conventional LWRs, the thermal hydraulics of HTGR cores show extra complexity from their multi-scale heat transfer mechanisms with varying importance during different operation modes or transient stages. The modeling of HTGR cores for the system analysis purpose therefore faces the challenge of reaching reasonable fidelity and accuracy while maintaining sufficient simplicity. Previous development and validation efforts have demonstrated that a “2-D ring model” with reasonable performance can be implemented in the System Analysis Module (SAM) for prismatic HTGR cores. In the current project, an alternative modeling methodology based on “representative cells” is proposed for the typical prismatic core of HTGRs. Different from the previous ring model approach, the proposed methodology first separately models and then combines the small- and large-scale thermal conductions, by connecting representative local heat transfer units (cells) with effective core- wise thermal resistance. A modeling practice for an integral high-temperature test facility (HTTF) elaborates the modeling details. Steady-state validation of the resultant candidate model is performed against a higher-resolution benchmark from a 3D-1D coupled simulation, which shows satisfactory prediction with reasonably captured global parameters and well-represented temperature fields. A postulated pressurized conduction cooldown (PCC) is also simulated and analyzed, demonstrating the model’s capability of transient prediction with physically captured phenomena resolved in both small and large temporal and spatial scales. In general, the current work achieves a preliminary success in proposing a methodology using representative cells to model prismatic HTGR cores in SAM. Future efforts are envisioned with more validation activities and with potential modeling improvements to eventually achieve the high confidence on a high-fidelity robust modeling methodology with reasonable accuracy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ligand-Controlled Energetics and Charge Transfer in Pure and Doped Nanocrystals

The research in this award period initially focused on the spectroscopy and dynamics of spherical CdSe nanocrystals and CdSe nanoplatelets. In the later part of the grant period we turned our attention to InP-based nanocrystals. It has long been believed that the transient absorption signal from approximately spherical CdSe nanocrystals is dominated by conduction band state filling. However, this long-held belief has recently been challenged, based on femtosecond absorption measurements. However, the recent studies challenging the conventional wisdom do not account of finite rates of spin-lattice relaxation. Our work showed that this is a crucial error and that the recent results are misinterpreted, that is, the previous prevailing wisdom is correct. Another study focused on CdSe nanocrystal photochemistry used transient absorption (TA) spectroscopy to determine the spatial extents of CdSe nanoplatelet (NPL) excitons. Our work shows that the spatial extents of the excitons in the NPLs are far less than the physical dimension of the NPL. Using a model developed to understand the transient absorption spectroscopy, we obtain an average excitonic area of 21.2 ± 2.5 nm 2 , independent of the nanoplatelet size. Our work on InP/ZnSe and InP/ZnS core/shell nanocrystals shows that when there is a small lattice mismatch (InP-ZnSe, 3.5%) a coherent core-shell interface is obtained. In contrast, the InP-ZnS lattice mismatch is much larger, 8.3%. In this case, the experimental results showed best agreement with calculations in which lattice strain is ignored, indicating that the interfaces in InP/ZnS nanocrystals are largely incoherent.

14 SOLAR ENERGY↗
DOE OSTI2024

Coupled hydro-thermal flow and radionuclide transport driven by spatial variation of heat-generating radioactive wastes in shale formations

Deep geologic disposal of multiple nuclear waste packages with various heat sources can induce nonuniform hydro-thermal behaviors in the near-field of the repository, consequently influencing the long-term radionuclide transport in the far-field once waste form breach initiates. Here, this study looks into three cases with variation in the spatial order of six groups of heat sources (10th, 50th, 75th, 90th, 95th, and 99th percentiles of heat outputs generated from 1,981 as-loaded dual-purpose canisters in the field site) in a shale-hosted repository with respect to the uni-directional groundwater flow (from west to east): (1) cooler waste packages from west to east, (2) hotter waste packages from west to east, and (3) hottest waste packages in the middle of the repository. Our field-scale PFLOTRAN simulation represents heat-driven multiphysics coupled mechanisms, including multiphase flow, heat transfer, and chemical/radioactive transport, and also, calculates the onset of waste form breach based on temperature-dependent canister vitality. The results from this sensitivity study will quantify the short- (less than 1 × 10 3 years) and long-term (up to 1 × 10 6 years) impacts of sporadic heat pulses from waste package on the spatio-temporal perturbation in hydro-thermal flow quantities and the rate of radionuclide transport in both near- and far-field of the repository system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

On the Interfacial Behavior of Catenated Poly(l-lactide) at the Air–Water Interface

Interfacial properties of polymeric materials are significantly influenced by their architectural structures and spatial features, while such a study of topologically interesting macromolecules is rarely reported. Here, we reported, for the first time, the interfacial behavior of catenated poly(l-lactide) (C-PLA) at the air–water interface and compared it with its linear analogue (L-PLA). The isotherms of surface pressure–area per repeating unit showed significant interfacial behavioral differences between the two polymers with different topologies. Isobaric creep experiments and compression–expansion cycles also showed that C-PLA demonstrated higher stability at the air–water interface. Interestingly, when the films at different surface pressures were transferred via the Langmuir–Blodgett method, successive atomic force microscopy imaging displayed distinct nanomorphologies, in which the surface of C-PLA exhibited nanofibrous structures, while that of the L-PLA revealed a smoother topology with less fiber-like structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Two-photon excited deep-red and near-infrared emissive organic co-crystals

Abstract Two-photon excited near-infrared fluorescence materials have garnered considerable attention because of their superior optical penetration, higher spatial resolution, and lower optical scattering compared with other optical materials. Herein, a convenient and efficient supramolecular approach is used to synthesize a two-photon excited near-infrared emissive co-crystalline material. A naphthalenediimide-based triangular macrocycle and coronene form selectively two co-crystals. The triangle-shaped co-crystal emits deep-red fluorescence, while the quadrangle-shaped co-crystal displays deep-red and near-infrared emission centered on 668 nm, which represents a 162 nm red-shift compared with its precursors. Benefiting from intermolecular charge transfer interactions, the two co-crystals possess higher calculated two-photon absorption cross-sections than those of their individual constituents. Their two-photon absorption bands reach into the NIR-II region of the electromagnetic spectrum. The quadrangle-shaped co-crystal constitutes a unique material that exhibits two-photon absorption and near-infrared emission simultaneously. This co-crystallization strategy holds considerable promise for the future design and synthesis of more advanced optical materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Windowed Fourier transform and cross-correlation algorithms for molecular tagging velocimetry

In this study, simulated and experimental molecular tagging velocimetry (MTV) images have been analyzed with a technique commonly used to process grid images on surfaces, the windowed Fourier transform with local spectrum analysis (WFT-LSA). A systematic synthetic image study of the modulation transfer function (MTF) and error tendencies of the WFT-LSA was performed and compared with a PIV-style cross-correlation algorithm to see if advanced strategies such as iterative image deformation can improve analysis of gridded images with high noise levels. Testing of single-pass algorithms showed that in typical MTV images, the WFT-LSA yields significantly lower bias errors than cross-correlation (CC) at displacements greater than 1 pixel but slightly higher random error at all displacements and image conditions. Analysis of the MTF shows that CC provided better resolution of spatial fluctuations than the WFT-LSA in many combinations of grid size and interrogation window. Tests of image deformation algorithms showed that the gap in performance between CC and WFT-LSA is maintained even as both methods improve. Additionally, WFT-LSA and CC methods are applied to real MTV experiments in high-speed gas jet flows. A preliminary analysis of phosphorescence lifetime provided by acetone vapor excited at 266 nm is used for assessing the required gas speeds for making MTV application feasible. The application of WFT-LSA to real MTV images demonstrates the ability of the algorithm to handle further real-world effects that could not be considered in the synthetic image analysis, like reduced signal-to-noise ratio and non-uniform intensity of the tagging grid across the image introduced by the actual laser beam energy distribution. With experimental images, CC is more accurate with shear flows but less robust to high noise levels than WFT-LSA, as predicted by the synthetic image analysis.

42 ENGINEERING↗
DOE OSTI2021

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↗