Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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,027 records · Page 57

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

LES↗

Flow Reconstruction in A Transonic Turbine Cascade Using Physics-Informed Neural Networks (PINNS)

This paper investigates the application of Physics-Informed Neural Networks (PINNs) for the analysis of turbine blades in a transonic cascade. The 2-D flow field in a transonic turbine cascade is reconstructed in three ways: the traditional forward approach (PINN not trained on experimental data), by training the PINN using discrete sets of experimentally measured pressure at midspan, and in the inverse sense where no inlet or outlet pressure boundary conditions are applied. Comparisons between the PINN solutions to measured data are made. This is repeated for three different turbine blades with distinct loading characteristics. Good agreement is shown between a CFD calculation of the CMC7 blade, and the PINN model trained with all data. The PINN is trained utilizing all available data, half the available data, data from only the leading edge region, and data from only the trailing edge region. The forward problem results deviate the most from experimental data but show promise. Solutions from the assisted training cases show that the PINN can reconstruct the flow field with acceptable accuracy when trained on measurements along the entire blade. In the inverse case, it is shown that to simultaneously achieve acceptable errors for inlet Mach number and outlet isentropic Mach number, the PINN must be trained on the static pressure data along the entire blade.

Machine Learning↗

The Atmospheric Grand Challenges in the Science of Wind Energy

Harvested by advanced technical systems honed over decades of research and development, wind energy has become a mainstream energy resource. However, continued innovation is needed to realize the potential of wind to serve the global demand for clean energy. Here, we outline three interdependent, cross-disciplinary grand challenges underpinning this research endeavor. The first is the need for a deeper understanding of the physics of atmospheric flow in the critical zone of plant operation. The second involves science and engineering of the largest dynamic, rotating machines in the world. The third encompasses optimization and control of fleets of wind plants working synergistically within the electricity grid. Addressing these challenges could enable wind power to provide as much as half of our global electricity needs and perhaps beyond.

17 WIND ENERGY↗

Transition to a Thorium-based Fuel Cycle - Choices, Challenges, and Options

The Charter for the Evaluation and Screening (E and S) of Nuclear Fuel Cycle Options requested by the US Department of Energy calls for the identification of promising fuel cycle options to help guide the future directions of research and development in this area. The final report for this study, which was produced by various national laboratory and industry experts, collected the expansive set of fuel cycle options into 40 evaluation groups (EGs). In support of this study, detailed steady-state physics and mass flow data for each of these groups were compiled into a set of Fuel Cycle Data Packages (FCDPs). The data in the FCDPs were then used to quantify the nine criteria detailed in the study to identify the most promising EGs. In addition to supporting the E and S study, data from these FCDPs can be used to benchmark fuel cycle simulation tools for equilibrium calculations. Once benchmarked, their dynamic modeling capabilities can support the continuing effort to determine the most effective approaches for transitioning to proposed advanced future fuel cycles. While steady-state conditions can be calculated directly in some cases, fuel cycle simulation tools allow for the transition scenarios to be simultaneously optimized on multiple parameters such as the fissile material availability, the amount of desirable material in storage, and the time required for transition. In addition to fuel cycle simulations where material availability is considered and is often a constraint to achieving transition, determining the technology needs to enable a transition to take place is also a key consideration (e.g., what technology is needed, when it is needed, and whether it is needed only during the transition and not for the equilibrium fuel cycle). This paper evaluates the material and technology issues and challenges associated with a transition from a uranium-based fuel cycle (such as that operating in the US today, represented by EG01) to a thorium-based fuel cycle, recognizing that thorium requires fissile material to start up. Based on the findings of the E and S report, the two most promising thorium fuel cycle options are considered in this report: (1) EG26, a thermal spectrum molten salt reactor (MSR), and (2) EG28, represented by a sodium fast reactor (SFR, fast spectrum) and a resource-renewable boiling water reactor (RBWR, intermediate spectrum). (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Novel Preprocessing Techniques for Enhancing Flowability of Miscanthus

Purpose-grown energy crops possess significant potential as promising feedstocks for biofuel and biochemical production. However, similar to other biomass feedstocks, their utilization has been hindered by handling and feeding challenges such as clogging and segregation, which prevent biorefineries from operating at full capacity. These challenges arise primarily due to particle attributes including irregular shape, high flexibility, and high compressibility, leading to interlocking, high compaction, and significant friction under stress consolidation. To enhance operational reliability, it is crucial to focus on preprocessing techniques to address these issues, alongside accurately characterizing the mechanical and physical properties and flow performance of these feedstocks. This study investigates novel preprocessing techniques, including pelletization and torrefaction, to fundamentally alter these particle attributes and improve flowability. Miscanthus, an herbaceous energy crop, was the primary focus of this study. Samples were processed using a Forest Concepts Crumbler rotary shear system and an oscillating multi-stack screen bed to achieve two nominal (baseline) particle sizes: 6mm and 2mm. The study evaluated the efficiency and energy consumption of size reduction, revealing that while 2mm particles required more energy, they exhibited better flowability and higher yield strength compared to 6mm particles. Shear and compression tests assessed the material's shear strength, internal friction, and elastic modulus, while wedge-shaped hopper tests measured the critical arching distance and mass flow rate. The results provide valuable insights into the mechanical properties and processing efficiency of Miscanthus, contributing to the development of a preprocessing framework that enhances energy efficiency and material flowability in bioenergy production. However, further preprocessing and development are required to optimize this framework fully.

09 - BIOMASS FUELS↗

Computational and Experimental Fluid-Structure Interaction Analysis of a High-Lift Wing with a Slat-Cove Filler for Noise Reduction

Airframe noise is a significant component of overall noise produced by transport aircraft during landing and approach (low speed maneuvers). A significant source for this noise is the cove of the leading-edge slat. The slat-cove filler (SCF) has been shown to be effective at mitigating slat noise. The objective of this work is to understand the fluid-structure interaction (FSI) behavior of a superelastic shape memory alloy (SMA) SCF in flow using both computational and physical models of a high-lift wing. Initial understanding of flow around the SCF and wing is obtained using computational fluid dynamics (CFD) analysis at various angles of attack. A framework compatible with an SMA constitutive model (implemented as a user material subroutine) is used to perform FSI analysis for multiple flow and configuration cases. A scaled physical model of the high-lift wing is constructed and tested in the Texas A&M 3 ft-by-4-foot wind tunnel. Initial validation of both CFD and FSI analysis is conducted by comparing lift, drag and pressure distributions with experimental results.

Scholten, William D.↗

Interpenetrating Plasma Simulations

This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the physical mechanisms that modify wind plant blockage in stable boundary layers

Abstract. Wind plants slow down the approaching wind, a phenomenon known as blockage. Wind plant blockage undermines turbine performance for front-row turbines and potentially for turbines deeper into the array. We use large-eddy simulations to characterize blockage upstream of a finite-size wind plant in flat terrain for different atmospheric stability conditions and investigate the physical mechanisms modifying the flow upstream of the turbines. To examine the influence of atmospheric stability, we compare simulations of two stably stratified boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. For a wind plant, a faster cooling rate at the surface, which produces stronger stably stratified flow in the boundary layer, amplifies blockage. As a novelty, we investigate the physical mechanisms amplifying blockage by evaluating the different terms in the momentum conservation equation within the turbine rotor layer. The velocity deceleration upstream of a wind plant is caused by an adverse pressure gradient and momentum advection out of the turbine rotor layer. The cumulative deceleration of the flow upstream of the front-row turbines instigates vertical motions. The horizontal flow is diverted vertically, reducing momentum availability in the turbine rotor layer. Although the adverse pressure gradient upstream of the wind plant remains unchanged with atmospheric stability, vertical advection of horizontal momentum is amplified in the more strongly stable boundary layer, mainly by larger shear of the horizontal velocity, thus increasing the blockage effect.

17 WIND ENERGY↗

Contributions to the understanding of large-scale coherent structures in developing free turbulent shear flows

Advances in the mechanics of boundary layer flow are reported. The physical problems of large scale coherent structures in real, developing free turbulent shear flows, from the nonlinear aspects of hydrodynamic stability are addressed. The presence of fine grained turbulence in the problem, and its absence, lacks a small parameter. The problem is presented on the basis of conservation principles, which are the dynamics of the problem directed towards extracting the most physical information, however, it is emphasized that it must also involve approximations.

Liu, J. T. C.↗

Relaxations of the steady optimal gas flow problem for a non-Ideal gas

Natural gas ranks second in U.S. primary energy consumption. Because most production sites are remote, gas must be transported through pipeline networks equipped with compressors, valves, and other components. For both economic efficiency and system reliability, it is desirable to operate these networks optimally. The governing physics across pipeline components entails nonlinear, non-convex equality and inequality constraints, and the most general steady-flow operations problem is a Mixed-Integer Nonlinear Program (MINLP).This work focuses on one such steady-flow problem-the Optimal Gas Flow (OGF) for a natural gas pipeline network-which minimizes production cost subject to the steady-flow physics. For day-to-day operations, the ability to quickly compute a globally optimal solution and a strong lower bound for varying demand profiles is crucial. A promising strategy is to build tight relaxations of the OGF’s nonlinear constraints. However, many nonlinearities arising from non-ideal equations of state either lack relaxations or have relaxations that do not scale to realistic network sizes. We address this gap by combining recent advances in polyhedral relaxations for univariate functions to construct tight, computationally efficient relaxations of the OGF with a non-ideal equation of state. These relaxations solve within seconds on a standard laptop. In conclusion, we demonstrate their quality through extensive numerical experiments on very large-scale test networks from the literature and find that the proposed approach proves optimality in 92% of tested instances.

03 NATURAL GAS↗

Numerical solutions of the complete Navier-Stokes equations

The physical phenomena within supersonic flows that sustain chemical reactions are investigated. An earlier study to develop accurate physical models for supersonic reacting flowfields focused on 2-D laminar shear layers. The objective is to examine the mixing and subsequent combustion within turbulent reacting shear layers. To conduct this study, a computer program has been written to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method uses a cell-centered finite volume approach and a Runge Kutta time stepping scheme. The Reynolds averaged equations are closed using the eddy viscosity concept. Several zero-equation models have been tested by making calculations for an H2-air nonreacting coaxial jet flow. Comparisons made with experimental data show that Cohen's eddy viscosity model provides best agreement. The finite rate chemistry model used in the study of 2-D laminar shear layers is incorporated into the computer program and data is compared from a recent experiment performed at NASA Langley.

Hassan, H. A.↗

MATCHED-INDEX-OF-REFRACTION FLOW FACILITY FOR FUNDAMENTAL AND APPLIED RESEARCH

Significant challenges face reactor designers with regard to thermal hydraulic design and associated modeling for advanced reactor concepts. Computational thermal hydraulic codes solve only a piece of the core. There is a need for a whole core dynamics system code with local resolution to investigate and understand flow behavior with all the relevant physics and thermo-mechanics. The matched index of refraction (MIR) flow facility at Idaho National Laboratory (INL) has a unique capability to contribute to the development of validated computational fluid dynamics (CFD) codes through the use of state-of-the-art optical measurement techniques, such as Laser Doppler Velocimetry (LDV) and Particle Image Velocimetry (PIV). PIV is a non-intrusive velocity measurement technique that tracks flow by imaging the movement of small tracer particles within a fluid. At the heart of a PIV calculation is the cross correlation algorithm, which is used to estimate the displacement of particles in some small part of the image over the time span between two images. Generally, the displacement is indicated by the location of the largest peak. To quantify these measurements accurately, sophisticated processing algorithms correlate the locations of particles within the image to estimate the velocity (Ref. 1). Prior to use with reactor deign, the CFD codes have to be experimentally validated, which requires rigorous experimental measurements to produce high quality, multi-dimensional flow field data with error quantification methodologies. Computational thermal hydraulic codes solve only a piece of the core. There is a need for a whole core dynamics system code with local resolution to investigate and understand flow behavior with all the relevant physics and thermo-mechanics. Computational techniques with supporting test data may be needed to address the heat transfer from the fuel to the coolant during the transition from turbulent to laminar flow, including the possibility of an early laminarization of the flow (Refs. 2 and 3) (laminarization is caused when the coolant velocity is theoretically in the turbulent regime, but the heat transfer properties are indicative of the coolant velocity being in the laminar regime). Such studies are complicated enough that computational fluid dynamics (CFD) models may not converge to the same conclusion. Thus, experimentally scaled thermal hydraulic data with uncertainties should be developed to support modeling and simulation for verification and validation activities. The fluid/solid index of refraction matching technique allows optical access in and around geometries that would otherwise be impossible while the large test section of the INL system provides better spatial and temporal resolution than comparable facilities. Benchmark data for assessing computational fluid dynamics can be acquired for external flows, internal flows, and coupled internal/external flows for better understanding of physical phenomena of interest. The core objective of this study is to describe MIR and its capabilities, and mention current development areas for uncertainty quantification, mainly the uncertainty surface method and cross-correlation method. Using these methods, it is anticipated to establish a suitable approach to quantify PIV uncertainty for experiments performed in the MIR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A physically consistent model for artificial dissipation in transonic potential flow computations

The effect that artificial dissipation has on numerical solutions of the transonic Full Potential Equation (FPE) are investigated by comparing the artificially dissipative FPE to a Physically Dissipative Potential (PDP) equation. Analytic expressions were derived from the variables C and M sub c that are used in the artificial density formulation. It was shown that these new values generate artificial dissipation which is equivalent to the physical dissipation existing in the PDP equation. The new expression for the variables C and M sub c can easily be incorporated into the existing full potential codes which are based either on the artificial density or on the artificial viscosity formulation. A comparison of Physically Dissipative Potential (PDP), Artificial Density or Viscosity (ADV), Artificial Mass Flux (AMF), and ADV with variable C and M sub c formulation (MCC) is also presented.

Dulikravich, George S.↗

Supersonic combustor modeling

The physical phenomena involved when a supersonic flow undergoes chemical reaction are discussed. Detailed physical models of convective and diffusive mixing, and finite rate chemical reaction in supersonic flow are presented. Numerical algorithms used to solve the equations governing these processes are introduced. Computer programs using these algorithms are used to analyze the structure of the reacting mixing layer. It is concluded that, as in subsonic flow, exothermic heat release in unconfined supersonic flows retards fuel/air mixing. Non mixing is shown to be a potential problem in reducing the efficiency of supersonic as well as subsonic combustion. Techniques for enhancing fuel/air mixing and combustion are described.

Riggins, David W.↗

Computation of oscillating airfoil flows with one- and two-equation turbulence models

The ability of one- and two-equation turbulence models to predict unsteady separated flows over airfoils is evaluated. An implicit, factorized, upwind-biased numerical scheme is used for the integration of the compressible, Reynolds-averaged Navier-Stokes equations. The turbulent eddy viscosity is obtained from the computed mean flowfield by integration of the turbulent field equations. One- and two-equation turbulence models are first tested for a separated airfoil flow at fixed angle of incidence. The same models are then applied to compute the unsteady flowfields about airfoils undergoing oscillatory motion at low subsonic Mach numbers. Experimental cases where the flow has been tripped at the leading-edge and where natural transition was allowed to occur naturally are considered. The more recently developed turbulence models capture the physics of unsteady separated flow significantly better than the standard kappa-epsilon and kappa-omega models. However, certain differences in the hysteresis effects are observed. For an untripped high-Reynolds-number flow, it was found necessary to take into account the leading-edge transitional flow region to capture the correct physical mechanism that leads to dynamic stall.

Ekaterinaris, J. A.↗

Investigation of mesoscale pre-convective forcing mechanisms

The installation and testing of a numerical mesoscale boundary layer model is discussed. A data analysis is first made using the rawinsonde network which shows a strong baroclinic zone developing in the region where convection was initiated. The observed wind field shows an easterly component developing in opposition to the synoptic flow producing convergence. The numerical boundary layer model is used to diagnose the physical development of the thermal and flow fields. Two physical mechanisms are examined; topography and cloud shading. The results shows that while heated sloping topography alone does create a baroclinic zone and convergence values, cloud shading is a stronger mechanism for producing a thermal field representative of observed conditions

Mcnider, R. T.↗

Direct Numerical Simulation of Turbulent Flow Inside a Differentially Heated Composite Cavity

In this work, turbulent flow inside a modified differentially heated cavity at high Rayleigh number (Ra ~ 10 9 ) has been studied through fully resolved direct numerical simulation (DNS) using the high-order spectral element method code Nek5000. The flow configuration includes two separate physical phenomena: the natural recirculation itself, and the flow inside a curved channel. Simulations have been carried out using both the Boussinesq approximation and the low-Mach compressible formulation. Significant discrepancies between the two methods inform of the extreme caution that should be exercised when using the Boussinesq approximation in the limits of its applicability. The DNS solutions are analyzed in terms of polynomial-order convergence and Reynolds stress budgets, and the turbulence quantities and velocity profiles are presented as a reference for the validation of turbulence models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗