Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

Expansion Valve With Temperature-Sensitive Flow Regulation

Joule-Thomson expansion valve designed so particles or condensed contaminants in vapor flowing through unlikely to clog orifice. In addition, new valve automatically adjusts flow of vapor. Allows high initial flow for rapid cooling, but gradually reduces flow as operating temperature reached. Flow adjusts manually when necessary. Developed for expansion of high-pressure vapors to condense them into low-pressure cryogenic liquids. Used to conserve helium coolant in superconducting system by condensing vapor boiling from liquid-helium bath.

Walker, Graham↗

Modeling of hypersonic reacting flows

The importance of thermochemical nonequilibrium in hypersonic flow is studied as well as the breakdown of a one-temperature kinetic model. Consideration is given to vibration relaxation, reaction rates, and conservation equations for chemical variables. Experimental data on the pitching moment characteristics of the raked-off blunted elliptic cone geometry proposed for use with an aeroassisted orbit transfer vehicle are presented.

Park, Chul↗

Particle migration in suspensions by thermocapillary or electrophoretic motion

Two problems of similar mathematical structure are presented: the thermocapillary motion of bubbles and the electrophoresis of colloidal particles. For both problems, it is shown that in a cloud of n particles surrounded by an infinite expanse of fluid, the velocity of each sphere under creeping flow conditions is equal to the velocity of an isolated particle, unchanged by interactions between the particles. However, when the cloud fills a container, conservation of mass shows that this result cannot continue to hold, and the average translational velocity must be calculated subject to a constraint on the mass flux. It is concluded that the average thermocapillary or electrophoretic translational velocity of a particle in the cloud is related to the effective conductivity of the cloud over the whole range of particle volume fractions, provided that the particles are identical, non-conducting and, for the thermocapillary problem, inviscid.

Acrivos, A.↗

Bulk properties and velocity distributions of water group ions at Comet Halley - Giotto measurements

In the region upstream of Comet Halley, pickup heavy ions of cometary origin were directly observed by the implanted ion spectrometer on Giotto. Diffusion of this population in pitch angle and in energy, during the approach to the comet and on the outbound leg is discussed. The two data sets are compared and qualitative ideas on scattering timescales are inferred. In addition the bulk parameters of these distributions have been computed and a comparison of the observed speed in the solar wind frame and the observed density with expectations is presented. Pitch angle scattering occurs more slowly than expected with filled shells appearing at 2,500,000 km, and significant energy diffusion does not occur until the bow shock region. Also the shell distributions downstream of the shock flow at the bispherical bulk speed (related to the Alfven speed) along the magnetic field with respect to the solar wind in accordance with conservation of energy between the pickup ions and the wave turbulence.

Coates, A. J.↗

Application of a Reynolds Stress turbulence model to a supersonic hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows undergoing hydrogen-air chemical reactions. An assumed Beta Probability Density Function is applied to account for the chemical source terms in the conservation equations. An algebraic Reynolds Flux model is used for the fluctuating density-velocity as well as the species mass fraction-velocity correlations. The variances of temperature and species fluctuations are also modelled using an algebraic flux technique. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion processes. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant.

Chandrasekhar, R.↗

Modelling the Centers of Galaxies

The key to studying central regions by means of nobody numerical experiments is to concentrate on the central few parsecs of a galaxy, replacing the remainder of the galaxy by a suitable boundary condition, rather after the manner in which stellar interiors can be studied without a detailed stellar atmosphere by replacing the atmosphere with a boundary condition. Replacements must be carefully designed because the long range gravitational force means that the core region is sensitive to mass outside that region and because particles can exchange between the outer galaxy and the core region. We use periodic boundary conditions, coupled with an iterative procedure to generate initial particle loads in isothermal equilibrium. Angular momentum conservation is ensured for problems including systematic rotation by a circular reflecting boundary and by integrating in a frame that rotates with the mean flow. Mass beyond the boundary contributes to the gravitational potential, but does not participate in the dynamics. A symplectic integration scheme has been developed for rotating coordinate systems. This combination works well, leading to robust configurations. Some preliminary results with this combination show that: (1) Rotating systems are extremely sensitive to non-axisymmetric external potentials, and (2) that a second core, orbiting near the main core (like the M31 second core system), shows extremely rapid orbital decay. The experimental setups will be discussed, along with preliminary results.

Smith, B. F.↗

The Collapse of Vapor Bubbles in a Spatially Non-Uniform Flow

Pressure gradients act differently on liquid particles and suspended bubbles and are, therefore, capable of inducing a relative motion between the phases even when no relative velocity initially exists. As a consequence of the enhanced heat transfer in the presence of convection, this fact may have a major impact on the evolution of a vapor bubble. The effect is particularly strong in the case of a collapsing bubble for which, due to the conservation of the system's impulse, the induced relative velocity tends to be magnified when the bubble volume shrinks. A practical application could be, for instance, the enhancement of the condensation rate of bubbles downstream of a heated region, thereby reducing the quality of a flowing liquid-vapor mixture. A simple model of the process, in which the bubble is assumed to be spherical and the flow potential, is developed in the paper.

Hao, Y.↗

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flow and pressure regulation in the cardiovascular system

Principles and descriptive fragments which may contribute to a model of the regulating chains in the cardiovascular system are presented. Attention is given to the strain sensitivity of blood vessels, the law of the autonomy of the heart beat oscillator, the law of the encapsulation of body fluids, the law of the conservation of protein, the law of minimum 'arterial' pressure, the design of the 'mammalian' kidney, questions of homeokinetic organization, and the development of self-regulatory chains. Details concerning the development program for the heart muscle are considered along with the speed of response of the breathing rate and the significance of the pulmonary vascular pressure-flow characteristics.

Iberall, A.↗

Design of supercritical swept wings

Computational fluid dynamics are used to discuss problems inherent to transonic three-dimensional flow past supercritical swept wings. The formulation for a boundary value problem for the flow past the wing is provided, including consideration of weak shock waves and the use of parabolic coordinates. A swept wing code is developed which requires a mesh of 152 x 10 x 12 points and 200 time cycles. A formula for wave drag is calculated, based on the idea that the conservation form of the momentum equation becomes an entropy inequality measuring the drag, expressible in terms of a small-disturbance equation for a potential function in two dimensions. The entropy inequality has been incorporated in a two-dimensional code for the analysis of transonic flow over airfoils. A method of artificial viscosity is explored for optimum pressure distributions with design, and involves a free boundary problem considering speed over only a portion of the wing.

Garabedian, P.↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transonic flow calculations using a flux vector splitting method for the Euler equations

A study of the flux vector splitting method of Steger and Warming for the solution of the time dependent Euler equations in strong conservation law form for arbitrary two-dimensional geometries is presented. The procedure employed here differs from that of Buning and Steger in that it uses a different algorithm and employs implicit boundary conditions. Moreover, the method, as implemented here, does not contain any explicit smoothing or any adjustable parameters. Calculations were carried out for an NACA 0012 airfoil at various Mach numbers and angles of attack, and cylinders. Steady symmetric solutions were obtained for the full cylinder at a freestream Mach number of .5 without imposing a symmetry condition. In general, good agreement with other methods was obtained.

Seaford, C. M.↗

Numerical solution of unsteady incompressible viscous flows in generalized moving coordinate systems

A solution method of the time-accurate, incompressible Navier-Stokes equations in generalized curvilinear moving coordinate systems is presented in this paper. Accuracy is achieved by a conservative finite-volume discretization which satisfies the geometric conservation laws in generalized moving coordinate systems. The solution method is second-order accurate in space and first-order accurate in time. A fractional step solution method is used to efficiently solve the discrete equations. The unknowns, namely the pressure and the volume-fluxes, are chosen to facilitate the formulation of a consistent Poisson equation and to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparisons to other numerical and experimental solutions. The comparisons show good agreement.

Rosenfeld, Moshe↗

Time-dependent solutions of viscous incompressible flows in moving co-ordinates

A time-accurate solution method for the incompressible Navier-Stokes equations in generalized moving coordinates is presented. A finite-volume discretization method that satisfies the geometric conservation laws for time-varying computational cells is used. The discrete equations are solved by a fractional-step solution procedure. The solution is second-order-accurate in space and first-order-accurate in time. The pressure and the volume fluxes are chosen as the unknowns to facilitate the formulation of a consistent Poisson equation and thus to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparing the solutions with other numerical and experimental results. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗

Advanced Small Perturbation Potential Flow Theory for Unsteady Aerodynamic and Aeroelastic Analyses

An advanced small perturbation (ASP) potential flow theory has been developed to improve upon the classical transonic small perturbation (TSP) theories that have been used in various computer codes. These computer codes are typically used for unsteady aerodynamic and aeroelastic analyses in the nonlinear transonic flight regime. The codes exploit the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP theory was developed methodically by first determining the essential elements required to produce full-potential-like solutions with a small perturbation approach on the requisite Cartesian grid. This level of accuracy required a higher-order streamwise mass flux and a mass conserving surface boundary condition. The ASP theory was further developed by determining the essential elements required to produce results that agreed well with Euler solutions. This level of accuracy required mass conserving entropy and vorticity effects, and second-order terms in the trailing wake boundary condition. Finally, an integral boundary layer procedure, applicable to both attached and shock-induced separated flows, was incorporated for viscous effects. The resulting ASP potential flow theory, including entropy, vorticity, and viscous effects, is shown to be mathematically more appropriate and computationally more accurate than the classical TSP theories. The formulaic details of the ASP theory are described fully and the improvements are demonstrated through careful comparisons with accepted alternative results and experimental data. The new theory has been used as the basis for a new computer code called ASP3D (Advanced Small Perturbation - 3D), which also is briefly described with representative results.

Batina, John T.↗

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗