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 415 records · Page 23

Nearly spherical constant-power detonation waves as driven by focused radiation.

Analysis is made of the flow within a 3-D spark created in a gas absorbing energy from a steady conical beam of radiation (typically obtained from an array of lasers) with nearly spherical symmetry. First-order deviations are represented by series of Legendre polynomials of the polar angle. Self-similarity is demonstrated. For harmonics 1 to 5 radial profiles of velocity components, density and pressure are presented. The linear perturbation scheme exhibits a singularity at the focus. This reflects the existence of a 'forbidden zone', as is shown by a study of the full non-linear flow equations.

George, Y. H.↗

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Approximate Turbulent Boundary-Layer Development in Plane Compressible Flow Along Thermally Insulated Surfaces with Application to Supersonic-Tunnel Contour Correction

Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

BOUNDARY LAYER - WIND TUNNEL WALLS↗

Survey of the status of finite element methods for partial differential equations

The finite element methods (FEM) have proved to be a powerful technique for the solution of boundary value problems associated with partial differential equations of either elliptic, parabolic, or hyperbolic type. They also have a good potential for utilization on parallel computers particularly in relation to the concept of domain decomposition. This report is intended as an introduction to the FEM for the nonspecialist. It contains a survey which is totally nonexhaustive, and it also contains as an illustration, a report on some new results concerning two specific applications, namely a free boundary fluid-structure interaction problem and the Euler equations for inviscid flows.

Temam, Roger↗

Cosmic-ray heating of cooling flows - A critical analysis

It is shown that a combination of MHD wave-mediated cosmic ray heating and thermal conduction could balance cooling in intracluster media and substantially reduce the rate of inflow. The appropriate system of steady state equations is solved, including a new self-consistent formulation for the cosmic-ray diffusivity. Models which can produce substantial positive temperature gradients in static configurations are found when conduction is reduced by a factor of 10 or more. These models have too-flat thermal pressure profiles compared with observations. It is found that cosmic-ray heating is unlikely either to stabilize positive density perturbations against condensation or to contribute appreciably to the powering of the optical filaments.

Loewenstein, Michael↗

Sound generated in a cascade by three-dimensional disturbances convected in a subsonic flow

Discrete tone sound generation in a subsonic fan subject to three-dimensional disturbances is investigated. The analytical model used treats the fan rotor and stator as linear cascades of thin airfoils in a rectangular duct subject to a three-dimensional gust for which a complete aerodynamic theory already exists. The sound pressure can then be cast as the sum of a finite number of discrete sound waves (modes) the magnitude of which depends on an unknown function satisfying a singular integral equation. Similarity rules are derived to reduce the problem to that of a two-dimensional gust. Three-dimensional effects on the cut-off condition, the sound pressure, and the acoustic power are first investigated for each mode. The theory is then applied to noise generated by typical rotor-wake-defect and rotor-tip-vortex disturbances interacting with a stator.

Atassi, H.↗

Adaptive Mesh Refinement Large Eddy Simulation of the Supercritical Carbon Dioxide Round Turbulent Jet

Supercritical carbon dioxide (sCO2) is of interest to a range of engineering problems, including carbon capture, utilization, and storage (CCUS) as well as advanced cycles for power generation. Non-ideal variations in physical properties of sCO2 impact the physics of these systems. In this study, we simulate turbulent sCO2 jets to gain a better understanding of these physics.We use a second order finite volume method with adaptive mesh refinement as implemented in the first-principles simulation code PeleC to perform a Large Eddy Simulation (LES) of three turbulent jets of sCO2. Additionally, we use the Soave-Redlich-Kwong equation of state to close the system and examine the impact of a cubic equation of state on the turbulent flow physics. We look at velocity and Reynolds stress profiles at different downstream locations for three cases in which the temperature of the jet andthat of the ambient fluid differ in order to capture the effects of widely varying thermal properties in the pseudocritical region. These results are then contrasted with established theory for ideal gas jets.

adaptive mesh refinement↗

Experimental and analytical sonic nozzle discharge coefficients for Reynolds numbers up to 8 x 10 to the 6th power

Sonic discharge coefficients are obtained for two different geometry flow nozzles using high-pressure nitrogen gas (100 atm) with significant real-gas flow corrections. Throat Reynolds number range extended up to 8 million. Discharge coefficients for both nozzles monotonically increase in value at the high throat Reynolds numbers. The 95-percent confidence band for each nozzle is shown. Analytical discharge coefficients for the continuous and finite radius of curvature nozzle are presented. These analytical results for the laminar and turbulent boundary-layer cases are compared to experimental values for sonic flow. Experimental values are also compared to values calculated from the best empirical curve fit equation for subsonic flow.

Szaniszlo, A. J.↗

Learning unknown physics of non-Newtonian fluids

Here, we use physics-informed neural networks (PINNs) to learn viscosity models of two non-Newtonian systems (polymer melts and suspensions of particles) using only velocity measurements. For synthetic velocity data generated with the power-law viscosity model, the PINN-inferred viscosity model agrees with the analytical model for shear rates with large absolute values but deviates for shear rates near zero where the analytical model has an unphysical singularity. Once the viscosity model is learned the PINN method can solve the momentum conservation equation using only the boundary conditions.

42 ENGINEERING↗

Passive Ultrasonic Deterrents to Reduce Bat Mortality in Wind Farms (Final Report)

The overall objective of this project was to develop active and passive ultrasonic deterrent concepts to mitigate bat mortality at wind turbines. The research was supported by the US Department of Energy under the contract DE-EE0008731. Novel ultrasonic bat deterrents are investigated in this study. The deterrents are based on the idea of aerodynamic whistles wherein flow-acoustic resonance is used to produce high-amplitude tonal sound at desired ultrasonic frequencies. The deterrents can be grouped into active and passive. The active deterrents/whistles are driven by pressurized air supplied by an external source (e.g., an air compressor). The passive deterrents/whistles are designed to be mounted on wind turbine blades and are “powered” by the energy in the air moving past the rotor blades. Multiple active and passive deterrent ideas are proposed and investigated numerically and experimentally. The computational analysis involves solving the compressible unsteady Reynolds-averaged Navier-Stokes (URANS) equations and coupling the near-field aerodynamic solution with an integral method based on the Ffowcs Williams-Hawkings acoustic analogy to predict ultrasound in the farfield. Experiments are performed in the anechoic chamber at ISU (for the active whistles) and in the Stability wind tunnel at Virginia Tech (for the passive whistles).

17 WIND ENERGY↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Turbulence kinetic energy equation for dilute suspensions

A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Abou-Arab, T. W.↗

Skylon Aerodynamics and SABRE Plumes

An independent partial assessment is provided of the technical viability of the Skylon aerospace plane concept, developed by Reaction Engines Limited (REL). The objectives are to verify REL's engineering estimates of airframe aerodynamics during powered flight and to assess the impact of Synergetic Air-Breathing Rocket Engine (SABRE) plumes on the aft fuselage. Pressure lift and drag coefficients derived from simulations conducted with Euler equations for unpowered flight compare very well with those REL computed with engineering methods. The REL coefficients for powered flight are increasingly less acceptable as the freestream Mach number is increased beyond 8.5, because the engineering estimates did not account for the increasing favorable (in terms of drag and lift coefficients) effect of underexpanded rocket engine plumes on the aft fuselage. At Mach numbers greater than 8.5, the thermal environment around the aft fuselage is a known unknown−a potential design and/or performance risk issue. The adverse effects of shock waves on the aft fuselage and plumeinduced flow separation are other potential risks. The development of an operational reusable launcher from the Skylon concept necessitates the judicious use of a combination of engineering methods, advanced methods based on required physics or analytical fidelity, test data, and independent assessments.

Mehta, Unmeel↗

Reinforcement Learning for Anomaly Detection in Nuclear Power Plant Operation and Maintenance

In nuclear power plants (NPPs), timely identification of sensor and human errors is critical to ensure safe and efficient plant operations. Anomaly detection models can be employed for this task. However, traditional anomaly detection approaches may have high dependency on labeled datasets and struggle with adaptability in complex, dynamic environments. Reinforcement learning (RL) has demonstrated significant potential in fault diagnosis and anomaly detection; however, its application to anomaly detection in NPPs remains a relatively underexplored research direction. Hence, to address this gap, in this study, we present a novel physics-informed reinforcement learning model, PIRL-AD: Physics-Informed Reinforcement Learning for Anomaly Detection, that integrates domain knowledge from calorimetric equations into the RL framework for enhanced sensor and human error anomaly detection. We evaluate the performance of PIRL-AD against a non-physics informed RL benchmark and a support vector machine (SVM) on data collected from a forced flow loop testbed. Experimental results suggest that PIRL-AD outperforms other baselines on a range of anomalous datasets that include both sensor and human-induced anomalies across key performance metrics, statistically outperforming the RL and SVM benchmarks with respect to geometric mean (respectively, 92.96% vs. 91.06% vs. 83.01%) and F1-score (respectively, 89.23% vs. 86.98% vs. 77.01%). Furthermore, the findings suggest the potential of physics-integrated reinforcement learning models for enhanced anomaly detection performance in NPPs.

Reinforcement learning↗

Forearm blood flow during body temperature transients produced by leg exercise

Subjects exercised for 30 min on a bicycle ergometer at 30, 50, and 70% of maximal aerobic power in ambient temperatures of 15, 25, and 35 C and vapor pressures of less than 18 torr. Exercise was used to vary internal temperature during an experiment, and different ambient temperatures were used to vary skin temperatures independently of internal temperature. Forearm skin temperature was fixed at about 36.5 C. Esophageal temperature was measured with a thermocouple at the level of the left atrium, and mean skin temperature was calculated from a weighted mean of thermocouple temperatures at eight skin sites. Forearm blood flow was measured by electrocapacitance plethysmography. Data are well accounted for by a linear equation independent of exercise intensity, although some subjects showed an equivocal vasodilator effect of exercise.

Wenger, C. B.↗

Shock development in sound transmitted through a nearly choked flow

A nonlinear quasi-one dimensional theory of sound transmitted through a converging-diverging duct section is extended to the case where the acoustical source is located well downstream of the throat, at a point where the flow Mach number is low. The development and subsequent effects of shocks in the acoustic quantities are of primary consideration. The analysis uses a method of matched asymptotic expansions, yielding a set of inner equations of motion and shock conditions valid in the near-sonic throat region. The analysis leads to a generalization of the 'equal area' relation of weak shock theory. The manner in which nonlinear effects increase with source strength, frequency, and throat Mach number is illustrated by the numerical results and corresponding graphs. The shock waves are shown to cause significant dissipation in acoustic power.

Myers, M. K.↗

Turbulent Fluid Motion 6: Turbulence, Nonlinear Dynamics, and Deterministic Chaos

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The turbulence is sustained by a nonrandom time-independent external force. The solutions, on the average, separate exponentially with time, having a positive Liapunov exponent. Thus, the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number (or strength of the forcing) is reduced. Several qualitatively different flows are noted. These are, respectively, fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, we classify only the fully chaotic flows as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows, although having positive Liapunov exponents, have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, as generally understood, is both time-dependent and aperiodic.

Deissler, Robert G.↗

Nuclear modules for space electric propulsion

Analysis of interplanetary cargo and piloted missions requires calculations of the performances and masses of subsystems to be integrated in a final design. In a preliminary and scoping stage the designer needs to evaluate options iteratively by using fast computer simulations. The Oak Ridge National Laboratory (ORNL) has been involved in the development of models and calculational procedures for the analysis (neutronic and thermal hydraulic) of power sources for nuclear electric propulsion. The nuclear modules will be integrated into the whole simulation of the nuclear electric propulsion system. The vehicles use either a Brayton direct-conversion cycle, using the heated helium from a NERVA-type reactor, or a potassium Rankine cycle, with the working fluid heated on the secondary side of a heat exchanger and lithium on the primary side coming from a fast reactor. Given a set of input conditions, the codes calculate composition. dimensions, volumes, and masses of the core, reflector, control system, pressure vessel, neutron and gamma shields, as well as the thermal hydraulic conditions of the coolant, clad and fuel. Input conditions are power, core life, pressure and temperature of the coolant at the inlet of the core, either the temperature of the coolant at the outlet of the core or the coolant mass flow and the fluences and integrated doses at the cargo area. Using state-of-the-art neutron cross sections and transport codes, a database was created for the neutronic performance of both reactor designs. The free parameters of the models are the moderator/fuel mass ratio for the NERVA reactor and the enrichment and the pitch of the lattice for the fast reactor. Reactivity and energy balance equations are simultaneously solved to find the reactor design. Thermalhydraulic conditions are calculated by solving the one-dimensional versions of the equations of conservation of mass, energy, and momentum with compressible flow.

Difilippo, F. C.↗