Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow analysis”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Calibration of reactive burn and Jones-Wilkins-Lee parameters for simulations of a detonation-driven flow experiment with uncertainty quantification

Here, uncertainties in the explosive-specific parameters of the Jones-Wilkins-Lee (JWL) equation of state (EOS) are carefully considered in hydrodynamic simulations of an explosive experiment to minimize the error in the flow prediction. Experimental data of the leading shock position in the transverse direction over time serves as the prediction metric for quantifying simulation prediction error. The uncertainty quantification technique, global sensitivity analysis, is utilized to determine the JWL parameters to which the transverse shock propagation is most sensitive. A polynomial response surface (PRS) is constructed in the space of the most influential JWL parameters, and the point of minimum error between the experimental data and the PRS yields calibrated JWL parameters for the experimental flow. The simulation results following the parameter calibration show good agreement with the experimental data. It was found that two significant parameters, the heat release per unit mass of reactant Q and JWL model exponent R 1 are strongly related, which makes it difficult to identify accurate values.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal Stress Modeling and Analysis of Packed-bed Thermocline Energy Storage Tank for INL Thermal Energy Distribution System (TEDS)

The Thermal Energy Distribution System (TEDS) at Idaho National Laboratory (INL) is a thermal-hydraulic flow loop to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit for TEDS. The packed-bed TES is adopted in TEDS because of its benefit as a low-cost single-tank storage option compared to the traditional two-tank storage. However, thermal ratcheting is one potential design concern which is caused by the rearrangement of granular filler inside a packed-bed tank during continuous thermal cycling operation of the packed-bed TES tank. If the thermally induced stress exceeds yield strength of the tank wall, it may cause catastrophic consequences like rupture of the thermal storage tank. Thus, it is crucial to understand the phenomenon to ensure the robust operation. Based on the temperature boundary conditions given by transient thermal analyses with computational fluid dynamics (CFD) simulations, the thermal ratcheting analysis is then conducted to evaluate the hoop stress and resultant thermal ratcheting potential of the TES tanks with two different modeling approaches: (1) infinite rigidity model and (2) Drucker-Prager (DP) model. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One Plant and evaluate the thermal ratcheting potential of the TEDS TES tank.

25 ENERGY STORAGE↗
DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrated system failure analysis software toolchain (IS-FAST)

Systems and methods are provided for generating faults and analyzing fault propagation and its effects. Starting from the ontologies of components, functions, flows, and faults, systems and methods are provided that describe, generate and track faults in a computer system across multiple domains throughout design and/or development. In order to construct the system and fault models, a series of concepts is introduced in the form of ontologies and their dependencies. An investigation is performed into the faults, including their type, cause, life-cycle aspects, and effect. Principles and rules are created to generate various faults based on system configurations. After the modeling process, a simulation engine is described to execute actions and simulate the process of fault generation and propagation. As a result, fault paths that impact components and functions can be obtained.

Diao, Xiaoxu↗
DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗
DOE OSTI2022

Analysis of solid uranium particulates on cotton swipes with an automated microextraction-ICP-MS system

Here, an automated microextraction method coupled to an inductively coupled plasma – mass spectrometer (ICP-MS) was developed for the direct analysis of solid uranium particulates on the surface of cotton swipes. The microextraction probe extracts particulates from the sample surface, in a flowing solvent, and directs the removed analyte to an ICP-MS for isotopic determination. The automated system utilizes a mechanical XY stage that is software controlled with the capability of saving and returning to specific locations and a camera focused to the swipe surface for optimal viewing of the extracted locations (i.e., material present). Here, particulates (n = 135) were extracted and measured by ICP-MS, including 35 depleted uranyl nitrate hexahydrate (UN) (used for mass bias corrections), 50 uranyl fluoride (UO 2 F 2 ), and 50 uranyl acetate (UAc) particulates. Blank extractions were performed on the cotton swipes between triplicate sample analyses. Between each swipe extraction, the probe was sent between two wells containing 10% and 5% HNO 3 to clean the probe head and to eliminate any analyte carryover between particulates. The measured 235 U/ 238 U and 234 U/ 238 U isotope ratios for the UO 2 F 2 particulates were 0.00725(8) and 0.000054(4), a percent relative difference (% RD) of –0.041% and –1.7% from the reference isotope ratios determined in-lab through multi-collector ICP-MS analysis of dissolved aliquots of the U material. The UAc samples had a measured 235 U/ 238 U isotope ratio of 0.00206(7), a –0.96% relative difference from the reference value of 0.00208(1). The 234 U/ 238 U and 236 U/ 238 U isotope ratios were 0.000008(1) and 0.000031(4), –5.1% RD and –4.3% RD, respectively. The automated sample stage enabled seamless and rapid particle analysis, leading to a significant increase in throughput versus what was previously possible. Additionally, the saved location capability reduced user sampling error as sampling locations were easily stored and recalled. Analysis of U particles on the swipe surface – including blanks, mass bias, and triplicate extractions – was completed in less than an hour without any sample preparation necessary.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Design of an Additively Manufactured Recuperator With 800 °C Inlet Temperature for sCO2 Power Cycle Application

Abstract Recuperators with design temperatures at and above 800 °C can further increase the thermal efficiency of supercritical CO2 power cycles by enabling higher turbine exhaust temperatures. Mar-M247 is a well-suited nickel-based superalloy for high temperature service due to its high creep strength that prevents excessive material thickness being required for pressure containment. Additive manufacturing using a high-speed laser-directed energy deposition (L-DED) process presents a promising solution, with build trials demonstrating the ability to produce nonconventional flow channels for enhanced heat transfer. A design process is presented that includes aerothermal and mechanical evaluation to maximize performance within the constraints of the manufacturing process. A 2-D heat transfer network and pressure drop code allows prediction of flow distribution and its effect on overall thermal performance. Established literature correlations, along with CFD simulation, inform the prediction of heat transfer coefficients and friction factors for the flowpaths and enhancement features in the heat exchanger core. Mechanical evaluation using finite element analysis (FEA) modeling with the intent of the ASME Boiler & Pressure Vessel Code (BPVC) Section VIII, Div. 2 assesses the operational safety of the design. The detailed design features annular finned passages that take advantage of helical flow paths to distribute the flow from separated headers to shared heat transfer surfaces. Performance predictions for the recuperator at a 50 kW scale provide insights into the feasibility of the additively manufactured (AM) process to produce recuperators on a commercial scale that extend existing operating envelopes.

Engineering↗
DOE OSTI2020

Pool Fire Analysis Methodology for Assessing Damage to Waste Containers

This document provides guidelines, whereby the physical parameters associated with postulated flammable/combustible liquid pool fires are determined (thermal analysis) in order to estimate the number of waste containers impacted. Pool fire damage to waste containers may occur due to a container being fully engulfed (within the periphery of the pool), partially engulfed, or within the critical heat flux distance. The methodology described here for estimating the size of the pool fire is developed using spill characteristics outlined in Chapter 65 of the SFPE Handbook (Ref. 7) and testing done by Hansen (Ref. 25). For exposures, the critical flux height is determined using the detailed method of Shokri and Beyler (Ref. 8) from the flame height, calculated as described by Heskestad (Ref. 7) and the average black body emissive power, determined using the method of Shokri and Beyler (Ref. 8). Three types of liquid pool fires should be considered; unconfined crash/instantaneous spill fires, unconfined continuously flowing (also called metered leak) spill fires, and confined spill fires. All fuel sources should be considered consistent with the facility Fire Hazard Analysis (FHA). However, the most likely liquid fuel sources present in TRU waste facilities are those present on TRU waste handling vehicles (forklifts, pallet movers, transport trucks, etc.). This method can also be used for assessing damage to items other than waste containers or other methods of analysis may be used with appropriate justification.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

A Modeling Toolkit for Comparing AC and DC Electrical Distribution Efficiency in Buildings

Recently, there has been considerable research interest in the potential for DC distribution systems in buildings instead of the traditional AC distribution systems. Due to the need for performing power conversions between DC and AC electricity, DC distribution may provide electrical efficiency advantages in some systems. To support comparative evaluations of AC-only, DC-only, and hybrid AC/DC distribution systems in buildings, a new modeling toolkit called the Building Electrical Efficiency Analysis Model (BEEAM) was developed and is described in this paper. To account for harmonics in currents or voltages arising from nonlinear devices, the toolkit implements harmonic power flow, along with nonlinear device behavioral descriptions derived from empirical measurements. This paper describes the framework, network equations, device representations, and an implementation of the toolkit in an open source software package, including a component library and graphical interface for creating circuits. Simulations of electrical behavior and device and system efficiencies using the toolkit are compared with experimental measurements of a small office environment in a variety of operating and load configurations. A detailed analysis of uncertainty estimation is also provided. Key findings were that a comparison of predicted versus measured efficiencies and power losses in the validation testbed using the initial toolkit implementation predicted device- and system-level efficiencies with reasonably good accuracy under both balanced and unbalanced AC scenarios. An uncertainty analysis also revealed that the maximum estimated error for system efficiency across all scenarios was 3%, and measured and modeled system efficiency agreed within the experimental uncertainty in approximately half of the scenarios. Based on the correspondence between simulation and measurement, the toolkit is proposed by the authors as a potentially useful tool for comparing efficiency in AC, DC, and hybrid AC/DC distribution systems in buildings.

DC distribution↗
DOE OSTI2024

The three-dimensional structure of black hole accretion flows within the plunging region

ABSTRACT We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e. at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry given by the standard low-magnetization ‘SANE’ configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the ‘thin disc’ limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two-dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling.

Astronomy & Astrophysics↗
DOE OSTI2023

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗
DOE OSTI2020

Combined Cr and S poisoning behaviors of La 1-x Sr x MnO 3±δ and La 1-x Sr x Co 1-y FeyO 3-δ cathodes in solid oxide fuel cells

Although the individual effects of airborne Cr and S contaminants on SOFC cathode performance degradation have been extensively studied, the combined effects of Cr and S contaminants remain largely unexplored. Under the real SOFC operating condition where the Cr and S species coexist, their effects may compete, affecting the poisoning behavior. Here, our investigation reveals that the combined Cr and S poisoning behavior of LSCF remains different from those of individual Cr and S poisonings, while the combined poisoning mechanism of LSM is equivalent to the sum of those of individual Cr and S effects. For LSCF electrode, gaseous Cr species are deposited mainly at the LSCF/GDC interface by electrochemical reduction, rather than forming SrCrO 4 on LSCF surfaces (as in the case of Cr-only poisoning), indicating no reaction between Cr vapors and SrO on the LSCF surface. Thermodynamic analysis demonstrates that the SrO on LSCF surface absorbs SO 2 (g) and thereby loses the Cr-gettering effect, allowing Cr vapors to flow through LSCF and to reach the LSCF/GDC interface where the Cr deposition occurs. Unlike LSCF, LSM electrode shows cumulative effects of Cr and S, as Cr accumulation occurs at the triple-phase boundary and S absorption takes place at localized Sr-rich regions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Using pollen in turbidites for vegetation reconstructions

Turbidites, deposited by sub-aqueous gravity flows, are common in sedimentary archives worldwide and present a unique challenge and opportunity when reconstructing past vegetation through pollen analysis. When sampling pollen from a sediment core for palaeovegetation records, it is common practice to target background sediments (i.e. pelagic sediment) and avoid sampling turbidites, as they are presumed to portray a misleading picture of past vegetation. This assumption stems from our limited understanding of pollen abundance and distribution through turbidites, meaning that palynologists overlook deposits that could potentially be used to reconstruct past vegetation and climate. We present pollen assemblage and sedimentological data from four recent (<150 years) deep marine turbidite deposits from the Hikurangi Subduction Margin, Aotearoa-New Zealand, with the aim of understanding the abundance and distribution of pollen in fine-grained turbidites. We find that pollen is diluted in the bases of turbidites, but despite this dilution, the proportions of different pollen taxa remain consistent through each turbidite. These results confirm that pollen can be sampled from turbidites for palaeovegetation reconstructions and that sampling the fine-grained upper parts of turbidites will provide the best pollen recovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Adaptable SEC‐SAXS data collection for higher quality structure analysis in solution

Abstract The two major challenges in synchrotron size‐exclusion chromatography coupled in‐line with small‐angle x‐ray scattering (SEC‐SAXS) experiments are the overlapping peaks in the elution profile and the fouling of radiation‐damaged materials on the walls of the sample cell. In recent years, many post‐experimental analyses techniques have been developed and applied to extract scattering profiles from these problematic SEC‐SAXS data. Here, we present three modes of data collection at the BioSAXS Beamline 4–2 of the Stanford Synchrotron Radiation Lightsource (SSRL BL4‐2). The first mode, the High‐Resolution mode, enables SEC‐SAXS data collection with excellent sample separation and virtually no additional peak broadening from the UHPLC UV detector to the x‐ray position by taking advantage of the low system dispersion of the UHPLC. The small bed volume of the analytical SEC column minimizes sample dilution in the column and facilitates data collection at higher sample concentrations with excellent sample economy equal to or even less than that of the conventional equilibrium SAXS method. Radiation damage problems during SEC‐SAXS data collection are evaded by additional cleaning of the sample cell after buffer data collection and avoidance of unnecessary exposures through the use of the x‐ray shutter control options, allowing sample data collection with a clean sample cell. Therefore, accurate background subtraction can be performed at a level equivalent to the conventional equilibrium SAXS method without requiring baseline correction, thereby leading to more reliable downstream structural analysis and quicker access to new science. The two other data collection modes, the High‐Throughput mode and the Co‐Flow mode, add agility to the planning and execution of experiments to efficiently achieve the user's scientific objectives at the SSRL BL4‐2.

Matsui, Tsutomu↗
DOE OSTI2025

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗