Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow rule”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

An improved method for the prediction of completely three-dimensional aerodynamic load distributions of configurations with leading edge vortex separation

The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.

Rubbert, p. E.↗
NASA NTRS1991

Flux-vector splitting algorithm for chain-rule conservation-law form

A flux-vector splitting algorithm with Newton-Raphson iteration was developed for the 'full compressible' Navier-Stokes equations cast in chain-rule conservation-law form. The algorithm is intended for problems with deforming spatial domains and for problems whose governing equations cannot be cast in strong conservation-law form. The usefulness of the algorithm for such problems was demonstrated by applying it to analyze the unsteady, two- and three-dimensional flows inside one combustion chamber of a Wankel engine under nonfiring conditions. Solutions were obtained to examine the algorithm in terms of conservation error, robustness, and ability to handle complex flows on time-dependent grid systems.

Shih, T. I.-P.↗
DOE OSTI2024

Simulation-Based Analysis of Feeder Operation with Different PV Smart Inverter Functions on an Actual Distribution System: Preprint

High penetration of photovoltaics (PV) in distribution feeders can cause problems, such as overvoltage, reverse power flow, and large net load changes. Traditional voltage regulation devices, such as capacitors and voltage regulators, can solve some of these problems but might have some delays. Today, smart inverters are gradually being used to provide voltage regulation and frequency support in distribution systems. Different smart inverter settings have been recommended in various rules and standards; however, the potential benefits and their impacts on distribution system operation are not well compared and studied. This paper presents a comparison of different smart inverter settings as applied to a distribution system. An actual feeder model from San Diego Gas & Electric Company is used to conduct the simulation. Additionally, a load disaggregation method is proposed to disaggregate the load and PV profile for each load location using advanced metering infrastructure net load measurements. Then, different smart inverter settings are applied to the PV systems in the feeder, and the simulation results are compared. The results show that the implementation of specific functions of smart inverters can reduce voltage exceedances, and the utility can determine the specific inverter setting based on its operational requirements.

distribution system↗
NASA NTRS1994

An initial investigation into methods of computing transonic aerodynamic sensitivity coefficients

The primary accomplishments of the project are as follows: (1) Using the transonic small perturbation equation as a flowfield model, the project demonstrated that the quasi-analytical method could be used to obtain aerodynamic sensitivity coefficients for airfoils at subsonic, transonic, and supersonic conditions for design variables such as Mach number, airfoil thickness, maximum camber, angle of attack, and location of maximum camber. It was established that the quasi-analytical approach was an accurate method for obtaining aerodynamic sensitivity derivatives for airfoils at transonic conditions and usually more efficient than the finite difference approach. (2) The usage of symbolic manipulation software to determine the appropriate expressions and computer coding associated with the quasi-analytical method for sensitivity derivatives was investigated. Using the three dimensional fully conservative full potential flowfield model, it was determined that symbolic manipulation along with a chain rule approach was extremely useful in developing a combined flowfield and quasi-analytical sensitivity derivative code capable of considering a large number of realistic design variables. (3) Using the three dimensional fully conservative full potential flowfield model, the quasi-analytical method was applied to swept wings (i.e. three dimensional) at transonic flow conditions. (4) The incremental iterative technique has been applied to the three dimensional transonic nonlinear small perturbation flowfield formulation, an equivalent plate deflection model, and the associated aerodynamic and structural discipline sensitivity equations; and coupled aeroelastic results for an aspect ratio three wing in transonic flow have been obtained.

Carlson, Leland A.↗
NASA NTRS1985

Continental volume and freeboard through geological time

The consequences of approximately constant freeboard for continental growth are explored using a model that relates the volumes of isostatically compensated continents and oceans to the secular decline in terrestrial heat flow. It is found that a post-Archean increase in freeboard by 200 m requires continental growth of only 10 percent, while a decrease in freeboard by 200 m during this same period necessitates a crustal growth of 40 percent. Shrinkage of the continental crust since the end of the Archean can be ruled out. Changes of more than 10 percent in post-Archean crustal thickness are highly unlikely.

Schubert, G.↗
NASA NTRS1957

On Slender-Body Theory at Transonic Speeds

The basic ideas of the slender-body approximation have been applied to the nonlinear transonic-flow equation for the velocity potential in order to obtain some of the essential features of slender-body theory at transonic speeds. The results of the investigation are presented from a unified point of view which demonstrates the similarity of slender-body solutions in the various Mach number ranges. The transonic area rule and some conditions concerning its validity follow from the analysis. (author)

Harder, Keith C↗
NASA NTRS1977

Liquid Rocket Lines, Bellows, Flexible Hoses, and Filters

Fluid-flow components in a liquid propellant rocket engine and the rocket vehicle which it propels are interconnected by lines, bellows, and flexible hoses. Elements involved in the successful design of these components are identified and current technologies pertaining to these elements are reviewed, assessed, and summarized to provide a technology base for a checklist of rules to be followed by project managers in guiding a design or assessing its adequacy. Recommended procedures for satisfying each of the design criteria are included.

Source record↗
NASA NTRS2011

Dynamic Simulation of Human Thermoregulation and Heat Transfer for Spaceflight Applications

Models of human thermoregulation and heat transfer date from the early 1970s and have been developed for applications ranging from evaluating thermal comfort in spacecraft and aircraft cabin environments to predicting heat stress during EVAs. Most lumped or compartment models represent the body as an assemblage cylindrical and spherical elements which may be subdivided into layers to describe tissue heterogeneity. Many existing models are of limited usefulness in asymmetric thermal environments, such as may be encountered during an EVA. Conventional whole-body clothing models also limit the ability to describe local surface thermal and evaporation effects in sufficient detail. A further limitation is that models based on a standard man model are not readily scalable to represent large or small subjects. This work describes development of a new human thermal model derived from the 41-node man model. Each segment is divided into four concentric, constant thickness cylinders made up of a central core surrounded by muscle, fat, and skin, respectively. These cylinders are connected by the flow of blood from a central blood pool to each part. The central blood pool is updated at each time step, based on a whole-body energy balance. Results show the model simulates core and surface temperature histories, sweat evaporation and metabolic rates which generally are consistent with controlled exposures of human subjects. Scaling rules are developed to enable simulation of small and large subjects (5th percentile and 95th percentile). Future refinements will include a clothing model that addresses local surface insulation and permeation effects and developing control equations to describe thermoregulatory effects such as may occur with prolonged weightlessness or with aging.

Miller, Thomas R.↗
DOE OSTI2024

ARPA-E Grid Optimization (GO) Competition Challenge 1

The ARPA-E Grid Optimization (GO) Competition Challenge 1, from 2018 to 2019, focused on the basic Security Constrained AC Optimal Power Flow problem (SCOPF) for a single time period. The Challenge utilized sets of unique datasets generated by the ARPA-E GRID DATA program. Each dataset consisted of a collection of power system network models of different sizes with associated operating scenarios (snapshots in time defining instantaneous power demand, renewable generation, generator and line availability, etc.). The datasets were of two types: Real-Time, which included starting-point information, and Online, which did not. Week-Ahead data is also provided for some cases but was not used in the Competition. Although most datasets were synthetic and generated by GRIDDATA, a few came from industry and were only used in the Final Event. All synthetic Input Data and Team Results for the GO Competition Challenge 1 for the Sandbox, Trial Events 1 to 3, and the Final Event along with problem, format, scoring and rules descriptions are available here. Data for industry scenarios will not be made public. Challenge 1, a minimization problem, required two computational steps. Solver 1 or Code 1 solved the base SCOPF problem under a strict wall clock time limit, as would be the case in industry, and reported the base case operating point as output, which was used to compute the Objective Function value that was used as the scenario score. The feasibility of the solution was provided by the Solver 2 or Code 2, which solves the power flow problem for all contingencies based on the results from Solver 1. This is not normally done in industry, so the time limits were relaxed. In fact, there were no time limits for Trial Event 1. This proved to be a mistake, with some codes running for more than 90 hours, and a time limit of 2 seconds per contingency was imposed for all other events. Entrants were free to use their own Solver 2 or use an open-source version provided by the Competition. Containers, such as Docker, were considered to improve the portability of codes, but none that could reliably support a multi-node parallel computing environment, e.g., MPI, could be found. For more information on the competition and challenge see the "GO Competition Challenge 1 Information" and "GO Competition Challenge 1 Additional Information" resources below.

ACOPF↗
DOE OSTI2024

Identifying Hydropower Operational Flexibilities in Presence of Streamflow and Net-load Uncertainty (Final Technical report)

In the existing operations, hydropower contributions to future system flexibility are generally modeled while maintaining traditional operating rules and constraints in supporting grid operation, such as the balancing of variable renewable energy production. Moreover, operation of large scale hydropower systems on major rivers has been investigated for decades, utilizing various systems engineering approaches, with the evolving electric grid, as the result of renewable resources integration, compounded by the changing climate (variability of river flows, intensification of hydrologic cycle resulting in more frequent extreme events) affecting water availability, the need for more advanced stochastic modeling and effective uncertainty analysis approaches have become necessary. The research results supported by this funding and presented in this report provide a new look at hydropower operational flexibility enforced by the changes identified above. Understanding how hydropower operates in response to the underlying uncertainties with respect to the system constraints is crucial in identifying its operational flexibility potentials. In this project, the flexibility of the operating hydropower facility is described by capturing uncertainties in both water and power system and formulating the operations as a multistage stochastic optimization problem. The proposed approach supports short- to seasonal-term operations and planning decision horizons.

13 HYDRO ENERGY↗
NASA NTRS1995

A test of a vortex method for the computation of flap side edge noise

Upon approach to landing, a major source location of airframe noise occurs at the side edges of the part span, trailing edge flaps. In the vicinity of these flaps, a complex arrangement of spanwise flow with primary and secondary tip vortices may form. Each of these vortices is observed to become fully three-dimensional. In the present study, a numerical model is developed to investigate the noise radiated from the side edge of a flap. The inherent three-dimensionality of this flow forces us to carefully consider a numerical scheme which will be both accurate in its prediction of the flow acoustics and also computationally efficient. Vortex methods have offered a fast and efficient means of simulating many two and three-dimensional, vortex dominated flows. In vortex methods, the time development of the flow is tracked by following exclusively the vorticity containing regions. Through the Biot-Savart law, knowledge of the vorticity field enables one to obtain flow quantities at any desired location during the flow evolution. In the present study, a numerical procedure has been developed which incorporates the Lagrangian approach of vortex methods into a calculation for the noise radiated by a flow-surface interaction. In particular, the noise generated by a vortex in the presence of a flat half plane is considered. This problem serves as a basic model of flap edge flow. It also permits the direct comparison between our computed results and previous acoustic analyses performed for this problem. In our numerical simulations, the mean flow is represented by the complex potential W(z) = Aiz(exp l/2), which is obtained through conformal mapping techniques. The magnitude of the mean flow is controlled by the parameter A. This mean flow has been used in the acoustic analysis by Hardin and is considered a reasonable model of the flow field in the vicinity of the edge and away from the leading and trailing edges of the flap. To represent the primary vortex which occurs near the flap, a point vortex is introduced just below the flat half plane. Using a technique from panel methods, boundary conditions on the flap surface are satisfied by the introduction of a row of stationary point vortices along the extent of the flap. At each time step in the calculation, the strength of these vortices is chosen to eliminate the normal velocity at intermediary collocation points. The time development of the overall flow field is then tracked using standard techniques from vortex methods. Vortex trajectories obtained through this computation are in good agreement with those predicted by the analytical solution given by Hardin, thus verifying the viability of this procedure for more complex flow arrangements. For the flow acoustics, the Ffowcs Williams-Hawkings equation is numerically integrated. This equation supplies the far field acoustic pressure based upon pressures occurring along the flap surface. With our vortex method solution, surface pressures may be obtained with exceptional resolution. The Ffowcs Williams-Hawkings equation is integrated using a spatially fourth order accurate Simpson's rule. Rational function interpolation is used to obtain the surface pressures at the appropriate retarded times. Comparisons between our numerical results for the acoustic pressure and those predicted by the Hardin analysis have been made. Preliminary results indicate the need for an improved integration technique. In the future, the numerical procedure developed in this study will be applied to the case of a rectangular flap of finite thickness and ultimately modified for application to the fully three-dimensional problem.

Martin, James E.↗
DOE OSTI2021

Learning viscoelasticity models from indirect data using deep neural networks

In this study, we propose a novel approach to model viscoelasticity materials, where rate-dependent and non-linear constitutive relationships are approximated with deep neural networks. We assume that inputs and outputs of the neural networks are not directly observable, and therefore common training techniques with input–output pairs for the neural networks are inapplicable. To that end, we develop a novel computational approach to both calibrate parametric and learn neural-network-based constitutive relations of viscoelasticity materials from indirect displacement data in the context of multiple-physics systems. We show that limited displacement data holds sufficient information to quantify the viscoelasticity behavior. We formulate the inverse computation – modeling viscoelasticity properties from observed displacement data – as a PDE-constrained optimization problem and minimize the error functional using a gradient-based optimization method. The gradients are computed by a combination of automatic differentiation and implicit function differentiation rules. The effectiveness of our method is demonstrated through numerous benchmark problems in geomechanics and porous media transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Correlations in the Initial Conditions of Heavy-Ion Collisions

Ultracentral collisions of heavy nuclei, in which the impact parameter is nearly zero, are especially sensitive to the details of the initial state model and the microscopic mechanism for collective flow. In a hydrodynamic “flow” picture, the final state momentum correlations are a direct response to the fluctuating initial geometry, although models of the initial geometry differ widely. Alternatively, dynamical mechanisms based in the color glass condensate (CGC) formalism can naturally lead to many-body correlations with very different systematics. Here we present a calculation of event-by-event elliptic flow in both the hydrodynamic and CGC paradigms and show that they can be qualitatively distinguished in ultracentral collisions of deformed nuclei. Specifically, the multiplicity dependence in such collisions is qualitatively opposite, with the CGC correlations increasing with multiplicity while the hydrodynamic correlations decrease. The consistency of the latter with experimental data on UU collisions appears to rule out a CGC-mediated explanation. We find that these qualitative features also persist in small deformed systems and can therefore be a valuable test of the microscopic physics in that regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

Growth of a diffusion flame in the field of a vortex

In the present study of the growth of a diffusion flame in the field of a vortex, the motion in the core is converted into a solid body rotation. The flame extension and distortion kinematics are presented, and the effect of the local flow field on local flame structure is analyzed in detail. The combustion field is found to consist of a totally reacted core region whose radius is time-dependent, and an external flame region which consists of a pair of spiral arms that extend at large radii toward their original positions on the horizontal axis. Two similarity rules are formulated which are independent of kinematic viscosity.

Marble, F. E.↗
NASA NTRS2002

CATS-based Air Traffic Controller Agents

This report describes intelligent agents that function as air traffic controllers. Each agent controls traffic in a single sector in real time; agents controlling traffic in adjoining sectors can coordinate to manage an arrival flow across a given meter fix. The purpose of this research is threefold. First, it seeks to study the design of agents for controlling complex systems. In particular, it investigates agent planning and reactive control functionality in a dynamic environment in which a variety perceptual and decision making skills play a central role. It examines how heuristic rules can be applied to model planning and decision making skills, rather than attempting to apply optimization methods. Thus, the research attempts to develop intelligent agents that provide an approximation of human air traffic controller behavior that, while not based on an explicit cognitive model, does produce task performance consistent with the way human air traffic controllers operate. Second, this research sought to extend previous research on using the Crew Activity Tracking System (CATS) as the basis for intelligent agents. The agents use a high-level model of air traffic controller activities to structure the control task. To execute an activity in the CATS model, according to the current task context, the agents reference a 'skill library' and 'control rules' that in turn execute the pattern recognition, planning, and decision-making required to perform the activity. Applying the skills enables the agents to modify their representation of the current control situation (i.e., the 'flick' or 'picture'). The updated representation supports the next activity in a cycle of action that, taken as a whole, simulates air traffic controller behavior. A third, practical motivation for this research is to use intelligent agents to support evaluation of new air traffic control (ATC) methods to support new Air Traffic Management (ATM) concepts. Current approaches that use large, human-in-the-loop simulations are unquestionably valuable for this purpose, but pose considerable logistical, fiscal, and experimental control problems. First, data analysis is extremely complicated, owing simply to the large number of participants and data sources in such simulations. In addition, experienced human air traffic controllers working adjacent sectors tend to flexibly adapt to the evolving control problem - potentially shifting to other strategies than those under investigation. In addition, their performance is tightly coupled to the control interface, which in the development phase may support some concepts and supporting strategies better than others. A simple shift in strategy by one controller can change the character of a particular traffic scenario dramatically, which makes experimental comparison of ATC performance under different traffic scenarios difficult. Training a given team of controllers on operations under a new ATM concept for a sufficient period of time could avert such difficulties, but instituting an adequate training program is expensive and logistically difficult.

Callantine, Todd J.↗
NASA NTRS1973

Relaxation solutions for inviscid axisymmetric transonic flow over blunt or pointed bodies.

A finite-difference relaxation method is presented for numerical solution of the full potential equation and exact boundary conditions for general axisymmetric bodies is inviscid, steady transonic flow. Body-normal coordinates are used in the nose region and sheared cylindrical coordinates are used on the afterbody to accommodate corners such as boattails and flares. An improved difference scheme is used which does not require that the flow be nearly alined with a coordinate direction in supersonic regions, and which treats either subsonic or supersonic free streams. Numerical results are illustrated for some simple classical shapes such as spheres and ellipsoids, and for more practical shapes like tangent-ogives with boattails. Special attention is given to bodies which have been studied for area-rule applications. Agreement with available experimental results is good in cases where viscous effects and wind-tunnel wall interference are not important.

South, J. C., Jr.↗
NASA NTRS2020

Inter-species molecular attraction effect in the development of a two-species mixing layer

Simulations of fluid mixing occurring at high-pressures mandate the use of real-fluid equations of state to capture the effect of boh intermolecular repulsive forces and attractive forces which are neglected in perfect gases. Dedicated mixing rules account for these forces by providing a generic mapping of the mixture state-space onto a pure fluid domain. However, the specific inter-species molecular forces depend on the particular characteristics of the molecular pairs under consideration, and typically require the knowledge of binary interaction coefficients k_{αβ}. The impact of k_{αβ} on a flow field is as of now unclear; further, k_{αβ} is unknown for many molecular pairs and when unknown, set to be null. This study addresses the impact of k_{αβ} on the temporal evolution of a real-fluid mixing layer and uses Direct Numerical Simulation (DNS) to explore this impact. The results show differences between augmented attraction (k_{αβ}<0) and diminished attraction (k_{αβ}>0), affecting both fluid dynamics and diffusion processes. Specifically, in a binary mixing layer, it is observed that augmented attraction leads to delayed transition and mixing layer growth, manifested in the momentum layer thickness. The fractal dimension, which is a manifestation of interface corrugation, is calculated. These results provide an hitherto undocumented mechanism affecting the development of mixing layers.

Banuti, Daniel↗
NASA NTRS1973

Theoretical prediction of airplane stability derivatives at subcritical speeds

The theoretical development and application is described of an analysis for predicting the major static and rotary stability derivatives for a complete airplane. The analysis utilizes potential flow theory to compute the surface flow fields and pressures on any configuration that can be synthesized from arbitrary lifting bodies and nonplanar thick lifting panels. The pressures are integrated to obtain section and total configuration loads and moments due side slip, angle of attack, pitching motion, rolling motion, yawing motion, and control surface deflection. Subcritical compressibility is accounted for by means of the Gothert similarity rule.

Tulinius, J.↗