SEARCH · Engineering Papers

Results for “Free surface flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Charge decay in the spatial afterglow of plasmas and its impact on diffusion regimes

Abstract The spatial afterglow is a region at the boundary of a non-equilibrium plasma where charged species relax into ambient equilibrium. In many applications, the spatial afterglow is the part of the plasma that interacts with surfaces, such as suspended particles or a material substrate. However, compared to the bulk plasma, there has been little effort devoted to studying the properties of the spatial afterglow, and a fundamental analysis has not yet been developed. Here, we apply double Langmuir probe measurements and develop an advection-diffusion-recombination model to provide a detailed description of charged species in the spatial afterglow over a wide range of pressures, temperatures, plasma dimensions, and flow rates. We find that the density of charged species in the spatial afterglow decays by orders of magnitude, which leads to a transition from ambipolar to free diffusion. These insights can be used to explain or predict experimental observations of phenomena, such as the charging of dust grains and the dose of charged species to a biomaterial.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

An interface-aware sub-scale dynamics multi-material cell model for solids with void closure and opening at all speeds

Here, we present a multi-material cell model (closure model) for demanding arbitrary Lagrangian-Eulerian (ALE) simulations of fluids and solids. It is based on the interface-aware sub-scale dynamics (IASSD) approach which utilizes the exact material interface geometry within the computational cell to calculate internal material interactions. Our formulation of the closure model also aims to improve the accuracy in low-speed impact events. Voids are used to represent ambient vacuum and internal free boundaries of the distinct materials. Void regions can close and open at contact surfaces, allowing a transition from contact physics to free motion in vacuum. The coupling of void closure and opening with a new formulation of the IASSD model for solids is tested on several one- and two-dimensional numerical examples, ranging from gas expansion in vacuum to planar and round object impacts at various speeds.

42 ENGINEERING↗

DOE OSTI2020

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗

DOE OSTI2026

Engineering Assembly Kinetics and Line Roughness in Solvent Vapor-Annealed Block Copolymer/Homopolymer Blends

Block copolymer (BCP) directed self-assembly (DSA) is a promising route to enhance lithography resolution by multiplying nanopattern density and reducing feature roughness. Eliminating kinetically trapped self-assembly defects requires fast self-assembly. However, acceleration strategies like solvent vapor annealing or homopolymer blending broaden domain interfaces, implying a trade-off in increased feature roughness. In this work, we experimentally investigate this apparent dilemma between self-assembly kinetics and line roughness for solvent vapor-annealed thin films of a lamellar poly(styrene-block-2-vinylpyridine) (PS-b-P2VP) BCP blended with PS and P2VP homopolymers. Binary blends with PS or P2VP homopolymers and ternary blends incorporating both in equal weight fractions were solvent vapor annealed using acetone, a near-neutral solvent for PS and P2VP, followed by P2VP-selective vapor-phase infiltration with alumina (AlOx) and polymer etching. Binary blends with P2VP exhibit a modest kinetic enhancement but also higher line-edge and -width roughness due to the increased frequency of P2VP protrusions and bridge defects in the alumina line patterns. In contrast, binary blends with PS self-assemble noticeably faster, while domain asymmetry from the added PS homopolymer reduces roughness by curbing the number of alumina protrusions and bridge defects. Ternary blends maintain DSA line patterns across a wider composition window and, at higher homopolymer loadings, reduce roughness at length scales near the lamellar period, consistent with a reduced impact of intradomain compositional fluctuations. These findings provide important insights for codesigning blend compositions and process flows to achieve high-resolution, defect free patterns with minimal roughness through BCP DSA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Microfluidic liquid sheets as large-area targets for high repetition XFELs

The high intensity of X-ray free electron lasers (XFELs) can damage solution-phase samples on every scale, ranging from the molecular or electronic structure of a sample to the macroscopic structure of a liquid microjet. By using a large surface area liquid sheet microjet as a sample target instead of a standard cylindrical microjet, the incident X-ray spot size can be increased such that the incident intensity falls below the damage threshold. This capability is becoming particularly important for high repetition rate XFELs, where destroying a target with each pulse would require prohibitively large volumes of sample. We present here a study of microfluidic liquid sheet dimensions as a function of liquid flow rate. Sheet lengths, widths and thickness gradients are shown for three styles of nozzles fabricated from isotropically etched glass. In-vacuum operation and sample recirculation using these nozzles is demonstrated. The effects of intense XFEL pulses on the structure of a liquid sheet are also briefly examined.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

A New Route Toward Atomically Flat and Defect-Free Ge/SiGe Planar Heterostructures

Germanium-based planar heterostructures are emerging as versatile platforms for realizing quantum devices. In particular, planar Ge/SiGe quantum wells (QWs) host hole states with exceptionally high mobility and strong, electrically tunable spin–orbit interactions, enabling full electrical control of quantum information. A key requirement for these systems is the growth of microscopic, defect-free, atomically flat Ge QW heterostructures on relaxed or reverse-graded SiGe buffer layers, as well as on commercial Ge substrates. While several physical deposition techniques have demonstrated high-quality planar Ge/SiGe QWs, a major challenge remains: minimizing defect density typically requires high growth temperatures, which are incompatible with standard CMOS process flows. Here, we present a convenient low-temperature process for realizing high-quality planar SiGe/Ge heterostructures using a combination of thermal and electron-beam evaporation. We systematically map the effects of ex-situ and in-situ substrate preparation protocols, growth temperature, and post-deposition annealing conditions, and correlate these parameters with surface roughness and defect density. We find that in-situ oxide desorption conditions and post-annealing parameters have the most pronounced impact on improving surface quality. Under optimized conditions, we achieve atomically smooth surfaces with root-mean-square roughness σrms ≤ 10 Å and negligible defect density. Interestingly, thermally evaporated Ge layers exhibit oriented triangular crystallites that elongate upon post-annealing in the presence of high Ge vapor pressure. These results demonstrate that this simple, low-temperature deposition approach is a viable and effective route for achieving high-quality Ge QW heterostructures, with strong potential for scalable quantum computing and sensing applications.

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

DOE OSTI2026

A New Route Toward Atomically Flat and Defect-Free Ge/SiGe Planar Heterostructures

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

DOE OSTI2024

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Equation‐Free Surrogate Modeling of Geophysical Flows at the Intersection of Machine Learning and Data Assimilation

Abstract There is a growing interest in developing data‐driven reduced‐order models for atmospheric and oceanic flows that are trained on data obtained either from high‐resolution simulations or satellite observations. The data‐driven models are non‐intrusive in nature and offer significant computational savings compared to large‐scale numerical models. These low‐dimensional models can be utilized to reduce the computational burden of generating forecasts and estimating model uncertainty without losing the key information needed for data assimilation (DA) to produce accurate state estimates. This paper aims at exploring an equation‐free surrogate modeling approach at the intersection of machine learning and DA in Earth system modeling. With this objective, we introduce an end‐to‐end non‐intrusive reduced‐order modeling (NIROM) framework equipped with contributions in modal decomposition, time series prediction, optimal sensor placement, and sequential DA. Specifically, we use proper orthogonal decomposition (POD) to identify the dominant structures of the flow, and a long short‐term memory network to model the dynamics of the POD modes. The NIROM is integrated within the deterministic ensemble Kalman filter (DEnKF) to incorporate sparse and noisy observations at optimal sensor locations obtained through QR pivoting. The feasibility and the benefit of the proposed framework are demonstrated for the NOAA Optimum Interpolation Sea Surface Temperature (SST) V2 data set. Our results indicate that the NIROM is stable for long‐term forecasting and can model dynamics of SST with a reasonable level of accuracy. Furthermore, the prediction accuracy of the NIROM gets improved by almost one order of magnitude by the DEnKF algorithm.

Pawar, Suraj↗

DOE OSTI2022

Preactivated zeolite nanosheet plate-tiled membrane on porous PVDF film: Synthesis and study of proton-selective ion conduction

Preactivated MFI zeolite nanosheet plates (ZNPs) with large areas (~2.0 × 2.0 μm 2 ) and nanometer thicknesses (~60 nm) were prepared directly in liquid-dispersed state. The individual ZNP was a stack of 4-nm-thick single crystalline zeolite nanosheets (ZNs) interlinked by Si–O–Si bonds between neighboring ZN surfaces. The dispersed open-pore ZNPs allowed formulating suspensions for tiling ZNP membranes on polymer substrates without requiring post-coating activation. A pinhole-free ZNP-tiled (ZNPT) membrane has been achieved on macroporous PVDF film by a self-repairing vacuum-assisted filtration coating method. The resultant ZNPT layer was ~500 nm-thick consisting of about 7 ZNP layers with inter-ZNP width and inter-ZNP entrance porosity around 9 nm and 2%, respectively. The ZNPT-PVDF membrane exhibited high selectivity to proton transport over vanadyl ion and low resistances to proton conduction in aqueous solutions. The membrane also demonstrated to function as an efficient ion separator for the vanadium redox flow battery. The reported synthesis of preactivated ZNP suspension may overcome the major hurdle to realizing polymer-supported ZN membranes, which is the inefficiency of existing methods in obtaining readily dispersible open-pore ZNs. Finally, the ZNPT-PVDF membrane can be a more affordable and sustainable alternative to the Nafion-based ion separation membranes.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Model and remote-sensing-guided experimental design and hypothesis generation for monitoring snow-soil–plant interactions

In this study, we develop a machine-learning (ML)-enabled strategy for selecting hillslope-scale ecohydrological monitoring sites within snow-dominated mountainous watersheds, with a particular focus on snow-soil–plant interactions. Data layers rely on spatial data layers from both remote sensing and hydrological model simulations. Specifically, a Landsat-based foresummer drought sensitivity index is used to define the dependency of the annual peak plant productivity on the Palmer drought severity index in the early growing season. Hydrological simulations provide the spatiotemporal dynamics of near-surface soil moisture and snow depth. In this framework, a regression analysis identifies the key hydrological variables relevant to the spatial heterogeneity of drought sensitivity. We then apply unsupervised clustering to these key variables, using the Gaussian mixture model, to group hillslopes into several zones that have divergent relationships regarding soil moisture, snow dynamics, and drought sensitivity. Using the datasets collected in the East River Watershed (Crested Butte, Colorado, United States), results show that drought sensitivity is significantly correlated with model-derived soil moisture and snow-free timing over space and time. The relationship is, however, non-linear, such that the correlation decreases above a threshold elevation and in a heavy snow year due to large snowpacks, lateral flow, and soil storage limitations. Clustering is then able to define the zones that have high or low sensitivity to drought, as well as the mid-elevation regions where sensitivity is associated with the topographic aspect and net potential radiation. In addition, the algorithm identifies the most representative hillslopes with road/trail access within each zone for installing monitoring sites. Our method also aims to significantly increase the use of ML and model-simulation results to guide critical zone and watershed monitoring activities.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Ring patterns generated by an expanding colloidal meniscus

The drop-and-dry is a common technique allowing for creation of periodic nanoparticle (NP) structures for sensing, photonics, catalysis, etc. However, the reproducibility and scalability of this approach for fabrication of NP-based structures faces serious challenges due to the complexity of the simple, at first glance, evaporation process. In this work we study the effect of the spatial confinement on the NP self-assembly under slow solvent evaporation, when the air-liquid-substrate contact line (CL) expands from the center towards the walls of a cylindrical cell, forming a toroid. Using in situ video monitoring of the stick-slip CL motion, we find regular hydrodynamic perturbations in the meniscus, and reveal fine details of the formation of quasiperiodic rings of close packed NP layers. We report that drying of the toroidal NP droplet has a number of important differences from drying of the classical hemispherical colloidal drops. In toroidal drops we observe linear-in-time average meniscus motion, in contrast to the hemispherical drops where the meniscus moves as a square root of time. While both droplet geometries produce NP ring patterns, the ring width for the toroidal drop decreases with increasing ring radius, while it decreases with decreasing the radius of the hemispherical drop. We suggest that free ligands are the main cause of the Marangoni instabilities driving the periodic vorticity in the meniscus. Finally, we show that the usually ignored contact line tension may yield a considerable contribution to the CL pinning causing the CL slip-stick motion and the ring formation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗

DOE OSTI2023

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Impurity transport in PISCES-RF

Linear plasma devices (LPD) utilizing a helicon plasma source, a high density light ion source, can generate impurities due to progressive erosion of the radio frequency (RF) transmission window caused by rectified sheath voltage. These source-born impurities can entrain and be transported by the plasma toward a target, affecting plasma-material interaction studies. Earlier work on material testing in Prototype-Materials Plasma Exposure eXperiment at ORNL revealed significant source impurity deposition on downstream targets. However, using a similar RF source, no target impurity deposition is observed in Plasma Interaction Surface Component Experimental Station (PISCES)-RF despite evidence of RF window erosion in the source region, thereby motivating the present work. Experimentally, using various magnetic field configurations upstream of the PISCES-RF plasma source and seeding titanium (Ti) impurities at various axial locations, impurity transport and deposition along the machine axis were investigated. It was found that Ti deposition was localized to the side of the plasma source where the Ti impurity was seeded. In contrast, aluminum (Al) deposition, originating from the sputtering of the helicon window, occurred predominantly upstream of the plasma source, suggesting an asymmetry in the axial transport of eroded RF window material. These observations suggest a stagnation of the parallel plasma flow immediately downstream of the plasma source, with impurity ions remaining unmagnetized near the source upstream. Al deposition in magnetic field-free regions in PISCES-RF indicates that sputtered Al impurities likely remained neutral due to their large ionization mean-free path under PISCES-RF conditions. Plasma modeling and simulation supported this, indicating that Al-neutrals transport toward the helicon source upstream for low electron density cases. It was found that the Larmor radius of the Al ions was greater than the plasma radius towards the source upstream and remained weakly magnetized in PISCES-RF, meaning that plasma source-born impurities are not efficiently entrained in the plasma flow. These findings provide critical insights into impurity transport in helicon plasma-based LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Roadmap on emerging concepts in the physical biology of bacterial biofilms: from surface sensing to community formation

Abstract Bacterial biofilms are communities of bacteria that exist as aggregates that can adhere to surfaces or be free-standing. This complex, social mode of cellular organization is fundamental to the physiology of microbes and often exhibits surprising behavior. Bacterial biofilms are more than the sum of their parts: single-cell behavior has a complex relation to collective community behavior, in a manner perhaps cognate to the complex relation between atomic physics and condensed matter physics. Biofilm microbiology is a relatively young field by biology standards, but it has already attracted intense attention from physicists. Sometimes, this attention takes the form of seeing biofilms as inspiration for new physics. In this roadmap, we highlight the work of those who have taken the opposite strategy: we highlight the work of physicists and physical scientists who use physics to engage fundamental concepts in bacterial biofilm microbiology, including adhesion, sensing, motility, signaling, memory, energy flow, community formation and cooperativity. These contributions are juxtaposed with microbiologists who have made recent important discoveries on bacterial biofilms using state-of-the-art physical methods. The contributions to this roadmap exemplify how well physics and biology can be combined to achieve a new synthesis, rather than just a division of labor.

Wong, Gerard C. L. (ORCID:0000000308936383)↗

DOE OSTI2021

Homotypic aggregates contribute to heterogeneity in B cell fates due to an intrinsic gradient of stimulant exposure

Highlights: • CD40-signaled B cells in cultures form three-dimensional homotypic aggregates. • These tight aggregates hinder the free diffusion of large molecule stimulants. • The resulting concentration gradient of stimulants leads to heterogeneous cell fates. • We describe a flow cytometric approach to quantify this positional information. Monocultures of several cell types result in the formation of robust clusters called homotypic aggregates (HAs). How this physical aggregation affects cell fates in immune cell cultures, is poorly understood. We studied anti-CD40-stimulated primary B cell cultures, where cells assembled into large three-dimensional LFA1-driven HAs by 72 h. The dense packing in these aggregates restricts the infiltration of stimulants, such as antibodies, to cells inside the clusters. This creates a concentration gradient of stimulant availability across the cross-section of HAs. We describe a method to retain this positional information even after the disruption of HAs, for analysis by flow cytometry. Comparison of stage-specific cell-surface markers showed that the extent of stimulant-binding affected multiple fates non-uniformly. While germinal center and lineage markers were moderately upregulated, immunoglobulins and markers associated with memory were more than doubled in the peripheral cells binding more anti-CD40. These cells also experienced a strong repression of the plasma cell regulator Prdm1 and an upregulation of the oncogene Myc. Thus, cells at different locations in HAs are subjected to unequal doses of stimulants, leading to a hitherto unreported source of heterogeneity in cell fates. These findings can be extrapolated to understand the dose-dependent effects of stimulants in other three-dimensional cell clusters.

60 APPLIED LIFE SCIENCES↗