Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid flow regimes”

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 505 records · Page 28

Turbulent entrainment in finite-length wind farms

In this article, we present an entrainment-based model for predicting the flow and power output of finite-length wind farms. The model is an extension of the three-layer approach of Luzzatto-Fegiz & Caulfield ( Phys. Rev. Fluids , vol. 3, 2018, 093802) for wind farms of infinite length, and assumes dependence of key flow quantities, such as the wind farm bulk velocity, on the streamwise distance from the farm entrance. To assist our analysis and validate the proposed model, we undertake a series of large-eddy simulations with different turbine spacing arrangements and layouts. Comparisons are also made with the top-down model with entrance effects of Meneveau ( J. Turbul. , vol. 13, 2012, N7) and data from the literature. The finite-length entrainment model is shown to be capable of capturing the power drop between contiguous rows of turbines as well as describing the advection and turbulent transport of kinetic energy in both the entrance and fully developed regions. The fully developed regime is approximated only deep in the wind farm, after approximately 15 rows of turbines. Our data suggest that for the cases considered in this study, the empirical coefficients that can be used to describe turbulent entrainment and transfers above the wind farm exhibit little dependence on the farm layout and may be considered constant for modelling purposes. However, the flow field within the wind farm layer can be strongly modulated by the turbine density (spacing) as well as the array layout, and to that extent it can be argued that they are both primary factors determining the wind farm power output.

17 WIND ENERGY↗

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗

Some aspects of unsteady separation

Unsteady separation can be forced in a variety of ways and in this presentation two fundamental means will be considered: (1) the introduction of convected vorticular disturbances into the flow; and (2) the influence of a specific type of three-dimensional geometry. In both situations a response of the viscous flow near the wall is provoked wherein the fluid near the surface abruptly focuses into a narrow region that erupts from the surface into the mainstream. In two-dimensional flows, the eruption takes the form of a narrow, explosively-growing spike, while in three-dimensional situations, examples are presented which indicate that the eruption is along a narrow zone in the shape of a crescent-shaped plume. The nature of the three-dimensional flow near a circular cylinder, which is mounted normal to a flat plate, is also examined in this study. Here the three-dimensional geometry induces complex three-dimensional separations periodically. The dynamics of the generation process is studied experimentally in a water channel using hydrogen bubble wires and a laser sheet, and the main features of the laminar regime through to transition are documented.

Smith, C. R.↗

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

Drag force to all orders in gradients

We study the energy loss of a heavy quark slowly moving through an evolving strongly coupled plasma. We use the linearized fluid/gravity correspondence to describe small perturbations of the medium flow with general spacetime dependence. This all order linearized hydrodynamics results in a drag force exerted on a heavy quark even when it is at rest with the fluid element. We show how the general contribution to the drag force can be derived order by order in the medium velocity gradients and provide explicit results valid up to the third order. We then obtain an approximate semi-analytic result for the drag force to all orders in the gradient expansion but linearized in the medium velocity. Thus, the effects of a class of hydrodynamic gradients on the drag force are re-summed, giving further insight into the dissipative properties of strongly coupled plasmas. The all order result allows us to study the drag force in the non-hydrodynamic regime of linear medium perturbations that vary rapidly in space and time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Updated magnetized transport coefficients: impact on laser-plasmas with self-generated or applied magnetic fields

Errors in the Epperlein and Haines (1986 Phys. Fluids ) transport coefficients were recently found at low electron magnetizations, with new magnetic transport coefficients proposed simultaneously by two teams (Sadler et al 2021 Phys. Rev. Lett. and Davies et al 2021 Phys. Plasmas ); these two separate sets of updated coefficients are shown in this paper to be in agreement. Here, the importance of these new coefficients in laser-plasmas with either self-generated or applied magnetic fields is demonstrated. When an external magnetic field is applied, the cross-gradient-Nernst term twists the field structure; this twisting is reduced by the new coefficients in the low magnetization regime. For plasmas where only self-generated magnetic fields are present, the new coefficients are found to result in the magnetic field moving with the Righi–Leduc heat-flow, enhancing the impact of MHD. Simulations of Biermann battery magnetic fields around ICF hot-spot perturbations are presented, with cross-gradient-Nernst transport increasing spike penetration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of barriers on wave propagation phenomena: With application for aircraft noise shielding

The frequency spectrum was divided into high and low frequency regimes and two separate methods were developed and applied to account for physical factors associated with flight conditions. For long wave propagation, the acoustic filed due to a point source near a solid obstacle was treated in terms of an inner region which where the fluid motion is essentially incompressible, and an outer region which is a linear acoustic field generated by hydrodynamic disturbances in the inner region. This method was applied to a case of a finite slotted plate modelled to represent a wing extended flap for both stationary and moving media. Ray acoustics, the Kirchhoff integral formulation, and the stationary phase approximation were combined to study short wave length propagation in many limiting cases as well as in the case of a semi-infinite plate in a uniform flow velocity with a point source above the plate and embedded in a different flow velocity to simulate an engine exhaust jet stream surrounding the source.

Mgana, C. V. M.↗

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

Evaluation of Transport-Equation-Based Transition Models for High-Speed Boundary Layers Using OVERFLOW

Accurate modeling of laminar-turbulent transition is crucial for the design of hypersonic flight systems. However, the current transition models used in production CFD codes are insufficient for high-speed flows. Many extensions to low-speed models have been suggested; however, a thorough verification and validation effort is needed before these models can be used in design settings. Challenges include potentially missing details of the model implementation requirements and/or a complete specification of the input parameters needed to replicate the test findings. A meaningful assessment of the generalization capability of these models is also hindered by a lack of information regarding the specific flow configurations and associated grids employed for model calibration. As a key first step toward model verification, we present an independent assessment of two recently proposed models for high-speed transition, namely, a model within the SST-$\gamma$ framework and a model based on the SST-$\gamma-\nu_L$ equations. These models are implemented in the NASA OVERFLOW 2.3e solver and their performance in predicting first mode, second mode, and crossflow transition has been evaluated for several test cases in the supersonic and hypersonic regimes. Besides the test cases employed by the model developers, which could have also been used for model calibration, the present assessment includes supplementary configurations that contribute to an unbiased assessment of the models. The outcomes presented in this study indicate the potential for the models to be applied to high-speed flight configurations. Key steps toward future improvements to these models are also outlined.

High-speed flow↗

Operational and geological controls of coupled poroelastic stressing and pore-pressure accumulation along faults: Induced earthquakes in Pohang, South Korea

Coupled poroelastic stressing and pore-pressure accumulation along pre-existing faults in deep basement contribute to recent occurrence of seismic events at subsurface energy exploration sites. Our coupled fluid-flow and geomechanical model describes the physical processes inducing seismicity corresponding to the sequential stimulation operations in Pohang, South Korea. Simulation results show that prolonged accumulation of poroelastic energy and pore pressure along a fault can nucleate seismic events larger than M w 3 even after terminating well operations. In particular the possibility of large seismic events can be increased by multiple-well operations with alternate injection and extraction that can enhance the degree of pore-pressure diffusion and subsequent stress transfer through a rigid and low-permeability rock to the fault. This study demonstrates that the proper mechanistic model and optimal well operations need to be accounted for to mitigate unexpected seismic hazards in the presence of the site-specific uncertainty such as hidden/undetected faults and stress regime.

58 GEOSCIENCES↗

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Nonlinear Modal Aeroelastic Solver for FUN3D

A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.

Goldman, Benjamin D.↗

A package for 3-D unstructured grid generation, finite-element flow solution and flow field visualization

A set of computer programs for 3-D unstructured grid generation, fluid flow calculations, and flow field visualization was developed. The grid generation program, called VGRID3D, generates grids over complex configurations using the advancing front method. In this method, the point and element generation is accomplished simultaneously, VPLOT3D is an interactive, menudriven pre- and post-processor graphics program for interpolation and display of unstructured grid data. The flow solver, VFLOW3D, is an Euler equation solver based on an explicit, two-step, Taylor-Galerkin algorithm which uses the Flux Corrected Transport (FCT) concept for a wriggle-free solution. Using these programs, increasingly complex 3-D configurations of interest to aerospace community were gridded including a complete Space Transportation System comprised of the space-shuttle orbitor, the solid-rocket boosters, and the external tank. Flow solutions were obtained on various configurations in subsonic, transonic, and supersonic flow regimes.

Parikh, Paresh↗

Transition to hard turbulence in thermal convection at infinite Prandtl number

Direct numerical simulations of two-dimensional high Rayleigh (Ra) number, base-heated thermal convection in large aspect-ratio boxes are presented for infinite Prandtl number fluids, as applied to the earth's mantle. A transition is characaterized in the flow structures in the neighborhood of Ra between 10 to the 7th and 10 to the 8th. These high Ra flows consist of large-scale cells with strong intermittent, boundary-layer instabilities. For Ra exceeding 10 to the 7th it is found that the heat-transfer mechanism changes from one characterized by mushroom-like plumes to one consisting of disconnected ascending instabilities, which do not carry with them all the thermal anomaly from the bottom boundary layer. Plume-plume collisions become much more prominent in high Ra situations and have a tendency of generating a pulse-like behavior in the fixed plume. This type of instability represents a distinct mode of heat transfer in the hard turbulent regime. Predictions of this model can be used to address certain issues concerning the mode of time-dependent convection in the earth's mantle.

Hansen, Ulrich↗

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗

EGS Collab Experiment 1: SIMFIP Notch-164 GRL Paper

Characterizing the stimulation mode of a fracture is critical to assess the hydraulic efficiency and the seismic risk related to deep fluid manipulations. We have monitored the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ~1.5 km depth in the crystalline rock of the Sanford Underground Research Facility (USA). The fracture initiates at 61% of the minimum horizontal stress by micro-shearing of the borehole on a foliation plane. As the fluid pressure increases further, borehole axial and radial displacements increase with injection time highlighting the opening and sliding of a new hydrofracture growing ~10 m away from the borehole, in accordance with the ambient normal stress regime and in alignment with the microseismicity. Our study reveals how fluid-driven fracture stimulation can be facilitated by a mixed-mode process controlled by the complex hydromechanical evolution of the growing fracture. The data presented in this submission refer to the SIMFIP measurements and analyses of the stimulation tests conducted on the 164 ft (50 m) notch of the Sanford Underground Research Facility (SURF), during the EGS-Collab test 1. In addition to the datafiles, there is the draft of a manuscript submitted to Geophysical Research Letters (GRL).

15 GEOTHERMAL ENERGY↗