Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗
DOE OSTI2020

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Thermal Optimization of a Silicon Carbide, Half-Bridge Power Module

This project describes the modeling process to design the packaging and heat exchanger for a half-bridge wide-bandgap (WBG) power semiconductor module. The module uses two silicon carbide, metal-oxide-semiconductor field-effect transistor (MOSFET) devices per switch position that are soldered to an aluminum nitride, direct-bond copper (DBC) substrate. A baseplate cooling configuration (e.g., no thermal grease) is used along with a water-ethylene glycol, jet-impingement-style heat exchanger. The heat exchanger was designed to be fabricated using prototyping equipment from the National Renewable Energy Laboratory, complies with automotive standards (for minimal channel sizes, flow rates, and coolant), and considers reliability aspects (i.e., erosion/corrosion). Device-scale computational fluid dynamics (CFD) is used first to design the slot jet impingement cooling configuration and compute the effective heat transfer coefficient (HTC) of the concept. The computed HTCs are then used as boundary conditions for a finite element study to optimize the package geometry (e.g., device layout and baseplate thickness) to minimize thermal resistance and minimize temperature variation between the module's four devices. Finally, a fluid manifold is designed to generate the slot jets and cool the devices. Module-scale CFD predicts a relatively low junction-to-fluid thermal resistance of 16.7 mm2 K/W, a 1.4 degrees C temperature variation between devices, and a total pressure drop of 5,860 Pa (0.85 psi) for the design. The thermal resistance of the module design is about 67% lower than the 2015 BMW i3 power electronics/modules thermal resistance.

DIRECT ENERGY CONVERSION↗
DOE OSTI2021

TRANSIENT THERMAL MODELING OF THE HIGH BURNUP DEMONSTRATION RESEARCH PROJECT CASK USING STAR-CCM+ AND COBRA-SFS

The Department of Energy in collaboration with the Electric Power Research Institute is in the process of conducting the High Burnup Demonstration Project. Where the objective is to characterize the performance of high-burnup fuel in long term storage. As part of this demonstration, a TN-32B dry storage cask was instrumented and loaded with spent nuclear fuel at North Anna Nuclear Generating Station in November 2017. The project cask provides a unique opportunity to gain information on spent fuel and cask performance from an in-service operational system. The cask was instrumented with thermocouples inside of the fuel assemblies, then loaded and dried using normal procedures. After the drying process a thermal soak period where the cask was left indoors was used to obtain steady temperatures for model comparison along with surface temperatures. This paper details thermal modeling validation work that was done to model both the steady state and transient cases. Two modeling tools were used to predict temperatures in the cask. The general purpose CFD and heat transfer code STAR-CCM+ was used with both a detailed pin-by-pin model was used along with a more efficient k-effective simplification. The other code is COBRA-SFS, a purpose-built detailed thermal modeling tool developed and maintained at PNNL. Results from all the modeling tools were compared blind to the test data. Each model utilized general design information and compared reasonably well to the blind steady state data using convection and surface temperature boundary conditions. Post-test adjustments were later made to better reflect the “as built” conditions of the cask. Based on lessons learned from the steady state analysis the vacuum drying process was modeled in both codes. Utilizing the measured loading conditions the transient models were able to compare very well with the measured data. Overall the project showed an ability to model spent fuel storage conditions very well and future work is planned to generalize the methodology used for vacuum drying.

Thermal Analysis, Spent Fuel, Dry Storage↗
DOE OSTI2022

Thermoacoustic Soret separation

In a fluid mixture in a channel with an axial time-averaged temperature gradient, high-amplitude oscillating flow can greatly increase the axial flux of thermal diffusion (Soret) separation of the components of the mixture. The enhancement occurs when the oscillating lateral temperature gradient greatly exceeds the axial gradient, causing a large oscillating concentration that can be favorably time-phased with the oscillating flow. This process can occur even with a negligible pressure oscillation or with a negligible temperature response to pressure, as is the case in most liquid solutions. The thermal boundary condition imposed by realistic solids on thermoacoustic liquids is imperfect, adding mathematical complications that are absent for typical gases, for which the solid surface is temporally isothermal. Compared with gas mixtures, the high Lewis number in typical liquid solutions reduces the separation flux associated with the time-averaged temperature gradient, but it also reduces the remixing associated with the time-averaged mole-fraction gradient. For large enough channels, the second-law separation efficiency is only slightly reduced from that of steady liquid Soret separation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enabling resonant ultrasound spectroscopy in high magnetic fields

Resonant ultrasound spectroscopy (RUS) is a powerful method to determine elastic constants with high accuracy and precision from a single measurement of the mechanical resonances of a sample. Conventionally, the quantitative extraction of elastic moduli with RUS assumes free boundary conditions which can often lead to the adoption of unstable sample positioning between ultrasonic transducers that is incompatible with extreme environments like high magnetic fields. We show that, under specific conditions, introducing a small amount of adhesive between a RUS sample and ultrasonic transducers introduces a perturbation to the free resonance condition which can be accounted for by a simple model. This means elastic constants can be determined to within the uncertainty of conventional RUS, but with significant improvements including sample stability and control of sample orientation. We demonstrate the efficacy of this approach with measurements on a range of materials including room temperature measurements on polycrystalline metals, temperature-dependent measurements of the structural phase transition in strontium titanate single crystals, and magnetic field-dependent measurements of magnetic phase transitions in gadolinium polycrystals up to 14 T.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
DOE OSTI2021

Learning the solution operator of parametric partial differential equations with physics-informed DeepONets

Partial differential equations (PDEs) play a central role in the mathematical analysis and modeling of complex dynamic processes across all corners of science and engineering. Their solution often requires laborious analytical or computational tools, associated with a cost that is markedly amplified when different scenarios need to be investigated, for example, corresponding to different initial or boundary conditions, different inputs, etc. In this work, we introduce physics-informed DeepONets, a deep learning framework for learning the solution operator of arbitrary PDEs, even in the absence of any paired input-output training data. We illustrate the effectiveness of the proposed framework in rapidly predicting the solution of various types of parametric PDEs up to three orders of magnitude faster compared to conventional PDE solvers, setting a previously unexplored paradigm for modeling and simulation of nonlinear and nonequilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Out-of-plane ferroelectricity and robust magnetoelectricity in quasi-two-dimensional materials

Thin-film ferroelectrics have been pursued for capacitive and nonvolatile memory devices. They rely on polarizations that are oriented in an out-of-plane direction to facilitate integration and addressability with complementary metal-oxide semiconductor architectures. The internal depolarization field, however, formed by surface charges can suppress the out-of-plane polarization in ultrathin ferroelectric films that could otherwise exhibit lower coercive fields and operate with lower power. Here, we unveil stabilization of a polar longitudinal optical (LO) mode in the n = 2 Ruddlesden–Popper family that produces out-of-plane ferroelectricity, persists under open-circuit boundary conditions, and is distinct from hyperferroelectricity. Our first-principles calculations show the stabilization of the LO mode is ubiquitous in chalcogenides and halides and relies on anharmonic trilinear mode coupling. We further show that the out-of-plane ferroelectricity can be predicted with a crystallographic tolerance factor, and we use these insights to design a room-temperature multiferroic with strong magnetoelectric coupling suitable for magneto-electric spin-orbit transistors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

AIVT: Inference of turbulent thermal convection from measured 3D velocity data by physics-informed Kolmogorov-Arnold networks

We propose the artificial intelligence velocimetry-thermometry (AIVT) method to reconstruct a continuous and differentiable representation of the temperature and velocity in turbulent convection from measured three-dimensional (3D) velocity data. AIVT is based on physics-informed Kolmogorov-Arnold networks and trained by optimizing a loss function that minimizes residuals of the velocity data, boundary conditions, and governing equations. We apply AIVT to a set of simultaneously measured 3D temperature and velocity data of Rayleigh-Bénard convection, obtained by combining particle image thermometry and Lagrangian particle tracking. This enables us to directly compare machine learning results to true volumetric, simultaneous temperature and velocity measurements. We demonstrate that AIVT can reconstruct and infer continuous, instantaneous velocity and temperature fields and their gradients from sparse experimental data at a high resolution, providing an additional approach for understanding thermal turbulence.

Science & Technology - Other Topics↗
DOE OSTI2020

Experimental characterization of fragile topology in an acoustic metamaterial

Symmetries crucially underlie the classification of topological phases of matter. Most materials, both natural as well as architectured, possess crystalline symmetries. Recent theoretical works unveiled that these crystalline symmetries can stabilize fragile Bloch bands that challenge our very notion of topology: Although answering to the most basic definition of topology, one can trivialize these bands through the addition of trivial Bloch bands. Here, we fully characterize the symmetry properties of the response of an acoustic metamaterial to establish the fragile nature of the low-lying Bloch bands. Additionally, we present a spectral signature in the form of spectral flow under twisted boundary conditions.

Science & Technology - Other Topics↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Computing Edge States without Hard Truncation

We present a numerical method which accurately computes the discrete spectrum and associated bound states of semi-infinite Hamiltonians which model electronic “edge” states localized at boundaries of one- and two-dimensional crystalline materials. The problem is nontrivial since arbitrarily large finite “hard” (Dirichlet) truncations of the Hamiltonian in the infinite bulk direction tend to produce spurious bound states partially supported at the truncation. Our method, which overcomes this difficulty, is to compute the Green's function of the semi-infinite Hamiltonian by imposing an appropriate boundary condition in the bulk direction; then, the spectral data is recovered via Riesz projection. We demonstrate our method's effectiveness by studies of edge states at a graphene zig-zag edge in the presence of defects modeled both by a discrete tight-binding model and a continuum PDE model under finite difference discretization. Our method may also be used to study states localized at domain wall-type edges in one- and two-dimensional materials where the edge Hamiltonian is infinite in both directions; we demonstrate this for the case of a tight-binding model of distinct honeycomb structures joined along a zig-zag edge. Here, we expect our method to be useful for designing novel devices based on precise wave-guiding by edge states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Explicit Runge–Kutta Methods that Alleviate Order Reduction

Explicit Runge–Kutta (RK) methods are susceptible to a reduction in the observed order of convergence when applied to an initial boundary value problem with time-dependent boundary conditions. We study conditions on explicit RK methods that guarantee high order convergence for linear problems; we refer to these conditions as weak stage order conditions. We prove a general relationship between the method’s order, weak stage order, and number of stages. Furthermore, we derive explicit RK methods with high weak stage order and demonstrate, through numerical tests, that they avoid the order reduction phenomenon up to any order for linear problems and up to order three for nonlinear problems.

explicit Runge–Kutta↗
DOE OSTI2023

Algebra of hyperbolic band theory under magnetic field

We explore algebras associated with the hyperbolic band theory under a magnetic field for the first time. We define the magnetic Fuchsian group associated with a higher genus Riemann surface. By imposing the magnetic boundary conditions for the hyperbolic Bloch states, we construct the hyperbolic magnetic Bloch states and investigate their energy spectrum. We give a connection between such magnetic Bloch states and automorphic forms. Our theory is a general extension of the conventional algebra associated with the band theory defined on a Euclidean lattice/space into that of the band theory on a general hyperbolic lattice/Riemann surface.

Physics↗
DOE OSTI2023

Inclusion of $K\Lambda $ electroproduction data in a coupled channel analysis

Exclusive electroproduction reactions provide an access to the structure of excited baryons. To extract electroproduction multipoles encoding this information, the Julich-Bonn-Washington (JBW) analysis framework is extended to the analysis of differential cross sections in KΛ electroproduction. Here this update enlarges the scope of previous coupled-channel analyses of pions and eta mesons, with photoproduction reactions as boundary condition in all analyzed electroproduction reactions. Polarization observables are predicted and compared to recent CLAS data. The comparison shows the relevance of these data to pin down baryon properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗