Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “network 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 217 records · Page 12

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

Filling a tank with cryogenic fluid is more challenging than filling a tank with water or any other fluid that is in liquid state at atmospheric condition. Filling a tank with cryogenic fluid is a two-step process. First the tank and the transfer line mustbe chilled. Liquid cryogens start flowing into the tank only after the tank and transfer lines are chilled to the fluid saturation temperature. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor, caused by heat transfer from the warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. Filling a cryogenic tank in the absence of gravity is more challenging because in a non-stratified environment liquid propellantmay not settle at the tank bottom as it does on earth. There is a strong possibility that liquid propellant may exit through the vent valve, which is typically located at the top of the tank to vent propellant vapor. There are several methodsof filling atank in space: Charge-Hold-Vent method [1], Vented-Chill / No-Vent-Fill (VCNVF)[2],and No Vent Fill with TVS assisted injector [3].This paper addressesnumerical modelingofthe No Vent Fill (NVF) processwhere a TVS (Thermodynamic Vent System) augmented injector was used with the vent valve closed during the entire filling process. NVF tests were conducted in the CRYOgenic Orbital Testbed (CRYOTE) tankat NASA/Marshall Space Flight Centerin 2018.The Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code [4] was used to simulate the chilland fill process.

Cryogenics↗

Nodal Numerical Modeling of Submerged Helium Injection in a Cryogenic Propellant Tank

Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper presents a mathematical model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. This mathematical model was incorporated in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. The numerical model developed with GFSSP has been validated against two sets of experimental data and has been shown to predict both propellant subcooling and helium consumption to within 30% in most cases. The test data used for the model validation were taken in 1-g, but the mass diffusion model was developed to be applied in both 1-g and micro-g environments.

Cryogenics↗

Evaluation of Heat Transfer Coefficient from Velocity Distribution in Boundary Layer

This paper describes a method of evaluating heat transfer coefficient from a velocity distribution in the boundary layer. The power law velocity profile and universal velocity distribution in a smooth pipe have been used to evaluate the wall shear stress, and the heat transfer coefficient was evaluated from the wall shear stress using the Chilton-Colburn analogy. The predicted pressure drop was compared with the pressure drop computed by the Colebrook equation, and the predicted heat transfer coefficient was compared with the heat transfer coefficient calculated by the Dittus-Boelter equation. The accuracy of pressure drop prediction was within 1%, and the accuracy of heat transfer coefficient prediction was within 8%. This method has the potential to use a Navier-Stokes based CFD solution to evaluate the heat transfer coefficient in nodal or network flow codes.

Heat Transfer Coefficient↗

Evaluation of Heat Transfer Coefficient from Velocity Distribution in Boundary Layer

This paper describes a method of evaluating heat transfer coefficient from a velocity distribution in the boundary layer. The power law velocity profile and universal velocity distribution in a smooth pipe have been used to evaluate the wall shear stress, and the heat transfer coefficient was evaluated from the wall shear stress using the Chilton-Colburn analogy. The predicted pressure drop was compared with the pressure drop computed by the Colebrook equation, and the predicted heat transfer coefficient was compared with the heat transfer coefficient calculated by the Dittus-Boelter equation. The accuracy of pressure drop prediction was within 1%, and the accuracy of heat transfer coefficient prediction was within 8%. This method has the potential to use a Navier-Stokes based CFD solution to evaluate the heat transfer coefficient in nodal or network flow codes.

Heat Transfer Coefficient↗

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

This paper presents preliminary analysis using a finite volume based computational method for network flow analysis for the simulation of chill and fill of a cryogenic tank. An axisymmetric two-dimensional model of a spherical tank was developed using a nodal approach. No Vent Filling (NVF) was accomplished using a Thermodynamic Vent System assisted injector where the injector was cooled by the cold gas generated by Joule-Thompson expansion of the liquid propellant. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around spray droplets, and c) condensation of vapor at the interface of the cooled injector and ullage. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data and one-dimensional (1D) model predictions. It was demonstrated that the axisymmetric 2D model provided a more realistic prediction than the 1D model.

Cryogenics↗

A Liquefaction Model for A Cryogenically Refrigerated Ellipsoidal Tank

Over the past few decades NASA has been developing technology to produce cryogenic propellants on the surfaces of the Moon or Mars. In-situ propellant production reduces the amount of propellants needed to be taken to Moon/Mars and ultimately reduces mission cost. Utilizing Lunar/Martian resources, the produced gases are liquefied and stored prior to use on the ascent vehicle. This paper presents a model for the liquefaction process of gaseous propellants in a cryogenically refrigerated tank. The tank is cylindrical with elliptical top and bottom domes. A transient model is developed based on the mass and energy conservation principles and wall-gas and liquid-gas interfacial mass and heat transfer correlations. The model is incorporated into the Generalized Fluid System Simulation Program (GFSSP), an MSFC in-house general-purpose computer program for flow network analysis. Description of the model and comparison of predicted results with available test data is presented.

Liquefaction↗

Nodal Modeling of Submerged Helium Injection Pressurization of a Cryogenic Propellant Tank

Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper describes a thermodynamic model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. The physics of helium dissolution during the pressurization process is also modeled, primarily for liquid hydrogen propellant where the dissolution is more significant. The analyses were performed in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. Hydrogen propellant testing of the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) conducted at NASA/Glenn Research Facility was also analyzed. The data used for the model validation were taken in 1-g, but the model was developed to be applicable in both multi-g and micro-g environments.

tank pressurization↗

Nodal Modeling of Liquid Propellant Feed and Pressurization System

Nodal or network flow modeling plays an important role in the design and operation of the feed and pressurization system of a liquid rocket engine. Model development and execution time is relatively short for nodal codes in comparison to the Navier-Stokes based CFD codes. Nodal models also allow the inclusion of several components into one model to predict the behavior of a larger system. Unlike CFD models, the nodal models do not need very fine discretization of the flow field because they use empirical correlations to model fluid friction and heat transfer. This paper presents several applications of nodal modeling of liquid propellant feed and pressurization systems using the Generalized Fluid System Simulation Program (GFSSP), a nodal code developed at NASA/ Marshall Space Flight Center. GFSSP discretizes the flow field into nodes which are connected by branches. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The momentum equations are solved at the branches to calculate flow rates. Applications include a) tank pressurization by inert gas as well as autogenous pressurization by gaseous propellant, b) submerged pressurization by helium, c) self-pressurization due to boil-off of cryogenic propellant, d) chilldown of a transfer line of a cryogenic propellant feed system, and e) chilldown and filling of a cryogenic tank. Each of the above-mentioned models was verified and validated by comparing with test data.

Nodal Model↗

Nodal Modeling of Liquid Propellant Feed and Pressurization System

Nodal or network flow modeling plays an important role in the design and operation of the feed and pressurization system of a liquid rocket engine. Model development and execution time is relatively short for nodal codes in comparison to the Navier-Stokes based CFD codes. Nodal models also allow the inclusion of several components into one model to predict the behavior of a larger system. Unlike CFD models, the nodal models do not need very fine discretization of the flow field because they use empirical correlations to model fluid friction and heat transfer. This paper presents several applications of nodal modeling of liquid propellant feed and pressurization systems using the Generalized Fluid System Simulation Program (GFSSP), a nodal code developed at NASA/ Marshall Space Flight Center. GFSSP discretizes the flow field into nodes which are connected by branches. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The momentum equations are solved at the branches to calculate flow rates. Applications include a) tank pressurization by inert gas as well as autogenous pressurization by gaseous propellant, b) submerged pressurization by helium, c) self-pressurization due to boil-off of cryogenic propellant, d) chilldown of a transfer line of a cryogenic propellant feed system, and e) chilldown and filling of a cryogenic tank. Each of the above-mentioned models was verified and validated by comparing with test data.

Nodal Model↗

Detection of anomalous events

A system is described for receiving a stream of events and scoring the events based on anomalousness and maliciousness (or other classification). The system can include a plurality of anomaly detectors that together implement an algorithm to identify low-probability events and detect atypical traffic patterns. The anomaly detector provides for comparability of disparate sources of data (e.g., network flow data and firewall logs.) Additionally, the anomaly detector allows for regulatability, meaning that the algorithm can be user configurable to adjust a number of false alerts. The anomaly detector can be used for a variety of probability density functions, including normal Gaussian distributions, irregular distributions, as well as functions associated with continuous or discrete variables.

97 MATHEMATICS AND COMPUTING↗

Modelling moisture sources due to air leakage with the HAMT extension for EnergyPlus

The moisture performance of building envelope components is one of the main challenges of energy efficient buildings. The combined Heat and Moisture Transfer (HAMT) model in EnergyPlus can provide moisture profiles through multiple material layers of an envelope assembly using a coupled, one-dimensional, finite volume model simulating the movement and storage of heat and moisture. Currently only a limited assessment of the effects of moisture sources inside an assembly to determine impacts on durability, indoor air quality, and energy performance is possible. This paper discusses different strategies to model moisture sources in building components due to air leakage. Models that enable air leakage to result in a moisture sink/source at an arbitrary location extend the existing mass transfer equation with a source term that includes water vapor deposit from air flow through an assembly. New sub-models with different levels of complexity to calculate the air flow through the assembly are developed. The air flow can be user defined, it can be computed as a result of the stack height and a user defined ventilation overpressure or completely derived from an air-flow network. We explain the model enhancements, discuss the required inputs for every option together with related limitations, showcase the model application and make the source code publicly available for external use.

Antretter, Florian↗

Heat Transfer Experiments of a 1st Stage Blade Cascade for Supercritical CO2 Oxy-Combustion Turbine Application

The results of internally cooled 1st stage blade (S1B) cascade testing in a supercritical CO2 environment is presented. The turbine blade design has been previously established for the end application of an oxy-combustion turbine operating in the Allam-Fetvedt cycle with turbine inlet conditions of 305 bar and 1150°C. The internally cooled blade features leading edge (LE) region impingement cooling, mid-section ribbed serpentine passages, and a pin-finned trailing edge (TE) region before cooling ejection holes. The geometry for the tested blade cascade has a cooled central blade with un-cooled blades on either side to match flowpath areas of the actual turbine. The flowpath reuses internal components previously employed for mid-section region ribbed serpentine passage experiments that established Nusselt number enhancement ratios over a range of Reynolds numbers from 100,000-400,000. New components include flow conditioning plates upstream and downstream of the blade cascade to adequately represent the flow field and blade external heat transfer coefficient profiles for the actual turbine. The cooled central blade utilizes uniform crystal temperature sensors (UCTS) with six sensors each on the blade pressure and suction surfaces distributed radially and from LE to TE. The post-processed UCTS quantified the maximum wall temperature seen at each installed sensor location. The test procedure consisted of establishing supercritical CO2 cooling flow temperature and flow rate and maintaining it throughout the test. The flow rate aims to match that for the actual in-service turbine blade design and is maintained through an orifice restriction to keep the pressure differential between internal cooling flow and external hot flow nearly constant. For the sCO2 flow path external to the blade, temperatures were ramped throughout the test via control of the test loop’s natural gas burner heater. The maximum temperature seen was 468°C and held constant for a duration of 10 minutes at which the blade metal temperature was predicted to be at its maximum before ramping down. For the turbine blade design for service inlet conditions, external flow path computational fluid dynamics (CFD) results and an internal cooling 1-D thermal and hydraulic flow network model using experimentally validated correlations served as thermal finite element (FE) boundary conditions to predict blade metal temperatures. These predicted temperatures were subsequently utilized in a structural FE model to predict blade life ratings dictated by Haynes 282 creep strength data, having a strong dependence on temperature. The boundary conditions experienced during testing are used in the same workflow and compared to the experimental results, with the goal of validating the analysis methodology and providing insight on the uncertainty in local metal temperature predictions.

20 FOSSIL-FUELED POWER PLANTS↗

Optimal Power Flow of Radial Networks and Its Variations: A Sequential Convex Optimization Approach

This paper proposes a sequential convex optimization method to solve broader classes of optimal power flow (OPF) problems over radial networks. The non-convex branch power flow equation is decomposed as a second-order cone inequality and a non-convex constraint involving the difference of two convex functions. Provided with an initial solution offered by an inexact second-order cone programming relaxation model, this approach solves a sequence of convexified penalization problems, where concave terms are approximated by linear functions and updated in each iteration. It could recover a feasible power flow solution, which usually appears to be very close, if not equal, to the global optimal one. Two variations of the OPF problem, in which non-cost related objectives are optimized subject to power flow constraints and the convex relaxation is generally inexact, are elaborated in detail. One is the maximum loadability problem, which is formulated as a special OPF problem that seeks the maximal distance to the boundary of power flow insolvability. The proposed method is shown to outperform commercial nonlinear solvers in terms of robustness and efficiency. The other is the bi-objective OPF problem. A non-parametric scalarization model is suggested, and is further reformulated as an extended OPF problem by convexifying the objective function. It provides a single trade-off solution without any subjective preference. The proposed computation framework also helps retrieve the Pareto front of the bi-objective OPF via the e-constraint method or the normal boundary intersection method. This paper also discusses extensions for OPF problems over meshed networks based on the semidefinite programming relaxation method.

bi-objective optimization↗

Measurement of a radial flow profile with eddy current flow meters and deep neural networks

Eddy current flow meters (ECFMs) measure flows of conductive fluids. Recent interest in ECFMs has increased due to applications in advanced nuclear reactors. ECFMs are well suited for such applications, as they can provide non-invasive measurements of flow in fluids that are often difficult to measure. Traditionally, ECFMs are operated using an alternating current at a single frequency, limiting ECFMs to measure average fluid velocities, blockages, or voids. Here, we expand the capabilities of ECFMs by measuring the fluid radial velocity profile of liquid mercury. To accomplish this, we made several ECFM sensitivity measurements at a range of frequencies. Different frequencies vary the electromagnetic skin depth of the device. By adjusting frequencies, we probed the fluid velocity at various radial locations and constructed a flow-velocity profile. The relationship between the ECFM measurements and velocity profile is nonlinear and requires solving an inverse problem. Using electromagnetic finite-element simulations to train a deep neural network (DNN), we created a model that provides a stable general relationship between the sensitivity measurements of an ECFM and the fluid velocity profile. Using ECFM measurements of liquid mercury, our DNN model calculates a flow profile that agrees well with computational fluid dynamics (CFD) simulations. This technique has potential to improve flow monitoring for optimization, safe operation of conductive fluid loops, and/or validating complex CFD models.

47 OTHER INSTRUMENTATION↗

Neural network–based pore flow field prediction in porous media using super resolution

Direct pore-scale simulations of fluid flow through porous media are computationally expensive to perform for realistic systems. Previous works have demonstrated using the geometry of the microstructure of porous media to predict the velocity fields therein based on neural networks. However, such trained neural networks do not perform well for unseen porous media with a large degree of heterogeneity. In this study we propose that incorporating a coarse velocity field in the input of neural networks is an effective way to improve the prediction performance. The coarse velocity field can be simulated with a low computational cost and provides global information to regularize the ill-posedness of the learning problem, which is usually caused by the use of local geometries due to the computational resource constraints. We show that incorporating the coarse-mesh velocity field significantly improves the prediction accuracy of the fine-mesh velocity field by comparison to the prediction that relies on geometric information alone, especially for the porous medium with a large interior vuggy pore space. We also show the flexibility of training the network in using coarse velocity fields with various resolutions. The results suggest that even using coarse velocity field with a very low resolution, the predictions are still enhanced and close to the ground truths. The feasibility of the method is further demonstrated by testing the trained network on real rocks. This study highlights the merits of incorporating a coarse-mesh velocity field into the input for neural networks, which provides global, physics-based information for the model, thereby improving the model's generalization capability.

42 ENGINEERING↗