Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Improvers”

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 775 records · Page 43

A comprehensive modeling of falling film evaporators subject to vapor flow, pass arrangements, and refrigerants

Improving the heat transfer performance of falling film evaporators is a crucial step for improving the energy efficiency of the heat pump or refrigeration systems. This study conducts a numerical investigation for the practical-scale falling film evaporator based on the epsilon–NTU method with an updated heat transfer correlation. The algorithm was validated with lab-scale and real-scale falling film evaporator experimental results, and the prediction reaches a mean absolute deviation of 12.5%. Here, the parametric study encompasses eight refrigerants: R-134a, R-410A, R-600a, R-717, R-1270, R-152a, R-1234yf, and R-1234ze(E). The results indicate that vaporization enthalpy of refrigerant is a key property in selecting an appropriate working fluid because it helps minimize severe heat transfer degradation caused by dry-out. Additionally, the vapor-flow–induced heat transfer degradation can be predicted using the critical Weber number. Furthermore, the trade-off between extending the tube length and increasing the number of tubes for heat transfer improvement is discussed. Finally, different two-pass arrangements show deviations of less than 4 %.

Dry-out↗

Investigation of the Impact of Flow on MHD Perturbations with the NIMROD Code (Final Report)

This work will enable improved understanding of tokamak edge flows and impurity-species dynamics and their impact on MHD stability as it relates to ELMs, RMPs and QH-mode by leveraging computational developments of the NIMROD code. Experimentally, it is established that the tokamak flows and impurities species can have a substantial effect on MHD stability. This is particularly true for the edge plasma where, for example, large flow shear is correlated with the occurrence of the Quiescent H-mode (QH- mode) state as opposed to operation with edge localized modes (ELMs). Understanding plasma flows is challenging because of the many physical effects that come into play: particle orbits loss, neutrals dynamics, and interaction between multiple ion species with decoupled motion. We will develop a model that incorporates multiple ion species into the MHD framework that includes self-consistent magnetic-field evolution. With this model we will understand how the presence of impurities impacts the tokamak edge-pedestal flows. Finally, we will quantify the impact of the new multiple species modeling on 3D QH-mode simulations. This work will have broader impacts on all MHD studies when multiple species are present (as is the case in modern tokamaks), as well as contributing to the edge and transport communities who also are interested in the studies of multiple-ion-species flows in the edge pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesized airfoil data method for prediction of dynamic stall and unsteady airloads

A detailed analysis of dynamic stall experiments has led to a set of relatively compact analytical expressions, called synthesized unsteady airfoil data, which accurately describe in the time-domain the unsteady aerodynamic characteristics of stalled airfoils. An analytical research program was conducted to expand and improve this synthesized unsteady airfoil data method using additional available sets of unsteady airfoil data. The primary objectives were to reduce these data to synthesized form for use in rotor airload prediction analyses and to generalize the results. Unsteady drag data were synthesized which provided the basis for successful expansion of the formulation to include computation of the unsteady pressure drag of airfoils and rotor blades. Also, an improved prediction model for airfoil flow reattachment was incorporated in the method. Application of this improved unsteady aerodynamics model has resulted in an improved correlation between analytic predictions and measured full scale helicopter blade loads and stress data.

Gangwani, S. T.↗

Higher order accuracy for upwind methods by using the compatibility equations

A new algorithm has been developed for obtaining higher order accuracy in upwind schemes for the Euler and Navier-Stokes equations. In this method, the compatibility relations for the Euler equations are used to construct formulas for the higher order interpolates. Using these formulas, computed results are obtained for steady, inviscid flow through a nozzle and also steady, inviscid and viscous flow over an airfoil. These results are compared with results obtained by using a minmod limiter function. The comparisons show improved accuracy throughout the flow field with the use of the compatibility equations.

Goorjian, Peter M.↗

FUN3D Grid Refinement and Adaptation Studies for the Ares Launch Vehicle

This paper presents grid refinement and adaptation studies performed in conjunction with computational aeroelastic analyses of the Ares crew launch vehicle (CLV). The unstructured grids used in this analysis were created with GridTool and VGRID while the adaptation was performed using the Computational Fluid Dynamic (CFD) code FUN3D with a feature based adaptation software tool. GridTool was developed by ViGYAN, Inc. while the last three software suites were developed by NASA Langley Research Center. The feature based adaptation software used here operates by aligning control volumes with shock and Mach line structures and by refining/de-refining where necessary. It does not redistribute node points on the surface. This paper assesses the sensitivity of the complex flow field about a launch vehicle to grid refinement. It also assesses the potential of feature based grid adaptation to improve the accuracy of CFD analysis for a complex launch vehicle configuration. The feature based adaptation shows the potential to improve the resolution of shocks and shear layers. Further development of the capability to adapt the boundary layer and surface grids of a tetrahedral grid is required for significant improvements in modeling the flow field.

Bartels, Robert E.↗

A dissipation-based method for improving the accuracy of computational fluid dynamics simulations of high level non-Newtonian wastes

Here this paper investigates a new viscosity model for laminar, transitional, and turbulent pipe flow simulations of a non-Newtonian slurry using the Star-CCM+ software. Herschel-Bulkley (HB) rheology equation was used in its original and modified forms to model the fluid viscosity. The realizable k-ε model was used for modeling of turbulence. It was observed that simulations using the HB model predicted inaccurate profiles of U (mean velocity) and U + (mean velocity in wall units). Therefore, we developed alternative values of viscosity obtained directly from non-linearity of stress–strain variation (α) in each iteration. The developed method, “Alpha method”, applied changes in places where dissipation rate of turbulent kinetic energy occurred and was designed based on local or global dissipation rate values. Inverse values of α were used to test effects of more viscosity reduction for all flow regimes. In prediction of the mean velocities, significant improvements were obtained through global adjustment of viscosity in inverse form during simulation of the laminar and turbulent flow cases. For the transitional flow, obtaining improvements in prediction of U and U + from the same model version faced challenges. While the global adjustments in direct form were found significantly effective in improving the prediction of U, they failed to predict the U + accurately. For this purpose, Local viscosity adjustments were found significantly effective.

42 ENGINEERING↗

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗

Distributed Aerodynamic Sensing and Processing Toolbox

A Distributed Aerodynamic Sensing and Processing (DASP) toolbox was designed and fabricated for flight test applications with an Aerostructures Test Wing (ATW) mounted under the fuselage of an F-15B on the Flight Test Fixture (FTF). DASP monitors and processes the aerodynamics with the structural dynamics using nonintrusive, surface-mounted, hot-film sensing. This aerodynamic measurement tool benefits programs devoted to static/dynamic load alleviation, body freedom flutter suppression, buffet control, improvement of aerodynamic efficiency through cruise control, supersonic wave drag reduction through shock control, etc. This DASP toolbox measures local and global unsteady aerodynamic load distribution with distributed sensing. It determines correlation between aerodynamic observables (aero forces) and structural dynamics, and allows control authority increase through aeroelastic shaping and active flow control. It offers improvements in flutter suppression and, in particular, body freedom flutter suppression, as well as aerodynamic performance of wings for increased range/endurance of manned/ unmanned flight vehicles. Other improvements include inlet performance with closed-loop active flow control, and development and validation of advanced analytical and computational tools for unsteady aerodynamics.

Brenner, Martin↗

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗

Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold

Understanding magnetohydrodynamic (MHD) phenomena associated with the flow of electrically conducting fluids in complex geometry ducts subject to a strong magnetic field is required to effectively design liquid metal (LM) blankets for fusion reactors. Specifically, accurately predicting the 3D MHD pressure drop and flow distribution is important. To explore these topics, we simulate a LM MHD flow through an electrically non-conducting prototypic manifold for a wide range of flow and geometry parameters using a 3D MHD solver, HyPerComp incompressible MHD solver for arbitrary geometry. The reference manifold geometry consists of a rectangular feeding duct which suddenly expands such that the duct thickness in the magnetic field direction abruptly increases by a factor r exp . Downstream of the sudden expansion, the LM is distributed into several parallel channels. As a first step in qualifying the flow, a magnitude of the curl of the induced Lorentz force was used to distinguish between inviscid, irrotational core flows and boundary and internal shear layers where inertia and/or viscous forces are important. Scaling laws have been obtained which characterize the 3D MHD pressure drop and flow balancing as a function of the flow parameters and the manifold geometry. Associated Hartmann and Reynolds numbers in the computations were ~10 3 and ~10 1 -10 3 , respectively, while r exp was varied from 4 to 12. An accurate model for the pressure drop was developed for the first time for inertial-electromagnetic and viscous-electromagnetic regimes based on 96 computed cases. Analysis reflects that flow balance can be improved by lengthening the distance between the manifold inlet and the entrances of the parallel channels by utilizing the effect of flow transitioning to a quasi-two-dimensional state in the expansion region of the manifold.

42 ENGINEERING↗

Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Separation of dilute plutonium(IV) from nitric acid solutions via Reillex HPQ anion exchange resin

Anion exchange is a common method for the separation of plutonium (Pu) in aqueous systems, providing high recovery and decontamination from impurities. Reillex HPQ anion exchange resin has been used around the Department of Energy complex in this application for decades but has typically been used for Pu feed concentrations greater than 1 g/L. It is desirable to recover Pu at more dilute concentrations from a dissolved used fuel or target element, for example. In this work, the performance of Reillex HPQ anion exchange resin was characterized with feed solutions of ca. 0.24 g/L Pu. A significant reduction in the performance of the resin was realized with this dilute concentration relative to a feed solution of 2.35 g/L Pu with similar feed flow rates and the same resin column, likely due to the increased driving force for mass transfer with a higher feed concentration. The maximum capacity of Reillex HPQ resin for Pu can be improved by reducing the flow rate and increasing the feed Pu concentration. Furthermore, the breakthrough point of Pu can be improved in the same manner and by increasing the resin bed height. Diffusion of Pu into the resin was illustrated to be the rate limiting step.

Plutonium separations↗

Detection of internal structure by scattered light intensity: Application to kidney cell sorting

Scattered light measurements in flow cytometry were sucessfully used to distinguish cells on the basis of differing morphology and internal structure. Differences in scattered light patterns due to changes in internal structure would be expected to occur at large scattering angles. Practically, the results of these calculations suggest that in experimental situations an array of detectors would be useful. Although in general the detection of the scattered light intensity at several intervals within the 10 to 60 region would be sufficient, there are many examples where increased sensitivity could be acheived at other angles. The ability to measure at many different angular intervals would allow the experimenter to empirically select the optimum intervals for the varying conditions of cell size, N/C ratio, granule size and internal structure from sample to sample. The feasibility of making scattered light measurements at many different intervals in flow cytometry was demonstrated. The implementation of simplified versions of these techniques in conjunction with independant measurements of cell size could potentially improve the usefulness of flow cytometry in the study of the internal structure of cells.

Goolsby, C. L.↗

Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation

This report gives an overview of physical modeling enhancements to the Wind-US flow solver which were made to improve the capabilities for simulation of hypersonic flows and the reliability of computations to complement hypersonic testing. The improvements include advanced turbulence models, a bypass transition model, a conjugate (or closely coupled to vehicle structure) conduction-convection heat transfer capability, and an upgraded high-speed combustion solver. A Mach 5 shock-wave boundary layer interaction problem is used to investigate the benefits of k- s and k-w based explicit algebraic stress turbulence models relative to linear two-equation models. The bypass transition model is validated using data from experiments for incompressible boundary layers and a Mach 7.9 cone flow. The conjugate heat transfer method is validated for a test case involving reacting H2-O2 rocket exhaust over cooled calorimeter panels. A dual-mode scramjet configuration is investigated using both a simplified 1-step kinetics mechanism and an 8-step mechanism. Additionally, variations in the turbulent Prandtl and Schmidt numbers are considered for this scramjet configuration.

Georgiadis, Nicholas J.↗

Effects of turbulence model selection on the prediction of complex aerodynamic flows

Numerical simulations of viscous transonic flow over a circular-arc airfoil and in a diffuser are described. The simulations are made with a new computer program designed to serve as a tool in the development of improved turbulence models for complex flows. The program incorporates zero-, one-, and two-equation eddy viscosity models and includes a variety of subsonic and supersonic boundary conditions. The airfoil flow contains a shock-separated boundary-layer interaction that has resisted previous attempts at simulation. The diffuser flow also contains a shock-boundary-layer interaction, which has not been simulated previously. Calculations using standard turbulence models, developed originally for incompressible unseparated flows, are described. Results indicate that although there are interesting differences in predictions between the various models, none of them predict the flows accurately. Suggestions for improved turbulence models are discussed.

Coakley, T. J.↗

Metric-discontinuous zonal grid calculations using the Osher scheme

Computations on zonal grids - in particular, grids with metric discontinuities resulting from the interspersion of highly clustered regions with coarse regions - are possible using a fully conservative form of the Osher upwind scheme. These zonal grids can result from an abrupt clustering of points near solution discontinuities or near other flow features that require improved resolution. The zonal approach is shown to capture shocks with almost 'shock-fitting' quality but with minimal effort. Results for inviscid flow, including quasi-one-dimensional nozzle flow, supersonic flow over a cylinder, and blast-wave diffraction by a ramp, are presented. These calculations demonstrate the powerful capabilities of the Osher scheme used in conjunction with zonal grids in simulating flow fields with complex shock patterns.

Rai, M. M.↗