Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional Normalizing Flow”

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 469 records · Page 26

Matching Inviscid/Boundary Layer Flowfields

New boundary-layer equations are developed and the solutions match all the boundary-layer properties, except for the normal velocity, exactly with the corresponding inviscid properties. The numerical procedure solves tridiagonal matrices at each marching station. As part of the solution an inviscid transpiration velocity at the surface is calculated from the boundary-layer solution. This transpiration velocity could be used as a boundary condition to calculate a new inviscid solution. If the inviscid/boundary layer solutions are iterated, then the normal velocity from the boundary layer solution will match the inviscid values exactly. Solutions are calculated for shear flows over a flat plate. Results from the present method compared well with Navier-Stokes solutions for incompressible constant shear flow and sinusoidal shear flow. Iterations of the inviscid/boundary layer solutions were not necessary. Compressible flow of an exponential shear at Mach 4 was calculated for cold, adiabatic, and hot walls. It was found that inviscid shear flows have a significant effect on skin friction and heating rates. The present method should yield accurate boundary-layer solutions at Reynolds numbers lower than the traditional boundary-layer equations.

DeJarnette, Fred R.↗

Influence of Source Propagation Direction and Shear Flow Profile in Impedance Eduction of Acoustic Liners

The acoustic impedance of liners is a key parameter for their design, and depends on the flow conditions, i.e., the sound pressure level and the presence of a grazing flow. The surface impedance of a locally reacting liner is defined as a local intrinsic property relating the acoustic pressure to the normal acoustic particle velocity at the liner surface. Impedance eduction techniques are now widely used to retrieve the impedance of liners in aeroacoustic facilities in the presence of a shear grazing flow. While surface impedance is intrinsic by definition, the educed impedance has recently been shown to depend on the direction of the incident waves relative to the mean flow. Different studies have investigated this issue by considering different acoustic propagation models used in the education process in the hope of matching the educed values. The purpose of the present work is to continue the previous investigations by evaluating the influence of the shear flow profile on the educed impedance, while considering a Bayesian inference process in order to evaluate the uncertainty on the educed values. The identified uncertainties were not able to totally account for the observed discrepancies between educed impedances.

Roncen, R.↗

BLSTA: A boundary layer code for stability analysis

A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.

Wie, Yong-Sun↗

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗

A direct numerical simulation study of flame structure and stabilization of an experimental high Ka CH 4 /air premixed jet flame

In the present work, a direct numerical simulation (DNS) of an experimental high Karlovitz number (Ka) CH 4 /air piloted premixed flame was analyzed to study the inner structure and the stabilization mechanism of the turbulent flame. A reduced chemical mechanism for premixed CH 4 /air combustion with NO x based on GRI-Mech3.0 was used, including 268 elementary reactions and 28 transported species. The evolution of the stretch factor, I 0 , indicates that the burning rate per unit flame surface area is considerably reduced in the near field and exhibits a minimum at x /D = 8. Downstream, the burning rate gradually increases. The stretch factor is different between different species, suggesting the quenching of some reactions but not others. Comparison between the turbulent flame and strained laminar flames indicates that certain aspects of the mean flame structure can be represented surprisingly well by flamelets if changes in boundary conditions are accounted for and the strain rate of the mean flow is employed; however, the thickening of the flame due to turbulence is not captured. The spatial development of displacement speeds is studied at higher Ka than previous DNS. In contrast to almost all previous studies, the mean displacement speed conditioned on the flame front is negative in the near field, and the dominant contribution to the displacement speed is normal diffusion with the reaction contribution being secondary. Further downstream, reaction overtakes normal diffusion, contributing to a positive displacement speed. The negative displacement speed in the near field implies that the flame front situates itself in the pilot region where the inner structure of the turbulent flame is affected significantly, and the flame stabilizes in balance with the inward flow. Notably, in the upstream region of the turbulent flame, the main reaction contributing to the production of OH, H+O 2 ⇌O+OH (R35), is weak. Moreover, oxidation reactions, H 2 +OH⇌H+H 2 O (R79) and CO+OH⇌CO 2 +H (R94), are influenced by H 2 O and CO 2 from the pilot and are completely quenched. Hence, the entire radical pool of OH, H and O is affected. Furthermore, the fuel consumption layer remains comparably active and generates heat, mainly via the reaction CH 4 +OH⇌CH 3 +H 2 O (R93).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The pre-conceptual design of the nuclear island of ASTRID

The CEA is involved in a substantial effort on the ASTRID (Advanced Sodium Technological Reactor for Industrial Demonstration) pre-conceptual design in cooperation with EDF, as experienced Sodium-cooled Fast Reactor (SFR) operator, AREVA, as experienced SFR Nuclear Island engineering company and components designer, ALSTOM POWER as energy conversion system designer and COMEX NUCLEAIRE as mechanical systems designer. The CEA is looking for other partnerships, in France and abroad. The ASTRID preliminary design is based on a sodium-cooled pool reactor of 1500 MWth generating about 600 MWe, which is required to guarantee the representativeness of the reactor core and the main components with regard to future commercial reactors. ASTRID lifetime target is 60 years. Two Energy Conversion Systems are studied in parallel until the end of 2012: Rankine steam cycle or Brayton gas based energy conversion cycle. ASTRID design is guided by the following major objectives: improved safety, simplification of structures, improved In Service Inspection and Repair (ISIR), improved manufacturing conditions for cost reduction and increased quality, reduction of risks related to sodium fires and water/sodium reaction, and improved robustness against external hazards. The core is supported by a diagrid, which lay on a strong back to transfer the weight to the main vessel. AREVA is involved in a substantial effort in order to improve the core support structure in particular regarding the ISIR and the connection to primary pump. In the preliminary design, the primary system is formed by the main vessel and the upper closure comprising the reactor roof, two rotating plugs - used for fuel handling - and the components plugs located in the roof penetrations. The Above Core Structure deflects the sodium flow in the hot pool and provides support to core instrumentation and guidance of the control rod drive mechanisms. The number of the major components in the main vessel, primary pumps, Intermediate Heat Exchangers, and Decay Heat Exchangers are now under consideration. Under normal conditions, power release is achieved using the steam/water plant (in case of Rankine steam cycle) or the gas plant (in case of Brayton gas cycle). The diverse design and operating modes of Decay Heat Removal systems provide protection against common cause failures. A Decay Heat Removal system through the reactor vault is in particular studied with the objective to complement Direct Reactor Cooling systems. At this stage of the studies, the secondary system comprises four independent sodium loops (two and three sodium loops configurations are also investigated). Each loop includes one mechanical pump (or a large capacity Annular Linear Induction Electromagnetic Pump), and three modular Steam Generator Units characterized by once through straight tube units with a ferritic tube bundle; nevertheless, helical coil steam generator with tubes made of Alloy 800, and inverted type steam generator with a ferritic tube bundle are also investigated. The limited power of each modular Steam Generator Unit allows the whole secondary loop to withstand a large water/sodium reaction consecutive to the postulated simultaneous rupture of all the heat exchange tubes of one module. The arrangement of the components is based on the 'Regain' concept, in which the secondary pump is situated at a low level in the circuit; conventional arrangement, as SUPERPHENIX type, is a back-up option. Alternative arrangements based on gas cycles are also studied together with Na-gas heat exchanger design. This paper presents a status of the ASTRID pre-conceptual design. The most promising options are highlighted as well as less risky and back-up options. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Wind tunnel tests of a blade subjected to midchord torsional oscillation at high subsonic stall flutter conditions

A mechanical drive system for oscillating blades in a wind tunnel at frequencies up to 767 hertz and amplitudes of + or - 1.2 deg is described. High-speed motion pictures of schlieren images of the flow over a double-circular arc blade oscillating in harmonic motion about the midchord revealed extensive shock patterns at a nominal free stream Mach number of 0.7, a mean angle of attack of 4 deg, and reduced frequency of about 0.7. A phase lag resulting from the slow response of the flow to the motion of the blade increased with increasing reduced frequency. This phase lag, based on the difference between the time the blade attained its maximum angle of attack and the time required for the normal shock to reach its extreme downstream position, was nominally 100 deg at the above conditions.

Boldman, D. R.↗

Langevin model of low-energy fission

Since the earliest days of fission, stochastic models have been used to describe and model the process. For a quarter century, numerical solutions of Langevin equations have been used to model fission of highly excited nuclei, where microscopic potential-energy effects have been neglected. In this paper I present a Langevin model for the fission of nuclei with low to medium excitation energies, for which microscopic effects in the potential energy cannot be ignored. I solve Langevin equations in a five-dimensional space of nuclear deformations. The macroscopic-microscopic potential energy from a global nuclear structure model well benchmarked to nuclear masses is tabulated on a mesh of approximately 10 7 points in this deformation space. The potential is defined continuously inside the mesh boundaries by use of a moving five-dimensional cubic spline approximation. Because of reflection symmetry, the effective mesh is nearly twice this size. For the inertia, I use a (possibly scaled) approximation to the inertia tensor defined by irrotational flow. A phenomenological dissipation tensor related to one-body dissipation is used. A normal-mode analysis of the dynamical system at the saddle point and the assumption of quasiequilibrium provide distributions of initial conditions appropriate to low excitation energies, and are extended to model spontaneous fission. A dynamical model of postscission fragment motion including dynamical deformations and separation allows the calculation of final mass and kinetic-energy distributions, along with other interesting quantities. The model makes quantitative predictions for fragment mass and kinetic-energy yields, some of which are very close to measured ones. Varying the energy of the incident neutron for induced fission allows the prediction of energy dependencies of fragment yields and average kinetic energies. With a simple approximation for spontaneous fission starting conditions, quantitative predictions are made for some observables which are close to measurements. In conclusion, this model is able to reproduce several mass and energy yield observables with a small number of physical parameters, some of which do not need to be varied after benchmarking to 235 U (n, f) to predict results for other fissioning isotopes.

233U↗

Investigations on the Thermal-hydraulic Behavior of Accident Tolerant Fuel Cladding Materials

Since the 2011 Fukushima accident, significant research has been devoted to developing accident-tolerant fuel (ATF) cladding materials. These investigations have mostly focused on the ability to resist runaway steam oxidation and retain mechanical strength and structural integrity under thermal shocks. However, much remains unknown about the materials’ thermal-hydraulic behavior, particularly under light-water reactor (LWR) operating conditions. Two phenomena that determine safety margins in both normal and off-normal operating conditions - i.e., the critical heat flux (CHF) and Leidenfrost point (LFP) temperature - have not been investigated thoroughly for ATF materials. Surface properties, such as wettability and roughness, are known to influence pool and flow boiling CHF, as well as the LFP temperature. However, little is known about the surface wettability of ATF materials, particularly at LWR pressure and temperature. Little has been done on flow-boiling CHF and quenching heat transfer, especially for droplet quenching. Nonetheless, it is known that CHF in transient conditions, e.g., an exponentially escalating power transient, can be significantly different from that expected in steady-state operation. Very few studies have investigated transient CHF, either on ATF materials or under LWR pressure and temperature. The thorough understanding of transient CHF under prototypical reactor conditions will benefit not only the deployment of ATF but also the upcoming national efforts to study transient ATF behavior in the Transient Reactor Test Facility (TREAT) at Idaho National Laboratory. In summary, as a part of ATF development efforts, there is an urgent need to understand how and how much these materials may affect two-phase heat transfer phenomena in nuclear reactor conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Critical flow and pressure ratio data for LOX flowing through nozzles

LOX and LN2 data for two-phase critical flow through nozzles were acquired with precision control. The principal measured parameters were inlet conditions, critical flow rate and critical flow pressure ratio. It is conclusively demonstrated that the principle of corresponding states can be applied to two-phase choked flow through nozzles. It is also shown that the proper normalizing parameters are developed, and that current theories can provide an adequate means for extrapolating data to other fluids.

Hendricks, R. C.↗

Critical flow and pressure ratio data for LOX flowing through nozzles

LOX and LN2 data for two-phase critical flow through nozzles have been acquired with precision control. The principal measured parameters were inlet conditions, critical flow rate and critical flow pressure ratio. The data conclusively demonstrate that the principle of corresponding states can be applied to two-phase choked flow through nozzles. These data also demonstrate that the proper normalizing parameters have been developed and current theories can provide an adequate means for extrapolating to other fluids.

Hendricks, R. C.↗

Multivalued velocity field model of turbulence

The considered study is based on new theoretical concepts regarding a post-instability model of a fluid discussed by Zak (1980). The model permits the completion of the governing equations of turbulence by introducing multivalued fields of velocities. Attention is given to the mechanism of energy dissipation, the characteristic wave propagation, a simplified model, the formation of turbulence around stagnation points, the formulation of boundary conditions, and the mechanism of turbulence formation. The mechanism of turbulence formation can be understood as propagation of initial discontinuities from the boundaries into a flow with the characteristic velocity which is defined by the normal (to the boundary) velocity components. These components emerge at the boundary as a result of jumps in the tangential components due to the continuity equation.

Zak, M.↗

Interaction of upstream flow distortions with high Mach number cascades

Features of the interaction of flow distortions, such as gusts and wakes with blade rows of advance type fans and compressors having high tip Mach numbers are modeled. A typical disturbance was assumed to have harmonic time dependence and was described, at a far upstream location, in three orthogonal spatial coordinates by a double Fourier series. It was convected at supersonic relative to a linear cascade described as an unrolled annulus. Conditions were selected so that the component of this velocity parallel to the axis of the turbomachine was subsonic, permitting interaction between blades through the upstream as well as downstream flow media. A strong, nearly normal shock was considered in the blade passages which was allowed curvature and displacement. The flows before and after the shock were linearized relative to uniform mean velocities in their respective regions. Solution of the descriptive equations was by adaption of the Wiener-Hopf technique, enabling a determination of distortion patterns through and downstream of the cascade as well as pressure distributions on the blade and surfaces. Details of interaction of the disturbance with the in-passage shock were discussed. Infuences of amplitude, wave length, and phase of the disturbance on lifts and moments of cascade configurations are presented. Numerical results are clarified by reference to an especially orderly pattern of upstream vertical motion in relation to the cascade parameters.

Englert, G. W.↗

Dynamic load measurements with delta wings undergoing self-induced roll-oscillations

The aerodynamic forces acting on a delta wing, mounted on a free-to-roll sting-balance apparatus, were measured. Two wing planforms having leading edge sweeps of 76 and 80 deg were tested, but only the wing with the 80 deg sweep would undergo periodic self-induced roll oscillation. The time dependent forces and roll angles for this wing were then recorded for various test conditions. In these tests a considerable drop in the average normal force of the free-to-roll wing was measured, relative to the normal force obtained in the static tests. Also, the helium-bubble flow visualization technique was used to gain some insight into the periodic motion of the separated leading edge vortices.

Levin, D.↗

SSME thrust chamber simulation using Navier-Stokes equations

The capability of the PHOENICS fluid dynamics code in predicting two-dimensional, compressible, and reacting flow in the combustion chamber and nozzle of the space shuttle main engine (SSME) was evaluated. A non-orthogonal body fitted coordinate system was used to represent the nozzle geometry. The Navier-Stokes equations were solved for the entire nozzle with a turbulence model. The wall boundary conditions were calculated based on the wall functions which account for pressure gradients. Results of the demonstration test case reveal all expected features of the transonic nozzle flows. Of particular interest are the locations of normal and barrel shocks, and regions of highest temperature gradients. Calculated performance (global) parameters such as thrust chamber flow rate, thrust, and specific impulse are also in good agreement with available data.

Przekwas, A. J.↗

Remote noncontacting measurements of heat transfer coefficients for detection of boundary layer transition in wind tunnel tests

An infrared measurement system is used that consists of a laser heating source, an infrared camera for data acquisition, and a video recorder for data storage. A laser beam is scanned over an airfoil, heating its surface to a few degrees above ambient. An infrared camera then measures the temperature of the airfoil over a two-dimensional field, and these temperatures are stored as a function of time on a video recorder. The resulting temperature pictures are digitized and an iterative approximation algorithm is used to extract the heat transfer coefficient. The resulting values are normalized to the natural convection condition. The technique has been applied in low-speed wind tunnel tests and compared to well-established hot-film measurements which were made simultaneously to confirm the flow conditions. Heat transfer coefficients were determined using a linear scanning pattern, to indicate the position of natural and of artificially induced transition on an airfoil, at various wind speeds. The technique is shown to be sensitive to transition at low Mach numbers. The advantages of the technique are discussed.

Heath, D. Michele↗

OPERATIONAL TOPOLOGY ESTIMATION OF THREE PHASE UNBALANCED POWER DISTRIBUTION SYSTEMS WITH OUTAGES

Distribution grid, in general, is constructionally loopy but a radially operated network. Opening and closing of suitable tie-switches allows transition from one radial topology to another. Tie-switches ensure reliability and provide restoration paths when faults (outages) occur. The estimation of the system's current operational topology is a precursor to implementing any optimal control actions (during normal operation) or restorative actions (during outage condition). The existing literature in topology detection suffers with several limitations including: (1) computationally intensive nonlinear formulations that are not scalable for large distribution systems, (2) methods requiring exhaustive search for the correct topology by enumerating the complete topological search space, and (3) approaches requiring large volumes of time-series measurements and are not applicable for estimating continuously changing operational topology. Furthermore, the existing literature on distribution grid topology estimation cannot be extended to estimate the operational topology during outage conditions. For operational topology estimation during outages, our objective is to estimate both energized radial feeder and the outaged sections. Starting with a radial topology, these methods identify outaged feeder sections and hence, are inapplicable when the current radial operational topology is not known. In this work, we present a mathematical programming approach to estimate the operational topology of a three-phase unbalanced power distribution system under both outage and normal operating conditions. The aforementioned limitations of existing literature were countered by modeling the problem as a mixed-integer linear programming (MILP) problem which minimizes the errors in available measurements while satisfying the power flow and structural constraints of the system. Specically, a minimum weighted least absolute value estimator is proposed that uses the instantaneous line flow measurements, historical/ forecasted load data, and ping measurements from a limited set of smart meters and solves an MILP to identify the operational topology based on a linearized three-phase unbalanced power flow model. The proposed approach is tested on IEEE 123-bus feeder and 1069 bus multi-feeder test system and the numerical results prove the applicability of the proposed method to large three phase unbalanced grids with high errors in measurements. It is observed that the approach performs well even at high percentages of measurement errors.

GANDLURU NAGA, ANANDINI↗

Investigation of Operating Characteristics of an Engine Equipped with Modifications to Eliminate Fuel-Evaporation Icing

Two modified fuel-injection systems, a drilled-inducer type and a spinner type, that prevent serious fuel-evaporation icing were installed on a V-type, liquid-cooled aircraft engine and a preliminary investigation was conducted to determine the effect on engine operating characteristics. The spinner system was also ground - and flight tested on a twin-engine fighter airplane. Flight measurements of cylinder-head temperature over a range of fuel-air ratios and engine power conditions were made at an altitude of approximately 10,000 feet. Starting and accelerating of the engine on the ground were unaffected by the fuel-injection modifications. During the flight investigation, no appreciable variation occurred between the maximum and minimum cylinder-head temperature with the standard and modified system for the same power condition and no irregularity of mixture distribution could be detected throughout the power range of the engine. Normal mixture distribution was also indicated by a similar response of cylinder-head temperature for variations of fuel-air ratio at manifold pressures of 25 and 35 inches of mercury absolute. Both modified fuel-injection systems required less fuel-nozzle pressure than the standard system to obtain the desired fuel-air ratio for given air-flow condition.

Mulholland, Donald R.↗