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,549 records · Page 86

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗

Dynamics of face and annular seals with two-phase flow

A detailed study was made of face and annular seals under conditions where boiling, i.e., phase change of the leaking fluid, occurs within the seal. Many seals operate in this mode because of flashing due to pressure drop and/or heat input from frictional heating. High pressure, water pumps, industrial chemical pumps, and cryogenic pumps are mentioned as a few of many applications. The initial motivation was the LOX-GOX seals for the space shuttle main engine, but the study was expanded to include any face or annular seal where boiling occurs. Some of the distinctive behavior characteristics of two-phase seals were discussed, particularly their axial stability. While two-phase seals probably exhibit instability to disturbances of other degrees of freedom such as wobble, etc., under certain conditions, such analyses are too complex to be treated at present. Since an all liquid seal (with parallel faces) has a neutral axial stiffness curve, and is stabilized axially by convergent coning, other degrees of freedom stability analyses are necessary. However, the axial stability behavior of the two-phase seal is always a consideration no matter how well the seal is aligned and regardless of the speed. Hence, axial stability is thought of as the primary design consideration for two-phase seals and indeed the stability behavior under sub-cooling variations probably overshadows other concerns. The main thrust was the dynamic analysis of axial motion of two-phase face seals, principally the determination of axial stiffness, and the steady behavior of two-phase annular seals. The main conclusions are that seals with two-phase flow may be unstable if improperly balanced. Detailed theoretical analyses of low (laminar) and high (turbulent) leakage seals are presented along with computer codes, parametric studies, and in particular a simplified PC based code that allows for rapid performance prediction. A simplified combined computer code for the performance prediction over the laminar and turbulent ranges of a two-phase seal is described and documented. The analyses, results, and computer codes are summarized.

Hughes, William F.↗

Arcjet load characteristics

Experiments were conducted to define the interface characteristics and constraints of 1 kW class arcjets run on simulated decomposition products of hydrazine and power processors. The impacts of power supply output current ripple on arcjet performance were assessed by variation of the ripple frequency from 100 Hz to 100 kHz with 10 percent peak-to-peak ripple amplitude at 1.2 kW. Ripple had no significant effects on thrust, specific impulse or efficiency. The impact of output ripple on thruster lifetime was not assessed. The static and dynamic impedances of the arcjet were quantified with two thrusters of nearly identical configuration. Superposition of an AC component on the DC arc current was used to characterize the dynamic impedance as a function of flow rate and DC current level. A mathematical model was formulated from these data. Both the static and dynamic impedance magnitude were found to be dependent on mass flow rate. The amplitude of the AC component was found to have little effect on the dynamic impedance. Reducing the DC level from 10 to 8 amps led to a large change in the magnitude of the dynamic impedance with no observable phase change. The impedance data compared favorably between the two thrusters.

Hamley, John A.↗

Arcjet load characteristics

Experiments were conducted to define the interface characteristics and constraints of 1 kW class arcjets run on simulated decomposition products of hydrazine and power processors. The impacts of power supply output current ripple on arcjet performance were assessed by variation of the ripple frequency from 100 Hz to 100 kHz with 10 percent peak-to-peak ripple amplitude at 1.2 kW. Ripple had no significant effects on thrust, specific impulse or efficiency. The impact of output ripple on thruster lifetime was not assessed. The static and dynamic impedances of the arcjet were quantitied with two thrusters of nearly identical configuration. Superposition of an AC component on the DC arc current was used to characterize the dynamic impedance as a function of flow rate and DC current level. A mathematical model was formulated from these data. Both the static and dynamic impedance magnitude were found to be dependent on mass flow rate. The amplitude of the AC component was found to have little effect on the dynamic impedance. Reducing the DC level from 10 to 8 amps led to a large change in the magnitude of the dynamic impedance with no observable phase change. The impedance data compared favorably between the two thrusters.

Hamley, John A.↗

Unsteady flow in cavitating turbopumps

Unsteady flow in a cavitating axial inducer pump is analyzed with the help of a simple two-dimensional cascade model. This problem was motivated by a desire to study the effect of unsteady cavitation on the so-called POGO instability in the operation of liquid rocket engines. Here, an important feature is a closed loop coupling between several different modes of oscillation, one of which is due to the basic unsteady characteristics of the cavitation itself. The approaching and leaving flow velocities up- and downstream of the inducer oscillate, and the cavity-blade system participates dynamically with the basic pulsating flow. In the present work, attention is focused on finding a transfer matrix that relates the set of upstream variables to those downstream. This quantity, which is essentially equivalent to cavitation compliance in the quasi-static analyses, is found to be complex and frequency dependent.

Kim, J. H.↗

Spherical accretion onto quasars

The problem of time-steady accretion of gas irradiated by a self-consistently generated quasar-like continuum is investigated. A controversy in the literature is resolved by a careful distinction between the existence of solutions and their stability. Specific criteria are found for when thermal instabilities may disrupt single-phase steady flow. As a function of initial temperature (1,000,000-100,000,000 K), the boundaries of stable flow in the luminosity are determined relative to the Eddington (L/LE)-efficiency plane. The thermal and dynamical history of these steady flows is analyzed. The observational data that define quasars and type I Seyfert galaxies are matched to find values of efficiency and L/LE required for consistent spherically symmetric steady-flow models. The efficiency is constrained to be less than 1/100, and L/LE resides between several times 1/1000-1/100.

Krolik, J. H.↗

Estimating unsteady aerodynamic forces on a cascade in a three-dimensional turbulence field

An analytical method has been developed to estimate the unsteady aerodynamic forces caused by flow field turbulence on a wind tunnel turning vane cascade system (vane set). This method approximates dynamic lift and drag by linearly perturbing the appropriate steady state force equations, assuming that the dynamic loads are due only to free stream turbulence and that this turbulence is homogeneous, isotropic, and Gaussian. Correlation and unsteady aerodynamic effects are also incorporated into the analytical model. Using these assumptions, equations relating dynamic lift and drag to flow turbulence, mean velocity, and vane set geometry are derived. From these equations, estimates for the power spectra and rms (root mean squared value, delta) loading of both lift and drag can be determined.

Norman, T.↗

Estimating unsteady aerodynamic forces on a cascade in a three-dimensional turbulence field

An analytical method has been developed to estimate tne unsteady aerodynamic forces caused by flow field turbulence on a wind tunnel turning vane cascade system (vane set). This method approximates dynamic lift and drag by linearly perturbing the appropriate steady state force equations, assuming that the dynamic loads are due only to free stream turbulence and that this turbulence is homogeneous, isotropic, and Gaussian. Correlation and unsteady aerodynamic effects are also incorporated into the analytical model. Using these assumptions, equations relating dynamic lift and drag to flow turbulence, mean velocity, and vane set geometry are derived. From these equations, estimates for the power spectra and rms (root mean squared value, delta) loading of both lift and drag can be determined.

Norman, T. R.↗

Preliminary Low-Speed Wind-Tunnel Investigation of Some Aspects of the Aerodynamic Problems Associated with Missiles Carried Externally in Positions Near Airplane Wings

A low-speed wind-tunnel investigation has been made of some aspects of the aerodynamic problems associated with the use of air-to-air missiles when carried externally on aircraft. Measurements of the forces and moments on a missile model for a range of positions under the mid-semispan location of a 45deg sweptback wing indicated longitudinal and lateral forces with regard to both carriage and release of the missiles. Surveys of the characteristics of the flow field in the region likely to be traversed by the missiles showed abrupt gradients in both flow angularity and in local dynamic pressure. Through the use of aerodynamic data on the isolated missile and the measured flow-field characteristics, the longitudinal forces and moments acting on the missile while in the presence of the wing-fuselage combination could be estimated with fair accuracy. Although the lateral forces and moments predicted were qualitatively correct, there existed some large discrepancies in absolute magnitude.

Alford, William J., Jr.↗

Comparison of CFD and DSMC Using Calibrated Transport Parameters

Hypersonic re-entry flows span a wide range of length scales where regions of both rarefied and continuum flow exist. Traditional computational fluid dynamics (CFD) techniques do not provide an accurate solution for the rarefied regions of such ‘mixed’ flow fields. Although direct simulation Monte Carlo (DSMC) can be used to accurately capture both the continuum and rarefied features of ‘mixed’ flow fields, they are computationally expensive when employed to simulate the low Knudsen number continuum regimes. Thus, a hybrid framework for seamlessly combining the two methodologies, CFD and DSMC, continues to be a topic of significant research effort. Ensuring consistency in the reaction kinetics and transport models employed within CFD and DSMC is a crucial requirement for obtaining a reliable solution from a hybrid framework for combined continuum/rarefied high speed flows. This paper represents one of the first studies to utilize the calibrated transport parameters developed to ensure consistency between CFD and DSMC solvers. The new variable soft sphere (VSS) parameters are compared to both previous “standard” variable hard sphere (VHS) parameters and also to solutions from the CFD transport properties that the new parameters were developed to reproduce.

Liechty, D. S.↗

Holistic Small-Signal Stability Analysis for Large-Scale Inverter-Intensive Power Systems with Coupled and Full-Order Dynamics from Control Systems and Power Networks

The increasing penetration of inverter-based resources (IBRs) into the existing power systems introduces tremendous benefits for enhanced sustainability but also poses inevitable challenges in terms of insufficient inertia, potential instability, and complex network dynamics, among others. However, the additional coupling introduced by the interactions among gridfollowing (GFL) and grid-forming (GFM) IBRs and the other components (i.e., synchronous generators [SGs], loads, and network, etc.) has not been clearly explored. A holistic, scalable, and quantitative stability analysis framework with the control systems and power networks is still missing. Here, in this paper, to fill in the technical gaps, a holistic small-signal model of the entire system with both rotating generation units and IBRs is established. An extended power flow model with operation dynamics from both generator control schemes and power networks is proposed to provide the varying steady-state operating points for small-signal modeling. The proposed method is compared with MATLAB solvers, and the results show that the proposed approach has a minimum calculation time, which can be less than 12 seconds for a large-scale power system with up to 2,000 buses. Furthermore, a quantitative method is developed to identify the impacts of IBRs on system performance with emphases on the potential stability issues with GFL IBRs, additional benefits of employing GFM IBRs, the feasibility of replacing SGs with GFM IBRs, and the impact of penetration level of different kinds of generation units. Finally, a field island power system is used to verify the proposed approach, and hardware-in-the-loop (HIL) tests are provided to further demonstrate the effectiveness of the proposed analysis.

14 SOLAR ENERGY↗

Unshrouded Centrifugal Turbopump Impeller

The ratio of rocket engine thrust to weight is a limiting constraint in placing more payload into orbit at a low cost. A key component of an engine's overall weight is the turbopump weight, Reducing the turbopump weight can result in significant engine weight reduction and hence, increased delivered payload. There are two main types of pumps: centrifugal and axial pumps. These types of pumps can be further sub-divided into those with shrouds and those without shrouds (unshrouded pumps). Centrifugal pumps can achieve the same pump discharge pressure as an axial pump and it requires fewer pump stages and lower pump weight than an axial pump. Also, with unshrouded centrifugal pumps (impeller), the number of stages and weight can be further reduced. However. there are several issues with regard to using an unshrouded impeller: 1) there is a pump performance penalty due to the front open face recirculation flow; 2) there is a potential pump axial thrust problem from the unbalanced front open face and the back shroud face; and, 3) since test data is very linu'ted for this configuration, there is uncertainty in the magnitude and phase of the rotordynamics forces due to the front impeller passage. The purpose of the paper is to discuss the design of an unshrouded impeller and to examine the design's hydrodynamic performance, axial thrust, and rotordynamics performance. The design methodology will also be discussed. This work will help provide some guidelines for unshrouded impeller design. In particular, the paper will discuss the design of three unshrouded impellers - one with 5 full and 5 partial blades (5+5). one with 6+6 blades and one with 8+8 blades. One of these designs will be selected for actual fabrication and flow test. Computational fluid dynamics (CFD) is used to help design and optimize the unshrouded impeller. The relative pump performance penalty is assessed by comparing the CFD results of the unshrouded impeller with the equivalent shrouded impeller for a particular design. Limited unshrouded - versus - shrouded impeller data from the J-2 pump is used to anchor the CFD. Since no detailed impeller blade force data is available, axial thrust and rotordynamic force predictions are based on the CFD model. For the axial thrust, the impeller front flow passage axial force is integrated from the CFD results and compared to the equivalent shrouded impeller axial force. For the rotordynamics forces, the fluid reaction forces are computed from unsteady flow CFD results using a moving boundary method; the rotor- shaft is moved at several whirl-to-speed frequency ratios to extract the rotordynamics coefficients.

Prueger, George↗

An Interactive Tool for Discrete Phase Analysis in Two-Phase Flows

Under a NASA MSFC SBIR Phase 1 effort an interactive software package has been developed for the analysis of discrete (particulate) phase dynamics in two-phase flows in which the discrete phase does not significantly affect the continuous phase. This package contains a Graphical User Interface (based on the X Window system and the Motif tool kit) coupled to a particle tracing program, which allows the user to interactively set up and run a case for which a continuous phase grid and flow field are available. The software has been applied to a solid rocket motor problem, to demonstrate its ease of use and its suitability for problems of engineering interest, and has been delivered to NASA Marshall Space Flight Center.

Dejong, Frederik J.↗

Focused Laser Differential Interferometry: Recent Developments and Applications for Flow Measurements

Focused laser differential interferometry (FLDI), first introduced in the 1970s, has developed into a powerful tool for the analysis of turbulent dynamics and fluctuations in fluid flows ranging from subsonic to hypersonic. This paper provides a review of recent FLDI developments and applications, including discussions on the basic theory of the instrument, various forms of FLDI, multi-point techniques, utilization in fluid flow measurements, quantitative measurement capabilities, and more.

Elizabeth K Benitez↗

Organizational Principles of Hyporheic Exchange Flow and Biogeochemical Cycling in River Networks Across Scales

Hyporheic zones increase freshwater ecosystem resilience to hydrological extremes and global environmental change. However, current conceptualizations of hyporheic exchange, residence time distributions, and the associated biogeochemical cycling in streambed sediments do not always accurately explain the hydrological and biogeochemical complexity observed in streams and rivers. Specifically, existing conceptual models insufficiently represent the coupled transport and reactivity along groundwater and surface water flow paths, the role of autochthonous organic matter in streambed biogeochemical functioning, and the feedbacks between surface-subsurface ecological processes, both within and across spatial and temporal scales. While simplified approaches to these issues are justifiable and necessary for transferability, the exclusion of important hyporheic processes from our conceptualizations can lead to erroneous conclusions and inadequate understanding and management of interconnected surface water and groundwater environments. This is particularly true at the landscape scale, where the organizational principles of spatio-temporal dynamics of hyporheic exchange flow (HEF) and biogeochemical processes remain largely uncharacterized. This article seeks to identify the most important drivers and controls of HEF and biogeochemical cycling based on a comprehensive synthesis of findings from a wide range of river systems. We use these observations to test current paradigms and conceptual models, discussing the interactions of local-to-regional hydrological, geomorphological, and ecological controls of hyporheic zone functioning. This improved conceptualization of the landscape organizational principles of drivers of HEF and biogeochemical processes from reach to catchment scales will inform future river research directions and watershed management strategies.

54 ENVIRONMENTAL SCIENCES↗

Tunable shear thickening, aging, and rejuvenation in suspensions of shape-memory-endowed liquid crystalline particles

The morphological features of particles, notably shape anisotropy, critically influence the rheological properties of dense suspensions, spanning both natural and engineered systems. This work explores the potential of using shape memory particles to dynamically regulate suspension fluid flow through controllable shape transformations. First, we synthesize shape-memory particles with programmable anisotropy from liquid crystal elastomers, such that the stiffness and shapes of the particles can be tuned by manipulating temperature. Our findings reveal that suspensions from such particles exhibit significant tunability in shear thickening behavior, transitioning from discontinuous shear thickening to a Newtonian-like response within a narrow temperature range of 60 ° C. This capability to modulate rheological responses in situ presents an approach for addressing processing challenges in many applications where control over flow behavior is paramount. Furthermore, we also show that suspensions composed of these anisotropic particles can undergo physical aging, and evolve into a glassy state. This state can be escaped upon activation of the shape memory effect. This reversibility underscores the potential for using such materials to engineer systems that can enter or come out of kinetic arrest by leveraging internal mechanical responses to external stimuli. The insights gained here not only broaden our understanding of the interplay between particle geometry and suspension dynamics but also pave the way for leveraging ensembles of stimuli-responsive objects to precisely control collective behaviors in many-body systems.

Science & Technology - Other Topics↗

Computationally efficient CFD prediction of bubbly flow using physics-guided deep learning

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale framework was proposed in this paper by applying a physics-guided data-driven approach. Instrumental to this framework, 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. In this work, 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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗