Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interface capturing”

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

Effect of Packing Nonuniformity at the Fiber Bundle–Case Interface on Performance of Hollow Fiber Membrane Gas Separation Modules

High-fidelity simulations of momentum and mass transfer within a hollow fiber gas separation membrane module are here reported. The simulations capture the potential detrimental effects of poor fiber packing at the bundle–case interface on fluid distribution and performance. Results are presented for both circular and planar fiber bundles. The length over which bundle–case gaps affects flow is determined. The length increases dramatically with increasing fiber packing fraction. As the packing fraction approaches 0.6, the impact extends over the entire bundle diameter for small modules (<1000 fibers). The results clearly demonstrate the detrimental effect of poor packing along the case and can be used to develop module manufacturing guidelines. To reduce computational costs, an equivalent planar bundle module approximation is developed. The approximate simulations agree well with results from full 3-D simulations and can reduce computational costs without sacrificing fidelity.

42 ENGINEERING↗
DOE OSTI2023

Synergizing Electron and Heat Flows in Photocatalyst for Direct Conversion of Captured CO 2

Abstract We report a ternary hybrid photocatalyst architecture with tailored interfaces that boost the utilization of solar energy for photochemical CO 2 reduction by synergizing electron and heat flows in the photocatalyst. The photocatalyst comprises cobalt phthalocyanine (CoPc) molecules assembled on multiwalled carbon nanotubes (CNTs) that are decorated with nearly monodispersed cadmium sulfide quantum dots (CdS QDs). The CdS QDs absorb visible light and generate electron‐hole pairs. The CNTs rapidly transfer the photogenerated electrons from CdS to CoPc. The CoPc molecules then selectively reduce CO 2 to CO. The interfacial dynamics and catalytic behavior are clearly revealed by time‐resolved and in situ vibrational spectroscopies. In addition to serving as electron highways, the black body property of the CNT component can create local photothermal heating to activate amine‐captured CO 2 , namely carbamates, for direct photochemical conversion without additional energy input.

Choi, Chungseok↗
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2023

The effect of terrain-influenced winds on fire spread in QUIC-Fire

Here in this manuscript, we describe the implementation of the terrain-following version of QUIC-URB into QUIC-Fire and a demonstration of the impacts of terrain-influenced winds on QUIC-Fire-simulated fire spread. No changes to the underlying QUIC-Fire fire spread algorithm were made other than what was required to correctly account for the inclusion of terrain. This paper summarizes simulations used to understand how the QUIC-URB terrain-influenced winds affect upslope fire behavior without additional changes to the QUIC-Fire fire spread algorithm. Previously published FIRETEC results are compared to simulation results from the modified QUIC-Fire (incorporating terrain-influenced winds) that use the same topographies and fuels. QUIC-Fire results showed overall similar behaviors in terms of how the topographies affected fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all non-flat topographies were over-predicted when compared to FIRETEC. Lateral fire growth showed similar trends with FIRETEC between topographies but did not capture the increase in spread due to a diagonal interface between grassland and forested fuel regions in the test domain. These simulations suggest possible refinements that are necessary to improve QUIC-Fire and thus guide ongoing efforts related to: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Data-driven approach to Determining the Fidelity in the Hardware-in-the-loop Systems using Subspace Identification Method

One of the major questions in any Hardware-in-the-loop (HiL) simulation is to understand the fidelity of the HiL simulation itself which is most often indicated qualitatively as either high or low instead of quantifying it. Being cognizant of the level of fidelity forms the crux to assess the validity and credibility of the HiL simulation. In this work, we address this issue by developing a systematic, data-driven approach to assess the fidelity of any HiL simulation, and more specifically, the fidelity of the interface between the simulator and hardware in an HiL simulation. Applying a subspace identification method, a linear system representing the interface is obtained from the time-series data captured between the hardware and simulator in the HiL simulation. Finally, the fidelity is defined based on this linear system properties. The proposed data-driven fidelity quantification framework is illustrated on the IEEE 123 node feeder system running on HYPERSIM and interacting with virtual protection relays.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-Fidelity Modeling of Curbside Driving Behavior in SUMO: Preprint

Recent advances in mobility technologies have considerably increased the need for curbside use and congestion at curb-lanes. At airports and in dense urban centers, the growing adoption of Transportation Network Companies (TNCs), such as Uber and Lyft, has led to an increase in curbside traffic due to the reduction in the use public transit and other shared modes. In dense city areas, the curbside demand and congestion were further augmented by the growing use of e-commerce, which requires home-delivery of goods. To model vehicles' curbside interactions realistically, there is a need for high-fidelity simulation models that can represent the complex interplay of vehicles as they complete pick-up or drop-off trips. Vehicles' curbside behaviors are usually not modeled with high-fidelity in microscopic simulators such Simulation of Urban of Mobility (SUMO). This paper presents such a high-resolution model for vehicle curbside behavior in SUMO, using SUMO's Traci interface. Simulations results demonstrate that the improved curbside behavior captures vehicles' curbside maneuvering more accurately than the default SUMO curbside model, and increases the curbside utilization and capacity.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Deformation profile and interface-mediated defect interaction in Cu/CuZr nanolaminates: An effective-temperature description

Both simulations and experiments have suggested that Cu/CuZr nanolaminates are stronger and more ductile than their individual constituents due to interface-mediated interactions between plasticity carriers. In this work, we use the effective-temperature theories of dislocation and amorphous shear-transformation-zone (STZ) plasticity to study amorphous-crystalline interface (ACI)-mediated plasticity in Cu/CuZr nanolaminates under mechanical straining. The model is shown to capture reasonably well the measured deformation response when strained either in tension parallel to or in compression normal to the amorphous-crystalline interface. Our analysis indicates that increasing CuZr or decreasing Cu layer thickness increases the maximum flow stress for both perpendicular and parallel loading cases. Furthermore, for the cases of parallel and perpendicular loading, the maximum flow stress values are 3.4 and 2.5GPa, respectively. Furthermore, increasing the strain rate for the parallel loading case decreases the slip strain in the amorphous and crystalline layers. Additionally, for the perpendicular loading case, an increase in strain rate decreases the amorphous layer slip but increases the crystalline layer slip. In all slip strain analyses, maximum slip strain occurs at the ACI, thus indicating that plasticity carriers accumulate at the interface and are absorbed there. These findings indicate a significant anisotropy in strength with greater sensitivity to layer thickness for the case of tensile loading parallel to the ACI. Further findings signify that slip strain is more sensitive when the nanolaminate is compressed perpendicular to the ACI.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Capturing thin structures in VOF simulations with two-plane reconstruction

A novel interface reconstruction strategy for volume of fluid (VOF) methods is introduced that represents the liquid-gas interface as two planes that co-exist within a single computational cell. In comparison to the piecewise linear interface calculation (PLIC), this new algorithm greatly improves the accuracy of the reconstruction, in particular when dealing with thin structures such as films. The placement of the two planes requires the solution of a non-linear optimization problem in six dimensions, which has the potential to be overly expensive. Further, an efficient solution to this optimization problem is presented here that exploits two key ideas: an algorithm for extracting multiple plane orientations from transported surface data, and an efficient and mass-conserving distance-finding algorithm that accounts for two planes with arbitrary orientation. Additionally, a simple and robust strategy is presented to accurately represent the surface tension forces produced at the interface of subgrid-thickness films. The performance of this new VOF reconstruction is demonstrated on several test cases that illustrate the capability to handle arbitrarily thin films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CO 2 adsorption mechanisms at the ZIF-8 interface in a Type 3 porous liquid

Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO 2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO 2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO 2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO 2 molecules were unable to diffuse into the ZIF-8 pore, with CO 2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO 2 adsorption should be based on the strength of gas binding at the solvated particle surface.

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

Study of Void Formation at the Lithium|Solid Electrolyte Interface

There is growing recognition of the critical role of void formation in lithium metal anodes in solid-state batteries and its impact on electrochemical performance. While experimental studies have demonstrated the challenges ensuing from void formation at the lithium metal interface with the solid electrolyte, there is a need to understand and quantify the role of intrinsic transport properties in lithium metal and the impact of external stimuli, such as temperature, pressure, and current density. Here, we develop this understanding by constructing a phase field-based model that captures the evolution of void domains at the lithium-solid electrolyte interface. Growth of the pores is driven by the fast removal of lithium from the interface during stripping at high current densities. Relative magnitudes of the bulk and surface lithium diffusivities, along with the applied current density, dictate the final pore morphology. Increasing the temperature results in faster diffusion, while external applied pressure causes creep flow of lithium, both of which help to mitigate the evolution of voids by quickly transporting metal from the bulk to the interface. Finally, a phase map as a function of temperature and pressure is developed as a guide to determine the regions that can lead to the stable cycling of lithium metal.

25 ENERGY STORAGE↗
DOE OSTI2020

An energy-based coupling approach to nonlocal interface problems

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, where traditional partial differential equations fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a rigorous nonlocal interface theory needed to support numerical developments. In this paper, we use an energy-based approach to develop a mathematically rigorous nonlocal interface theory which provides a physically consistent extension of the classical perfect interface PDE formulation. Numerical examples validate the proposed framework and demonstrate the scope of our theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Ab Initio Study of the Interface of the Solid-State Electrolyte Li 9 N 2 Cl 3 with a Li-Metal Electrode

Ab initio molecular dynamics study of an electrochemical interface between a solid-state-electrolyte Li 9 N 2 Cl 3 and a Li-metal is performed to analyze interphase formation due to electrochemical instability. An external electric field is applied to emulate the interface under charging conditions. At open circuit, we observed electrochemical stability at both sides of the interface, without electrolyte decomposition, or new phase formation, in good agreement with reported experiments. However, a structural re-accommodation takes place at both sides of the interface. Some oxidation of the surface Li-atoms takes place when Cl and N anions capture Li from the metal to complete octahedral structures, a Bader charge and electron density analyses explain the difference between Cl and N behavior. At the interface with the Li 9 N 2 Cl 3 (002) facet, the larger negative charge on low-coordinated N atoms is able to retain more Li atoms than the Cl atoms do. At charging conditions applying an external electric field (0.5 V/Å), no major changes are observed in atomic diffusion or structural changes at the interphase; however, the main difference is the change of the Li-metal crystal structure. Due to migration of Li from the metal anode to the electrolyte, vacancies are created on the anode, causing a structural re-arrangement at the anode surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-speed synchrotron X-ray phase-contrast imaging for evaluating microscale damage mechanisms and tracking cracking behaviors inside cross-ply GFRCs

In this study, we integrated the high-speed synchrotron X-ray phase-contrast imaging (PCI) with a modified Kolsky compression bar loading platform to visualize the dynamic failure processes of cross-ply glass fiber reinforced composites (GFRCs). Four S-2 glass/SC-15 composite specimens were prepared, having similar thicknesses but different stacking sequences, namely as [0 12 /90 12 ], [90 12 /0 12 ], [0 8 /90 8 /0 8 ], and [90 8 /0 8 /90 8 ]. Three-dimensional synchrotron X-ray computed tomography and scanning electron microscopy (SEM) were employed to examine the microstructures and quantify the fiber volume fractions. Each specimen was notched and subjected to a dynamic three-point flexural loading. The onset of cracking close to the notch tip, crack propagation in 0° or 90° plies and their interface, crack opening, and ultimately failure of the specimen were captured by high-speed synchrotron X-ray PCI. Additional dynamic experiments were performed to determine the average time when the stress wave propagated through the specimen and correlate the X-ray images with the specimen&#x27;s force-deflection response. Finally, the surface morphology of each specimen after the dynamic loading was imaged by SEM. Comparison between real-time X-ray images and post-fracture SEM images demonstrated the capability of the X-ray method to record damage evolution inside composites. Furthermore, the high spatial and temporal resolutions of the X-ray setup and edge enhancement by PCI enabled the identification of microscale damage features within 1 μs. Two damaging processes were identified, crack growth was quantified, and fracture toughness of the composites was evaluated. The method is deemed useful to reveal microscale damage mechanisms and track cracking behaviors inside cross-ply GFRCs under dynamic loading in real time.

42 ENGINEERING↗
DOE OSTI2025

Autonomous Synthesis of Metastable Materials Using a Modular Mixed-Flow Reactor

Understanding and controlling atomic-level processes at solid-liquid interfaces is key to advancing technologies in energy storage, carbon capture, critical element recovery, and materials synthesis. Many of these processes are dominated by the formation of short-lived intermediate precipitates that determine the final properties of synthesized materials. However, studying these intermediates is challenging due to their sensitivity and the reliance on trial-and-error methods. To address this, we developed an automated variable-volume mixed-flow reactor (MFR) to optimize metastable material synthesis and investigate rapid kinetic processes. This state-of-the-art MFR system, paired with an automated modeling framework, enables efficient synthesis and real-time analysis of transient phases. Benchmarking with advanced capabilities, such as wide-/small-angle X-ray scattering, allows us to resolve fast nucleation and growth dynamics that were previously inaccessible. By combining automation, ML-guided optimization, and tailored kinetic modeling, this approach provides a robust platform for improving material design and achieving precise control over solid-liquid reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗