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,153 records · Page 64

Data-driven wind turbine wake modeling via probabilistic machine learning

Wind farm design primarily depends on the variability of the wind turbine wake flows to the atmospheric wind conditions and the interaction between wakes. Physics-based models that capture the wake flow field with high-fidelity are computationally very expensive to perform layout optimization of wind farms, and, thus, data-driven reduced-order models can represent an efficient alternative for simulating wind farms. In this work, we use real-world light detection and ranging (LiDAR) measurements of wind-turbine wakes to construct predictive surrogate models using machine learning. Specifically, we first demonstrate the use of deep autoencoders to find a low-dimensional latent space that gives a computationally tractable approximation of the wake LiDAR measurements. Then, we learn the mapping between the parameter space and the (latent space) wake flow fields using a deep neural network. Additionally, we also demonstrate the use of a probabilistic machine learning technique, namely, Gaussian process modeling, to learn the parameter-space-latent-space mapping in addition to the epistemic and aleatoric uncertainty in the data. Finally, to cope with training large datasets, we demonstrate the use of variational Gaussian process models that provide a tractable alternative to the conventional Gaussian process models for large datasets. Furthermore, we introduce the use of active learning to adaptively build and improve a conventional Gaussian process model predictive capability. Overall, we find that our approach provides accurate approximations of the wind-turbine wake flow field that can be queried at an orders-of-magnitude cheaper cost than those generated with high-fidelity physics-based simulations.

Deep neural networks↗

UV/IR symmetries of the S -matrix and RG flow

The low energy S-matrix which describes non-relativistic scattering arising from finite-range forces has UV/IR symmetries that are hidden in the corresponding effective field theory (EFT) action. Furthermore, it is shown that the S-matrix symmetries are manifest as geometric symmetries of the RG flow of coupling constants in the EFT action in both three and two spatial dimensions. An example is given demonstrating that UV/IR symmetry breaking in the S-matrix implies strong constraints on the RG flow of the coefficients of the corresponding symmetry-breaking operators in the EFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sensitivity of D meson azimuthal anisotropies to system size and nuclear structure

Recent experimental data at the Large Hadron Collider (LHC) confirmed that the soft sector azimuthal anisotropies in central XeXe collisions are sensitive to deformations in the wave function of the Xenon nucleus. Additionally, the CMS experiment found that D meson flow is slightly suppressed compared to other particle species when considering quark number scaling in small systems of pPb compared to PbPb. In this talk we used the D and B Mesons Modular (DAB-MOD) code coupled to Trento+v-USPhydro to calculate the R AA and v n ’s of D mesons using different energy loss models and Langevin techniques. Comparing D mesons R AA and v n in PbPb to XeXe collisions it is found that in central collisions D meson azimuthal anisotropies are sensitive to details in the nuclear structure wave function of Xenon. Additionally, we find that in mid-central collisions the smaller system size of XeXe suppresses D meson flow.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗

Two-Phase Flow Technology Developed and Demonstrated for the Vision for Exploration

NASA s vision for exploration will once again expand the bounds of human presence in the universe with planned missions to the Moon and Mars. To attain the numerous goals of this vision, NASA will need to develop technologies in several areas, including advanced power-generation and thermal-control systems for spacecraft and life support. The development of these systems will have to be demonstrated prior to implementation to ensure safe and reliable operation in reduced-gravity environments. The Two-Phase Flow Facility (T(PHI) FFy) Project will provide the path to these enabling technologies for critical multiphase fluid products. The safety and reliability of future systems will be enhanced by addressing focused microgravity fluid physics issues associated with flow boiling, condensation, phase separation, and system stability, all of which are essential to exploration technology. The project--a multiyear effort initiated in 2004--will include concept development, normal-gravity testing (laboratories), reduced gravity aircraft flight campaigns (NASA s KC-135 and C-9 aircraft), space-flight experimentation (International Space Station), and model development. This project will be implemented by a team from the NASA Glenn Research Center, QSS Group, Inc., ZIN Technologies, Inc., and the Extramural Strategic Research Team composed of experts from academia.

Sankovic, John M.↗

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

D 0 -Meson R AA in PbPb Collisions at s N N = 5.02 TeV and Elliptic Flow in pPb Collisions at s N N = 8.16 TeV with CMS

The study of charm production in heavy-ion collisions is considered an excellent probe for the properties of the hot and dense medium created in heavy-ion collisions. Measurements of D 0 -meson nuclear modification factor can provide strong constraints into the mechanisms of in-medium energy loss and charm flow in the medium. The measurement of D0-meson elliptic flow in pPb collisions helps us understand the strength of charm quarks coupling to significantly reduced systems which demonstrate hydrodynamic properties. In this paper, the measurements of the D 0 -meson nuclear modification factor in PbPb collisions at 5.02 TeV together with the new measurement of D 0 -meson elliptic flow in high multiplicity pPb collisions at 5.02 TeV using the two-particle correlation method will be presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Flame dynamics in a micro-channeled combustor

The increasing use of Micro-Electro-Mechanical Systems (MEMS) has generated a significant interest in combustion-based power generation technologies, as a replacement of traditional electrochemical batteries which are plagued by low energy densities, short operational lives and low power-to-size and power-to-weight ratios. Moreover, the versatility of integrated combustion-based systems provides added scope for combined heat and power generation. This paper describes a study into the dynamics of premixed flames in a micro-channeled combustor. The details of the design and the geometry of the combustor are presented in the work by Kariuki and Balachandran [1]. This work showed that there were different modes of operation (periodic, a-periodic and stable), and that in the periodic mode the flame accelerated towards the injection manifold after entering the channels. The current study investigates these flames further. We will show that the flame enters the channel and propagates towards the injection manifold as a planar flame for a short distance, after which the flame shape and propagation is found to be chaotic in the middle section of the channel. Finally, the flame quenches when it reaches the injector slots. The glow plug position in the exhaust side ignites another flame, and the process repeats. It is found that an increase in air flow rate results in a considerable increase in the length (and associated time) over which the planar flame travels once it has entered a micro-channel, and a significant decrease in the time between its conversion into a chaotic flame and its extinction. It is well known from the literature that inside small channels the flame propagation is strongly influenced by the flow conditions and thermal management. An increase of the combustor block temperature at high flow rates has little effect on the flame lengths and times, whereas at low flow rates the time over which the planar flame front can be observed decreases and the time of existence of the chaotic flame increases. The frequency of re-ignition of successive flames decreases at higher flow rates and increases at higher temperatures. The data and results from this study will not only help the development of new micro-power generation devices, but they will also serve as a validation case for combustion models capable of predicting flame behavior in the presence of strong thermal and flow boundary layers, a situation common to many industrial applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Framework for Integrated Component and System Analyses of Instabilities

Instabilities associated with fluid handling and operation in liquid rocket propulsion systems and test facilities usually manifest themselves as structural vibrations or some form of structural damage. While the source of the instability is directly related to the performance of a component such as a turbopump, valve or a flow control element, the associated pressure fluctuations as they propagate through the system have the potential to amplify and resonate with natural modes of the structural elements and components of the system. In this paper, the authors have developed an innovative multi-level approach that involves analysis at the component and systems level. The primary source of the unsteadiness is modeled with a high-fidelity hybrid RANS/LES based CFD methodology that has been previously used to study instabilities in feed systems. This high fidelity approach is used to quantify the instability and understand the physics associated with the instability. System response to the driving instability is determined through a transfer matrix approach wherein the incoming and outgoing pressure and velocity fluctuations are related through a transfer (or transmission) matrix. The coefficients of the transfer matrix for each component (i.e. valve, pipe, orifice etc.) are individually derived from the flow physics associated with the component. A demonstration case representing a test loop/test facility comprised of a network of elements is constructed with the transfer matrix approach and the amplification of modes analyzed as the instability propagates through the test loop.

Ahuja, Vineet↗

Representing lateral groundwater flow from land to river in Earth system models

Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.

54 ENVIRONMENTAL SCIENCES↗

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ice velocity at the ice front of the Filchner-Ronne Ice Shelf, Antarctica, as observed with ERS interferometry

ERS images of the two ends of the ice front of the Filchner-Ronne ice shelf (Antarctica) were utilized interferometrically to study the deformation rate of the ice shelf in response to viscous creep. On the western bank of Berkner Island (BI), near Hemmen Ice Rise (HIR), a time series of ERS data acquired in February 1992, in both ascending and descending mode, and with a three-day time interval, were utilized to map the ice velocity in two dimensions. Finite-element ice-shelf flow simulations are compared with the ERS interferograms to interpret the ice motion in terms of the physical constraints on ice-shelf flow. The efforts to fit artificial interferograms generated with model velocity output suggest that the flow regime is strongly influenced by three processes. First, a void-creation process responsible for rifts at coastal margins tends to uncouple the ice shelf from the ice rise and neighboring coast of BI. Secondly, sea ice within the void space appears to act as a binding agent between discrete ice-shelf fragments, allowing rigid-body rotations of these fragments. Third, strain rates appear to be enhanced in a narrow zone adjacent to HIR, implying significant strain softening along the boundary. It is believed that synthetic aperture radar (SAR) interferogram/model intercomparison represents a powerful impetus towards the development of better, more physically realistic ice-shelf flow models.

Rignot, Eric↗

Charge-dependent flow induced by electromagnetic fields in heavy ion collisions

The colliding heavy ions create extremely strong magnetic and electric fields that significantly affect the evolution of the produced quark-gluon plasma (QGP). The knowledge of these fields is essential for establishing the role of topological fluctuations in the QGP through the chiral magnetic effect and related anomaly-induced phenomena. In this talk, we describe our work on the evolution of the QGP in electric and magnetic fields in the framework of hydrodynamics supplemented, in a perturbative fashion, by the dynamical electromagnetism. The evolution of the QGP fluid is described within the iEBE-VISHNU framework. We find that the electromagnetically induced currents result in a charge-odd directed flow Δv 1 and a charge-odd Δv 3 flow both of which are odd in rapidity. Furthermore, while the predicted magnitude of these charge-odd flows agrees with the data from RHIC and LHC, the sign of the predicted asymmetry between the flows of positive and negative hadrons is opposite to the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding and controlling capsule symmetry in near vacuum hohlraums at the National Ignition Facility

The near vacuum hohlraum platform is an inertial confinement fusion design at the National Ignition Facility (NIF) that uses the lowest practical density of helium gas of 30 μg/cc to fill the hohlraum, which is ten times lower than now used routinely. This has several advantages, such as high laser coupling; however, the inability to understand and simulate the symmetry of the imploded capsule has limited the use of this platform. This work presents the first simulations that are able to accurately capture the highly prolate implosion seen experimentally without unphysical, ad hoc model changes. While previous investigations attributed this asymmetry to multi-species interpenetration in the hohlraum, we find that this alone has little effect on symmetry. Instead, it is the presence of crossed-beam energy transfer (CBET), occurring with no applied wavelength shift between the laser beams, that increases the laser power to the inner cones and causes a more prolate implosion. The effect of CBET is increased in the simulation model when the hohlraum laser entrance hole hardware is included. Using this understanding, CBET is exploited by shifting the inner-beam wavelength by –0.75 Å (at 1ω) with respect to the outer-beams. This transfers laser power to the outer-beams in contrast to positive wavelength shifts as done routinely on NIF and produces a round capsule implosion in our simulations. Furthermore, this work shows the possibility of the near vacuum hohlraum as a viable experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An assessment of some vorticity field-leading edge interactions

A review is undertaken of some recent advances in the following areas: the encounter of distributed and concentrated vorticity at a single frequency with a leading-edge, the interaction of multiple concentrations of vorticity with an edge, the distortion due to finite scale and angle of attack of the leading-edge, and the penetration of the leading-edge interaction into the near-wake region. Time-sequence flow visualization and phase-coordinated pressure measurements provide physical insight into the mechanisms of flow distortion and loading. Promising paths for future research are suggested on the basis of recent progress.

Rockwell, D.↗