Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow cell”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Assessment and validation of NEAMS tools for high-fidelity multiphysics transient modeling of microreactors: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The NEAMS Multiphysics Applications team aims at providing assessment of code useability and functionality for microreactor design and analyses, together with demonstration of their capabilities to properly capture the steady-state and time-dependent behavior of different microreactor concepts. In FY-24, significant progress was achieved in improving multi-physics models of several microreactors systems: HP-MR, GC-MR and KRUSTY. These efforts focused on solving more complex multiphysics problems enabled by enhanced tools capability, verifying and validating results obtained, providing feedback to developers for suggested improvements, and sharing these models to facilitate user training. A series of new multiphysics transients were completed on the HP-MR (using Griffin/BISON/Sockeye) with core startup transient, control drum inadvertent rotation accident, and hydrogen leakage from hydride moderator (also including SWIFT). On the GC-MR, a new full-core model was developed and analyzed through a series of new multiphysics (Griffin/BISON/SAM) transients to simulate moderator leakage (also including SWIFT), flow blockage and coolant depressurization. Additional and updated TRISO failure analyses were completed on the HP-MR unit-cell and GC-MR assembly models leveraging improved TRISO modeling capabilities. The amount of SiC failure following accidental transients at end-of-life was null. However, GC-MR assembly TRISO analysis highlighted Pd penetration rate can be problematic and may require design changes on the studied microreactor concept. The neutronics discrepancies observed on the KRUSTY model in previous years were resolved using hybrid set of Monte Carlo/Deterministic cross-sections. The multiphysics (Griffin neutronics / BISON thermal-mechanics) 15₵ insertion transient simulation displayed good agreement when comparing with experimental data. Initial modeling of the 30 ₵ reactivity insertion also displays promising results. Such close agreement provides important validation data that can be leveraged by the NEAMS program and by microreactor vendors to support licensing of their technology. Finally, important experience was gathered with the NEAMS tools leading to several user feedback shared with tools developers, especially with regards to MOOSE mesh generator and Griffin. This project led to many publications demonstrating modeling capabilities, and to three models shared on the Virtual Test Bed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A perfused multi-well bioreactor platform to assess tumor organoid response to a chemotherapeutic gradient

There is an urgent need to develop new therapies for colorectal cancer that has metastasized to the liver and, more fundamentally, to develop improved preclinical platforms of colorectal cancer liver metastases (CRCLM) to screen therapies for efficacy. To this end, we developed a multi-well perfusable bioreactor capable of monitoring CRCLM patient-derived organoid response to a chemotherapeutic gradient. CRCLM patient-derived organoids were cultured in the multi-well bioreactor for 7 days and the subsequently established gradient in 5-fluorouracil (5-FU) concentration resulted in a lower IC 50 in the region near the perfusion channel versus the region far from the channel. We compared behaviour of organoids in this platform to two commonly used PDO culture models: organoids in media and organoids in a static (no perfusion) hydrogel. The bioreactor IC 50 values were significantly higher than IC 50 values for organoids cultured in media whereas only the IC 50 for organoids far from the channel were significantly different than organoids cultured in the static hydrogel condition. Using finite element simulations, we showed that the total dose delivered, calculated using area under the curve (AUC) was similar between platforms, however normalized viability was lower for the organoid in media condition than in the static gel and bioreactor. Our results highlight the utility of our multi-well bioreactor for studying organoid response to chemical gradients and demonstrate that comparing drug response across these different platforms is nontrivial.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗
DOE OSTI2025

Evaluating inkjet dispenser/liquid vortex capture-mass spectrometry for single-cell metabolomics in Hep G2 steatosis caused by tamoxifen

Single-cell mass spectrometry (MS) is advancing our understanding of metabolic pathways in heterogeneous cell populations; however, many techniques are slow or require disruptive sample preparations. This study evaluated coupling a modified HP D100 single-cell inkjet dispenser with liquid vortex capture-mass spectrometry (D100/LVC-MS). The D100 is a single-cell inkjet dispenser capable of titrating solutions and isolating single cells via disposable cassettes equipped with microfluidic channels and an impedance sensor. The LVC-MS enables high-throughput capture, lysis, and ionization of analytes for mass spectrometric analysis. The D100/LVC-MS system was characterized through titration and single-cell experiments. Propranolol titration demonstrated linearity across a broad concentration range using the D100/LVC-MS system. Additionally, Hep G2 hepatocarcinoma cells and Chlamydomonas reinhardtii algae were used to showcase the D100’s high-throughput or low-buffer-volume single-cell dispensing strategies. The D100/LVC-MS system’s performance was validated by evaluating tamoxifen-induced steatosis in Hep G2 cells. Tamoxifen, associated with nonalcoholic fatty liver disease in breast cancer patients following long-term use, was tested in Hep G2 cells at 20 µM for 72 h against DMSO-treated controls. High-throughput analysis of 500 cells per condition, completed in 25 min per run, demonstrated the system’s efficiency. The D100/LVC-MS system simultaneously quantified tamoxifen and measured triglycerides and phosphatidylcholines. Triglycerides were upregulated in the tamoxifen-treated cells and the results indicated two distinct cell populations, differing in tamoxifen and phosphatidylcholines levels, suggesting heterogeneity within the treated population. In conclusion, these findings highlight the D100/LVC-MS system as a cost-effective, high-throughput platform for single-cell metabolomics and lipidomics, with significant potential for evaluating metabolic alterations.

Chromatography↗
DOE OSTI2020

In Situ Structural and Electrical Conductivity Characterization of Sr2MMoO6-δ Double Perovskite Solid Oxide Fuel Cell Anode Materials

In this study, a series of Sr 2 MMoO 6-δ double perovskite solid oxide fuel cell (SOFC) anode materials (where 0 <= δ <= 0.041 and M = Fe, Co, or Ni) were synthesized, and the changes in their morphologies under relevant SOFC operating conditions were explored. Ultra-small-angle X-ray scattering (USAXS), small-angle X-ray scattering (SAXS), and wide-angle X-ray scattering (WAXS) were used to determine changes in the microstructures and phase compositions of the anode materials at high temperatures and under a reducing atmosphere. The stability of the double perovskite structure was found to be highly dependent on the identity of the cation M, such that when M = Fe, the material remained stable over the course of the experiments. However, when M = Co or Ni, significant changes in the microstructure and phase composition were observed. An in situ study of the M = Co sample, in which electrical conductivity and USAXS/SAXS/WAXS measurements were conducted simultaneously, revealed the structural degradation mechanisms and electrical conductivity changes over a range of temperatures. To conduct these measurements, a cell was developed that allowed for a sample mounted on Pt wires to be placed in the X-ray beam and heated under gas flow comprising 4% mass H 2 and 96% mass N 2 . The resulting measurements allowed for the direct comparison of the electrical and morphological changes occurring in the material under operating conditions, such that increases in conductivity could be attributed to the growth of new phases.

25 ENERGY STORAGE↗
DOE OSTI2023

Characterization of Lactobacilli Phage Endolysins and Their Functional Domains–Potential Live Biotherapeutic Testing Reagents

Phage endolysin-specific binding characteristics and killing activity support their potential use in biotechnological applications, including potency and purity testing of live biotherapeutic products (LBPs). LBPs contain live organisms, such as lactic acid bacteria (LAB), and are intended for use as drugs. Our approach uses the endolysin cell wall binding domains (CBD) for LBP potency assays and the endolysin killing activity for purity assays. CBDs of the following five lactobacilli phage lysins were characterized: CL1, Jlb1, Lj965, LL-H, and ΦJB. They exhibited different bindings to 27 LAB strains and were found to bind peptidoglycan or surface polymers. Flow cytometry based on CBD binding was used to enumerate viable counts of two strains in the mixture. CL1-lys, jlb1-lys, and ΦJB-lys and their enzymatic domains (EADs) exhibited cell wall digestive activity and lytic activity against LAB. Jlb1-EAD and ΦJB-EAD were more sensitive than their respective hololysins to buffer pH and NaCl changes. The ΦJB-EAD exhibited stronger lytic activity than ΦJB-lys, possibly due to ΦJB-CBD-mediated sequestration of ΦJB-lys by cell debris. CBD multiplex assays indicate that these proteins may be useful LBP potency reagents, and the lytic activity suggests that CL1-lys, jlb1-lys, and ΦJB-lys and their EADs are good candidates for LBP purity reagent development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Holistic design principles for flow batteries: Cation dependent membrane resistance and active species solubility

Cation dependent resistance of a commercial cation exchange membrane, Nafion™ 212, as well as the solubility of select active materials are investigated, demonstrating practical consequences of cation choice on redox flow battery (RFB) performance. The conductivity of alkali and methyl- or -alcohol substituted ammonium cations through pretreated Nafion™ 212 reveal that cation size is the dominant factor affecting ionic transport. Among the alkali series the resistance increases with larger ionic radii and similarly, with increasing steric bulk through substitution on ammonia. The solubility of ferrocyanide, a commonly used redox flow battery posolyte (or catholyte), as well as iron, cobalt, and chromium metal-organic compounds are determined for the same series of cations. Full cells are assembled using various cation mixtures and membrane pre-treatments to demonstrate practical effects of each variable. Here, we highlight the need to carefully select the correct cations for efficient electrolyte design, since certain cations may lead to high conductivity but low solubility, or vice versa, and provide design principles for RFBs using Nafion™ 212.

25 ENERGY STORAGE↗
DOE OSTI2022

A weighted state redistribution algorithm for embedded boundary grids

State redistribution is an algorithm that stabilizes cut cells for embedded boundary grid methods. This work extends the earlier algorithm in several important ways. First, state redistribution is extended to three spatial dimensions. Second, we discuss several algorithmic changes and improvements motivated by the more complicated cut cell geometries that can occur in higher dimensions. In particular, we introduce a weighted version with less dissipation in an easily generalizable framework. Third, we demonstrate that state redistribution can also stabilize a solution update that includes both advective and diffusive contributions. Notably, the stabilization algorithm is shown to be effective for incompressible as well as compressible reacting flows. Finally, we discuss the implementation of the algorithm for several exascale-ready simulation codes based on AMReX, demonstrating ease of use in combination with domain decomposition, hybrid parallelism and complex physics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗
DOE OSTI2022

Spatially resolved performance and degradation in a perfluorinated anion exchange membrane fuel cell

Anion exchange membrane fuel cells may enable future operation with non-precious metal-based catalysts. These systems have a delicate sensitivity to operating conditions such as humidification levels and the presence of CO2 in the air oxidant stream. We present spatially resolved in-situ performance results that shed light on phenomena that are unique to anion exchange membrane fuel cells. For cell construction, a highly conductive perfluorinated anion exchange polymer was used as the membrane and the material in powder form as the ionomer. Experiments were conducted to investigate the effects of humidification, fuel/oxidant concentration, and carbonation effects on the performance and its distribution in the cell. The results indicated that (i) dry conditions at the cathode have a stronger effect than at the anode on overall cell performance, (ii) performance significantly suffered when humidification was below 90%, (iii) fuel and oxidant dilution effects lead mass-transport losses and were stronger than flow rate effects, (iv) CO2 in the cathode feed stream creates an equilibration disparity between the inlet and outlet sections and CO2 purging is affected by flooding conditions, and (v) after >500 h of operation, performance deteriorates predominantly at the inlet.

2+1D fuel cell model↗
DOE OSTI2021

Towards performance portability in the Spark astrophysical magnetohydrodynamics solver in the Flash-X simulation framework

Simulations of core-collapse supernovae, and other astrophysical phenomena, are quintessential extreme-scale computing challenges. For core-collapse supernova simulations to be carried out by the ExaStar project under the Exascale Computing Project umbrella, a robust, efficient, and state-of-the-art magnetohydrodynamics solver is a critical requirement. In Flash-X, the primary software instrument for ExaStar, a new magnetohydrodynamics solver has been designed and implemented from the ground up to achieve accuracy and efficiency for simulations of complex astrophysical flows. This new solver, dubbed Spark, uses high-order spatial reconstruction, Runge-Kutta time integration, and an efficient cell-centered approach to satisfying the divergence-free condition for the magnetic fields. Spark was written to be optimized for data locality in cache hierarchy of CPUs. Since data locality optimizations for cache hierarchy are not directly compatible with those of accelerators, we have taken the approach of using program synthesis to avoid massive amounts of code replication that would be necessary if we were to maintain two different versions of the solver. Our program synthesis relies on a simple key-dictionary approach, implemented in python, that enables us to assemble the version of the solver suitable for the target hardware from code fragments identified by specific keys. In this work, we describe the data locality optimizations of the solver for CPUs and accelerators and the program synthesis tools that enable this portability. We also detail the parallel performance of Spark for both CPUs and accelerators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The effect of blood viscosity on shear‐induced hemolysis using a magnetically levitated shearing device

Abstract Introduction Blood contacting medical devices, including rotary blood pumps, can cause shear‐induced blood damage that may lead to adverse effects in patients. Due in part to an inadequate understanding of how cell‐scale fluid mechanics impact red blood cell membrane deformation and damage, there is currently not a uniformly accepted engineering model for predicting blood damage caused by complex flow fields within ventricular assist devices (VADs). Methods We empirically investigated hemolysis in a magnetically levitated axial Couette flow device typical of a rotary VAD. The device is able to accurately control the shear rate and exposure time experienced by blood and to minimize the effects of other uncharacterized stresses. Using this device, we explored the effects of both hematocrit and plasma viscosity on shear‐induced hemolysis to characterize blood damage based on the viscosity‐independent shear rate, rather than on shear stress. Results Over a shear rate range of 20 000 – 80 000 1/s, the Index of Hemolysis (IH) was found to be dependent upon and well‐predicted by the shear rate alone. IH was independent of hematocrit, bulk viscosity, or the suspension media viscosity and less correlated to shear stress (MSE = 0.46–0.75) than to shear rate (MSE = 0.06–0.09). Conclusion This study recommends that future investigations of shear‐induced blood damage report findings with respect to the viscosity‐neutral term of shear rate, in addition to the bulk whole blood viscosity measured at an appropriate shear rate relevant to the flow conditions of the device.

Krisher, James A.↗
DOE OSTI2023

As a Matter of Dynamical Range – Scale Dependent Energy Dynamics in MHD Turbulence

Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Heat Transfer Characteristics of Particle and Air Flow Through Additively Manufactured Lattice Frame Material Based on Octet-Shape Topology

Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchanger is a critical component in next-generation concentrating solar power (CSP) plants. The inherently low heat transfer between falling particles and sCO 2 imposes a challenge toward economic justification of levelized cost of electricity produced through solar energy. Introduction of integrated porous media with the walls bounding particle flow has the potential to enhance the overall particle-to-sCO 2 heat exchanger performance. Here, this paper presents an experimental study on heat transfer characterization of additively manufactured lattice frame material based on Octet-shaped unit cell with particles and air as working fluids. The lattice structures were additively manufactured in stainless steel (SS) 316L and SS420 (with 40% bronze infiltration) via Binder jetting process, where the lattice porosities were varied between 0.75 and 0.9. The mean particle diameters were varied from 266 μm to 966 μm. The effective thermal conductivity and averaged heat transfer coefficient were determined through steady-state experiments. It was found that the presence of lattice enhances the effective thermal conductivity by 2–4 times when compared to packed bed of particles alone. Furthermore, for gravity-assisted particle flow through lattice panel, significantly high convective heat transfer coefficients ranging from 200 W/m2K to 400 W/m2K were obtained for the range of particle diameters tested. The superior thermal transport properties of Octet-shape-based lattice frame for particle flow makes it a very promising candidate for particle-to-sCO 2 heat exchanger for CSP application.

14 SOLAR ENERGY↗
DOE OSTI2020

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

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

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.

42 ENGINEERING↗
DOE OSTI2022

A method for generating moving, orthogonal, area preserving polygonal meshes

A new method for generating locally orthogonal polygonal meshes from a set of generator points is presented in which polygon areas are a constraint. The area constraint property is particularly useful for particle methods where moving polygons track a discrete portion of material. Because Voronoi polygon meshes have some very attractive mathematical and numerical properties for numerical computation, a generalization of Voronoi polygon meshes was formulated that enforces a polygon area constraint. Area constrained moving polygonal meshes allow one to develop hybrid particle-mesh numerical methods that display some of the most attractive features of each approach. It is shown that this mesh construction method can continuously reconnect a moving, unstructured polygonal mesh in a pseudo-Lagrangian fashion without change in cell area/volume, and the method's ability to simulate various physical scenarios is shown. Overall, the advantages are identified for incompressible fluid flow calculations, with demonstration cases that include material discontinuities of all three phases of matter and large density jumps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗