Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow dynamics”

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

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

Detonation↗

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

detonation↗

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗

Validation of blade-resolved computational fluid dynamics for a MW-scale turbine rotor in atmospheric flow

In this study, simulation results of two different computational fluid dynamics codes, Nalu-Wind and EllipSys3D, are presented for a wind turbine rotor in complex yawed and sheared inflow. The results are compared to measurements from the DanAero experiments, to validate computed pressures and azimuthal trends. Despite different code methodologies and grid setups, the codes agree well in computed pressures and integrated forces along the blade for all blade azimuthal positions, however with some discrepancy in the very yawed case. Additionally, both codes capture well the azimuthal trends and force levels seen in measurements. Investigation into discrepancies shows that expanding grids before the rotor, lead to smearing of the wind profiles, which is likely the main cause of the differences in the results between the codes. Additionally, omission of the ground constraint cause discrepancies in relative velocity seen by the passing blade, due to an over speeding beneath the rotor.

17 WIND ENERGY↗

Power flows and Mechanical Intensities in structural finite element analysis

The identification of power flow paths in dynamically loaded structures is an important, but currently unavailable, capability for the finite element analyst. For this reason, methods for calculating power flows and mechanical intensities in finite element models are developed here. Formulations for calculating input and output powers, power flows, mechanical intensities, and power dissipations for beam, plate, and solid element types are derived. NASTRAN is used to calculate the required velocity, force, and stress results of an analysis, which a post-processor then uses to calculate power flow quantities. The SDRC I-deas Supertab module is used to view the final results. Test models include a simple truss and a beam-stiffened cantilever plate. Both test cases showed reasonable power flow fields over low to medium frequencies, with accurate power balances. Future work will include testing with more complex models, developing an interactive graphics program to view easily and efficiently the analysis results, applying shape optimization methods to the problem with power flow variables as design constraints, and adding the power flow capability to NASTRAN.

Hambric, Stephen A.↗

A coupled fluid-structure dynamic analysis for water impact loads

In this paper a nonlinear transient, hydroelastic method is developed for response analysis of the Space Shuttle recoverable booster on water impact. At each time step, coupled equations of fluid flow and structure dynamics are solved with an iterative, self-correcting process. Rigid body motions are represented by large angle equations. Local deformations are represented by a NASTRAN-generated model of the three-dimensional structure. The fluid is represented by series solutions of the potential flow equation which include the effects of the local structure motions. The resulting water loads and structural stresses obtained will aid in qualification of the final structural design. Comparisons are made to experimental impact data to validate the method.

Herting, D. N.↗

The source of the Earth's long wavelength geoid anomalies: Implications for mantle and core dynamics

The long wavelength components of the Earth's gravity field result mainly from density contrasts associated with convection in the mantle. Direct interpretation of the geoid for mantle convection is complicated by the fact that convective flow results in dynamically maintained deformation of the surface of the Earth, the core mantle boundary (CMB), and any interior chemical boundaries which might exist. These boundary deformations effect the geoid opposite in sign and are comparable in magnitude to those of the interior density contrasts driving the flow. The total difference of two relatively large quantities.

Hager, B. H.↗

Experimental study of vortex breakdown in a cylindrical, swirling flow

The stability of a steady, vortical flow in a cylindrical container with one rotating endwall has been experimentally examined to gain insight into the process of vortex breakdowwn. The dynamics of the flow are governed by the Reynolds number (Re) and the aspect ratio of the cylinder. Re is given by Omega R(sup 2)/nu, where Omega is the speed of rotation of the endwall, R is the cylinder radius, and nu is the kinematic viscosity of the fluid filling the cylinder. The aspect ratio is H/R, where H is the height of the cylinder. Numerical simulation studies disagree whether or not the steady breakdown is stable beyond a critical Reynolds number, Re(sub c). Previous experimental researches have considered the steady and unsteady flows near Re(sub c), but have not explored the stability of the steady breakdown structures beyond this value. In this investigation, laser induced fluorescence was utilized to observe both steady and unsteady vortex breakdown at a fixed H/R of 2.5 with Re varying around Re(sub c). When the Re of a steady flow was slowly increased beyond Re(sub c), the breakdown structure remained steady even though unsteadiness was possible. In addition, a number of hysteresis events involving the oscillation periods of the unsteady flow were noted. The results show that both steady and unsteady vortex breakdown occur for a limited range of Re above Re(sub c). Also, with increasing Re, complex flow transformations take place that alter the period at which the unsteady flow oscillates.

Stevens, J. L.↗

Preliminary efforts toward development of data handling and analysis software for unsteady flow measurements: An application for aeroelastic transonic flow configurations

A few years ago the Structural Dynamics Division at LaRC started ambitious experimental research efforts known as the Benchmark Models Program. The primary objective of this program was to provide experimental data that may serve as a calibration source for computational fluid dynamics (CFD) efforts that deal with aeroelastic unsteady flow configurations. It also focuses on the understanding of complex flow phenomenon associated with unsteady flow developments. The overall plan for the program has been described by Bennett, including a presentation of initial test results of flutter of a rigid wing mounted on flexible supports. An example of a test model employed to measure the dynamic response along with corresponding pressure distributions is shown. This model incorporates eighty pressure transducers distributed along two spanwise stations. In addition, the models are equipped with four accelerometers and two strain gages. The data handling system for the Benchmark Model Program is under development. Several interactive computer routines designed for the user interface, dynamic memory allocation, unsteady flow measurements data extraction, and further data processing were developed. To present a few examples of measured data, the unsteady pressure distributions and the wing model dynamic response were plotted.

Finaish, Fathi↗

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗

Comparative analysis of internal energy excitation and dissociation of nitrogen predicted by independently developed ab initio potential energy surfaces

In this article we present a comparative atomic level study analyzing the vibrational excitation and dissociation of molecular nitrogen due to N 2 ⁢( 1 Σ g ⁢ + )+N⁢( 4 S u ) and N 2 ⁢( 1 Σ g ⁢ + )+N 2 ⁢( 1 Σ g ⁢+ ) interactions governed by independently developed potential energy surfaces at the University of Minnesota and NASA Ames Research Center. Further, vibrational excitation was studied for N 2 +N 2 interactions from T=10000 to 25000 K and for N 2 +N from T=5000 to 30000 K. Nonequilibrium dissociation is studied from T=10000 to 30000 K under the quasi-steady-state condition for N 2 +N 2 and N 2 +N interactions. Finally, an inviscid Mach 20 dissociating nitrogen flow over a cylinder with a Knudsen number of 0.015 is carried out to study the impact of molecular interactions predicted by independently developed potential energy surfaces on a canonical hypersonic flow.

74 ATOMIC AND MOLECULAR PHYSICS↗

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

Global Qualitative Flow-Path Modeling for Local State Determination in Simulation and Analysis

For qualitative modeling and analysis, a general qualitative abstraction of power transmission variables (flow and effort) for elements of flow paths includes information on resistance, net flow, permissible directions of flow, and qualitative potential is discussed. Each type of component model has flow-related variables and an associated internal flow map, connected into an overall flow network of the system. For storage devices, the implicit power transfer to the environment is represented by "virtual" circuits that include an environmental junction. A heterogeneous aggregation method simplifies the path structure. A method determines global flow-path changes during dynamic simulation and analysis, and identifies corresponding local flow state changes that are effects of global configuration changes. Flow-path determination is triggered by any change in a flow-related device variable in a simulation or analysis. Components (path elements) that may be affected are identified, and flow-related attributes favoring flow in the two possible directions are collected for each of them. Next, flow-related attributes are determined for each affected path element, based on possibly conflicting indications of flow direction. Spurious qualitative ambiguities are minimized by using relative magnitudes and permissible directions of flow, and by favoring flow sources over effort sources when comparing flow tendencies. The results are output to local flow states of affected components.

Malin, Jane T.↗

On the Pathways of the Return Flow of the Meridional Overturning Circulation in the Tropical Atlantic

A numerical model of the tropical Atlantic ocean is used to investigate the upper layer pathways of the Meridional Overturning Circulation (MOC) in the tropical Atlantic. The main focus of this thesis is on those parts of the tropical circulation that are thought to be important for the MOC return flow, but whose dynamics have not been understood yet. It is shown how the particular structure of the tropical gyre and the MOO act to inhibit the flow of North Atlantic water into the equatorial thermocline. As a result, the upper layers of the tropical Atlantic are mainly fed by water from the South Atlantic. The processes that carry the South Atlantic water across the tropical Atlantic into the North Atlantic as part of the MOO are described here, and three processes that were hitherto not understood are explained as follows: The North Brazil Current rings are created as the result of the reflection of Rossby waves at the South American coast. These Rossby waves are generated by the barotropically unstable North Equatorial Countercurrent. The deep structure of the rings can be explained by merger of the wave's anticyclones with the deeper intermediate eddies that are generated as the intermediate western boundary current crosses the equator. The bands of strong zonal velocity in intermediate depths along the equator have hitherto been explained as intermediate currents. Here, an alternative interpretation of the observations is offered: The Eulerian mean flow along the equator is negligible and the observations are the signature of strong seasonal Rossby waves. The previous interpretation of the observations can then be explained as aliasing of the tropical wave field. The Tsuchyia Jets are driven by the Eliassen-Palm flux of the tropical instability waves. The equatorial current system with its strong shears is unstable and generates tropical instability waves.

Jochum, Markus↗

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗

Controls From Above and Below: Snow, Soil, and Steepness Drive Diverging Trends of Subsurface Water and Streamflow Dynamics

ABSTRACT The importance of subsurface water dynamics, such as water storage and flow partitioning, is well recognised. Yet, our understanding of their drivers and links to streamflow generation has remained elusive, especially in small headwater streams that are often data‐limited but crucial for downstream water quantity and quality. Large‐scale analyses have focused on streamflow characteristics across rivers with varying drainage areas, often overlooking the subsurface water dynamics that shape streamflow behaviour. Here we ask the question: What are the climate and landscape characteristics that regulate subsurface dynamic storage, flow path partitioning, and dynamics of streamflow generation in headwater streams? To answer this question, we used streamflow data and a widely‐used hydrological model (HBV) for 15 headwater catchments across the contiguous United States. Results show that climate characteristics such as aridity and precipitation phase (snow or rain) and land attributes such as topography and soil texture are key drivers of streamflow generation dynamics. In particular, steeper slopes generally promoted more streamflow, regardless of aridity. Streams in flat, rainy sites (< 30% precipitation as snow) with finer soils exhibited flashier regimes than those in snowy sites (> 30% precipitation as snow) or sites with coarse soils and deeper flow paths. In snowy sites, less weathered, thinner soils promoted shallower flow paths such that discharge was more sensitive to changes in storage, but snow dampened streamflow flashiness overall. Results here indicate that land characteristics such as steepness and soil texture modify subsurface water storage and shallow and deep flow partitioning, ultimately regulating streamflow response to climate forcing. As climate change increases uncertainty in water availability, understanding the interacting climate and landscape features that regulate streamflow will be essential to predict hydrological shifts in headwater catchments and improve water resources management.

Kerins, Devon [Department of Civil and Environment↗

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynmcal hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothes pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at lentgh-scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto-geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core- mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order l20 years is pseudo-scale-invarient. Magnetic spectra of other planets may differ; however, if a transition to non-conductmg fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗