Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power flow 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 307 records · Page 17

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical modeling of an arcjet thruster

A numerical model of an arcjet thruster was developed and the results of calculations using this model are compared to the experimental data from a 30 kW arcjet using ammonia as a propellant. The model contains equations for the conservation of mass, radial, axial, and azimuthal momentum, energy, and the radial component of the magnetic field and provides a solution for three components of velocity, temperature, and the axial and radial current as a function of position in an axisymmetric flowfield. The model predicts 12 percent more specific impulse than the experiment when the mass flow and power input are constrained to the experimental values. A major conclusion of this study is that improved procedures for calculating transport properties are necessary if the accuracy of the model is to be improved.

Rhodes, Robert P.↗

Users' Manual for Computer Code SPIRALI Incompressible, Turbulent Spiral Grooved Cylindrical and Face Seals

The SPIRALI code predicts the performance characteristics of incompressible cylindrical and face seals with or without the inclusion of spiral grooves. Performance characteristics include load capacity (for face seals), leakage flow, power requirements and dynamic characteristics in the form of stiffness, damping and apparent mass coefficients in 4 degrees of freedom for cylindrical seals and 3 degrees of freedom for face seals. These performance characteristics are computed as functions of seal and groove geometry, load or film thickness, running and disturbance speeds, fluid viscosity, and boundary pressures. A derivation of the equations governing the performance of turbulent, incompressible, spiral groove cylindrical and face seals along with a description of their solution is given. The computer codes are described, including an input description, sample cases, and comparisons with results of other codes.

Walowit, Jed A.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models have been implemented into the FUN3D code to evaluate their ability to improve the computation of challenging flows encountered in aerospace propulsion, including mixing flows. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the mixing length. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. characteristics, no significant changes in the downstream flow behavior compared to the baseline case are observed. Furthermore, the total power consumed by the fans for different incoming flow conditions also remain marginally the same. It is hoped that the results, albeit obtained at very low speeds. would serve as a database for this technologically interesting flow field that has not been explored adequately before.

computational fluid dynamics↗

Developing VSC-HVDC Oscillation Damping Control Constraints in Unit Commitment

High-voltage direct current (HVDC) systems within the existing power grids will play a pivotal role in enabling high share of renewable power integration and transportation electrification. There is a need for transmission planning models to accurately capture the impact of VSC-HVDC oscillation damping control on the system planning and model the constraints in the security constrained unit commitment (SCUC) problem. In this paper, the SCUC problem in a hybrid DC-AC grid considering damping on electromechanical inter-area oscillation is studied. An electromechanical oscillation-driven transmission capability constraints are derived from system swing equation and incorporated into the SCUC. Multiple system damping control scenarios are set based on the extracted plausible bounds of damping coefficient and inertia constant to linearize the nonlinear expressions. The SCUC problem with proposed constraints is tested on the IEEE 73-bus reliability test system (RTS) case via PLEXOS. The impact of the proposed constraints on regional generation, inertia, power flow, annual system cost, and renewable energy curtailment are fully assessed.

damping control↗

The free-wake prediction of rotor hover performance using a vortex embedding method

A method is developed to predict the rotor hover performance. This method solves the compressible mass conservation equation much like current full potential codes and can therefore predict the transonic flows on a rotor. However, the newly developed approach also allows for the free convection of shed vorticity and permits the computation of the entire hover wake system. The method uses a vortex embedding scheme in potential flow and has been implemented in a computer code, HELIX -I. To predict power we implement a simple boundary layer and two different induced-drag integration schemes. The induced-drag is obtained from surface pressure integration and an energy flux integral. Comparisons between computations and experiment show good agreement for the prediction of power polars, surface pressure distribution, and tip vortex geometry.

Ramachandran, K.↗

De-Aliasing Through Over-Integration Applied to the Flux Reconstruction and Discontinuous Galerkin Methods

High-order methods are quickly becoming popular for turbulent flows as the amount of computer processing power increases. The flux reconstruction (FR) method presents a unifying framework for a wide class of high-order methods including discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV). It offers a simple, efficient, and easy way to implement nodal-based methods that are derived via the differential form of the governing equations. Whereas high-order methods have enjoyed recent success, they have been known to introduce numerical instabilities due to polynomial aliasing when applied to under-resolved nonlinear problems. Aliasing errors have been extensively studied in reference to DG methods; however, their study regarding FR methods has mostly been limited to the selection of the nodal points used within each cell. Here, we extend some of the de-aliasing techniques used for DG methods, primarily over-integration, to the FR framework. Our results show that over-integration does remove aliasing errors but may not remove all instabilities caused by insufficient resolution (for FR as well as DG).

High-Order↗

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗

Developing VSC-HVDC Oscillation Damping Control Constraints in Unit Commitment

High-voltage direct current (HVDC) systems within the existing power grids will play a pivotal role in enabling high share of renewable power integration and transportation electrification. There is a need for transmission planning models to accurately capture the impact of VSC-HVDC oscillation damping control on the system planning and model the constraints in the security constrained unit commitment (SCUC) problem. In this paper, the SCUC problem in a hybrid DC-AC grid considering damping on electromechanical inter-area oscillation is studied. An electromechanical oscillation-driven transmission capability constraints are derived from system swing equation and incorporated into the SCUC. Multiple system damping control scenarios are set based on the extracted plausible bounds of damping coefficient and inertia constant to linearize the nonlinear expressions. The SCUC problem with proposed constraints is tested on the Reliability Test System Grid Modernization Lab Consortium (RTS-GMLC) via PLEXOS. The impact of the proposed constraints on regional generation, inertia, power flow, annual system cost, and renewable energy curtailment are fully assessed.

damping control↗

Thermal interaction between an impinging hot jet and a conducting solid surface

Powered-lift aircraft may produce severe high-temperature environments which are potentially damaging to a landing surface or the aircraft. The interaction betweean the high temperature flow field and a nonadiabatic landing surface is analyzed with a coupled computational fluid dynamics/solid thermal conduction computer code, HOTJET. The HOTJET code couples time-accurate, implicit, factored solution schemes for the governing fluid dynamics equations (Reynolds-averaged Navier-Stokes equations) to the unsteady thermal conduction equation, which governs heat flux within a solid. HOTJET is validated against exact solutions to the thermal conduction and Navier-Stokes equations. First-of-a-kind results are included which show the impact of surface material properties on the fluid physics and the coupled fluid/material thermal fields.

Abeloff, P. A.↗

Numerical Experiments in Error Control for Sound Propagation Using a Damping Layer Boundary Treatment

This paper presents results from numerical experiments for controlling the error caused by a damping layer boundary treatment when simulating the propagation of an acoustic signal from a continuous pressure source. The computations are with the 2D Linearized Euler Equations (LEE) for both a uniform mean flow and a steady parallel jet. The numerical experiments are with algorithms that are third, fifth, seventh and ninth order accurate in space and time. The numerical domain is enclosed in a damping layer boundary treatment. The damping is implemented in a time accurate manner, with simple polynomial damping profiles of second, fourth, sixth and eighth power. At the outer boundaries of the damping layer the propagating solution is uniformly set to zero. The complete boundary treatment is remarkably simple and intrinsically independant from the dimension of the spatial domain. The reported results show the relative effect on the error from the boundary treatment by varying the damping layer width, damping profile power, damping amplitude, propagtion time, grid resolution and algorithm order. The issue that is being addressed is not the accuracy of the numerical solution when compared to a mathematical solution, but the effect of the complete boundary treatment on the numerical solution, and to what degree the error in the numerical solution from the complete boundary treatment can be controlled. We report maximum relative absolute errors from just the boundary treatment that range from O[10-2] to O[10-7].

Algorithm↗

Numerical Experiments in Error Control for Sound Propagation Using a Damping Layer Boundary Treatment

This paper presents results from numerical experiments for controlling the error caused by a damping layer boundary treatment when simulating the propagation of an acoustic signal from a continuous pressure source. The computations are with the 2D Linearized Euler Equations (LEE) for both a uniform mean flow and a steady parallel jet. The numerical experiments are with algorithms that are third, fifth, seventh and ninth order accurate in space and time. The numerical domain is enclosed in a damping layer boundary treatment. The damping is implemented in a time accurate manner, with simple polynomial damping profiles of second, fourth, sixth and eighth power. At the outer boundaries of the damping layer the propagating solution is uniformly set to zero. The complete boundary treatment is remarkably simple and intrinsically independant from the dimension of the spatial domain. The reported results show the relative effect on the error from the boundary treatment by varying the damping layer width, damping profile power, damping amplitude, propagtion time, grid resolution and algorithm order. The issue that is being addressed is not the accuracy of the numerical solution when compared to a mathematical solution, but the effect of the complete boundary treatment on the numerical solution, and to what degree the error in the numerical solution from the complete boundary treatment can be controlled. We report maximum relative absolute errors from just the boundary treatment that range from O[10-2] to O[10-7].

Computational Aeroacoustics↗

TRITON: A Multi-GPU open source 2D hydrodynamic flood model

A new open source multi-GPU 2D flood model called TRITON is presented in this work. The model solves the 2D shallow water equations with source terms using a time-explicit first order upwind scheme based on an Augmented Roe's solver that incorporates a careful estimation of bed strengths and a local implicit formulation of friction terms. Here, the scheme is demonstrated to be first order accurate, robust and able to solve for flows under various conditions. TRITON is implemented such that the model effectively utilizes heterogeneous architectures, from single to multiple CPUs and GPUs. Different test cases are shown to illustrate the capabilities and performance of the model, showing promising runtimes for large spatial and temporal scales when leveraging the computer power of GPUs. Under this hardware configuration, communication and input/output subroutines may impact the scalability. The code is developed under an open source license and can be freely downloaded in https://code.ornl.gov/hydro/triton.

2D flood model↗

Characterization and Simulation of the Thermoacoustic Instability Behavior of an Advanced, Low Emissions Combustor Prototype

Extensive research is being done toward the development of ultra-low-emissions combustors for aircraft gas turbine engines. However, these combustors have an increased susceptibility to thermoacoustic instabilities. This type of instability was recently observed in an advanced, low emissions combustor prototype installed in a NASA Glenn Research Center test stand. The instability produces pressure oscillations that grow with increasing fuel/air ratio, preventing full power operation. The instability behavior makes the combustor a potentially useful test bed for research into active control methods for combustion instability suppression. The instability behavior was characterized by operating the combustor at various pressures, temperatures, and fuel and air flows representative of operation within an aircraft gas turbine engine. Trends in instability behavior versus operating condition have been identified and documented, and possible explanations for the trends provided. A simulation developed at NASA Glenn captures the observed instability behavior. The physics-based simulation includes the relevant physical features of the combustor and test rig, employs a Sectored 1-D approach, includes simplified reaction equations, and provides time-accurate results. A computationally efficient method is used for area transitions, which decreases run times and allows the simulation to be used for parametric studies, including control method investigations. Simulation results show that the simulation exhibits a self-starting, self-sustained combustion instability and also replicates the experimentally observed instability trends versus operating condition. Future plans are to use the simulation to investigate active control strategies to suppress combustion instabilities and then to experimentally demonstrate active instability suppression with the low emissions combustor prototype, enabling full power, stable operation.

DeLaat, John C.↗

Topological Control of Triply Periodic Minimal Surfaces for Thermal Design and Advanced Manufacturing: A Gyroid Case Study

Recently, there has been a heightened interest in using triply periodic minimal surfaces (TPMSs) in the design of compact process engineering components. The benefits of high surface area per unit volume, modular form, and inherent periodicity provide a holistic self-supporting network and flow-conducive features. Applications of importance include thermal power management, biomimetic scaffolds and structures, and feasibility of advanced manufacturing. This study presents a novel approach to the manipulation of the characteristic Schwarz-G, or gyroid TPMS, for thermal design in the context of advanced manufacturing. The study presents relationships between design parameters and resulting surface area as a target response using the characteristic equation of a gyroid. Through parametric control, the characteristic equation is manipulated to produce a 20-fold increase in achievable area over a baseline design characteristic of 25.4 mm through controlled combinations of design parameters. A second relationship is presented as a function of the maximum area achieved and manipulated design parameters. Through the analysis, the study presents a framework to identify and maximize the achievable area of TPMSs for advanced manufacturing and thermal management applications.

gyroid↗

Magnetic shielding of large high-power-satellite solar arrays using internal currents

Present concepts for solar power satellites involve dimensions up to tens of kilometers and operating internal currents up to hundreds of kiloamperes. A question addressed is whether the local magnetic fields generated by these strong currents during normal operation can shield the array against impacts by plasma ions and electrons (and from thruster plasmas) which can cause possible losses such as power leakage and surface erosion. One of several prototype concepts was modeled by a long narrow rectangular panel 2 km wide and 20 km long. The currents flow in a parallel across the narrow dimension (sheet current) and along the edge (wire currents). The wire currents accumulate from zero to 100 kiloamp and are the dominant sources. The magnetic field is approximated analytically. The equations of motion for charged particles in this magnetic field are analyzed. The ion and electron fluxes at points on the surface are represented analytically for monoenergetic distributions and are evaluated.

Parker, L. W.↗

Novel Integrated Radial and Axial Magnetic Bearing

Typically, fully active magnetically suspended systems require one axial and two radial magnetic bearings. Combining radial and axial functions into a single device allows for more compact and elegant packaging. Furthermore, in the case of high-speed devices such as energy storage flywheels, it is beneficial to minimize shaft length to keep rotor mode frequencies as high as possible. Attempts have been made to combine radial and axial functionality, but with certain drawbacks. One approach requires magnetic control flux to flow through a bias magnet reducing control effectiveness, thus resulting in increased resistive losses. This approach also requires axial force producing magnetic flux to flow in a direction into the rotor laminate that is undesirable for minimizing eddy-current losses resulting in rotational losses. Another approach applies a conical rotor shape to what otherwise would be a radial heteropolar magnetic bearing configuration. However, positional non-linear effects are introduced with this scheme and the same windings are used for bias, radial, and axial control adding complexity to the controller and electronics. For this approach, the amount of axial capability must be limited. It would be desirable for an integrated radial and axial magnetic bearing to have the following characteristics, separate inputs for radial and axial control for electronics and control simplicity, all magnetic control fluxes should only flow through their respective air gaps and should not flow through any bias magnets for minimal resistive losses, be of a homopolar design to minimize rotational losses, position related non-linear effects should be minimized, and dependent upon the design parameters, be able to achieve any radial/axial force or power ratio as desired. The integrated radial and axial magnetic bearing described in this paper exhibits all these characteristics. Magnetic circuit design, design equations, and analysis results will be presented.

Blumenstock, Kenneth A.↗

Novel Integration Radial and Axial Magnetic Bearing

Typically, fully active magnetically suspended systems require one axial and two radial magnetic bearings. Combining radial and axial functions into a single device allows for more compact and elegant packaging. Furthermore, in the case of high-speed devices such as energy storage flywheels, it is beneficial to minimize shaft length to keep rotor mode frequencies as high as possible. Attempts have been made to combine radial and axial functionality, but with certain drawbacks. One approach requires magnetic control flux to flow through a bias magnet reducing control effectiveness, thus resulting in increased resistive losses. This approach also requires axial force producing magnetic flux to flow in a direction into the rotor laminate that is undesirable for minimizing eddy-current losses resulting in rotational losses. Another approach applies a conical rotor shape to what otherwise would be a radial heteropolar magnetic bearing configuration. However, positional non-linear effects are introduced with this scheme and the same windings are used for bias, radial, and axial control adding complexity to the controller and electronics. For this approach, the amount of axial capability must be limited. It would be desirable for an integrated radial and axial magnetic bearing to have the following characteristics; separate inputs for radial and axial control for electronics and control simplicity, all magnetic control fluxes should only flow through their respective air gaps and should not flow through any bias magnets for minimal resistive losses, be of a homopolar design to minimize rotational losses, position related non-linear effects should be minimized, and dependent upon the design parameters, be able to achieve any radial/axial force or power ratio as desired. The integrated radial and axial magnetic bearing described in this paper exhibits all these characteristics. Magnetic circuit design, design equations, and magnetic field modeling results will be presented.

Blumenstock, Kenneth↗