Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous 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 1,333 records · Page 74

HLPW-5: Overview and Workshop Summary

The Fifth AIAA CFD High-Lift Prediction Workshop was held with the goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this endeavor was the use of Technology Focus Groups (TFG), which included both mesh generation and flow-solver experts working together to accelerate advancements for their particular CFD methodology, by addressing key questions of importance prior to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop, and was used for three unique test cases. Wind-tunnel data were available for comparison for one of the test cases. Altogether, 365 datasets of CFD results were submitted by 47 teams, with 41 teams contributing to the multiple configurations of Case 1, 40 to Case 2, and 18 to Case 3. This paper provides a high-level summary of the results and conclusions from the workshop. As concluded from past workshops, fixed-grid Reynolds-averaged Navier-Stokes methods continued to be inaccurate and inconsistent for high-lift flows near maximum lift. However, application of mesh-adaptation technology helped to achieve improved consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics, particularly at maximum lift. Best practices for these methods were refined over the course of the workshop and new challenges were identified.

Adam M Clark↗

Computation of two-phase shear-layer flow using an Eulerian-Lagrangian analysis

A combined Eulerian-Lagrangian analysis which combines a linearized block implicit Navier-Stokes analysis for the continuous phase with a Lagrangian analysis for the discrete phase has been utilized to simulate a two-phase shear-layer flow in a two-dimensional channel flow configuration. The motion of the particles is tracked in computational coordinate space, resulting in improved computational efficiency, and interphase coupling terms for the Eulerian analysis are computed from the instantaneous distribution of the particles rather than the trajectory information. The use of the implicit Navier-Stokes analysis for the continuous phase has made it possible to efficiently use a highly stretched mesh. A low Reynolds number form of the k-epsilon turbulence model was used with near-wall resolution, thus eliminating the need to use the 'wall-function' approach. The streamwise velocity profiles for the continuous and the discrete phases have been compared with experimental data for two test cases.

Sabnis, Jayant S.↗

Monthly Sewer Monitoring Report for LLNL Livermore Site, April 2020

Lawrence Livermore National Laboratory (LLNL) collects effluent samples from its sewer outfall at the B196 Sewer Monitoring Station (SMS). Effluent flow-proportional composite samples are collected at the SMS on a daily (midnight-to-midnight), weekly (Thursday through Wednesday), and “monthly” (composited from daily) basis; effluent grab samples are also collected each month at that same location. Certified contract laboratories analyzes these compliance samples, and results (Tables 1–3) are used to establish LLNL compliance with the 2019–2020 Wastewater Discharge Permit (Permit 1250) granted by the City of Livermore Water Resources Division (WRD). Supplemental analyses for biochemical oxygen demand (BOD) and total suspended solids (TSS) are performed onsite and reported (Table 4) for the calculation and assessment of sewer service charges. Quarterly, effluent is sampled for metals (24-hour composite) and cyanide (grab sample) concentrations (mg/L). The results are presented in Table 5. In addition to the sampling noted above, the SMS continuously monitors LLNL sewage effluent in real-time for flow rate, pH, metals, and radioactivity. Standard measurement methods are used to monitor flow and pH to determine permit compliance. Unique nonstandard analytical methods monitor for the presence of metals and radioactivity using x-ray fluorescence (XRF) and gamma spectroscopy, respectively. If an anomalous condition is detected by the monitoring system, an alarm activates LLNL’s Sewer Diversion Facility (SDF) and Livermore Water Reclamation Plant (LWRP) operators are notified that a potential release has occurred.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Effects of Transient Purge Flow on Stator-Rotor Seal Performance

Abstract As modern engine designs target higher efficiencies through increased turbine inlet temperatures, critical turbine components are at increased risk of damage from conditions exceeding material melting temperatures. In particular, improperly designed underplatform hardware components are susceptible to damage when hot main gas path flow is ingested into the stator-rotor cavity. While all turbines inherently experience transients during operation, a majority of publicly available turbine studies have been executed using steady operating conditions or inherently transient “blowdown” rigs. For this reason, routine transient events are not well understood. To address this need, the present study utilized a combined experimental and computational approach. The test article is a continuous-duration, one-stage test turbine operating with true-scale engine hardware and seal geometries at engine-representative flow conditions. The nature of the continuous-duration facility uniquely supports the direct assessment of transient events through its ability to transition between steady-state operating conditions. The effects of a transient purge flow were investigated in this study to identify general trends for transient events in a full-scale engine. Results from multiple measurement techniques in the wheelspace region show an interdependence of transient purge flow with a thermal lag of the underplatform hardware. Through experiments conducted at different coolant-to-main gas path temperature ratios, the use of pressure measurements as an indicator of fully purged behavior was introduced, and a thermally driven influence on rim seal performance was quantified. The computational results show good agreement with experimental pressure measurements and provide insight into the physical mechanisms that drive the relationship between pressure and sealing effectiveness measurements observed for the tested geometry.

Engineering↗

A rapid method for the approximate determination of nonlinear solutions Application to aerodynamic flows

A method for determining highly accurate approximations to families of strongly nonlinear solutions which are either continuous or discontinuous, and which represent variations in some arbitrary parameters, is developed and evaluated. The procedure consists of defining a unit perturbation by employing two or more nonlinear solutions which differ from one another by a nominal change in some geometric or flow parameter, and then using that unit perturbation to predict a family of related nonlinear solutions over a range of parameter variation. Coordinate straining is used in determining the unit perturbation to account for the movement of discontinuities and maxima of high-gradient regions due to the perturbation. Although the procedure is generally applicable, results are presented here for nonlinear aerodynamic applications. Attention is focused in particular on transonic flows which are strongly supercritical and exhibit large surface shock movement over the parametric range studied; and on subsonic flows which display large pressure variations in the stagnation and peak suction pressure regions. Flows past both isolated airfoils and compressor cascades involving a variety of flow and geometry parameter changes are considered. Comparisons with the corresponding 'exact' nonlinear solutions indicate a remarkable accuracy and range of validity of such a procedure. Computational time is trivial.

Stahara, S. S.↗

A Hydrogen Iron Flow Battery with High Current Density and Long Cyclability Enabled Through Circular Water Management

The hydrogen‐iron (HyFe) flow cell has great potential for long‐duration energy storage by capitalizing on the advantages of both electrolyzers and flow batteries. However, its operation at high current density (high power) and over continuous cycling testing has yet to be demonstrated. In this paper, we discuss our design and demonstration of a water management strategy that supports high current and long cycling performance of a HyFe flow cell. Water molecules associated with the movement of protons from the iron electrode to the hydrogen electrode are sufficient to hydrate the membrane and electrode at a low current density of 100 mA cm −2 during the charge process. At higher charge current density, more aggressive measures must be taken to counter back‐diffusion driven by the acid concentration gradient between the iron and hydrogen electrodes. Our water management approach is based on water vapor feeding in the hydrogen electrode, and water evaporation in the iron electrode, thus enabling the high current density operation of 300 mA cm −2 .

25 ENERGY STORAGE↗

A Hydrogen Iron Flow Battery with High Current Density and Long Cyclability Enabled Through Circular Water Management

The hydrogen‐iron (HyFe) flow cell has great potential for long‐duration energy storage by capitalizing on the advantages of both electrolyzers and flow batteries. However, its operation at high current density (high power) and over continuous cycling testing has yet to be demonstrated. In this article, we discuss our design and demonstration of a water‐management strategy that supports high current and long‐cycling performance of a HyFe flow cell. Water molecules associated with the movement of protons from the iron electrode to the hydrogen electrode are sufficient to hydrate the membrane and electrode at a low current density of 100 mA cm −2 during the charge process. At higher charge current density, more aggressive measures must be taken to counter back‐diffusion driven by the acid concentration gradient between the iron and hydrogen electrodes. Our water‐management approach is based on water vapor feeding in the hydrogen electrode and water evaporation in the iron electrode, thus enabling high current density operation of 300 mA cm −2 .

Yan, Litao↗

Characteristics of Nozzle Clogging and Evolution of Oxide Inclusion for Al-Killed Ti-Stabilized 18Cr Stainless Steel

Characteristics of nozzle clogging and evolution of oxide inclusion were studied by sampling the deposits of submerged entry nozzle (SEN) and the molten steel in tundish during Al-killed Ti-stabilized 18Cr stainless steelmaking process. The deposits found in the clogging of SEN were mainly composed of frozen steel and inclusions. The inclusions in the deposits and tundish were both mainly (MgO-Al{sub 2}O{sub 3}){sub rich}-CaO-TiO{sub x}. The compositions of most inclusions were Liquid+Spinel or Liquid+TiSp phase. The (MgO-Al{sub 2}O{sub 3}){sub rich}-CaO inclusions that were not effectively modified by calcium treatment could be a source of clogging. By combining both the experimental results and thermodynamic calculation, (MgO-Al{sub 2}O{sub 3}){sub rich}-TiO{sub x} inclusions (TiSp inclusions) would be formed with the decreasing temperature in steel during continuous casting, which could be transported and accumulated on the inner wall of the SEN. Increasing the calcium content in the molten steel could decrease the formation of TiSp inclusions. During continuous casting, these inclusions accumulated on the inner wall of the SEN and gradually formed porous network of clog. Molten steel would gradually solidify in the porous space of the clog due to the increased residence time of the molten steel when flowing through it. The clogging gradually increased and caused the shutdown of continuous casting.

36 MATERIALS SCIENCE↗

SAM Two-Phase Flow Model Development and Applications for Operational Transients in Advanced Reactors

As advanced nuclear technologies continue to develop, the need for the flexible operation and generation of these advanced reactors becomes necessary to maximize economic potential. As large-scale experiments are not always feasible, modeling and simulations of advanced reactors play a crucial role in design optimization and analysis. The SAM (System Analysis Module) code developed at Argonne National Laboratory is a state-of-the-art system-level thermal-hydraulic code aimed at simulating advanced reactor systems. Recent code developments have implemented two-phase flow modeling using the homogeneous equilibrium model, and a new steam generator component has been developed to utilize the two-phase flow implementation. In addition to verification tests, a load-following simulation was performed to model a realistic load-following transient in a proposed integrated system consisting of a conceptual advanced reactor known as the Advanced Burner Test Reactor (ABTR) and thermal energy storage (TES) tanks. The integrated system model uses two large TES tanks designed for sodium and a model helical coil steam generator to simulate the operational load-following transient. The flow rates of the feedwater and secondary loops are regulated to meet a prescribed steam generator load consistent with the electricity demand over a 24-h period. In conclusion, the results found the ABTR system was able to maintain stable reactor conditions and primary- and secondary-side characteristics over the course of the load-following transient.

Advanced Burner Test Reactor (ABTR)↗

Forced Convection and Sedimentation Past a Flat Plate

The steady laminar flow of a well-mixed suspension of monodisperse solid spheres, convected steadily past a horizontal flat plate and sedimenting under the action of gravity, is examined. It is shown that, in the limit as Re approaches infinity and epsilon approaches 0, where Re is the bulk Reynolds number and epsilon is the ratio of the particle radius a to the characteristic length scale L, the analysis for determining the particle concentration profile has several aspects in common with that of obtaining the temperature profile in forced-convection heat transfer from a wall to a fluid stream moving at high Reynolds and Prandtl numbers. Specifically, it is found that the particle concentration remains uniform throughout the O(Re(exp -1/2)) thick Blasius boundary layer except for two O(epsilon(exp 2/3)) thin regions on either side of the plate, where the concentration profile becomes non-uniform owing to the presence of shear-induced particle diffusion which balances the particle flux due to convection and sedimentation. The system of equations within this concentration boundary layer admits a similarity solution near the leading edge of the plate, according to which the particle concentration along the top surface of the plate increases from its value in the free stream by an amount proportional to X(exp 5/6), with X measuring the distance along the plate, and decreases in a similar fashion along the underside. But, unlike the case of gravity settling on an inclined plate in the absence of a bulk flow at infinity considered earlier, here the concentration profile remains continuous everywhere. For values of X beyond the region near the leading edge, the particle concentration profile is obtained through the numerical solution of the relevant equations. It is found that, as predicted from the similarity solution, there exists a value of X at which the particle concentration along the top side of the plate attains its maximum value phi(sub m) and that, beyond this point, a stagnant sediment layer will form that grows steadily in time. This critical value of X is computed as a function of phi(sub s), the particle volume fraction in the free stream. In contrast, but again in conformity with the similarity solution, for values of X sufficiently far removed from the leading edge along the underside of the plate, a particle-free region is predicted to form adjacent to the plate. This model, with minor modifications, can be used to describe particle migration in other shear flows, as, for example, in the case of crossflow microfiltration.

Pelekasis, Nikolaos A.↗

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗

Fractional-Step, Finite-Volume Computation Of Flow

Method of solving Navier-Stokes equations of incompressible flow in general nonorthogonal curvilinear coordinates incorporates fractional-step and finite-volume approaches. Developed in continuing effort to achieve accuracy without need for excessive computation time in numerical simulation of time-dependent, three-dimensional flows bounded by surfaces of complicated shape.

Kwak, Dochan↗

Paralinear Oxidation of Silicon Nitride in a Water Vapor/Oxygen Environment

Three silicon nitride materials were exposed to dry oxygen flowing at 0.44 cm/s at temperatures between 1200 and 1400 C. Reaction kinetics were measured with a continuously recording microbalance. Parabolic kinetics were observed. When the same materials were exposed to a 50% H2O - 50% O2 gas mixture flowing at 4.4 cm/s, all three types exhibited paralinear kinetics. The material is oxidized by water vapor to form solid silica. The protective silica is in turn volatilized by water vapor to form primarily gaseous Si(OH)4. Nonlinear least squares analysis and a paralinear kinetic model were used to determine both parabolic and linear rate constants from the kinetic data. Volatilization of the protective silica scale can result in accelerated consumption of Si3N4. Recession rates under conditions more representative of actual combustors are compared to the furnace data.

Fox, Dennis S.↗

Local convection characteristics of inline arrangement of Kagome-shaped unit cells in a square duct

Transient liquid crystal thermography experiments have been conducted to determine detailed convective heat transfer coefficients at the endwalls of lattice-frame configurations based on Kagome-shaped unit cells. Kagome unit cells with porosity of 0.88 were arranged in an inline manner, where two such arrangements were studied. In the first arrangement, a total of ten unit cells were placed next to each other along the streamwise direction resulting in a continuous configuration. In the 2nd arrangement, alternate unit cells from the continuous configuration were dropped, resulting in a discrete configuration which featured a total of 5 Kagome unit cells. Due to the asymmetric nature of the strut connections within the unit cell, the convective heat transfer coefficient maps were determined for the two opposite walls where the struts meet the endwalls. Transient liquid crystal experiments were conducted for Reynolds number ranging between 10,000 and 30,000. Here, the study was focused on the developing nature of flow along the streamwise direction and the local convection characteristics for continuous and discrete arrangement of unit cells. For the continuous configuration, the Nusselt number ratios (N u /Nu 0 ) varied between 3.15–3.46 and 2.75–2.89 for the two walls A and B, respectively. For the discrete configuration, convective heat transfer coefficients varied between 2.56–2.71 and 2.09–2.40 for the two walls. Kagome unit cells have the potential to be fabricated through inexpensive manufacturing routes such as wire-woven method and these unit cells find their applications in the areas which require different heat transfer levels on opposite walls.

Kagome↗

The role of rifting in the generation of melt: Implications for the origin and evolution of the Lada Terra-Lavinia Planitia region of Venus

A rift over 6000 km in total length runs along the border of Lada Terra, a highland in the southern hemisphere of Venus, and Lavinia Planitia, a basin that has been interpreted as a site of early-stage mantle downwelling. Along the length of the rift are a number of volcanic centers of widely varying morphology and volcanic output. These include coronae, radially fractured domes, and large flow fields similar in scale to terrestrial flood basalts. We develop a model for the origin of extension related to passive rifting in response to stresses created by the adjacent downwelling. Volcanism and extension at other rifts on Venus, such as Devana Chasma, have been attributed to deep-seated mantle plume activity. In contrast, we interpret the origin of extension and volcanism along the Lada rift to be linked to upwelling and decompression melting of mantle material due to rifting and, possibly, to counterflow associated with downwelling. Extension occurred generally prior to the formation of volcanic centers and the eruption of large-scale flow fields, although most of the volcanic centers have been fractured by continued extension along the rift. Current debate over the formation of terrestrial flood basalts centers on the necessity of preexisting extension and stretched and thinned lithosphere to produce enhanced decompression melting within a large plume head or mantle thermal anomaly. Our studies of large-scale flow fields associated with the Lada rift and coronae on Venus indicate that extension is a prerequisite for the formation of the majority of large-scale flow units on Venus.

Magee, Kari P.↗

A three-dimensional axisymmetric calculation procedure for turbulent flows in a radial vaneless diffuser

An analytical model is proposed to calculate the three-dimensional axisymmetric turbulent flowfield in a radial vneless diffuser. The model assumes that the radial and tangential boundary layer profiles are approximated by power law profiles. Then, using the integrated radial and tangential momentum and continuity equations for the boundary layer and corresponding inviscid equations for the core flow, there results six ordinary differential equation in six unknowns which are easily solved using a Runge-Kutta technique. A model is also proposed for fully developed flow. The results using this technique were compared with the results from a three-dimensional viscous, axisymmetric duct code and with experimental data and good quantitative agreement was obtained.

Schumann, L. F.↗

Code Development for Unsteady Inlet Flows using Parallel Processing, Iced Airfoils

As a part of an effort for the development of an "Integrated Solution Process for Iced Airfoils" which combines a CFD (Computational Fluid Dynamics) code and an ice accretion code for accurate prediction of ice growth and performance degradation, a study on the effect of iced-geometry-smoothing was initiated. As a first step in this study, the degree of smoothing was defined by the number of control points for the given iced airfoil geometry. Then, reducing these number of control points in a systematical way provided various degrees of grid generation. This study will be continued by comparing CFD computed data such as pressure, lift, drag, flow separation, and wake flow patterns between the iced airfoils against any existing experimental data.

Chung, Joongkee↗

Time-dependent numerical simulation of hot ion outflow from the polar ionosphere

The time-dependent hydrodynamic model of Gombosi et al. (1985) was used to investigate the effects of bulk plasma heating on ion flows in the polar ionosphere. In the analysis, coupled time-dependent continuity, momentum, and energy equations of a two-ion (H/+/ and O/+/) quasi-neutral plasma were solved between 200 and 8000 km for polar wind conditions. In spite of its several limitations, the model can serve to indicate the gross behavior of polar region plasma flows. Two specific cases, ion heating and electron heating, were considered. The calculated temperature and flow profiles suggest that ion heating, rather than electron heating, generates the upwelling ion events, for which large transient O(+) outflows accompanied by elevated ion temperatures were observed.

Gombosi, T. I.↗