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 451 records · Page 25

Shuttle Continuous Open Wind Tunnel (SCOWT)

The aerodynamic and heat transfer coefficients are studied within the range of thermo-fluid-dynamic conditions experienced by a satellite during Tethered Satellite System (TSS) atmospheric flights. The gasdynamic processes occuring downstream of the bow wave standing in front of the satellite are also studied. The knowledge of the chemistry and physics of the upper atmosphere related to satellite aerothermodynamics is furthered. The existence of an overshooting of the air drag coefficient of the sphere in the transition regime is validated. A complete set of measurements is performed in order to provide the data base to develop and validate theoretical models of free molecule transition flow fields.

Carlomagno, G. M.↗

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗

A High-Fidelity Modeling Tool to Support the Design of Oxy-Combustors for Direct-Fired Supercritical CO2 Cycles

The challenge in the design of oxy-combustors for direct-fired supercritical CO 2 (sCO 2 ) cycles is in addressing disparate performance metrics and objectives. Key design parameters to consider include, among others, injector design for mixing and flame stability, split of recycled CO 2 diluent between injectors and cooling films, target flame temperatures to control noncondensable products, and strategies to inject the diluent CO 2 for film cooling and thermal control. In order to support novel oxy-combustor designs, a high-fidelity yet numerically efficient modeling framework based on the CRUNCH CFD® flow solver is presented, featuring key physics-based submodels relevant in this regime. For computational efficiency in modeling large kinetic sets, a flamelet/progress variable (FPV) based tabulated-chemistry approach is utilized featuring a three-stream extension to allow for the simulation of the CO 2 film cooling stream in addition to the fuel and oxidizer streams. Finite rate chemistry effects are modeled in terms of multiple progress variables for the primary flame as well as for slower-evolving chemical species such as NO x and SO x contaminants. Real fluid effects are modeled using advanced equations of states. The predictive capabilities of this computationally tractable design support tool are demonstrated on a conceptual injector design for an oxy-combustor operating near 30 MPa. Simulations results provide quantitative feedback on the effectiveness of the film cooling as well as the level of contaminants (CO, NO, and N) in the exhaust due to impurities entering from the injectors. These results indicate that this framework would be a useful tool for refining and optimizing the oxy-combustor designs as well as risk mitigation analyses.

Engineering↗

The kinematics of buoyancy induced mixing

A model problem to study mixing driven by buoyancy flow fields due to a time dependent body force which consists of a steady and an oscillatory component is considered. Flow fields generated by a time dependent body force, which consists of a steady and an oscillatory component from vibrations or g-jitter, mix two fluids inside a cavity by stretching and folding the interface into various morphological patterns. These patterns show two basic structures depending on the dominant components of the body force. These structures are whirls and tendrils, corresponding to the dominance of either the steady or oscillatory component, respectively. A combination of these structures also occurs in the interaction region. In this region the effect of the steady component opposes the oscillatory component and the result is to smooth the interface, and prevent breakup of the interface. For the parametric range considered, two basic flow regimes occur: convective and chaotic. In the convective regime, the morphological patterns are topological. In the chaotic regime, mixing occurs by a repetitive sequence of bubble formation, necking, and breakup which spreads the concentration field throughout the cavity. The length stretch of the interface increases exponentially with time. The dynamical system characteristics show that in the convective regime, at a point in space, the flow field is oscillatory with elliptical phase space trajectories, and its power spectrum indicates that the flow field responds to the corresponding input frequency. Whereas, in the chaotic regime, the flow field is aperiodic, the phase space trajectories show irregular patterns, and the power spectrum shows a broadband distribution.

Duval, Walter M. B.↗

Parametric Analysis of Cyclic Phase Change and Energy Storage in Solar Heat Receivers

A parametric study on cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, has been performed. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results in the two-phase regime. Results indicate that parametric changes in receiver gas inlet temperature and receiver heat input effects higher sensitivity to changes in receiver gas exit temperatures.

Hall, Carsie A., III↗

Convergence Analysis of Turbulent Flow Solutions

Data from the "Turbulence Modeling Resource" website for turbulent flow over an NACA-0012 airfoil is analyzed to determine the convergence behavior of three second-order CFD (Computational Fluid Dynamics) codes: CFL3D (Computational Fluids Lab 3 Dimensional flow solver), FUN3D (Fully Unstructured Navier-stokes flow solver), and TAU (German Aerospace Center (DLR) 2 dimensional code for unstructured hybrid grids solving the Reynolds-Averaged Navier-Stokes equations or the Euler equations). The convergence of both integrated properties and pointwise data are examined. Several different methods for estimating errors and computing convergence rates are compared. A high-order extension to the Richardson extrapolation is developed that improves the accuracy of the mesh limit values and provides a quantitative estimate of the threshold of the asymptotic regime. The coefficient of total drag exhibits second-order convergence for all three codes, and convergence is monotone over a sequence of 7 grids. Other force coefficients are not so well behaved. The convergence rates of the viscous component of drag on the three nest grids ranges from 3:0 for CFL3D to 1:0 for FUN3D. The three codes are converging to similar but not identical solutions. The largest differences between the codes are in the coefficient of lift for which the difference between CFL3D and FUN3D is greater than 10 (sup minus 4). The best agreement occurs in the viscous component of drag, which is the only force component for which all three codes are converging towards each other at a rate of second-order. The agreement between the two unstructured grid codes is good with all properties except lift converging towards common values at a rate of second-order. No one code was universally better than the other. The TAU code has the lowest error in total drag, FUN3D has the lowest error in lift, and CFL3D has the lowest error in the viscous component of drag.

Atkins, Harold L.↗

Cosmic-ray transport in accelerating flows

The quasi-linear transport equation of energetic charged particles is derived, including scattering by Alfven waves propagating parallel and antiparallel to a uniform magnetic field, losses by synchrotron radiation, and acceleration of the cold background fluid supporting the waves. As in the comoving frame equations of radiative transfer, scattering, and loss terms are evaluated in the frame in which the fluid is locally at rest. The diffusion approximation is applied to the resulting equation, yielding a transport equation for the isotropic part of the distribution function and a first-order approximation to the anisotropy. The conditions under which this approximation may be applied are discussed. In the nonrelativistic regime, the standard equation of the diffusion approximation is recovered. As contributions previously investigated, new effects appear from the spatial and temporal variations of the fluid velocity as well as the synchrotron radiation terms. These are discussed in detail.

Kirk, John G.↗

A pumped, two-phase flow heat transport system for orbiting instrument payloads

A pumped two-phase (heat absorption/heat rejection) thermal transport system for orbiting instrument payloads is investigated. The thermofluid characteristics necessary for the system design are discussed. A preliminary design with a series arrangement of four instrument heat stations and six radiators in a single loop is described in detail, and the total mass is estimated to be 134 kg, with the radiators, instrument heat stations, and fluid reservoir accounting for approximately 86, 24, and 12 kg, respectively. The evaluation of preliminary test results shows that the system has potential advantages; however, further research is necessary in the areas of one-g and zero-g heat transfer coefficients/fluid regimes, fluid by-pass temperature control, and reliability of small pumps.

Fowle, A. A.↗

Observations of Stably Stratified Flow through a Microscale Gap

This paper reports the findings of a comprehensive field investigation on flow through a mountain gap subject to a range of stably stratified environmental conditions. This study was embedded within the Perdigão field campaign, which was conducted in a region of parallel double-ridge topography with ridge-normal wind climatology. One of the ridges has a well-defined gap (col) at the top, and an array of in situ and remote sensors, including a novel triple Doppler lidar system, was deployed around it. The experimental design was mostly guided by previous numerical and theoretical studies conducted with an idealized configuration where a flow (with characteristic velocity U 0 and buoyancy frequency N ) approaches normal to a mountain of height h with a gap at its crest, for which the governing parameters are the dimensionless mountain height G= (Nh)/U 0 and various gap aspect ratios. Modified forms of G were proposed to account for real-world atmospheric variability, and the results are presented in terms of a gap-averaged value G c . The nature of gap flow was highly dependent on G c , wherein a nearly neutral flow regime (G c <1), a transitional mountain wave regime (G c ~O(1)), and a gap-jetting regime ( G c >O(1)) were identified. The measurements were in broad agreement with previous numerical and theoretical studies on a single ridge with a gap or double ridge topography, although details vary. This is the first ever detailed field study reported on microscale (O(100 m)) gap flows, and provides useful data and insights for future theoretical and numerical studies.

Microscale, Gap, environmental fluid dynamics, win↗

Xcompact3D: An open-source framework for solving turbulence problems on a Cartesian mesh

Xcompact3D is a Fortran 90–95 open-source framework designed for fast and accurate simulations of turbulent flows, targeting CPU-based supercomputers. It is an evolution of the flow solver Incompact3D which was initially designed in France in the mid-90’s for serial processors to solve the incompressible Navier–Stokes equations. Incompact3D was then ported to parallel High Performance Computing (HPC) systems in the early 2010’s. Very recently the capabilities of Incompact3D have been extended so that it can now tackle more flow regimes (from incompressible flows to compressible flows at low Mach numbers), resulting in the design of a new user-friendly framework called Xcompact3D. The present manuscript presents an overview of Xcompact3D with a particular focus on its functionalities, its ready-to-run simulations and a few case studies to demonstrate its impact.

17 WIND ENERGY↗

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗

Block-Structured Operator Inference for Coupled Multiphysics Model Reduction

This work presents a block-structured formulation of Operator Inference as a way to learn structured reduced-order models for multiphysics systems. The approach specifies the governing equation structure for each physics component and the structure of the coupling terms. Once the multiphysics structure is specified, the reduced-order model is learned from snapshot data following the nonintrusive Operator Inference methodology. In addition to preserving physical system structure, which in turn permits preservation of system properties such as stability and second-order structure, the block-structured approach has the advantages of reducing the overall dimensionality of the learning problem and admitting tailored regularization for each physics component. The numerical advantages of the block-structured formulation over a monolithic Operator Inference formulation are demonstrated for aeroelastic analysis, which couples aerodynamic and structural models. For the benchmark test case of the AGARD 445.6 wing, block-structured Operator Inference provides an average 20% online prediction speedup over monolithic Operator Inference across subsonic and supersonic flow conditions in both the stable and fluttering parameter regimes while preserving the accuracy achieved with monolithic Operator Inference.

42 ENGINEERING↗

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transition prediction in hypersonic regime on complex geometries with RANS-based models

In the near future, RANS computations will continue to play a significant role in the design of hypersonic vehicles with complex geometries. Therefore, it is imperative to continue testing, benchmarking, and refining the RANS models. In the present study, we evaluate RANS-like, transport equations-based models for predicting laminar-turbulent transition over a full-scale scale model of the BOLT flight configuration that was tested in the CUBRC LENS-II wind tunnel test facility. Based on the availability of the experimental results, comparisons are made between computations from several pre-existing transition models and Computational Fluid Dynamics (CFD) codes, with an emphasis on using the same computational meshes and flow conditions for all computations. The analysis covers the sensitivity of the transition predictions to the input parameters for five different transition models (four mainstream models designed for low-speed flows and a uniquely high-speed model), grid resolution, and the details of model implementation across three different flow solvers. The results show that the phenomenological models can describe significant aspects of the measured transition front. However, a number of additional improvements are required before these models can offer more reliable estimates of transition in high-speed flows.

CFD↗

Transition Prediction in Hypersonic Regime on Complex Geometries With Rans-Based Models

In the near future, RANS computations will continue to play a significant role in the design of hypersonic vehicles with complex geometries. Therefore, it is imperative to continue testing, benchmarking, and refining the RANS models. In the present study, we evaluate RANS-like, transport equations-based models for predicting laminar-turbulent transition over a full-scale scale model of the BOLT flight configuration that was tested in the CUBRC LENS-II wind tunnel test facility. Based on the availability of the experimental results, comparisons are made between computations from several pre-existing transition models and Computational Fluid Dynamics (CFD) codes, with an emphasis on using the same computational meshes and flow conditions for all computations. The analysis covers the sensitivity of the transition predictions to the input parameters for five different transition models (four mainstream models designed for low-speed flows and a uniquely high-speed model), grid resolution, and the details of model implementation across three different flow solvers. The results show that the phenomenological models can describe significant aspects of the measured transition front. However, a number of additional improvements are required before these models can offer more reliable estimates of transition in high-speed flows.

CFD↗

Simulations of divertor plasmas with inverse sheaths

The effect of strong electron emission from material surfaces has been proposed to form an “inverse sheath”: a region with a positive potential relative to the near-wall plasma which prevents the flow of ions to the wall [M. D. Campanell, “Negative plasma potential relative to electronemitting surfaces,” Phys. Rev. E. 88, 033103 (2013); M. D. Campanell and M. V. Umansky, “Strongly emitting surfaces unable to float below plasma potential,” Phys. Rev. Lett. 116, 1–5 (2016); M. D. Campanell and G. R. Johnson, “Thermionic cooling of the target plasma to a sub-ev temperature,” Phys. Rev. Lett. 122, 1–5 (2019)]. We assess the viability of this regime in a tokamak device using the 2D edge plasma transport code UEDGE [T. Rognlien et al ., “A fully implicit, time dependent 2-D fluid code for modeling tokamak edge plasmas,” J. Nucl. Mater. 196–198, 347–351 (1992)]. Since the UEDGE code does not consider the sheath region directly, we apply boundary conditions at the divertor targets which emulate the physics of both “standard” and “inverse” sheath regimes [R. Masline et al ., “Influence of the inverse sheath on divertor plasma performance in tokamak edge plasma simulations,” Contrib. Plasma Phys. 60, e201900097 (2020)]. Using these boundary conditions, we perform scoping studies to assess plasma parameters near the target by varying the density at the core-edge interface. We observe a smooth transition in the resultant profiles of plasma parameters for the standard sheath, and a bifurcation across the simulation set for plasmas with an inverse sheath. The cause of this bifurcation is assessed by performing the parameter scan both with and without impurity radiation; we observe that the bifurcation persists in both cases, indicating that this bifurcation is caused by plasma recombination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Cyclic Phase Change and Energy Storage in Solar Heat Receivers

Numerical results pertaining to cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, have been reported. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed and results compared with available experimental data. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results for comparisons with GTD data for both the subcooled and two-phase regimes. While qualitative trends were in close agreement for the balanced orbit modes, excellent quantitative agreement was observed for steady-state modes.

Hall, Carsie A., III↗

Decoupling the Roles of Inertia and Gravity on Particle Dispersion

Inertial and gravitational forces determine a particle's motion in a turbulent flow field. Gravity plays the dominant role in this motion by pulling the particles through adjacent regions of fluid turbulence. To better understand and model how a particle's inertia effects its displacement, one must examine the dispersion in a turbulent flow in the absence of gravity. In this paper, we present the particle experiments planned for NASA's KC-135 Reduced-Gravity Aircraft, which generates microgravity conditions for about 20 seconds. We also predict the particle behavior using simulation and ground-based experiments. We will release particles with Stokes numbers of 0.1, 1, and 10 into an enclosed tank of near-isotropic, stationary, and homogenous turbulence. These particle Stoke numbers cover a broad range of flow regimes of interest. Two opposed grids oscillating back and forth generate the turbulent field in the tank with a range of turbulence scales that covers about three orders of magnitude and with turbulence intensities of about ten times the mean velocity. The motion of the particles will be tracked using a stereo image velocimetry technique.

Groszmann, D. E.↗