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 163 records · Page 9

Miniature Heat Transport System for Spacecraft Thermal Control

Loop heat pipes (LHP) are efficient devices for heat transfer and use the basic principle of a closed evaporation-condensation cycle. The advantage of using a loop heat pipe over other conventional methods is that large quantities of heat can be transported through a small cross-sectional area over a considerable distance with no additional power input to the system. By using LHPs, it seems possible to meet the growing demand for high-power cooling devices. Although they are somewhat similar to conventional heat pipes, LHPs have a whole set of unique properties, such as low pressure drops and flexible lines between condenser and evaporator, that make them rather promising. LHPs are capable of providing a means of transporting heat over long distances with no input power other than the heat being transported because of the specially designed evaporator and the separation of liquid and vapor lines. For LHP design and fabrication, preliminary analysis on the basis of dimensionless criteria is necessary because of certain complicated phenomena that take place in the heat pipe. Modeling the performance of the LHP and miniaturizing its size are tasks and objectives of current research. In the course of h s work, the LHP and its components, including the evaporator (the most critical and complex part of the LHP), were modeled with the corresponding dimensionless groups also being investigated. Next, analysis of heat and mass transfer processes in the LHP, selection of the most weighted criteria from known dimensionless groups (thermal-fluid sciences), heat transfer rate limits, (heat pipe theory), and experimental ratios which are unique to a given heat pipe class are discussed. In the third part of the report, two-phase flow heat and mass transfer performances inside the LHP condenser are analyzed and calculated for Earth-normal gravity and microgravity conditions. On the basis of recent models and experimental databanks, an analysis for condensing two-phase flow regimes, pressure gradients, and local heat transfer coefficients using ammonia, propylene, and R134, are carried out.

Ochterbeck, Jay M.↗

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

Towards Robust Designs Via Multiple-Objective Optimization Methods

Fabricating and operating complex systems involves dealing with uncertainty in the relevant variables. In the case of aircraft, flow conditions are subject to change during operation. Efficiency and engine noise may be different from the expected values because of manufacturing tolerances and normal wear and tear. Engine components may have a shorter life than expected because of manufacturing tolerances. In spite of the important effect of operating- and manufacturing-uncertainty on the performance and expected life of the component or system, traditional aerodynamic shape optimization has focused on obtaining the best design given a set of deterministic flow conditions. Clearly it is important to both maintain near-optimal performance levels at off-design operating conditions, and, ensure that performance does not degrade appreciably when the component shape differs from the optimal shape due to manufacturing tolerances and normal wear and tear. These requirements naturally lead to the idea of robust optimal design wherein the concept of robustness to various perturbations is built into the design optimization procedure. The basic ideas involved in robust optimal design will be included in this lecture. The imposition of the additional requirement of robustness results in a multiple-objective optimization problem requiring appropriate solution procedures. Typically the costs associated with multiple-objective optimization are substantial. Therefore efficient multiple-objective optimization procedures are crucial to the rapid deployment of the principles of robust design in industry. Hence the companion set of lecture notes (Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks ) deals with methodology for solving multiple-objective Optimization problems efficiently, reliably and with little user intervention. Applications of the methodologies presented in the companion lecture to robust design will be included here. The evolutionary method (DE) is first used to solve a relatively difficult problem in extended surface heat transfer wherein optimal fin geometries are obtained for different safe operating base temperatures. The objective of maximizing the safe operating base temperature range is in direct conflict with the objective of maximizing fin heat transfer. This problem is a good example of achieving robustness in the context of changing operating conditions. The evolutionary method is then used to design a turbine airfoil; the two objectives being reduced sensitivity of the pressure distribution to small changes in the airfoil shape and the maximization of the trailing edge wedge angle with the consequent increase in airfoil thickness and strength. This is a relevant example of achieving robustness to manufacturing tolerances and wear and tear in the presence of other objectives.

Man Mohan, Rai↗

Reduction of wind tunnel wall interference by controlled wall flow

An alternate method of testing was developed in which flow through the porous walls of the tunnel was actively controlled so as to approximate free air conditions in the neighborhood of the model during the test. The amount and distribution of the controlled flow through the walls is computed using a potential flow representation of the model based on the measured lift. Theoretical analysis is presented to prove the convergence of the method to free air conditions and to substantiate the general three-dimensional theory of operation when the normal flow distribution is continuous. A two-dimensional tunnel was constructed to evaluate the concept. Results show that substantial reduction of wall interference may be achieved with relatively low values of porosity of actively controlled walls.

Bernstein, S.↗

Time-Accurate Unsteady Flow Simulations Supporting the SRM T+68-Second Pressure Spike Anomaly Investigation (STS-54B)

Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The objectives of the present analysis are to: (1) capture the dynamic responses of the motor combustion gas flow to correlate with available low-frequency (less than 12.5 sample/sec) data and (2) observe the high-frequency (up to 50 Hz) characteristics of the response to determine any potentials for dynamic coupling.

Dougherty, N. S.↗

CFD Modeling of Cryogenic Chilldown in a Complex Channel under Normal and Low Gravity Conditions

Future NASA architectures have baselined cryogenic propulsion systems as well as cryogenic fluid management to support lunar missions and ultimately to support future missions to Mars. These missions will require chilling hardware down prior to engine restart as well as chilling lines and tanks prior to transferring and refueling these propulsion elements in orbit. In lieu of expensive tests conducted on-orbit, accurate predictive computational models of these chilldown processes can be used to reduce system and propellant mass as well as mission risk. To gain confidence in these computational models, appropriate anchoring and validation to experimental data in a relevant environment needs to be performed. Recent ground and sub-orbital flight experiments conducted by the Japan Aerospace Exploration Agency (JAXA) investigated chilldown of a complex channel resembling a turbopump bearing cavity at low flow rates. This work presents Computational Fluid Dynamics (CFD) model development of the chilldown experiment employing two-phase flow boiling models available in commercial CFD software STAR-CCM+ using Volume of Fluid (VOF) and the traditional Euler-Euler multiphase flow solvers. Comparisons of the numerical and experimental results under normal and low-gravity conditions are presented. An assessment of solid wall temperatures and phase distribution yielded important insights into multiphase solver choice, dependence on gravity environment, and challenges associated with cryogenic flow boiling prediction and validation.

CFD↗

Multiphysics modeling of the FW/Blanket of the U.S. fusion nuclear science facility (FNSF)

The dual coolant lead-lithium (DCLL) blanket concept, which is utilized in the Fusion Nuclear Science Facility (FNSF) conceptual design, is based on a helium-cooled first wall and blanket structure with RAFS (Reduced Activation Ferritic Steel) and a self-cooled LiPb breeding zone. The objective of this work is to develop a multiphysics modeling process in order to optimize the design and achieve long lifetime, maintainability, and high reliability. 3D finite element multiphysics modeling of the DCLL first wall and blanket (midplane of one sector) has been performed using COMSOL 5.2. The multiphysics aspect of the design is demonstrated via coupling of Computational Fluid Dynamics (CFD), conjugate heat transfer and solid mechanics. Both normal and off-normal loading conditions have been analyzed. The results of velocity, pressure, and temperature distributions of helium flow, as well as the primary and thermal stress of the structure were obtained. This was followed by determination of the factors of safety along three critical paths based on the ITER Structural Design Criteria for In-vessel Components (ISDC-IC). We show here that the structural design meets the ITER-ISDC design rules under both normal and off-normal operating conditions, though the safety factors under off-normal condition with 8 MPa helium pressure are marginal. Thus simple design optimization was conducted based on a parametric study on first wall dimensions to improve the design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of a low order panel method to complex three-dimensional internal flow problems

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

Ashby, D. L.↗

Axisymmetric internal flows with shocks

Internal flow problems with supersonic entrance conditions and subsonic exit conditions are studied. The suitability of applying the transonic-small disturbance theory to internal flows is examined and found to be very limited. The full inviscid equations in nonconservative form are solved by relaxation. Jump conditions which conserve mass and normal momentum are applied explicitly at shock waves. This method is suitable for solving flows which include supersonic-subsonic shock waves nearly normal to the streamlines.

Rizk, M. H.↗

High-Pressure, High-Temperature Thermal Hydraulic Test Facility for Nuclear-Renewable Hybrid Energy System Studies; Facility Design Description and Status Report

A high-pressure, high-temperature water flow loop has been designed at Idaho National Laboratory (INL) for deployment in the INL Dynamic Energy Transport and Integration Laboratory (DETAIL) within the Energy Systems Laboratory (ESL) D100 northwest high bay. This loop is the first of three thermally coupled flow loops that will comprise the Advanced Reactor Technology Integral System Test (ARTIST) facility. The ARTIST facility will ultimately include a high temperature helium loop and a liquid salt loop, in addition to the PWR loop. The high-pressure water loop will operate at Pressurized Water Reactor (PWR) conditions. It will support experimental research on nuclear-renewable hybrid energy systems as well as advanced reactor technology topics. Within DETAIL, the water flow loop will also serve as a thermal energy source, emulating a reactor system. The loop will be thermally integrated with co-located energy systems including a thermal energy transport loop, thermal energy storage system, and a 25 kW high-temperature electrolysis system for hydrogen production. DETAIL will include additional energy systems such as an integrated microgrid and renewable energy systems including solar photovoltaics. The electrically heated flow loop will be dynamically controlled to simulate nuclear fuel behavior under normal and off-normal operating conditions. DETAIL will be designed for characterization of complex system dynamics. The various components will be interfaced with Digital Real-Time Simulators for supervisory control and dynamic simulation of complete systems. The design and assembly of the high-pressure water flow loop and supporting systems has been subdivided into two Phases. Phase 1, which will be completed by the end of FY17, includes installation of cable trays and supports, electrical conduit, DI water system, drain lift station and the concrete pad for the chiller located outside of the laboratory. Phase 2 includes the detailed design and final assembly of the flow loop, supporting structures and personnel platforms, the chiller and associated piping. As of September, 2017, the final detailed design of the flow loop is complete. Funding for Phase 2 procurement, assembly, system operability checkout, shakedown testing and initial operation of the flow loop is requested for FY18.

42 ENGINEERING↗

Radiant extinction of gaseous diffusion flames

The absence of buoyancy-induced flows in microgravity significantly alters the fundamentals of many combustion processes. Substantial differences between normal-gravity and microgravity flames have been reported during droplet combustion, flame spread over solids, candle flames, and others. These differences are more basic than just in the visible flame shape. Longer residence time and higher concentration of combustion products create a thermochemical environment which changes the flame chemistry. Processes such as flame radiation, that are often ignored under normal gravity, become very important and sometimes even controlling. This is particularly true for conditions at extinction of a microgravity diffusion flame. Under normal-gravity, the buoyant flow, which may be characterized by the strain rate, assists the diffusion process to transport the fuel and oxidizer to the combustion zone and remove the hot combustion products from it. These are essential functions for the survival of the flame which needs fuel and oxidizer. Thus, as the strain rate is increased, the diffusion flame which is 'weak' (reduced burning rate per unit flame area) at low strain rates is initially 'strengthened' and eventually it may be 'blown-out'. Most of the previous research on diffusion flame extinction has been conducted at the high strain rate 'blow-off' limit. The literature substantially lacks information on low strain rate, radiation-induced, extinction of diffusion flames. At the low strain rates encountered in microgravity, flame radiation is enhanced due to: (1) build-up of combustion products in the flame zone which increases the gas radiation, and (2) low strain rates provide sufficient residence time for substantial amounts of soot to form which further increases the flame radiation. It is expected that this radiative heat loss will extinguish the already 'weak' diffusion flame under certain conditions. Identifying these conditions (ambient atmosphere, fuel flow rate, fuel type, etc.) is important for spacecraft fire safety. Thus, the objective is to experimentally and theoretically investigate the radiation-induced extinction of diffusion flames in microgravity and determine the effect of flame radiation on the 'weak' microgravity diffusion flame.

Atreya, Arvind↗

Development Plan for Ultrasonic Pressure Sensors [Slides]

Current light water reactors (LWRs) utilize differential pressure cells, based on strain gauges, which are remotely located from the reactor vessel via impulse lines. These pressure sensors are used to measure liquid level within the vessel and to give a general indication of system health with regard to leaks and flow system fouling. The advanced reactor types currently under development generally operate at higher temperatures than light-water reactors and utilize coolants that may be highly corrosive (i.e., liquid metals). They will use pressure measurements for varied measurands, such as level and flow rate, as well as indications of accident conditions or off-normal operating conditions. More robust pressure sensors are needed to address these more challenging conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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.↗

Efficient Fluid Dynamic Design Optimization Using Cartesian Grids

This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.

Dadone, A.↗

The instability of sheared liquid layers

The critical values of the Reynolds number (Rc) above which unstable travelling waves exist in a steady shear flow are calculated with the inclusion of a term in the normal stress boundary condition omitted by Miles (1960). Results show a value of Rc = 34.2, compared to Miles' value of 203, for the case of a zero-surface tension interface. The critical wavenumber, critical phase speed, and Rc are calculated as functions of a nondimensional surface tension, and the mechanism of the instability is investigated by examining the disturbance energy equation. A long wave instability is found at small values of the Reynolds number when the shear flow has a parabolic velocity profile. In addition, these results are extended to larger values of the wavenumber by numerical methods, and the relation between this long wave instability and the profile curvature is examined.

Smith, M. K.↗

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↗

Spent Crystalline Silicotitanate Storage Study

Washington River Protection Solutions is working to support initial production of immobilized low-activity waste (LAW) by feeding Hanford tank supernate from tank farms to the Hanford Waste Treatment and Immobilization Plant (WTP) LAW Facility. This goal incorporates the design of a Tank-Side Cesium Removal (TSCR) system, which filters tank waste supernate to remove suspended solids and then removes Cs by processing it through crystalline silicotitanate (CST) ion exchange media manufactured by Honeywell UOP, LLC. The 137 Cs-depleted product is intended to be sent to the WTP for vitrification. The Cs-loaded CST columns will be stored indefinitely, with a goal of eventual CST removal and treatment. Thus, the spent CST needs to be recoverable. The testing described herein looks to potential upset process conditions where CST storage may be required before various rinse steps are completed. This study evaluated upset conditions at three sequential processing steps envisioned for TSCR (feed, 0.1 M NaOH rinse, water rinse), and in-column drying with compressed air (normal end step) to assess impacts on the nature of the CST bed. Testing was conducted at the small scale (12-mL bed volume); simulated AP-105 tank waste was used as the feed. Following process disruption, the CST bed was dried in place at 70 °C. Post-dried CST bed physical properties (angle of repose, penetration depth, particle morphological changes) were measured to evaluate how CST moved and flowed. The testing is intended to provide a preliminary assessment of issues that may arise from desiccation of CST with the indicated salt solutions in place. Since these were small-scale tests, the processing conditions will not match full scale conditions exactly; however, the tests do provide insight into the impact of stopping processing at an earlier step than normal. Except for the feed that was dried in-place, all other process stop-conditions showed the CST bed flowed well after drying. At this small scale, CST beds would not present an issue for retrievability. The feed that was dried in place had solidified into a rock-hard monolith with no movement possible. Samples had to be chipped from the surface.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Further developments in the prediction of oscillatory aerodynamics in mixed transonic flow

This paper presents improvements made to the oscillatory transonic aerodynamic kernel function method given in AIAA Paper No. 74-359. The improved method is fully developed for flows with continuously varying local Mach number and oscillating normal shocks. The boundary conditions developed for oscillating shocks require that a line doublet be added at the foot of the shock which represents the loads induced by shock movement. Numerical results show good agreement with experiment and significant deviation from linear aerodynamic theory. The method is also used to predict wind-tunnel wall interference effects.-

Cunningham, A. M., Jr.↗