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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Development of a Trajectory-Centric CFD-RBD Framework for Advanced Multidisciplinary/Multiphysics Simulation

The desire to model increasingly complex unsteady flow phenomena drives coupling of physics-based disciplinary analysis tools, such as coupled aerodynamics-rigid body dynamics simulations. This paper documents the creation of a framework linking the six-degree-of-freedom trajectory propagator POST2 with NASA’s FUN3D computational fluid dynamics flow solver. Cross-code verification between the framework and a CFD-centric 6DOF code is performed using the Army-Navy Finner projectile experiencing unsteady accelerating flow. Free-flight simulations of an entry vehicle ballistic range test are validated against physical and computational experiments.

Zachary J Ernst↗
NASA NTRSDate not supplied

Development of a Trajectory-Centric CFD-RBD Framework for Advanced Multidisciplinary/Multiphysics Simulation

The desire to model increasingly complex unsteady flow phenomena drives coupling of physics-based disciplinary analysis tools, such as coupled aerodynamics-rigid body dynamics simulations. This paper documents the creation of a framework linking the six-degree-of-freedom trajectory propagator POST2 with NASA’s FUN3D computational fluid dynamics flow solver. Cross-code verification between the framework and a CFD-centric 6DOF code is performed using the Army-Navy Finner projectile experiencing unsteady accelerating flow. Free-flight simulations of an entry vehicle ballistic range test are validated against physical and computational experiments.

Zachary J Ernst↗
NASA NTRS2018

Measurements to Understand the Flow Mechanisms Contributing to Tandem-Rotor Outwash

Downwash and outwash characteristics of a model-scale tandem-rotor system in the presence of the ground were analyzed by identifying and understanding the physical mechanisms contributing to the observed flow field behavior. A building block approach was followed in simplifying the problem, separating the effects of the fuselage, effects of one rotor on the other, etc. Flow field velocities were acquired in a vertical plane at four aircraft azimuths of a small-scale tandem rotor system using the particle image velocimetry (PIV) technique for radial distances up to 4 times the rotor diameter. Results were compared against full-scale CH-47D measurements. Excellent correlation was found between the small- and full-scale mean flow fields (after appropriate normalization using rotor and wall jet parameters). Following the scalability analysis, the effect of rotor height on the outwash was also studied. Close to the aircraft, an increase in rotor height above ground decreased the outwash velocity at all aircraft azimuths. However, farther away, the longitudinal and lateral axes of the aircraft showed increasing and decreasing outwash velocities, respectively, with increasing rotor height. Measurements also indicated the presence of large-scale (of the size of the rotor height) shear-layer vortical structures along the ground that could be the source of low-frequency (approximately 1 Hz) flow variation observed in the full-scale measurements. Flow visualization studies and PIV measurements were also made on jets of different sizes to complement the observations made on rotors wherever possible. Baseline rotor measurements were made out-of-ground effect to understand the nature of inflow distribution for realistic rotor configurations and their modified characteristics in the presence of ground. Lastly, a feasibility study on applying high-fidelity CFD simulations for outwash study was conducted using Helios to model an isolated rotor configuration IGE at full scale Reynolds number. The results were encouraging and demonstrated the practical challenges associated with predicting rotor outwash.

Measurements↗
NASA NTRS2015

Measurements to Understand the Flow Mechanisms Contributing to Tandem-Rotor Outwash

Downwash and outwash characteristics of a model-scale tandem-rotor system in the presence of the ground were analyzed by identifying and understanding the physical mechanisms contributing to the observed flow field behavior. A building block approach was followed in simplifying the problem, separating the effects of the fuselage, effects of one rotor on the other, etc. Flow field velocities were acquired in a vertical plane at four aircraft azimuths of a small-scale tandem rotor system using the particle image velocimetry (PIV) technique for radial distances up to 4 times the rotor diameter. Results were compared against full-scale CH-47D measurements. Excellent correlation was found between the small- and full-scale mean flow fields (after appropriate normalization using rotor and wall jet parameters). Following the scalability analysis, the effect of rotor height on the outwash was also studied. Close to the aircraft, an increase in rotor height above ground decreased the outwash velocity at all aircraft azimuths. However, farther away, the longitudinal and lateral axes of the aircraft showed increasing and decreasing outwash velocities, respectively, with increasing rotor height. Measurements also indicated the presence of large-scale (of the size of the rotor height) shear-layer vortical structures along the ground that could be the source of low-frequency (approximately 1 Hz) flow variation observed in the full-scale measurements. Flow visualization studies and PIV measurements were also made on jets of different sizes to complement the observations made on rotors wherever possible. Baseline rotor measurements were made out-of-ground effect to understand the nature of inflow distribution for realistic rotor configurations and their modified characteristics in the presence of ground. Lastly, a feasibility study on applying high-fidelity CFD simulations for outwash study was conducted using Helios to model an isolated rotor configuration IGE at full-scale Reynolds number. The results were encouraging and demonstrated the practical challenges associated with predicting rotor outwash.

Ramasamy, Manikandan↗
DOE OSTI2020

Shear-induced lift force on spheres in a viscous linear shear flow at finite volume fractions

Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.

42 ENGINEERING↗
DOE OSTI2018

Integral Full Core Multi-Physics PWR Benchmark with Measured Data

In recent years, the importance of modeling and simulation has been highlighted extensively in the DOE research portfolio with concrete examples in nuclear engineering with the CASL and NEAMS programs. These research efforts and similar efforts worldwide aim at the development of high-fidelity multi-physics analysis tools for the simulation of current and next-generation nuclear power reactors. Like all analysis tools, verification and validation is essential to guarantee proper functioning of the software and methods employed. The current approach relies mainly on the validation of single physic phenomena (e.g. critical experiment, flow loops, etc.) and there is a lack of relevant multiphysics benchmark measurements that are necessary to validate high-fidelity methods being developed today. This work introduces a new multi-cycle full-core Pressurized Water Reactor (PWR) depletion benchmark based on two operational cycles of a commercial nuclear power plant that provides a detailed description of fuel assemblies, burnable absorbers, in-core fission detectors, core loading and re-loading patterns. This benchmark enables analysts to develop extremely detailed reactor core models that can be used for testing and validation of coupled neutron transport, thermal-hydraulics, and fuel isotopic depletion. The benchmark also provides measured reactor data for Hot Zero Power (HZP) physics tests, boron letdown curves, and three-dimensional in-core flux maps from 58 instrumented assemblies. The benchmark description is now available online and has been used by many groups. However, much work remains to be done on the quantification of uncertainties and modeling sensitivities. This work aims to address these deficiencies and make this benchmark a true non-proprietary international benchmark for the validation of high-fidelity tools. This report details the BEAVRS uncertainty quantification for the first two cycle of operations and serves as the final report of the project.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Interfacial stress balances in structured continua and free surface flows in ferrofluids

Interfacial linear and internal angular momentum balances are obtained for a structured continuum and for the special case of a ferrofluid, a suspension of magnetic nanoparticles in a Newtonian fluid. The interfacial balance equations account for the effects of surface tension and surface tension gradient, magnetic surface excess forces, antisymmetric stresses, and couple stresses in driving interfacial flows in ferrofluids. Application of the interfacial balance equations is illustrated by obtaining analytical expressions for the translational and spin velocity profiles in a thin film of ferrofluid on an infinite flat plate when a rotating magnetic field is applied with axis of rotation parallel to the ferrofluid/air interface. The cases of zero and non-zero spin viscosity are considered for small applied magnetic field amplitude. Expressions for the maximum translational velocity, slope of the translational velocity profile at the ferrofluid/air interface, and volumetric flow rate are obtained and their use to test the relevance of spin viscosity and couple stresses in the flow situation under consideration is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A multiclass vehicular dynamic traffic flow model for main roads and dedicated lanes/roads of multimodal transport network

We present in this paper a model of vehicular traffic flow for a multimodal transportation road network. We introduce the notion of class of vehicles to refer to vehicles of different transport modes. Our model describes the traffic on highways (which may contain several lanes) and network transit for pubic transportation. The model is drafted with Eulerian and Lagrangian coordinates and uses a Logit model to describe the traffic assignment of our multiclass vehicular flow description on shared roads. The paper also discusses traffic streams on dedicated lanes for specific class of vehicles with event-based traffic laws. An Euler-Lagrangian-remap scheme is introduced to numerically approximate the model’s flow equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Physical properties of the inner shocks in hot, tilted black hole accretion flows

Simulations of hot, pressure-supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two nonaxisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10° and 15°. By integrating fluid trajectories in the simulations through the shocks and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction (≅ 3%-12%) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

TU-G-204-01: BEST IN PHYSICS (IMAGING): Dynamic CT Myocardial Perfusion Measurement and Its Comparison to Fractional Flow Reserve

Purpose: To evaluate a first pass analysis (FPA) technique for CT perfusion measurement in a swine animal and its validation using fractional flow reserve (FFR) as a reference standard. Methods: Swine were placed under anesthesia and relevant physiologic parameters were continuously recorded. Intra-coronary adenosine was administered to induce maximum hyperemia. A pressure wire was advanced distal to the first diagonal branch of the left anterior descending (LAD) artery for FFR measurements and a balloon dilation catheter was inserted over the pressure wire into the proximal LAD to create varying levels of stenosis. Images were acquired with a 320-row wide volume CT scanner. Three main coronary perfusion beds were delineated in the myocardium using arteries extracted from CT angiography images using a minimum energy hypothesis. The integrated density in the perfusion bed was used to calculate perfusion using the FPA technique. The perfusion in the LAD bed over a range of stenosis severity was measured. The measured fractional perfusion was compared to FFR and linear regression was performed. Results: The measured fractional perfusion using the FPA technique (P-FPA) and FFR were related as P-FPA = 1.06FFR – 0.06 (r{sup 2} = 0.86). The perfusion measurements were calculated with only three to five total CT volume scans, which drastically reduces the radiation dose as compared with the existing techniques requiring 15–20 volume scans. Conclusion: The measured perfusion using the first pass analysis technique showed good correlation with FFR measurements as a reference standard. The technique for perfusion measurement can potentially make a substantial reduction in radiation dose as compared with the existing techniques.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Multiscale Characterization of Structural Compositional and Textural Heterogeneity of Nano-porous Geomaterials

The purpose of the project was to perform multiscale characterization of low permeability rocks to determine the effect of physical and chemical heterogeneity on the poromechanical and flow responses of shales and carbonate rocks with a broad range of physical and chemical heterogeneity . An integrated multiscale imaging of shale and carbonate rocks from nanometer to centimeter scales include s dual focused ion beam - scanning electron microscopy (FIB - SEM) , micro computed tomography (micro - CT) , optical and confocal microscopy, and 2D and 3D energy dispersive spectroscopy (EDS). In addition, mineralogical mapping and backscattered imaging with nanoindentation testing advanced the quantitative evaluat ion of the relationship between material heterogeneity and mechanical behavior. T he spatial distribution of compositional heterogeneity, anisotropic bedding patterns, and mechanical anisotropy were employed as inputs for brittle fracture simulations using a phase field model . Comparison of experimental and numerical simulations reveal ed that proper incorporation of additional material information, such as bedding layer thickness and other geometrical attributes of the microstructures, can yield improvements on the numerical prediction of the mesoscale fracture patterns and hence the macroscopic effective toughness. Overall, a comprehensive framework to evaluate the relationship between mechanical response and micro-lithofacial features can allow us to make more accurate prediction of reservoir performance by developing a multi - scale understanding of poromechanical response to coupled chemical and mechanical interactions for subsurface energy related activities.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Error and Uncertainty Quantification in the Numerical Simulation of Complex Fluid Flows

The failure of numerical simulation to predict physical reality is often a direct consequence of the compounding effects of numerical error arising from finite-dimensional approximation and physical model uncertainty resulting from inexact knowledge and/or statistical representation. In this topical lecture, we briefly review systematic theories for quantifying numerical errors and restricted forms of model uncertainty occurring in simulations of fluid flow. A goal of this lecture is to elucidate both positive and negative aspects of applying these theories to practical fluid flow problems. Finite-element and finite-volume calculations of subsonic and hypersonic fluid flow are presented to contrast the differing roles of numerical error and model uncertainty. for these problems.

Barth, Timothy J.↗
NASA NTRS2020

A strongly changing accretion morphology during the outburst decay of the neutron star x-ray binary 4U 1608-52

It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below ∼10−2 of the Eddington limit (LEdd). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work, we use NuSTAR and NICER observations obtained during the 2018 accretion outburst of the neutron star LMXB 4U 1608−52, to study changes in the reflection spectrum. We find that the broad Fe–Kα line and Compton hump, clearly seen during the peak of the outburst when the X-ray luminosity is ∼1037 erg s−1 (∼0.05 LEdd), disappear during the decay of the outburst when the source luminosity drops to ∼4.5 × 1035 erg s−1 (∼0.002 LEdd). We show that this non-detection of the reflection features cannot be explained by the lower signal-to-noise ratio at lower flux, but is instead caused by physical changes in the accretion flow. Simulating synthetic NuSTAR observations on a grid of inner disc radius, disc ionization, and reflection fraction, we find that the disappearance of the reflection features can be explained by either increased disc ionization (log ξ  4.1) or a much decreased reflection fraction. A changing disc truncation alone, however, cannot account for the lack of reprocessed Fe–Kα emission. The required increase in ionization parameter could occur if the inner accretion flow evaporates from a thin disc into a geometrically thicker flow, such as the commonly assumed formation of a radiatively inefficient accretion flow at lower mass accretion rates.

J. van den Eijnden↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2018

Steady Flow Regimes in the Problem of Intense Wind-Driven Circulation in a Thin Layer of Viscous Rotating Fluid

The problem of wind-driven circulation in a thin layer of homogeneous viscous rotating fluid is considered. In the limit of small Froude and Burger numbers it is described by two parameters—the Rossby number proportional to the forcing and the Ekman number. For small Ekman numbers and moderate Rossby numbers, the regimes of steady fluid circulation are investigated through numerical simulations of a Kármántype flow. Two different types of steady flows are shown to exist, and their physical interpretation based on a simple analytical model is proposed. A regime diagram on the plane of problem parameters is provided. The simultaneous existence of steady flows of two different types turns out to be possible in some domain of parameters, i.e., there is a dependence of the stationary solution to the problem on initial data. The asymptotic dependence of steady flow characteristics is investigated in the limit of small Ekman numbers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experimental and numerical investigation of reactive shock-accelerated flows

The main goal of this program was to establish a qualitative and quantitative connection, based on the appropriate dimensionless parameters and scaling laws, between shock-induced distortion of astrophysical plasma density clumps and their earthbound analog in a shock tube. These objectives were pursued by carrying out laboratory experiments and numerical simulations to study the evolution of two gas bubbles accelerated by planar shock waves and compare the results to available astrophysical observations. The experiments were carried out in an vertical, downward-firing shock tube, 9.2 m long, with square internal cross section (25×25 cm 2 ). Specific goals were to quantify the effect of the shock strength (Mach number, M) and the density contrast between the bubble gas and its surroundings (usually quantified by the Atwood number, i.e. the dimensionless density difference between the two gases) upon some of the most important flow features (e.g. macroscopic properties; turbulence and mixing rates). The computational component of the work performed through this program was aimed at (a) studying the physics of multi-phase compressible flows in the context of astrophysics plasmas and (b) providing a computational connection between laboratory experiments and the astrophysical application of shock-bubble interactions. Throughout the study, we used the FLASH4.2 code to run hydrodynamical and magnetohydrodynamical simulations of shock bubble interactions on an adaptive mesh.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

Physical Mechanisms of Glaze Ice Scallop Formations on Swept Wings

An experiment was conducted to understand the physical mechanisms that lead to the formation of scallops on swept wings. Icing runs were performed on a NACA 0012 swept wing tip at 45 deg, 30 deg, and 15 deg sweep angles. A baseline case was chosen and direct measurements of scallop height and spacing, castings, video data and close-up photographic data were obtained. The results showed the scallops are made of glaze ice feathers that grow from roughness elements that have reached a minimum height and are located beyond a given distance from the attachment line. This distance depends on tunnel conditions and sweep angle, and is the critical parameter in the formation of scallops. It determines if complete scallops, incomplete scallops or no scallops are going to be formed. The mechanisms of growth for complete and incomplete scallops were identified. The effect of velocity, temperature and LWC on scallop formation was studied. The possibility that cross flow instability may be the physical mechanism that triggers the growth of roughness elements into glaze ice feathers is examined.

Vargas, Mario↗
DOE OSTI2020

Mean-field control for efficient mixing of energy loads

In this work, we pose an engineering challenge of controlling an ensemble of energy devices via coordinated, implementation-light, and randomized on/off switching as a problem in nonequilibrium statistical mechanics. We show that mean-field control with nonlinear feedback on the cumulative consumption, assumed available to the aggregator via direct physical measurements of the energy flow, allows the ensemble to recover from its use in the demand response regime, i.e., transition to a statistical steady state, significantly faster than in the case of the fixed feedback. Moreover when the nonlinearity is sufficiently strong, one observes the phenomenon of “super-relaxation,” where the total instantaneous energy consumption of the ensemble transitions to the steady state much faster than the underlying probability distribution of the devices over their state space, while also leaving almost no devices outside of the comfort zone.

25 ENERGY STORAGE↗