Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On the measurement of S diff splitting caused by lowermost mantle anisotropy

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (S diff ) are potentially strongly influenced by lowermost mantle anisotropy. S diff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. S diff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SH diff –SV diff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of S diff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which S diff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on S diff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SH diff waves. In particular, we investigate how SH diff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to S diff splitting that can be clearly distinguished from other waveform effects. Here, these results enable us to lay out a strategy for the analysis of S diff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of S diff and the influence of source-side anisotropy. We demonstrate our S diff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained S diff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Experimental approach for development of a powder spreading metric in additive manufacturing

The powder spreading is a vital step of powder-based additive manufacturing (AM) processes. The quality of spread powder can considerably influence the properties of fabricated parts. Poorly packed powder beds with high surface roughness result in printed part layers with large porosity and low dimensional accuracy, leading to poor mechanical properties. Therefore, the powder spreadability and its dependence on process parameters and powder characteristics should be quantified to improve the efficiency of powder-based AM methods. This study proposes a novel dimensionless powder spreadability metric that can be commonly used in different powder-based AM processes. The quality of spread powder in terms of powder bed density and surface roughness was evaluated by adjusting the process parameters including recoating velocity and layer thickness, and powder characteristics including particle size distribution. In addition, the dynamic repose angle was proposed and examined as another powder spreadability metric. Furthermore, the results showed that these two proposed metrics were strongly correlated and lower recoating velocity and larger layer thickness led to higher spreadability and lower dynamic repose angle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A RELAP5-3D Model of the Lobo Lead Loop

The goal of our research is to build upon the capability of RELAP5-3D to model molten lead systems. Molten lead has several potential uses in future advanced reactors like the lead fast reactor or fusion reactors that utilize dual-coolant lead lithium blankets. This potential for use in future generations of reactors highlights the necessity to develop molten lead models to ensure that they can accurately predict the thermohydraulic behavior. We have developed a RELAP5-3D model of the Lobo Lead Loop facility, located at the University of New Mexico, to verify the accuracy of RELAP5-3D via comparison to existing computational fluid dynamics results and analytical calculations. It was found that RELAP5-3D accurately calculated radiative heat transfer (within < 1%) when compared to theoretical calculations. In addition, pressure drop calculations done in RELAP5-3D demonstrated reasonable agreement within 20 kPa, mostly within ~7-15%, when compared to the computational fluid dynamics model of the facility developed by the University of New Mexico, and captured the dependence of pressure drop on flow velocity accurately. Finally, a hypothetical loss of flow transient was imposed on the RELAP5-3D model to determine the feasibility of performing a similar experiment with the Lobo Lead Loop. It was found that such an experiment could be possible as the RELAP5-3D model indicated that the temperatures of the fluid would not exceed limiting temperatures of the structure (1658 K) nor the maximum temperature of the electromagnetic pump inlet (823 K). Although there is not experimental data to begin validation, the model will be readily available for future validation studies when the experimental data is generated, especially as the model continues to evolve over time. Furthermore, the results so far demonstrate a promising first step in the verification/validation of the RELAP5-3D model of the Lobo Lead Loop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Co-Optimization of Velocity and Charge-Depletion for Plug-In Hybrid Electric Vehicles: Accounting for Acceleration and Jerk Constraints

Abstract Recent advances in vehicle connectivity and automation technologies promote advanced control algorithms that co-optimize the longitudinal dynamics and powertrain operation of hybrid electric vehicles. Typically, a sequential optimization with the vehicle dynamics optimized followed by powertrain optimization is adopted to manage a number of complexities such as the inherent mixed-integer nature of the hybrid powertrain, the numerous state and control variables, the differing time scales of vehicle and powertrain subsystems, time-varying state constraints, and large horizon lengths. Instead, we solve the offline optimization problem in a centralize manner assuming exact knowledge of the lead vehicle's position over the entire trip by applying a discrete-time single shooting-based numerical approach, Discrete Mixed-Integer Shooting (DMIS), including a linearly increasing computational complexity to the problem horizon. In particular, the hierarchical problem structure is exploited to decompose the computationally intensive Hamiltonian minimization step into a set of low-dimensional optimizations. DMIS allows us to compute the direct fuel minimization problem including the vehicle and powertrain dynamics in a centralized manner to its full horizon while systematically tuning weighting factors that penalize passenger discomfort. For the first time, this study reveals that practically implemented sequential optimization exhibits similar fuel optimality as co-optimization when a certain level of passenger comfort is required.

Automation & Control Systems↗
DOE OSTI2022

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗
DOE OSTI2022

Predicting Asteroid Material Properties from a DART-like Kinetic Impact

Abstract NASA’s Double Asteroid Redirection Test (DART) mission is the first full-scale test of the kinetic impactor method for asteroid deflection, in which a spacecraft intentionally impacts an asteroid to change its trajectory. DART represents an important first step for planetary defense technology demonstration, providing a realistic assessment of the effectiveness of the kinetic impact approach on a near-Earth asteroid. The momentum imparted to the asteroid is transferred from the impacting spacecraft and enhanced by the momentum of material ejected from the impact site. However, the magnitude of the ejecta contribution is dependent on the material properties of the target. These properties, such as strength and shear modulus, are unknown for the DART target asteroid, Dimorphos, as well as most asteroids since such properties are difficult to characterize remotely. This study examines how hydrocode simulations can be used to estimate material properties from information available post-impact, specifically the asteroid size and shape, the velocity and properties of the impacting spacecraft, and the final velocity change imparted to the asteroid. Across >300 three-dimensional simulations varying seven material parameters describing the asteroid, we found many combinations of properties could reproduce a particular asteroid velocity. Additional observations, such as asteroid mass or crater size, are required to further constrain properties like asteroid strength or outcomes like the momentum enhancement provided by impact ejecta. Our results demonstrate the vital importance of having as much knowledge as possible prior to an impact mission, with key material parameters being the asteroid’s mass, porosity, strength, and elastic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Sensitivity Analysis of an Electrospray Dehumidification System

Electrically enhanced condensation achieved by an electrospray system can enhance dehumidification in HVAC applications if successfully scaled up. The first step towards scale-up is identifying critical design and operating parameters that significantly impact dehumidification for the electrospray system. Sensitivity analysis is one approach for identifying critical parameters. This paper presents a methodology that can efficiently perform quantitative sensitivity analysis on a computationally expensive high fidelity simulation model, which is a computational fluid dynamics (CFD) model of the electrospray system. The methodology employs surrogate modeling and Sobol sensitivity analysis. We estimate the main and total effect Sobol sensitivity indices of four operating parameters of the electrospray system. The results reveal that the most important parameter is the size of the sprayed droplets, followed by the volumetric flow rate of the bulk air and injected water in the system. The water injection velocity of the spray was found to have little to no impact on dehumidification.

Young, David↗
DOE OSTI2024

Solid–liquid separation of lignocellulosic sugars from biomass by rotating ceramic disc filtration

In many biomass conversion processes, the separation of cellulosic sugars from residual, lignin-rich solids is a critical step, and achieving high recovery yields of sugars by conventional tangential crossflow and vacuum filtration is challenged by the presence of biomass solids, which rapidly foul filters, resulting in decreased throughput. Considering the performance limitations of traditional filtration methods, dynamic filtration, which generates high shear at the membrane surface to decrease fouling, is emerging as a viable alternative for demanding solid–liquid separations. For high solids separations, there is little available information regarding the performance, limitations, and energy consumption of dynamic filtration. To that end, here we characterized the performance of a dynamic filtration module, specifically a rotating ceramic disc (RCD) filter, for the aseptic recovery of cellulosic sugars from biomass solids following pretreatment and enzymatic hydrolysis. We show how RCD rotational velocity and percent biomass solids impact the filter throughput. Additionally, we used computational fluid dynamics (CFD) simulations to estimate the shear rate at the membrane surface and to visualize hydrodynamic profiles within the module. With the combined CFD simulations and experimental results, we estimated the energy demand and operating expenses for a viable dynamic filtration system operating with a lignocellulosic feed slurry. Our results indicate that an RCD filter can achieve ≥95% recovery of sugars and produce a retentate slurry containing 12 wt% insoluble solids with low energy consumption (a 2.2-fold improvement over cross-flow filtration) and low operating costs ($\$$0.06 per kg sugars). These results show a viable path towards operationally reliable, energy efficient, and cost-effective separations of sterilized cellulosic sugars from biomass solids and highlight the potential of dynamic filtration systems for challenging solid–liquid separations.

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

A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different -space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. Here, the new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA / mc 2 , exceeding for typical cell sizes and time steps.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Transient behavior of drift and ionization in atmospheric pressure nitrogen discharge

The fluid models are frequently used to describe a non-thermal plasma such as a streamer discharge. The required electron transport data and rate coefficients for the fluid model are parametrized using the local field approximation (LFA) in first order models and the local-mean-energy approximation (LMEA) in second order models. We performed Monte Carlo simulations in Nitrogen gas with step changes in the E/N (reduced electric field) to study the behavior of the transport properties in the transient phase. During the transient phase of the simulation, we extract the instantaneous electron mean energy, which is different from the steady state mean electron energy, and the corresponding transport parameters and rate coefficients. Our results indicate that the mean electron energy is not a suitable parameter for mobility/drift of electrons due to big difference in momentum relaxation and energy relaxation. Furthermore, the high energy threshold rates such as ionization show a strong correlation to mean electron energy. In second order models where the energy-balance equation is solved, we suggest that it would rather be appropriate to use the local electric field to find electron drift velocity in gases such as Nitrogen and the local mean electron energy to determine the ionization and excitation rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗
DOE OSTI2022

Thermite and intermetallic projectiles examined experimentally in air and inert gas environments

Intermetallic (aluminum and zirconium) and thermite (aluminum and molybdenum trioxide) projectiles were launched using a high velocity impact ignition testing system. The experiments were designed to simulate reactivity in high (argon) and low (air) altitude environments. The projectiles were launched into a chamber that included a steel target plate for projectile penetration before impacting a rear witness plate. The chamber was semi-sealed and instrumented for quasi-static pressure data. The results provide an understanding of energy release from the projectile materials and of the environmental influence on performance. The transient pressure traces provide insight into reaction kinetics. A bifurcation in transient pressure rise was an indication of a shift in reaction kinetics from the inherent reactive material to metal oxidation with the environment. The bifurcation was delayed by about 0.15 ms for the intermetallic relative to the thermite, evidence that the thermite reaction proceeded faster upon impact than the intermetallic. The two-step process (impact ignition of the reactive material followed by metal oxidation) was shown to produce higher energy conversion efficiencies than projectiles composed of pure fuel (i.e., aluminum) reported previously. Both reactive materials showed energy conversion efficiencies greater than 30% (for air) and 50% (for argon), and an explanation of underestimated efficiency and energy losses is provided. These results have implications for advancing formulations for ballistic applications. Structural reactive materials can be used to modify the effective reactivity of metal-containing formulations in varied atmospheric environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reduced recombination via tunable surface fields in perovskite thin films

Here, the ability to reduce energy loss at semiconductor surfaces through passivation or surface field engineering is an essential step in the manufacturing of efficient photovoltaic (PV) and optoelectronic devices. Similarly, surface modification of emerging halide perovskites with quasi-two-dimensional (2D) heterostructures is now ubiquitous to achieve PV power conversion efficiencies (PCEs) >25%, yet a fundamental understanding to how these treatments function is still generally lacking. Here we use a unique combination of depth-sensitive nanoscale characterization techniques to uncover a tunable passivation strategy and mechanism found in perovskite PV devices that were the first to reach the >25% PCE milestone. Namely, treatment with hexylammonium bromide leads to the simultaneous formation of an iodide-rich 2D layer along with a Br halide gradient that extends from defective surfaces and grain boundaries into the bulk three-dimensional (3D) layer. This interface can be optimized to extend the charge carrier lifetime to record values >30 μs and to reduce interfacial recombination velocities to values as low as <7 cm s −1 .

Devices for energy harvesting↗
DOE OSTI2020

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗
DOE OSTI2021

Validation of Computational Fluid Dynamics Simulations for Biological Performance Assessment in Hydropower units (Final Report)

The biological performance assessment (BioPA) toolset developed by Pacific Northwest National Laboratory (PNNL) estimates the relative biological performance of fish passage at a hydroelectric power turbine unit. The tool is based on the use of computational fluid dynamics (CFD) and fish biological response relationships. The recent release, BioPA-v3, is based on directly computed trajectory and collision of material Lagrangian particles using CFD simulation codes rather than the prior version that relies on Tecplot to compute streamtrace trajectories. Before modifying the toolset, a series of validation tests were performed at the various steps of modification in the toolset. Validation is a critical step of any numerical investigation that reflects the accuracy and reliability of the predicted results. It raises the confidence level of the user to use the modified version of the BioPA toolset. Several test cases were simulated and compared, where available, to observed data. The trajectory and collision of the small spherical and cylindrical particles in a water flume were compared to in-house experiments. The CFD predicted collision rate and flow field compared well with experimental observation for vane array and large cylinder as target bodies. Next, the CFD-predicted flow field and hydraulic performance of a laboratory-scale model of a Francis turbine was also successfully validated. Note that the trajectory of the particles is significantly affected by the flow field in such extreme conditions. In addition to the particle trajectories and flow field, the collision detection method employed in the CFD simulations was also successfully validated. The CFD predicted impact velocity, collision time, velocity, and trajectory of a sphere excellently matched with analytical value for a bouncing ball in the elastic collision. A similar approach was also tested and successfully validated for a collision of sphere with a 45° inclined plane. After successfully validating different cases, the BioPA toolset was modified to use direct output of the CFD prediction and the new version can be used in evaluating biological performance at hydroelectric turbines.

13 HYDRO ENERGY↗