Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

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 199 records · Page 11

kL-based Linear and Nonlinear Two-Equation Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) two-equation turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important differences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The present report documents the formulation of two model levels of turbulence models as implemented in the computational fluid dynamics FUN3D code. The levels are the two-equation linear k-kL and the two-equation algebraic Reynolds stress model (ARSM). Free shear, separated and corner flow cases are documented and compared with experimental, and other turbulence model data. The results show generally very good comparisons with experimental data. The results from this formulation are similar or better than results using the SST two-equation turbulence model. ARSM shows great promise with similar level of computational resources as basic two-equation turbulence models.

Abdol-Hamid, K. S.↗

Analysis and Modeling of a Two-Phase Jet Pump of a Thermal Management System for Aerospace Applications

Jet pumps are devices capable of pumping fluids to a higher pressure by inducing the motion of a secondary fluid employing a high speed primary fluid. The main components of a jet pump are a primary nozzle, secondary fluid injectors, a mixing chamber, a throat, and a diffuser. The work described in this paper models the flow of a two-phase primary fluid inducing a secondary liquid (saturated or subcooled) injected into the jet pump mixing chamber. The model is capable of accounting for phase transformations due to compression, expansion, and mixing. The model is also capable of incorporating the effects of the temperature and pressure dependency in the analysis. The approach adopted utilizes an isentropic constant pressure mixing in the mixing chamber and at times employs iterative techniques to determine the flow conditions in the different parts of the jet pump.

Sherif, S.A.↗

Real-Time Measurements of Aft Dome Insulation Erosion on Space Shuttle Reusable Solid Rocket Motor

Real-time erosion of aft dome internal insulation was measured with internal instrumentation on a static test of a lengthened version of the Space Shuffle Reusable Solid Rocket Motor (RSRM). This effort marks the first time that real-time aft dome insulation erosion (Le., erosion due to the combined effects of thermochemical ablation and mechanical abrasion) was measured in this kind of large motor static test [designated as Engineering Test Motor number 3 (ETM3)I. This paper presents data plots of the erosion depth versus time. The data indicates general erosion versus time behavior that is in contrast to what would be expected from earlier analyses. Engineers have long known that the thermal environment in the aft dome is severe and that the resulting aft dome insulation erosion is significant. Models of aft dome erosion involve a two-step process of computational fluid dynamics (CFD) modeling and material ablation modeling. This modeling effort is complex. The time- dependent effects are difficult to verify with only prefire and postfire insulation measurements. Nozzle vectoring, slag accumulation, and changing boundary conditions will affect the time dependence of aft dome erosion. Further study of this data and continued measurements on future motors will increase our understanding of the aft dome flow and erosion environment.

McWhorter, Bruce↗

Advancements in engineering turbulence modeling

Some new developments in two-equation models and second order closure models are presented. Two-equation models (k-epsilon models) have been widely used in computational fluid dynamics (CFD) for engineering problems. Most of low-Reynolds number two-equation models contain some wall-distance damping functions to account for the effect of wall on turbulence. However, this often causes the confusion and difficulties in computing flows with complex geometry and also needs an ad hoc treatment near the separation and reattachment points. A set of modified two-equation models is proposed to remove the aforementioned shortcomings. The calculations using various two-equation models are compared with direct numerical simulations of channel flow and flat boundary layers. Development of a second order closure model is also discussed with emphasis on the modeling of pressure related correlation terms and dissipation rates in the second moment equations. All the existing models poorly predict the normal stresses near the wall and fail to predict the 3-D effect of mean flow on the turbulence (e.g. decrease in the shear stress caused by the cross flow in the boundary layer). The newly developed second order near-wall turbulence model is described and is capable of capturing the near-wall behavior of turbulence as well as the effect of 3-D mean flow on the turbulence.

Shih, T.-H.↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗

Progress in turbulence modeling for complex flow fields including effects of compressibility

Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.

Wilcox, D. C.↗

The AFDD International Dynamic Stall Workshop on Correlation of Dynamic Stall Models with 3-D Dynamic Stall Data

A variety of empirical and computational fluid dynamics two-dimensional (2-D) dynamic stall models were compared to recently obtained three-dimensional (3-D) dynamic stall data in a workshop on modeling of 3-D dynamic stall of an unswept, rectangular wing, of aspect ratio 10. Dynamic stall test data both below and above the static stall angle-of-attack were supplied to the participants, along with a 'blind' case where only the test conditions were supplied in advance, with results being compared to experimental data at the workshop itself. Detailed graphical comparisons are presented in the report, which also includes discussion of the methods and the results. The primary conclusion of the workshop was that the 3-D effects of dynamic stall on the oscillating wing studied in the workshop can be reasonably reproduced by existing semi-empirical models once 2-D dynamic stall data have been obtained. The participants also emphasized the need for improved quantification of 2-D dynamic stall.

Tan, C. M.↗

Computing Incompressible Flows With Free Surfaces

RIPPLE computer program models transient, two-dimensional flows of incompressible fluids with surface tension on free surfaces of general shape. Surface tension modeled as volume force derived from continuum-surface-force model, giving RIPPLE both robustness and accuracy in modeling surface-tension effects at free surface. Also models wall adhesion effects. Written in FORTRAN 77.

Kothe, D.↗

The nonlinear dynamics of a spacecraft coupled to the vibration of a contained fluid

The dynamics of a linear spacecraft mode coupled to a nonlinear low gravity slosh of a fluid in a cylindrical tank is investigated. Coupled, nonlinear equations of motion for the fluid-spacecraft dynamics are derived through an assumed mode Lagrangian method. Unlike linear fluid slosh models, this nonlinear slosh model retains two fundamental slosh modes and three secondary modes. An approximate perturbation solution of the equations of motion indicates that the nonlinear coupled system response involves fluid-spacecraft modal resonances not predicted by either a linear, or a nonlinear, uncoupled slosh analysis. Experimental results substantiate the analytical predictions.

Peterson, Lee D.↗

Observations of Keyhole Porosity and Comparisons to Analytical Models for Ti-6Al-4V Powder Bed Fusion

Keyhole porosity defects are a common concern in powder bed fusion (PBF) processing. Keyhole porosity prediction models have generally fallen into two categories – high fidelity computational fluid dynamics simulations and low fidelity analytical or empirical models based on keyhole vibration dynamics. The low computational cost of low fidelity models allows them to better approach part-scale predictions. However, studies on low fidelity modeling techniques are limited by the lack of experimental data to assess the validity of the calibration over a wide range of processing conditions. This work extracts and quantifies keyhole porosity across 14 laser velocities, two laser powers, and six+ repetitions for a total of 176 independent trials. The measured porosity data are compared to low fidelity keyhole models to assess their success in predictive porosity occurrence. This work impacts the field by providing independent validation of keyhole porosity models for PBF for use in part-scale defect prediction.

powder bed fusion↗

Structure of ion acoustic solitons and shock waves in a two-component plasma.

Time-independent solitary waves and shocks are investigated in a two-component plasma using a fluid model and kinetic theory. It is found that very small concentrations of a light ion can drastically alter the structure, changing the potential maximum by an order of magnitude. For a fixed Mach number, a critical density ratio of light to heavy ions is found at which the potential maximum changes discontinuously from a value large enough to reflect the light ions to one which allows them to traverse the shock front and enter the downstream flow. The downstream oscillatory structure normally seen in a shock is completely quenched by dissipation due to light ion reflection at concentrations of 3-8% He in an Ar plasma for typical electron to ion temperature ratios and Mach number values.

White, R. B.↗

Geophysical fluid flow model experiments in spherical geometry

An experimental and theoretical program was undertaken to assist in the design of geophysical fluid flow model experiments for Spacelab. Two new nonintrusive measurement techniques were developed. A theoretical calculation was carried out to guide the design of a proposed atmospheric general circulation model experiment.

Fowlis, W. W.↗

Particle dispersion in confined turbulent swirling flows

This paper reports a numerical investigation of confined swirling flows of gas-particle mixtures. A recently developed two-fluid, multiple-scale mixing model is applied to study the influence of particles on the intensity of the reverse flows of the gas phase and the effects of swirl on the particle dispersion in an annular expanding chamber under isothermal condition. The calculations were made for different swirl strength of injection of the annular jet into the mixing chamber. Results agree well qualitatively with experimental information available. It is also found that the calculated flow fields depend heavily on the prescription of the inlet flow conditions.

Chen, C. P.↗

A mapped finite difference study of noise propagation in nonuniform ducts with mean flow

The primary objective of this work is to study noise propagation in acoustically lined variable area ducts with mean fluid flow. The method of study is numerical in nature and involves a body-fitted grid mapping procedure in conjunction with a factored-implicit finite difference technique. The mean fluid flow model used is two-dimensional, inviscid, irrotational, incompressible, and nonheat conducting. Fully-coupled solutions of the linearized gasdynamic equations are obtained for both positive and negative Mach numbers as well as for hard and soft wall conditions. The factored-implicit finite difference technique used did give rise to short wavelength perturbations, but these were dampened by the introduction of higher order artificial dissipation terms into the scheme. Results compared favorably with available numerical and experimental data.

Raad, Peter E.↗

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗

Modeling of an Atmospheric-Boundary-Layer Profile in Support of Experiments in the NASA Langley Transonic Dynamics Tunnel

This paper presents the results of the Computational Fluid Dynamics (CFD) analyses to model two atmospheric-boundary-layer (ABL) profiles inside the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The CFD models include tunnel walls and all additional hardware needed to simulate the ABL profiles. This hardware includes Irwin spires and floor-mounted roughness elements. The numerical simulations show that the application of higher-fidelity numerical methods is necessary to compute boundary-layer and turbulenceintensity profiles that match experimental data. The ABL profiles are computed both inside the numerical model of the TDT and in a classical free-air model. Both Unsteady Reynolds Averaged Navier-Stokes (URANS) with Spalart-Allmaras (SA) turbulence model and Modified Delayed Detached Eddy (MDDES) simulation methods are used. The results show that the MDDES-simulation results match the experimental data very well while URANS-SA does not. The results also show that the wind-tunnel walls have a significant effect on ABL prediction.

Chwalowski, Pawel↗

Noise produced by turbulent flow into a rotor: Theory manual for atmospheric turbulence prediction and mean flow and turbulence contraction prediction

Prediction of helicopter main rotor noise due to ingestion of atmospheric turbulence was analyzed. The analysis combines several different models that describe the fluid mechanics of the turbulence and the ingestion process. Two models, atmospheric turbulence, and mean flow and turbulence contraction were covered. The third model, covered in a separate report, describes the rotor acoustic mode. The method incorporates the atmospheric turbulence model and a rapid distortion turbulence contraction description to determine the statistics of the anisotropic turbulence at the rotor plane. The analytical basis for a module was provided which was incorporated in NASA's ROTONET helicopter noise prediction program. The mean flow and turbulence statistics associated with the atmospheric boundary layer were modeled including effects of atmospheric stability length, wind speed, and altitude. The turbulence distortion process is modeled as a deformation of vortex filaments (which represent the turbulence field) by a mean flow field due to the rotor inflow.

Simonich, J. C.↗

Computer-aided-engineering system for modeling and analysis of ECLSS integration testing

The accurate modeling and analysis of two-phase fluid networks found in environmental control and life support systems is presently undertaken by computer-aided engineering (CAE) techniques whose generalized fluid dynamics package can solve arbitrary flow networks. The CAE system for integrated test bed modeling and analysis will also furnish interfaces and subsystem/test-article mathematical models. Three-dimensional diagrams of the test bed are generated by the system after performing the requisite simulation and analysis.

Sepahban, Sonbol↗