Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid solid interfaces”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Detached Melt and Vapor Growth of InI in SUBSA hardware

Indium iodide (InI) single crystals, grown by the Bridgman process, have shown a significant promise as room temperature detector material for X-rays and γ-rays due to the large, close to ideal, bandgap of Eg = 2 eV, and high electrical resistivity. Yet the transport properties of the charge and the detector performance remain inadequate. The Bridgman growth of semiconductor crystals in microgravity have typically resulted in partial dewetting and Marangoni convection, which disturbed the diffusion-controlled growth and produced nonuniform dopant distribution. However, unconfined or partially confined solidification drastically reduced the number of stress-induced defects compared to the confined growth in crucibles on Earth.The goal of the present investigation was to explore melt growth of a heavy metal halide in microgravity, focusing on dewetting, in order to improve the crystalline perfection, i.e., to reduce the number of stress-related defects. InI was chosen because it is non-toxic and has a low melting point of only 365 C (~150 C below the melting point of InSb, used in the SUBSA experiments in 2002). Furthermore, InI melts and evaporates congruently, making it suitable for both growth by directional solidification and growth from the vapor phase. In preparation for the flight experiments, ground-based experiments were conducted in the SUBSA ground unit, consisting of 4 melt growth and 2 vapor growth experiments. The crystals were grown by a gradient freeze method by lowering the setpoint of the furnace. The transparent section of the SUBSA furnace did not allow observations of the solid-liquid interface. Contrary to our expectations, the possible occurrence of dewetting during growth could not be observed. We will present i) the ampoule design and preparation, ii) the experiments conducted in the SUBSA ground unit, iii) the images obtained during the microgravity experiments and iv) the results of the characterization of the grown crystals and detector fabrication.

Crystal Growth↗
DOE OSTI2020

Efficient Mesofluidic Separation of Large Particles in nuclear slurries - 20408

Efficient and effective particle separation is essential to cleanup of many nuclear wastes. For example, particle separation may be used to accelerate settle-decant bottlenecks that limit the throughput necessary to achieve the U.S. Department of Energy waste cleanup mission in a timely manner. Particle separation may be used to protect downstream processing equipment from pulses of solids that may be encountered near (within a few feet of) liquid-solid interfaces during waste transfer operations, thereby preventing work stoppages. Effective particle separation may permit efficient sludge washing. These are only a few of the ways in which particle separation is important. Yet, efficient and effective particle separation in nuclear processing environments remains challenging. For example, settle-decant operations permit larger and heavier particles to settle, leaving smaller, lighter particles suspended. However, settle-decant operations are slow, convection currents may resuspend solids, hindered settling of multicomponent slurries remains incompletely understood, and predictive models for settling in graduated cylinders fail to match limited observations of settling in large waste tanks. Additionally, pumping operations can cause turbulent resuspension of particles when the pump intake is close to the settled solids layer. Filtration techniques that use membranes or partially permeable barriers retain larger particles, permitting smaller particles to permeate. However, filters, including dead-end filtration, are prone to clogging and caking, operate at elevated pressures due to minimal void volume, and increase in pressure during operations, which require more control systems. Therefore, the need for high throughput particle separation techniques that operate with modest pressure drops persists. A novel mesofluidic separator presents the opportunity to effectively and efficiently accelerate the waste cleanup mission. The separator separates large particles from process streams across a broad range of particle sizes and has no moving parts or media to replace, regenerate, or clean. This separator design has an unusually large void volume, permitting operation at much higher flow rates (and lower pressures) than traditional filtration (e.g., dead-end filtration). Industrial-scale flow rates have been demonstrated. In performance testing, mesofluidic separators operate at might flow rates (>90 gpm (0.0006 m{sup 3}/s) in piping 3 inches (0.08 m) in diameter; Re>10{sup 5}), with modest pressure drop (∼25 psi (170 KPa) in testing). In complex, aggregating waste simulants, the separator loses <25% of flow projected over months without back pulsing or chemical cleaning. Mesofluidic separation presents infrastructure advantages, reduces risk, and provides mission impact. Infrastructure advantages include implementation within existing transfer systems, plug and play without facility modifications to safety systems, and operation at low pressures. The separator is flexible in location, may be positioned within or outside of waste tanks, and may replace or augment dead-end and cross-flow filters. Mesofluidic separation reduces risk by minimizing waste (no media to replace or regenerate and no cleaning chemicals to handle and dispose) and by minimizing filter change-out consistent with as low as reasonably achievable (ALARA) exposure to workers. The potential mission impact of these separators is substantial. The separator has potential to support, simplify, and accelerate in-farm transfers and waste feed delivery. Furthermore, the separator can free up settling and holding tanks in the Direct Feed Low-Activity Waste (DFLAW) mission by reducing or eliminating post-transfer settling and wait times with the potential to free up 1 Mgal (4000 m{sup 3}) of double-shell tank space. This paper quantitatively compares mesofluidic separation to dead-end filtration, discusses scale-up results, and considers the separator's potential to efficiently and effectively reduce the long-term environmental ability of particulate-rich nuclear wastes. (authors)

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

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗
DOE OSTI2024

Characterization of compositional variation and solidification condition on rapid solidification behavior of beta stabilized titanium

Rapid solidification (RS) is a field of study that has increased in importance due to the rise of additive manufacturing (AM) techniques. The breakdown of equilibrium at the solid-liquid interface results in an inability to use traditional solidification models that have been applied to great success in casting. As such, there is a current need for a high-throughput testing method to screen alloys at rapid solidification conditions without the significant cost for powder atomization. In this work, twin hammer splat quenching (SQ) is utilized and theoretically developed as a technique uniquely placed to test alloy systems for AM applications. Using analytical thermal spray models and numerical heat transfer simulations, the cooling rate and solidification conditions of splat quenching are compared across several parameters to evaluate the key parameters for generating cooling rates relevant to AM process. Once a framework for the cooling rates and relevancy of SQ are developed, SQ was used to process alloy compositions in the Ti5553 alloy space (AMS 7026) to study the formation of detrimental secondary phases occurring upon rapid quenching from the melt. A martensitic phase is found to occur at the lower end of the standard, but the composition listed in the standard were found to be safe. Finally, the rapid solidification behavior of Ti-5553 is evaluated by comparing splat quenching, laser welding, and plasma atomization that offer a wide window of solidification conditions. The segregation and solidification structure are analyzed and compared to models to provide a deeper understanding of the solidification behavior of beta stabilized titanium alloys.

Williamson, Christopher J.↗
NASA NTRS1993

Effect of the baffle on the spacecraft fluid propellant viscous stress and moment fluctuations

The dynamical behavior of fluids affected by asymmetric gravity jitter oscillations, in particular the effect of surface tension on rotating fluids in a partially filled dewar tank imposed by time-dependent directions of background gravity have been investigated. Results show that a lower frequency gravity jitter imposed on the time-dependent variations of the direction of background gravity induced a greater amplitude of oscillations and a stronger degree of asymmetry in liquid-vapor interface geometry than that caused by higher-frequency gravity jitter. As viscous force between liquid and solid interface greatly contributes to the damping effect of slosh wave excitation, a rotating dewar with installed baffle provides more areas of liquid-solid interface than a rotating dewar without baffle. Damping by baffle reduces the amplitude of slosh wave excitation, angular momentum, and fluid moment fluctuation and lowers the degree of asymmetry in the liquid-vapor distribution.

Hung, R. J.↗
NASA NTRS1995

Effect of Baffle on Gravity-Gradient-Excited Slosh Waves and Spacecraft Moment and Angular-Momentum Fluctuations in Microgravity

The dynamical behavior of fluids affected by the asymmetric gravity gradient acceleration has been investigated. In particular, the effects of surface tension on partially filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank with and without baffles are studied. Results of slosh wave excitation along the liquid-vapor interface induced by gravity gradient acceleration indicate that the gravity gradient acceleration is equivalent to the combined effect of a twisting force and a torsional moment acting on the spacecraft. The results are clearly seen from one-up one-down and one-down one-up oscillations in the cross-section profiles of two bubbles in the vertical (r, z)-plane of the rotating dewar, and from the eccentric contour of the bubble rotating around the axis of the dewar in a horizontal (r, theta)-plane. As the viscous force, between liquid and solid interface, greatly contributes to the damping of slosh wave excitation, a rotating dewar with baffles provides more areas of liquid-solid interface than that of a rotating dewar without baffles. Results show that the damping effect provided by the baffles reduces the amplitude of slosh wave excitation and lowers the degree of asymmetry in liquid-vapor distribution. Fluctuations of angular momentum and fluid moment caused by the slosh wave excited by gravity gradient acceleration with and without baffle boards are also investigated. It is also shown that the damping effect provided by the baffles greatly reduces the amplitudes of angular momentum and fluid moment fluctuations.

Hung, R. J.↗
DOE OSTI2022

Manipulating Interfacial Dissolution–Redeposition Dynamics to Resynthesize Electrode Surface Chemistry

The solid-liquid electrochemical interface offers a two-dimensional environment for geometrically confined interfacial reactions to tailor electrode surface chemistry under operating conditions. In this study, we demonstrate that the dissolution and redeposition kinetics of transition metal cations, a ubiquitous phenomenon at the electrochemical interface, can be manipulated to regulate the chemical composition and crystal structure of the electrode surface as well as the overall electrochemical performance. Foreign cations (e.g., Ti 4+ ), either added as electrolyte additives or dissolved from surface coatings, can rapidly participate in the electrode dissolution-redeposition process, and facilitate the establishment of the dissolution-redeposition equilibrium. Our work expands the control over the electrochemical reactions at the solid-liquid interface and provides new insights into interfacial studies in electrochemistry and surface science.

25 ENERGY STORAGE↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

In situ x-ray imaging to understand subsurface behavior during continuous wave laser drilling

A limited understanding regarding the underlying dynamics and mechanisms of material removal during continuous wave laser drilling has presented significant challenges in achieving precision and process control. Here, to address this, we employed high-fidelity, in situ synchrotron x-ray imaging to reveal previously unknown material behaviors during continuous wave laser drilling with power modulation. Our findings highlight that high-aspect ratio drill holes are achieved when the laser modulation frequency falls within the range of 8–12 kHz, provided that the laser average power and modulation amplitude levels meet the specified limits. Under these conditions, we identified a material removal mechanism driven by incremental accumulation of recoil pressure that gradually pushes material upward from deep within the substrate to the surface. This mechanism manifested as a low-frequency fluctuation in the vapor depression depth, resulting in periodic instances of material ejection. Furthermore, our study underscores that rapid expansion of the melt pool and the widening of the drill hole opening can impede effective material removal by redirecting energy from material ejection to increasing the melt pool size. This investigation contributes essential insights into the subsurface dynamics involved in the drilling of high-aspect ratio holes, furthering our fundamental understanding of this intricate process.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
DOE OSTI2021

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D 2 O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

Interfacial Force Field Characterization in a Constrained Vapor Bubble Thermosyphon

Isothermal profiles of the extended meniscus in a quartz cuvette were measured in the earth's gravitational field using an image-analyzing interferometer that is based on computer-enhanced video microscopy of the naturally occurring interference fringes. These profiles are a function of the stress field. Experimentally, the augmented Young-Laplace equation is an excellent model for the force field at the solid-liquid-vapor interfaces for heptane and pentane menisci on quartz and tetradecane on SFL6. The effects of refractive indices of the solid and liquid on the measurement techniques were demonstrated. Experimentally obtained values of the disjoining pressure and dispersion constants were compared to those predicted from the Dzyaloshinskii - Lifshitz - Pilaevskii theory for an ideal surface and reasonable agreements were obtained. A parameter introduced gives a quantitative measurement of the closeness of the system to equilibrium. The nonequilibrium behavior of this parameter is also presented

DasGupta, Sunando↗
DOE OSTI2020

Liquid–Solid Boundaries Dominate Activity of CO 2 Reduction on Gas-Diffusion Electrodes

Electrochemical CO 2 electrolysis to produce hydrocarbon fuels or material feedstocks offers a renewable alternative to fossilized carbon sources. Gas-diffusion electrodes (GDEs), composed of solid electrocatalysts on porous supports positioned near the interface of a conducting electrolyte and CO 2 gas, have been able to demonstrate the substantial current densities needed for future commercialization. These higher reaction rates have often been ascribed to the presence of a three-phase interface, where solid, liquid, and gas provide electrons, water, and CO 2 , respectively. Conversely, mechanistic work on electrochemical reactions implicates a fully two-phase reaction interface, where gas molecules reach the electrocatalyst’s surface by dissolution and diffusion through the electrolyte. Because the discrepancy between an atomistic three-phase versus two-phase reaction has substantial implications for the design of catalysts, gas-diffusion layers, and cell architectures, the nuances of nomenclatures and governing phenomena surrounding the three-phase-region require clarification. In this paper we outline the macro, micro, and atomistic phenomena occurring within a gas-diffusion electrode to provide a focused discussion on the architecture of the often-discussed three-phase region for CO 2 electrolysis. From this information, we comment on the outlook for the broader CO 2 electroreduction GDE cell architecture.

CO2 reduction↗
DOE OSTI2021

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. Furthermore, we also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Parallel Grand Canonical Monte Carlo (ParaGrandMC) Simulation Code

This report provides an overview of the Parallel Grand Canonical Monte Carlo (ParaGrandMC) simulation code. This is a highly scalable parallel FORTRAN code for simulating the thermodynamic evolution of metal alloy systems at the atomic level, and predicting the thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces.

Vesselin I Yamakov↗
NASA NTRS2022

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗
DOE OSTI2020

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice Boltzmann model for conjugate heat transfer across thin walls

A lattice Boltzmann (LB) model with an efficient and accurate interface treatment for conjugate heat transfer across a thin wall between two different media is developed. The proposed interface treatment avoids fine meshing and computation within the thin layer; instead, the energy balance within the thin layer and the conjugate conditions on each interface are utilized to construct explicit updating schemes for the microscopic distribution functions of the LB model at the interior lattice nodes of the two media next to the thin layer. The proposed interface schemes reduce to the standard interface scheme for conjugate conditions in the literature in the limit of zero thickness of the thin layer, and thus it can be considered a more general interface treatment. A simplified version of the interface treatment is also proposed when the heat flux variation along the tangential direction of the thin layer is negligible. Three representative numerical tests are conducted to verify the applicability and accuracy of the proposed interface schemes. The results demonstrate that the intrinsic second-order accuracy of the LB model is preserved with the proposed interface schemes for thin layers with constant tangential fluxes, while for general situations with varying tangential fluxes, first-order accuracy is obtained. Finally, this interface treatment within the LB framework is attractive in conjugate heat transfer modeling involving thin layers for its simplicity, accuracy, and significant reduction in computational resources.

42 ENGINEERING↗