Search NASA⌕ Search

SEARCH · Search NASA

Results for “CHANNEL 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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Additively manufactured cryogenic microchannel distillation device for air separation

Abstract The efficiency of air separation is tested using three different small‐scale cryogenic distillation columns. The performance of a random packed column is compared to the performance of two microchannel distillation (MCD) columns that use thin wicking structures and gas flow channels to achieve process intensification. The MCD columns tested include a plate‐type layered column and an additively manufactured porous honeycomb (AMPH) column. For columns with 25.4 cm of active height and run under similar conditions, the packed, plate‐type layering, and AMPH columns achieved approximate height equivalent of a theoretical plate (HETP) values of 5.5, 3.7, and 3.2 cm for nitrogen, and 5.9, 4.9, and 3.3 cm for argon. The AMPH column can produce up to 0.4 SLM of more than 90% purity oxygen with 12 W of cooling lift. These results demonstrate the feasibility of using additive manufacturing to construct MCD devices and pave a way for constructing novel MCD designs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Coupled Thermo-Hydro-Mechanical Processes in Fractured Rocks: Some Past Scientific Highlights and Future Research Directions

Abstract Coupled thermo-hydro-mechanical (THM) processes in fractured rocks have been a topic of intense scientific research for more than 30 years. The present paper takes a look into the past and highlights some scientific advances which are of an unusual “out-of-the-box” nature, and then looks forward and discusses possible directions of future research in this interesting field of study. Concerning future research directions, we see a trend from a focus on coupled THM processes in single fractures or a few interacting fractures, to the study of coupled THM behavior in complex fracture network systems where the fractures act collectively giving rise to local stress concentration points and points of large pressure gradients. Three examples of future research directions are presented. First is an effort towards identifying characterizing parameters of a fracture network that play a direct controlling role in major coupled THM phenomena (such as induced seismicity and flow channeling), rather than parameters of stochastic distributions of fractures in the network. The second example of research direction is accounting for the heterogeneity and hierarchy of fractures in a fault or fracture zone which has been associated with major THM events in a number of geo-energy projects. The third example is at the opposite end of the first; here it is recognized that in some cases, the coupled THM processes in fractured rocks may be controlled dominantly by only a few key bridges. Identification, characterization, and evaluation of these key bridges should be one of the important research directions in the coming days.

58 GEOSCIENCES↗
DOE OSTI2024

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗
DOE OSTI2024

Reynolds-number scaling of wall-pressure–velocity correlations in wall-bounded turbulence

Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure–velocity coupling requires investigation to properly design a controller with such input data so that it can actuate upon the desired turbulent structures. A comprehensive database from direct numerical simulations (DNS) of turbulent channel flow is used for this purpose, spanning a Reynolds-number range$Re_\tau \approx 550\unicode{x2013}5200$. Spectral analysis reveals that the streamwise velocity is most strongly coupled to the linear term of the wall pressure, at a Reynolds-number invariant distance-from-the-wall scaling of$\lambda _x/y \approx 14$(and$\lambda _x/y \approx 8$for the wall-normal velocity). When extending the analysis to both homogeneous directions in$x$and$y$, the peak coherence is centred at$\lambda _x/\lambda _z \approx 2$and$\lambda _x/\lambda _z \approx 1$for$p_w$and$u$, and$p_w$and$v$, respectively. A stronger coherence is retrieved when the quadratic term of the wall pressure is concerned, but there is only little evidence for a wall-attached-eddy type of scaling. An experimental dataset comprising simultaneous measurements of wall pressure and velocity complements the DNS-based findings at one value of$Re_\tau \approx 2$k, with ample evidence that the DNS-inferred correlations can be replicated with experimental pressure data subject to significant levels of (acoustic) facility noise. It is furthermore shown that velocity-state estimations can be achieved with good accuracy by including both the linear and quadratic terms of the wall pressure. An accuracy of up to 72 % in the binary state of the streamwise velocity fluctuations in the logarithmic region is achieved; this corresponds to a correlation coefficient of$\approx$0.6. This thus demonstrates that wall-pressure sensing for velocity-state estimation – e.g. for use in real-time control of wall-bounded turbulence – has merit in terms of its realization at a range of Reynolds numbers.

Mechanics↗
DOE OSTI2020

The Sensitivity of Hyporheic Exchange to Fractal Properties of Riverbeds

Hyporheic exchange in riverbeds is driven by current-bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large-scale models. Further, they can also be used to design artificial topographies that target hyporheic ecosystem services.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Permeability Hysteresis From Microchannels Opening During Dissolution/Reprecipitation Cycle

Abstract Permeability is a critical parameter for geological resources characterization. Its evolution with respect to porosity is particularly interesting and many research initiatives focus on deriving such relationships, to understand some hydraulic impacts of microstructure alteration. Permeability evolution from chemical reactions, for instance, can become complex as flow channels may open during rock dissolution. In this contribution, we show that this phenomenon can lead to irregular porosity‐permeability curves and permeability hysteresis after reprecipitation. Current approaches describing permeability as a simple function of porosity can therefore not capture this behavior, and we advocate instead the use of dynamic modeling for such scenarios. We demonstrate our approach by modeling a dissolution/precipitation cycle for a unit cell pore channel and quantify the process at larger scale on three different rock samples, whose microstructures are reconstructed from segmented micro‐Computerised Tomography scans.

Lesueur, Martin↗
DOE OSTI2021

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗
DOE OSTI2024

Mixing of the Connecticut River Plume During Ambient Flood Tides: Spatial Heterogeneity and Contributions of Bottom–Generated and Interfacial Mixing

The Connecticut River plume is influenced by energetic ambient tides in the Long Island Sound receiving waters. The objectives of this modeling study are (a) characterizing the spatial heterogeneity of turbulent buoyancy fluxes, (b) partitioning turbulent buoyancy fluxes into bottom-generated and interfacial shear contributions, and (c) quantifying contributions to plume-integrated mixing within the tidal plume. The plume formed during ambient flood tides under low river discharge, spring tides, and no winds is analyzed. Turbulent buoyancy fluxes (B) and depth-integrated B through the plume (Bd) are characterized by pronounced spatial heterogeneity. Strong mixing (Bd ~ 10 –5 -10 –4 m 3 /s 3 ) occurs near the mouth, in the nearfield plume turning region, over shoals, and nearshore shallow areas. Low to moderate mixing (Bd ~ 10 –8 -10 –6 m 3 /s 3 ) occupies half the plume. Buoyancy fluxes are first partitioned based on the depth of the shear stress minimum between plume-generated and bottom-generated shear maxima. Four other tested partitioning methods are based on open channel flow and stratified shear flow parameterizations. Interfacial and bottom-generated shear contribute to different areas of intense and moderate mixing. All methods indicate a significant plume mixing role for bottom-generated mixing, but interfacial mixing is a bigger contributor. Plume-integrated total and interfacial mixing peak at max ambient flood and the timing of peak bottom-generated mixing varies among partitioning methods. Two-thirds of the mixing occurs in concentrated intense mixing areas. A parameter space with the ambient tidal Froude number and plume thickness to depth ratio as axes indicates many tidally modulated plumes are moderately to dominantly influenced by bottom-generated tidal mixing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Significance of Low‐Velocity Zones on Solute Retention in Rough Fractures

Natural fractures are characterized by high internal heterogeneity. This internal variability is the cause of flow channeling, which in turn leads to contaminant transport taking place primarily along the high-velocity channels. Mass exchange between the high-velocity channels and the low-velocity zones has the potential to enhance contaminant retention, due to solute diffusion into the low-velocity zones and subsequent exposure to additional surface area for diffusion into the bordering rock matrix. Here, we derive a random walk particle tracking method for heterogeneous fractures, which includes an additional term to account for the aperture gradient. The method takes into account advection, diffusion in the fracture and matrix diffusion. The developed numerical framework is applied to assess the effect of low-velocity zones in rough self-affine fractures. The results show that diffusion into low-velocity zones has a visible but modest impact on contaminant retention. The magnitude of this impact does not change considerably, regardless of whether diffusion into the rock matrix is considered in the model, and increases for a decreasing average Péclet number of the fracture.

58 GEOSCIENCES↗
DOE OSTI2025

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Relevance of Tritium Breeder Irradiation Testing in a Fusion Prototypic Neutron Source

A fusion prototypic neutron source (FPNS) will play a fundamental role in testing and developing fusion materials that can withstand the harsh fusion environment. We conducted a thorough assessment to determine whether blanket materials would benefit from FPNS testing. Because of the cross sections for the various blanket materials, lithium-containing materials can generally be tested adequately using existing fission irradiation capabilities instead of waiting for or relying solely on future FPNS experiments. However, non-lithium blanket materials, such as beryllium multipliers and flow channel inserts, should be tested under a fusion neutron spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Exploration of non-resonant divertor features on the Compact Toroidal Hybrid

Abstract Non-resonant divertors (NRDs) separate the confined plasma from the surrounding plasma facing components (PFCs). The resulting striking field line intersection pattern on these PFCs is insensitive to plasma equilibrium effects. However, a complex scrape-off layer (SOL), created by chaotic magnetic topology in the plasma edge, connects the core plasma to the PFCs through varying magnetic flux tubes. The Compact Toroidal Hybrid (CTH) serves as a test-bed to study this by scanning across its inductive current. Simulations observe a significant change of the chaotic edge structure and an effective distance between the confined plasma and the instrumented wall targets. The intersection pattern is observed to be a narrow helical band, which we claim is a resilient strike line pattern. However, signatures of finger-like structures, defined as heteroclinic tangles in chaotic domains, within the plasma edge connect the island chains to this resilient pattern. The dominant connection length field lines intersecting the targets are observed via heat flux modeling with EMC3-EIRENE. At low inductive current levels, the excursion of the field lines resembles a limited plasma wall scenario. At high currents, a private flux region is created in the area where the helical strike line pattern splits into two bands. These bands are divertor legs with distinct SOL parallel particle flow channels. The results demonstrate the NRD strike line pattern resiliency within CTH, but also show the underlying chaotic edge structure determining if the configuration is diverted or limited. This work supports future design efforts for a mechanical structure for the NRD.

Physics↗
DOE OSTI2022

Marginally stable current sheets in collisionless magnetic reconnection

Noncollisional current sheets that form during the nonlinear development of spontaneous magnetic reconnection are characterized by a small thickness, of the order of the electron skin depth. They can become unstable to the formation of plasmoids, which allows the magnetic reconnection process to reach high reconnection rates. In this work, we investigate the marginal stability conditions for the development of plasmoids when the forming current sheet is purely collisionless and in the presence of a strong guide field. We analyze the geometry that characterizes the reconnecting current sheet, and what promotes its elongation. Once the reconnecting current sheet is formed, we identify the regimes for which it is plasmoid unstable. Finally, our study shows that plasmoids can be obtained, in this context, from current sheets with an aspect ratio much smaller than in the collisional regime, and that the plasma flow channel of the marginally stable current layers maintains an inverse aspect ratio of 0.1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparing Experimental Results for Large Particle Separation from Non-Newtonian Slurries Using Full and Tapered Bump Arrays

To separate particles from carrier fluids arrays of staggered posts inserted into the flow stream can be employed. The smaller particles follow the flow stream while the larger particles move laterally to one side. During flow, the concentration of particulate increases in one direction and decreases in the other crosswise direction. This technique separates the particles from the carrier fluid. Experiments using a slurry with non-Newtonian rheology were conducted to evaluate large particle separation from the non-Newtonian carrier fluid using bump arrays. The bentonite kaolin clay slurry with non-Newtonian rheology was spiked with large diameter inert glass particles. Experiments with flow through lateral displacement arrays evaluated the performance of particle separation using two array configurations: a full array with all posts and a tapered array with posts removed in a triangular portion of the array past the rows of anticipated particle of separation. The posts are staggered posts to promote particle segregation to one side of the flow channel. Using this type of array for removing particles from non-Newtonian yield stress slurries is novel and we present unique results. In this paper, the performance of full and tapered arrays are compared. These results have many industrial applications including removing particles from slurries of nuclear waste. Large particle removal is an important step in waste processing.

slurry, solids removal, deterministic lateral disp↗
DOE OSTI2022

Image transformers for classifying acute lymphoblastic leukemia

Cancer is the leading cause of death by disease in American children. Each year, nearly 16,000 children in the United States and over 300,000 children globally are diagnosed with cancer. Leukemia is a form of blood cancer that originates in the bone marrow and accounts for one-third of pediatric cancers. This disease occurs when the bone marrow contains 20% or more immature white blood cell blasts. Acute lymphoblastic leukemia is the most prevalent leukemia type found in children, with half of all annual cases in the U.S. diagnosed for subjects under 20 years of age. To diagnose acute lymphoblastic leukemia, pathologists often conduct a morphological bone marrow assessment. This assessment determines whether the immature white blood cell blasts in bone marrow display the correct morphological characteristics, such as size and appearance of nuclei. Pathologists also use immunophenotyping via multi-channel flow cytometry to test whether certain antigens are present on the surface of blast cells; the antigens are used to identify the cell lineage of acute lymphoblastic leukemia. These manual processes require well-trained personnel and medical professionals, thus being costly in time and expenses. Computerized decision support via machine learning can accelerate the diagnosis process and reduce the cost. Training a reliable classification model to distinguish between mature and immature white blood cells is essential to the decision support system. Here, we adopted the Vision Transformer model to classify white blood cells. The Vision Transformer achieved superb classification performance compared to state-of-the-art convolutional neural networks while requiring less computational resources for training. Additionally, the latent self-attention architecture provided attention maps for a given image, providing clues as to which portion(s) of the image were significant in decision-making. We applied the Vision Transformer model and a convolutional neural network model to an acute lymphoblastic leukemia classification dataset of 12,528 samples and achieved accuracies of 88.4% and 86.2%.

Cho, Priscilla↗