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 1,747 records · Page 97

An accuracy assessment of Cartesian-mesh approaches for the Euler equations

A critical assessment of the accuracy of Cartesian-mesh approaches for solving the Euler equations is made. An exact solution of the Euler equations (Ringleb's flow) is used not only to infer the order of error of the Cartesian mesh approaches, but also to compare the magnitude of the error directly to that obtained with a structured mesh approach. The effect of cell merging is investigated as well as the use of two different K-exact reconstruction procedures. The solution methodology of the schemes is explained and tabulated results are presented to compare the solution accuracies. Adaptive and uniform mesh refinement is evaluated for Ringleb's flow and the supersonic flow through an axisymmetric inlet.

Coirier, William J.↗

Numerical simulation of unsteady incompressible viscous flows in generalized coordinate systems

Several numerical solutions of the three dimensional unsteady incompressible Navier-Stokes equations in generalized coordinate systems are presented. The governing equations are discretized by finite volumes with special care to the accurate approximation of the geometric quantities. The unknowns are the pressure and the volume fluxes over the computational cell faces. This formulation results in a robust fractional step solution method for solving discrete equations. Although this method is formulated for the three dimensional case, only two dimensional unsteady results are given. Results are presented for a lid driven two dimensional cavity flow at Reynolds number of 10,000, and for the flow over a circular cylinder with vortex shedding for several Reynolds numbers in the range 100 less than Re less than 1000.

Rosenfeld, Moshe↗

FTIR Spectroscopy of HNO3 and NO2 Relevant to Stratospheric Wake Analysis

The Fourier Transform Infrared (FTIR) technique has been employed to measure absolute concentrations of nitric acid (HNO3) and nitrogen dioxide (NO2) with 1/cm resolution and an absorption pathlength of 4 m under quasi-static and flow conditions at atmospheric pressure and room temperature. Water features seen under quasi-static conditions diminished in intensity under flowing conditions. Nitric acid was observed in the 1660-1760/cm range, while nitrogen dioxide was detected both in the 1536-1660 and 1213-1400/cm ranges. Concentrations of nitrogen dioxide and nitric acid were determined to be 11.9 and 4.35 parts per million (ppm), respectively, with an uncertainty of 0.2 ppm. Experiments are underway with a 10 m cell to measure the absorption of nitric acid, water, sulfur dioxide, hydrochloric acid and ammonia on various materials such as glass, teflon, stainless steel and aluminum used for implementation of the flow system. Such materials will be used for the measurements of stratospheric trace gases by the Quartz Crystal Microbalance (QCM) and Surface Acoustic Wave (SAW) devices.

Abina, Rafiu A.↗

Solar Supergranulation Revealed as a Superposition of Traveling Waves

40 years ago two new solar phenomena were described: supergranulation and the five-minute solar oscillations. While the oscillations have since been explained and exploited to determine the properties of the solar interior, the supergranulation has remained unexplained. The supergranules, appearing as convective-like cellular patterns of horizontal outward flow with a characteristic diameter of 30 Mm and an apparent lifetime of 1 day, have puzzling properties, including their apparent superrotation and the minute temperature variations over the cells. Using a 60-day sequence of data from the MDI (Michelson-Doppler Imager) instrument onboard the SOHO (Solar and Heliospheric Observatory) spacecraft, we show that the supergranulation pattern is formed by a superposition of traveling waves with periods of 5-10 days. The wave power is anisotropic with excess power in the direction of rotation and toward the equator, leading to spurious rotation rates and north-south flows as derived from correlation analyses. These newly discovered waves could play an important role in maintaining differential rotation in the upper convection zone by transporting angular momentum towards the equator.

Gizon, L.↗

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↗

Relativistic particle acceleration in plerions

We discuss recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is either pure electrons and positrons or primarily electrons and positrons with an admixture of heavy ions. Particle-in-cell simulation techniques as well as analytic theory have been used to show that such shocks in pure pair plasmas are fully thermalized -- the downstream particle spectra are relativistic Maxwellians at the temperature expected from the jump conditions. On the other hand, shocks containing heavy ions which are a minority constituent by number but which carry most of the energy density in the upstream medium do put approximately 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) varies as E(exp -2), where N(E)dE is the number of particles with energy between E and E+dE. The mechanism of thermalization and particle acceleration is found to be synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front. The synchrotron maser modes radiated by the heavy ions are absorbed by the pairs at their (relativistic) cyclotron frequencies, allowing the maximum energy achievable by the pairs to be gamma(sub +/-)m(sub +/-)c squared = m(sub i)c squared gamma(sub 1)/Z(sub i), where gamma(sub 1) is the Lorentz factor of the upstream flow and Z(sub i) is the atomic number of the ions. The shock's spatial structure is shown to contain a series of 'overshoots' in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. This shock model is applied to an interpretation of the structure of the inner regions of the Crab Nebula, in particular to the 'wisps,' surface brightness enhancements near the pulsar. We argue that these surface brightness enhancements are the regions of magnetic overshoot, which appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure gamma(sub 1) and a number of other properties of the pulsar's wind. We also discuss applications of the shock theory to the termination shocks of the winds from rotation-powered pulsars embedded in compact binaries. We show that this model adequately accounts for (and indeed predicted) the recently discovered X-ray flux from PSR 1957+20, and we discuss several other applications to other examples of these systems.

Arons, Jonathan↗

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↗

NASA Tech Briefs, August 2004

Topics covered include: Data Relay Board with Protocol for High-Speed, Free-Space Optical Communications; Software and Algorithms for Biomedical Image Data Processing and Visualization; Rapid Chemometric Filtering of Spectral Data; Prioritizing Scientific Data for Transmission; Determining Sizes of Particles in a Flow from DPIV Data; Faster Processing for Inverting GPS Occultation Data; FPGA-Based, Self-Checking, Fault-Tolerant Computers; Ultralow-Power Digital Correlator for Microwave Polarimetry; Grounding Headphones for Protection Against ESD; Lightweight Stacks of Direct Methanol Fuel Cells; Highly Efficient Vector-Inversion Pulse Generators; Estimating Basic Preliminary Design Performances of Aerospace Vehicles; Framework for Development of Object-Oriented Software; Analyzing Spacecraft Telecommunication Systems; Collaborative Planning of Robotic Exploration; Tools for Administration of a UNIX-Based Network; Preparing and Analyzing Iced Airfoils; Evaluating Performance of Components; Fuels Containing Methane of Natural Gas in Solution; Direct Electrolytic Deposition of Mats of MnxOy Nanowires; Bubble Eliminator Based on Centrifugal Flow; Inflatable Emergency Atmospheric-Entry Vehicles; Lightweight Deployable Mirrors with Tensegrity Supports; Centrifugal Adsorption Cartridge System; Ultrasonic Apparatus for Pulverizing Brittle Material; Transplanting Retinal Cells using Bucky Paper for Support; Using an Ultrasonic Instrument to Size Extravascular Bubbles; Coronagraphic Notch Filter for Raman Spectroscopy; On-the-Fly Mapping for Calibrating Directional Antennas; Working Fluids for Increasing Capacities of Heat Pipes; Computationally-Efficient Minimum-Time Aircraft Routes in the Presence of Winds; Liquid-Metal-Fed Pulsed Plasma Thrusters; Personal Radiation Protection System; and Attitude Control for a Solar-Sail Spacecraft.

Source record↗

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↗

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↗

Flow resolution and domain of influence in rarefied hypersonic blunt-body flows

The study assesses the effects of the upstream domain size and grid resolution upon flow properties and body aerodynamics computed for rarefied flows over cold blunt bodies with a direct simulation Monte Carlo (DSMC) particle method. Empirical correlations are suggested for aerodynamic coefficients for two-dimensional flows past a perpendicular flat plate. Free-stream parameters which were varied in the study include the Mach number, Knudsen number, surface temperature, and intermolecular potential. Insufficient grid resolution leads to overprediction of aerodynamic heating and forces in the DSMC method. Solution accuracy correlates well with the Reynolds number defined at the wall temperature and the stagnation mean free path relative to the cell dimension. Insufficient upstream domain size in the DSMC method leads to overprediction of heating and drag. Errors in aerodynamic coefficients correlate well with the distance ahead of the body where flow temperature reaches half of its peak value. Simulation of a hard-sphere gas is more sensitive to grid resolution, while simulation of a Maxwell gas is more sensitive to upstream domain size.

Haas, Brian L.↗

A Comparison of NO Laser-Induced Fluorescence Models at Conditions Relevant to Supersonic and Hypersonic Flows

Planar laser-induced fluorescence (PLIF) of the nitric oxide (NO) molecule has been widely used in wind tunnel facilities for flow visualization, velocity, and temperature measurements. The experimental PLIF measurements are often compared with synthetic PLIF images using computationally derived temperatures, pressures, velocities, and species mole fractions. This approach is commonly referred to as computational flow imaging (CFI). In the present work, we compare signal intensity from PLIF models with experimental PLIF measurements obtained within a low pressure gas cell system at pressures and NO mole fractions relevant to supersonic and hypersonic flowfields. Experimental measurements were compared to several different laser induced-fluorescence models reported in the literature including LIFBASE, LINUS, and a NASA two-level model. The experimental measurements agreed well with all of the models at lower pressures and lower NO mole fractions; the fluorescence there is linear with both of these parameters. However, at higher pressures and mole fractions, the signal becomes nonlinear with respect to these parameters as self-quenching limits the signal and absorption further limits the signal. In fact, for the experimental path length of the experiment, the combination of high pressure and high NO mole fraction causes the experimental results to deviate significantly from the predicted results that neglect absorption of the incident laser sheet. The LINUS model, which allows absorption to be calculated, provided results that agreed better with the experimental measurements. Since supersonic and hypersonic flowfields may contain a region of the flow with high pressures and measurements in large-scale facilities often include a long path length, neglecting absorption may have a significantly negative effect on the CFI comparison to experimental PLIF images. As a result, PLIF models that account for absorption should be included in computational flow imaging approaches for laser induced fluorescence.

laser-induced fluorescence↗

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↗

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↗

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↗

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.↗

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↗

Plasma engineering

Flow of conducting liquid in annular gap, effects of high fields on surface charge distributions, and materials study of silicon solar cells

Pollack, S.↗