Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BERNOULLI EQUATION”

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 37 records · Page 2

The Beginner's Guide to Wind Tunnels with TunnelSim and TunnelSys

The Beginner's Guide to Wind Tunnels is a Web-based, on-line textbook that explains and demonstrates the history, physics, and mathematics involved with wind tunnels and wind tunnel testing. The Web site contains several interactive computer programs to demonstrate scientific principles. TunnelSim is an interactive, educational computer program that demonstrates basic wind tunnel design and operation. TunnelSim is a Java (Sun Microsystems Inc.) applet that solves the continuity and Bernoulli equations to determine the velocity and pressure throughout a tunnel design. TunnelSys is a group of Java applications that mimic wind tunnel testing techniques. Using TunnelSys, a team of students designs, tests, and post-processes the data for a virtual, low speed, and aircraft wing.

Benson, Thomas J.↗

Numerical Modeling of Cavitating Venturi: A Flow Control Element of Propulsion System

In a propulsion system, the propellant flow and mixture ratio could be controlled either by variable area flow control valves or by passive flow control elements such as cavitating venturies. Cavitating venturies maintain constant propellant flowrate for fixed inlet conditions (pressure and temperature) and wide range of outlet pressures, thereby maintain constant, engine thrust and mixture ratio. The flowrate through the venturi reaches a constant value and becomes independent of outlet pressure when the pressure at throat becomes equal to vapor pressure. In order to develop a numerical model of propulsion system, it is necessary to model cavitating venturies in propellant feed systems. This paper presents a finite volume model of flow network of a cavitating venturi. The venturi was discretized into a number of control volumes and mass, momentum and energy conservation equations in each control volume are simultaneously solved to calculate one-dimensional pressure, density, and flowrate and temperature distribution. The numerical model predicts cavitations at the throat when outlet pressure was gradually reduced. Once cavitation starts, with further reduction of downstream pressure, no change in flowrate is found. The numerical predictions have been compared with test data and empirical equation based on Bernoulli's equation.

Majumdar, Alok↗

Variations in the modal characteristics of a telescopically deploying beam

The equations of motion for a two-segment deploying telescopic beam are derived through application of Lagrange's equation. The outer tube of the beam is fixed at one end and the inner tube slides freely relative to the fixed segment. The resulting nonlinear, non-autonomous set of equations is linearized and simplified to the standard Euler-Bernoulli partial differential equations for an elastic beam by freezing the deployment process at various stages of deployment, and examining the small amplitude and natural modes of vibration of the resulting configuration. Application of the natural boundary conditions and compatibility of motion relations for the two segments in their common region of overlap leads to a transcendental characteristic equation in the frequency parameter Beta(L). Numerical solution of the equation for the characteristic roots determines the modal frequencies, and the corresponding mode shapes are obtained from the general solution of the Euler-Bernoulli equation tailored to the natural boundary conditions. Sample results of modal frequencies and shapes are presented for various stages of deployment and discussed. It is shown that for all intermediate stages of deployment (between 0 and 100 percent) the spectral distribution is drastically altered by the appearance of regions of very closely spaced modal frequencies. The sources of this modal agglomeration are explored.

Amos, Anthony K.↗

Calculation of the Pressure Distribution on Bodies of Revolution in the Subsonic Flow of a Gas: Axially Symmetrical Flow - Part 1

The present report concerns a method of computing the velocity and pressure distributions on bodies of revolution in axially symmetrical flow in the subsonic range. The differential equation for the velocity potential Phi of a compressible fluid motion is linearized tn the conventional manner, and then put in the form Delta(Phi) = 0 by affine transformation. The quantity Phi represents the velocity potential of a fictitious incompressible flow, for which a constant superposition of sources by sections is secured by a method patterned after von Karman which must comply with the boundary condition delta(phi)/delta(n) = 0 at the originally specified contour. This requirement yields for the "pseudo-stream function" psi a differential equation which must be fulfilled for as many points on the contour as source lengths are assumed. In this manner, the problem of defining the still unknown source intensities is reduced to the solution of an inhomogeneous equation system. The pressure distribution is then determined with the aid of Bernoulli's equation and adiabatic equation of state. Lastly, the pressure distributions in compressible and incompressible medium are compared on a model problem.

Bilharz, Herbert↗

Relativistic hadrons and the origin of relativistic outflows in active galactic nuclei

We examine the hydrodynamic origin of relativistic outflows in active galactic nuclei (AGN). Specifically, we propose that the presence of a population of relativistic hadrons in the AGN 'central engine' and the associated neutron production suffices to produce outflows which under rather general conditions could be relativistic. The main such condition is that the size of the neutron production region be larger than the neutron flight path tau(sub n) approximately 3 x 10(exp 13) cm. This condition guarantees that the mean energy per particle in the proton fluid, resulting from the decay of the neutrons outside their production region, be greater than the proton rest mass. The expansion of this fluid can then lead naturally to a relativistic outflow by conversion of its internal energy to directed motion. We follow the development of such flows by solving the mass, energy as well as the kinetic equation for the proton gas in steady state, taking into account the source terms due to compute accurately the adiabatic index of the expanding gas, and in conjunction with Bernoulli's equation the detailed evolution of the bulk Lorentz factor. We further examine the role of large-scale magnetic fields in confining these outflows to produce the jets observed at larger scales.

Contopoulos, John↗

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗

Multi-fidelity modeling to predict the rheological properties of fiber suspensions

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications like biofuel production. The 3D numerical simulations of the suspension of fibers are often computationally expensive and time-consuming. Machine learning methods such as neural networks can simplify the prediction of rheological behavior; however, they require a relatively large training data set. Multi -fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified physical equations, can lead to optimized predictions. Here, we focus on a neural network with two levels of fidelity, i.e., high and low fidelity networks. To produce high-fidelity data, we perform direct numerical simulations to model the fibers as one-dimensional inextensible slender bodies that obey the Euler- Bernoulli beam equation. The Navier-Stokes equations govern the suspended fluid, and an immersed boundary method is used to couple the fluid and solid motion. The low-fidelity data is produced by using constitutive equations. Noticeable improvements have been observed in the accuracy of predicting the rheological behavior when a multi-fidelity network is used compared to the single-fidelity network.

Boodaghidizaji, Miad↗

Development of MCAERO wing design panel method with interactive graphics module

A reliable and efficient iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical pressure distribution. The design process is initialized by using MCAERO (MCAIR 3-D Subsonic Potential Flow Analysis Code) to analyze a baseline configuration. A second program DMCAERO is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter by applying a first-order expansion to the baseline equations in MCAERO. This matrix is calculated only once but is used in each iteration cycle to calculate the geometry perturbation and to analyze the perturbed geometry. The potential on the new geometry is calculated by linear extrapolation from the baseline solution. This extrapolated potential is converted to velocity by numerical differentiation, and velocity is converted to pressure by using Bernoulli's equation. There is an interactive graphics option which allows the user to graphically display the results of the design process and to interactively change either the geometry or the prescribed pressure distribution.

Hawk, J. D.↗

Subsonic panel method for designing wing surfaces from pressure distribution

An iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical distribution of pressure. The calculations are initialized by using a surface panel method to analyze a baseline wing or wing-fuselage configuration. A first-order expansion to the baseline panel method equations is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter. In every iteration cycle, the matrix is used both to calculate the geometry perturbation and to analyze the perturbed geometry. The distribution of potential on the perturbed geometry is established by simple linear extrapolation from the baseline solution. The extrapolated potential is converted to pressure by Bernoulli's equation. Not only is the accuracy of the approach good for very large perturbations, but the computing cost of each complete iteration cycle is substantially less than one analysis solution by a conventional panel method.

Bristow, D. R.↗

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

Fluidmagnetic buoyancy.

Fluidmagnetic buoyancy effects in ferrohydrodynamic /FHD/ flows and implications of Ernshaw levitation theorem

FLUID MECHANICS↗

Physics associated with cavitation induced material damage

The role of cavitation in mechanical failure is discussed. Some of the most common types of material damage associated with the presence of cavitation are surface material removal, delamination and structural vibration. This occurs in external flows such as on propellers, hydrofoils, and high speed non-lifting surfaces. In internal flows, pipe bends, inlets, constructions, pumps and turbines are typical. Nominally nonflowing liquids are also susceptible in, for example, strong acoustic fields and high energy particle detectors. For flowing systems, Bernoulli's equation shows how a local pressure is reduced as the fluid's velocity is increased. At sufficiently high velocities, a tension can actually develop and this has, in fact, been demonstrated experimentally. Once the pressure is reduced below the fluid vapor pressure a vapor cavity can be nucleated. Various aspects of this process are simply shown by considering the flow over a lifting surface.

Peterson, F. B.↗

A simplified theory of oscillating airfoils in transonic flow - Review and extension

Significant old and new results are presented to show to what extent a simplified theory for transonic flow may be used. Solutions are obtained by classical techniques and compared with experiment. Results are given for two-dimensional, steady and unsteady flow and three-dimensional, steady flow. The effects of flow separation and improvements in Bernoulli's equations and the surface boundary condition are also briefly discussed.

Dowell, E. H.↗

Implementation of nonreflective boundary condition at the outflow boundary

Numerical techniques for the solution of free surface flows about a disturbance on or near a free surface were developed. Any disturbance moving on or near a liquid surface generates surface waves that propagate away from the disturbance. The potential flow problems in which the evolution of free surface waves is the principal phenomenon of interest are considered. The distinct feature of these problems is that the position of the surface is also unknown and is to be determined as part of the solution. Two boundary conditions are required at the free surface since there are two unknown functions there, the free surface height and the potential function. The required boundary conditions are satisfied by the kinematic condition which states that no flow can cross the surface and the dynamic condition which matches the fluid pressure to the pressure acting on the surface through Bernoulli's equation.

Hall, D. R.↗

A simplified theory of oscillating aerofoils in transonic flow - Review and extension

A simplified theory of the dynamic motion of aerofoils of finite thickness in transonic flow is presented which excludes the effect of shock waves on the aerofoil itself and, thus, is restricted to free stream Mach numbers equal to unity or above. Numerical examples are analyzed for two-dimensional steady and unsteady (including transient) aerofoil motion and three-dimensional steady and unsteady flow over delta wings. The effects of flow separation and improvements in Bernoulli's equation and the surface boundary condition are briefly discussed.

Dowell, E. H.↗

A design procedure for active control of beam vibrations

The transverse vibrations of beams is discussed and a methodology for the design of an active damping device is given. The Bernoulli-Euler equation is used to derive a transcendental transfer function, which relates a torque applied at one end of the beam to the rotational position and velocity at that point. The active damping device consists of a wire, a linear actuator and a short torque arm attached to one end of the beam. The action of the actuator varies a tension in the wire and creates a torque which opposes the rotation of the beam and thus damps vibration. A design procedure for such an active damper is given. This procedure shows the relationships and trade-offs between the actuator stroke, power required, stress levels in the wire and beam and the geometry of the beam and wire. It is shown that by consideration of the frequency response at the beam natural frequencies, the aforementioned relationships can be greatly simplified. Similarly, a simple way of estimating the effective damping ratios and eigenvalue locations of actively controlled beams is presented.

Dickerson, S. L.↗

Parameter estimation for the Euler-Bernoulli-beam

An approximation involving cubic spline functions for parameter estimation problems in the Euler-Bernoulli-beam equation (phrased as an optimization problem with respect to the parameters) is described and convergence is proved. The resulting algorithm was implemented and several of the test examples are documented. It is observed that the use of penalty terms in the cost functional can improve the rate of convergence.

Graif, E.↗