Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Parametric Decay of Alfvénic Wave Packets in Nonperiodic Low-beta Plasmas

The parametric decay of finite-size Alfvén waves in nonperiodic low-beta plasmas is investigated using one-dimensional (1D) hybrid simulations. Compared with the usual small periodic system, a wave packet in a large system under the absorption boundary condition shows different decay dynamics, including reduced energy transfer, localized density cavitation, and ion heating. The resulting Alfvén wave dynamics are influenced by several factors relating to this instability, including the growth rate, central wave frequency, and unstable bandwidth. A final steady state of the wave packet may be achieved when the instability does not have enough time to develop within the residual packet, and the packet size shows well-defined scaling dependencies on the growth rate, wave amplitude, and plasma beta. Under the proper conditions, enhanced secondary decay can also be excited in the form of a narrow, amplified wave packet. These results may help to interpret laboratory and spacecraft observations of Alfvén waves, and to refine our understanding of the associated energy transport and ion heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reduced-Order Modeling Method for Phase-Change Thermal Energy Storage Heat Exchangers

Thermal energy storage can facilitate the effective utilization of renewable energy. To speed up the design process of thermal energy storage devices, it is critical to develop fast and accurate modeling methods for phase change material embedded heat exchangers (PCM HXs). This study developed and compared two approximation-assisted reduced-order PCM HX models for the simulation of thermal storage components and systems, which were verified against a validated finite-volume model. They are a pure black-box model, and a grey-box model based on the Number of Transfer Units (or effectiveness-NTU) approach. We used the reduced-order models to predict the performance of standalone PCM HX for 1000 cases with different fluid inlet conditions and HX designs. We also integrated the models with a vapor compression system to predict the compressor energy consumption and total charging time under various conditions. The results show that overall, the black-box model gave more accurate results than the grey-box model. On average, the mean absolute deviation in the PCM HX fluid outlet temperature was 0.05 K and 0.1 K for the black-box model and grey-box model, respectively. The grey-box model gave larger temperature deviation toward the end of the phase change process, due to the simplified two-node PCM representation. On the system level, the mean absolute deviation in compressor energy consumption was 0.2% and 0.3% for the black-box model and grey-box model, respectively. Their mean absolute deviation in total charging time was 1.1% and 2.6%, respectively. In terms of computation efficiency, the system simulation speedup ratio gained by using the reduced-order models was 11 to 57. System simulation time was decreased from an average of 1465 s to an average of 59 s. This shows that the proposed reduced-order modeling methods can be used to predict system performance of PCM HXs with less than 3% accuracy penalty, and 25 times less computational time than finite-volume models, enabling faster design and evaluation of PCM thermal storage devices.

25 ENERGY STORAGE↗
DOE OSTI2023

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Large-eddy simulations of turbulent flows in arrays of helical- and straight-bladed vertical-axis wind turbines

Effects of helical-shaped blades on the flow characteristics and power production of finite-length wind farms composed of vertical-axis wind turbines (VAWTs) are studied numerically using large-eddy simulation (LES). Two helical-bladed VAWTs (with opposite blade twist angles) are studied against one straight-bladed VAWT in different array configurations with coarse, intermediate, and tight spacings. Statistical analysis of the LES data shows that the helical-bladed VAWTs can improve the mean power production in the fully developed region of the array by about 4.94%–7.33% compared with the corresponding straight-bladed VAWT cases. The helical-bladed VAWTs also cover the azimuth angle more smoothly during the rotation, resulting in about 47.6%–60.1% reduction in the temporal fluctuation of the VAWT power output. Using the helical-bladed VAWTs also reduces the fatigue load on the structure by significantly reducing the spanwise bending moment (relative to the bottom base), which may improve the longevity of the VAWT system to reduce the long-term maintenance cost.

17 WIND ENERGY↗
DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗
DOE OSTI2024

Nanoindentation Stress Relaxation to Quantify Dislocation Velocity–Stress Exponent

This work reports a new methodology using indentation stress relaxation to characterize the dislocation velocity–stress exponent. Through the indentation stress relaxation process, the dislocation structure builds up at the rate governed by dislocation velocity, which is a function of the externally applied stress. The relationship between the dislocation velocity and stress can thus be derived from the indentation stress relaxation data of the stress as a function of time. In this study, instrumented nanoindentation stress relaxation experiments were performed on pure aluminum samples, following three different initial displacement rates of 100, 400, and 800 nm/s. Based on the scaling properties of dislocation kinetics, the data were interpreted to derive a dislocation velocity–stress exponent of 2.5 ± 0.5 for room-temperature aluminum. Crystal plasticity finite-element simulations were performed to illustrate the sensitivity of the proposed nanoindentation stress relaxation methodology to the dislocation velocity–stress exponent value.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Experimentally Validated Computational Modeling of Creep and Creep-Cracking for Nuclear Concrete Structures

In a Nuclear Power Plant, one of the most important components is the concrete nuclear reactor cavity, which serves both a structural and protective function as the biological radiation shield. Given that creep has been identified as a major knowledge gap in the assessment of nuclear structures (NUREG/CR-7153), this work helps to further the understanding of creep behavior of massive concrete containment structures for decades to enable safe and long-term operation of these facilities. This project has developed a robust, experimentally validated model to predict creep in nuclear concrete structures for up 60 years using short-term creep data thereby enabling a longer service life of critical facilities and early detection of structural failure. The work presented in this report is a pairing of computational and experimental methods. For the first time, the time temperature superposition (TTS) principle was successfully used to generate a uniaxial creep compliance master curve to predict mortar creep response for up to 22,500 days (nearly 60 years) at a reference temperature of 20°C. These data were used as input into finite element analysis (FEA) codes that use highly realistic random, 3D concrete microstructures from reconstructed coarse limestone aggregates. Finite element analysis performed provides the ability to quickly upscale mortar viscoelastic behavior to long-term concrete creep/relaxation data. A master creep compliance curve, constructed from the TTS principle, spanning 27 years, was used to validate two and a half decades of simulated concrete creep. Concurrently, three different simulated wall specimens were designed to mimic the behavior of post-tensioned concrete nuclear containment facility vessel walls over time as a result of concrete creep. The specimens were designed with different thicknesses, transverse and longitudinal reinforcement ratios, and level of post-tensioning stress. Each specimen contained various instrumentation to measure internal concrete temperature, concrete strain, and post-tensioning strain hourly for over 3 years. The concrete creep model developed in this project, based on the FEA concrete simulations, was applied to simulate the structural-scale experiments of prestressed concrete walls conducted in this project using the Grizzly code. These models can represent the effects of reinforcing and prestressing. Although there are some discrepancies with the experimental data, the model can predict the overall trends of the creep response in these experiments. One of these experimental models was also applied to an extended time to demonstrate how the findings from this study can be used to predict the behavior of actual structures of interest that have been in service for extended periods of time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Effects of finite sampling on fatigue damage estimation of wind turbine components: A statistical study

The variability of the wind turbine loads complicates fatigue assessment in the design phase, as performing simulations covering the entire lifetime is computationally expensive. The current work provides important information for assessing the uncertainty in fatigue damage estimation due to finite data. Here we study the sample size effect on mean, variance, and skewness of damage in each wind bin, identify the important wind bins, and study the uncertainty propagation from each wind bin to the lifetime damage using 3600 aeroelastic simulations and bootstrapping. To achieve less than 1% error in the damage estimation across all load channels in the current case study, at least 100 turbulence seeds are needed. Damage in different wind bins follows a lognormal distribution when using the conventional approach of six seeds. The provided insights and information allow the designer to achieve a specific level of accuracy for a given computational cost using strategic bin sampling.

17 WIND ENERGY↗
DOE OSTI2023

Stress and Curvature Effects in Layered 2D Ferroelectric CuInP 2 S 6

Nanoscale ferroelectric 2D materials offer the opportunity to investigate curvature and strain effects on materials functionalities. Among these, CuInP 2 S 6 (CIPS) has attracted tremendous research interest in recent years due to combination of room temperature ferroelectricity, scalability to a few layers thickness, and ferrielectric properties due to coexistence of 2 polar sublattices. Here, we explore the local curvature and strain effect on polarization in CIPS via piezoresponse force microscopy and spectroscopy. To explain the observed behaviors and decouple the curvature and strain effects in 2D CIPS, we introduce the finite element Landau–Ginzburg–Devonshire model, revealing strong changes in hysteresis characteristics in regions subjected to tensile and compressive strain. The piezoresponse force microscopy (PFM) results show that bending induces ferrielectric domains in CIPS, and the polarization-voltage hysteresis loops differ in bending and nonbending regions. Furthermore, these studies offer insights into the fabrication of curvature-engineered nanoelectronic devices.

2D materials↗
DOE OSTI2021

Reentrant superconductivity through a quantum Lifshitz transition in twisted trilayer graphene

A series of recent experiments has demonstrated robust superconductivity in magic-angle twisted trilayer graphene (TTG). In particular, a recent work by Cao et al. [Nature (London) 595, 526 (2021)] studies the behavior of the superconductor in an in-plane magnetic field and an out-of-plane displacement field, finding that the superconductor is unlikely to have purely spin-singlet pairing. Herein, this work also finds that at high magnetic fields and a smaller range of dopings and displacement fields, the superconductor undergoes a transition to a distinct field-induced superconducting state. Inspired by these results, we develop an understanding of the superconductivity in TTG using a combination of phenomenological reasoning and microscopic theory. We describe the role that an in-plane field plays in TTG, and we use this understanding to argue that the reentrant transition may be associated with a quantum Lifshitz phase transition, with the high-field phase possessing finite-momentum pairing. We argue that the superconductor is likely to involve a superposition of singlet and triplet pairing, and we describe the structure of the normal state. We also draw lessons for twisted bilayer graphene (TBG), and we explain the differences in the phenomenology with TTG despite their close microscopic relationship. We propose that a singlet-triplet superposition is realized in the TBG superconductor as well, and that the ν = –2 correlated insulator may be a time-reversal protected Z 2 topological insulator obtained through spontaneous spin symmetry breaking.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Many-body Green's function approaches to the doped Fröhlich solid: Exact solutions and anomalous mass enhancement

In polar semiconductors and insulators, the Fröhlich interaction between electrons and long-wavelength longitudinal optical phonons induces a many-body renormalization of the carrier effective masses and the appearance of characteristic phonon sidebands in the spectral function, commonly dubbed “polaron satellites.” The simplest model that captures these effects is the Fröhlich model, whereby electrons in a parabolic band interact with a dispersionless longitudinal optical phonon. The Fröhlich model has been employed in a number of seminal papers, from early perturbation-theory approaches to modern diagrammatic Monte Carlo calculations. One limitation of this model is that it focuses on undoped systems, thus ignoring carrier screening and Pauli blocking effects that are present in real experiments on doped samples. Here, to overcome this limitation, we extend the Fröhlich model to the case of doped systems, and we provide exact solutions for the electron spectral function, mass enhancement, and polaron satellites. We perform the analysis using two approaches, namely, Dyson’s equation with the Fan-Migdal self-energy, and the second-order cumulant expansion. We find that these two approaches provide qualitatively different results. In particular, Dyson’s approach yields better quasiparticle masses and worse satellites, while the cumulant approach provides better satellite structures, at the price of worse quasiparticle masses. Both approaches yield an anomalous enhancement of the electron effective mass at finite doping levels, which in turn leads to a breakdown of the quasiparticle picture in a significant portion of the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Many-Body Signatures of Collective Decay in Atomic Chains

Fully inverted atoms placed at exactly the same location synchronize as they deexcite, and light is emitted in a burst (known as “Dicke’s superradiance”). We investigate the role of finite interatomic separation on correlated decay in mesoscopic chains and provide an understanding in terms of collective jump operators. We show that the superradiant burst survives at small distances, despite Hamiltonian dipole-dipole interactions. However, for larger separations, competition between different jump operators leads to dephasing, suppressing superradiance. Collective effects are still significant for arrays with lattice constants of the order of a wavelength, and lead to a photon emission rate that decays nonexponentially in time. We calculate the two-photon correlation function and demonstrate that emission is correlated and directional, as well as sensitive to small changes in the interatomic distance. Furthermore, these features can be measured in current experimental setups, and are robust to realistic imperfections.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Simulation for the fault current limiting operation of REBCO CORC superconducting cables with different core materials

With high power density and low loss, CORC superconducting cables composed of RE- Ba 2 Cu 3 O 7–δ (REBCO) coated conductors are of great interest for power transmission applications. They can also effectively be designed to serve as fault current limiters, thanks to their sharp superconducting-to-normal transition. Coupled electromagnetic-thermal finite element (FE) simulations implemented in COMSOL Multiphysics were developed and validated against the published results for the fault current limiting performance of two CORC cables of REBCO coasted conductors wound on either copper or stainless steel cores. For improved accuracy, temperature dependencies of electrical and thermal properties of all component materials were considered in the simulations. The simulations were performed in the cross-sections of the cables for every individual layer of each REBCO tape to deliver comprehensive understanding of evolution of the current distribution and temperature rise in those layers. Furthermore, the studies can suggest approaches to optimize cable design for fault current limiting (FCL) applications by assessing the role of individual components. In the paper, possible computational errors, challenges and improvements are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effects of pillar design on the thermal performance of vacuum-insulated glazing

Vacuum-insulated glazing (VIG) with a low-emittance coating has a great market potential as an effective transparent insulator. The thermal insulating performance of VIG is determined by its design, including material selection and configuration of different components. Thermal conductance of the vacuum gap as a transport bottleneck is one of the primary factors controlling the thermal transport across VIG. In particular, because support pillars provide the main thermal transport channels across the vacuum gap, increasing the pillar thermal resistance is a key strategy for creating effective thermal insulation while maintaining the vacuum space. Here, the effects of various pillar design parameters, such as thermal conductivity, geometry, and arrangement, on the VIG thermal performance were comprehensively investigated via the finite element method. In addition, analytical models for thermal transport were examined and thermal conductance across the VIG unit was experimentally measured for validation. The pillar design parameters, especially the height, shape, spacing, and arrangement of the pillars, showed significant effects on the thermal performance of VIG. This research also shows that the smaller contact area for horizontal pillars can effectively decrease the heat loss by more than 30%. Because current VIG analytical equations of thermal performance are only applicable to cylindrical pillars, an analytical equation that can better describe the thermal performance of rectangular parallelepiped pillars is presented, along with a discussion about the mechanism of thermal transfer for different pillar shapes. Thermo-mechanical analyses based on 3D FEM simulations can provide valuable insights into the effect of various design parameters on the overall performance the VIG, allowing for the development of an optimal VIG design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗