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 397 records · Page 22

Electrohydrodynamically coupled Kelvin-Helmholtz instability

An inviscid model of the surface dynamics at the interface between two parallel-flowing insulating dielectric fluids, subject to a tangential electric field with a weak gradient, is developed. The two flows are oppositely directed, resulting in an electrohydrodynamically coupled Kelvin-Hemholtz instability. A simple stability criterion is used successfully to correlate experimental data obtained with a test configuration.

Perry, M. P.↗

Fundamental Study of Material Flow in Friction Stir Welds

The presented research project consists of two major parts. First, the material flow in solid-state, friction stir, butt-welds as been investigated using a marker insert technique. Changes in material flow due to welding parameter as well as tool geometry variations have been examined for different materials. The method provides a semi-quantitative, three-dimensional view of the material transport in the welded zone. Second, a FSW process model has been developed. The fully coupled model is based on fluid mechanics; the solid-state material transport during welding is treated as a laminar, viscous flow of a non-Newtonian fluid past a rotating circular cylinder. The heat necessary for the material softening is generated by deformation of the material. As a first step, a two-dimensional model, which contains only the pin of the FSW tool, has been created to test the suitability of the modeling approach and to perform parametric studies of the boundary conditions. The material flow visualization experiments agree very well with the predicted flow field. Accordingly, material within the pin diameter is transported only in the rotation direction around the pin. Due to the simplifying assumptions inherent in the 2-D model, other experimental data such as forces on the pin, torque, and weld energy cannot be directly used for validation. However, the 2-D model predicts the same trends as shown in the experiments. The model also predicts a deviation from the "normal" material flow at certain combinations of welding parameters, suggesting a possible mechanism for the occurrence of some typical FSW defects. The next step has been the development of a three-dimensional process model. The simplified FSW tool has been designed as a flat shoulder rotating on the top of the workpiece and a rotating, cylindrical pin, which extends throughout the total height of the flow domain. The thermal boundary conditions at the tool and at the contact area to the backing plate have been varied to fit experimental data such as temperature profiles, torque and tool forces. General aspects of the experimentally visualized material flow pattern are confirmed by the 3-D model.

Reynolds, Anthony P.↗

A Reduced-Order Model For Zero-Mass Synthetic Jet Actuators

Accurate details of the general performance of fluid actuators is desirable over a range of flow conditions, within some predetermined error tolerance. Designers typically model actuators with different levels of fidelity depending on the acceptable level of error in each circumstance. Crude properties of the actuator (e.g., peak mass rate and frequency) may be sufficient for some designs, while detailed information is needed for other applications (e.g., multiple actuator interactions). This work attempts to address two primary objectives. The first objective is to develop a systematic methodology for approximating realistic 3-D fluid actuators, using quasi-1-D reduced-order models. Near full fidelity can be achieved with this approach at a fraction of the cost of full simulation and only a modest increase in cost relative to most actuator models used today. The second objective, which is a direct consequence of the first, is to determine the approximate magnitude of errors committed by actuator model approximations of various fidelities. This objective attempts to identify which model (ranging from simple orifice exit boundary conditions to full numerical simulations of the actuator) is appropriate for a given error tolerance.

Yamaleev, Nail K.↗

CFD Multiphysics Tool

The recent bold initiatives to expand the human presence in space require innovative approaches to the design of propulsion systems whose underlying technology is not yet mature. The space propulsion community has identified a number of candidate concepts. A short list includes solar sails, high-energy-density chemical propellants, electric and electromagnetic accelerators, solar-thermal and nuclear-thermal expanders. For each of these, the underlying physics are relatively well understood. One could easily cite authoritative texts, addressing both the governing equations, and practical solution methods for, e.g. electromagnetic fields, heat transfer, radiation, thermophysics, structural dynamics, particulate kinematics, nuclear energy, power conversion, and fluid dynamics. One could also easily cite scholarly works in which complete equation sets for any one of these physical processes have been accurately solved relative to complex engineered systems. The Advanced Concepts and Analysis Office (ACAO), Space Transportation Directorate, NASA Marshall Space Flight Center, has recently released the first alpha version of a set of computer utilities for performing the applicable physical analyses relative to candidate deep-space propulsion systems such as those listed above. PARSEC, Preliminary Analysis of Revolutionary in-Space Engineering Concepts, enables rapid iterative calculations using several physics tools developed in-house. A complete cycle of the entire tool set takes about twenty minutes. PARSEC is a level-zero/level-one design tool. For PARSEC s proof-of-concept, and preliminary design decision-making, assumptions that significantly simplify the governing equation sets are necessary. To proceed to level-two, one wishes to retain modeling of the underlying physics as close as practical to known applicable first principles. This report describes results of collaboration between ACAO, and Embry-Riddle Aeronautical University (ERAU), to begin building a set of level-two design tools for PARSEC. The "CFD Multiphysics Tool" will be the propulsive element of the tool set. The name acknowledges that space propulsion performance assessment is primarily a fluid mechanics problem. At the core of the CFD Multiphysics Tool is an open-source CFD code, HYP, under development at ERAU. ERAU is renowned for its undergraduate degree program in Aerospace Engineering the largest in the nation. The strength of the program is its applications-oriented curriculum, which culminates in one of three two-course Engineering Design sequences: Aerospace Propulsion, Spacecraft, or Aircraft. This same philosophy applies to the HYP Project, albeit with fluid physics modeling commensurate with graduate research. HYP s purpose, like the Multiphysics Tool s, is to enable calculations of real (three-dimensional; geometrically complex; intended for hardware development) applications of high speed and propulsive fluid flows.

Perrell, Eric R.↗

New cellular automaton model for magnetohydrodynamics

A new type of two-dimensional cellular automation method is introduced for computation of magnetohydrodynamic fluid systems. Particle population is described by a 36-component tensor referred to a hexagonal lattice. By appropriate choice of the coefficients that control the modified streaming algorithm and the definition of the macroscopic fields, it is possible to compute both Lorentz-force and magnetic-induction effects. The method is local in the microscopic space and therefore suited to massively parallel computations.

Chen, Hudong↗

Simulation and Modeling of Flow Generated by Gust Generator in a Wind Tunnel

A computational study of a gust field generated by a gust generator in a low-speed wind tunnel. The gust generator is designed for the University of Washington Aeronautical Laboratory (UWAL) Kirsten wind tunnel for a gust load alleviation (GLA) control experiment of a Common Research Model (CRM) flexible wing utilizing the Variable Camber Continuous Trailing Edge Flap (VCCTEF). The gust generator comprises four horizontal NACA 0015 gust generator vanes placed upstream of the test section. Computational fluid dynamics simulations using a two-dimensional (2D) Unsteady-Reynolds-Averaged-Navier-Stokes (URANS) with k-ω Shear Stress Transport (SST) turbulence model provide detailed time-resolved information about the generated flow by the gust generator under prescribed sinusoidal motion. The characteristics of the induced flow by the gust generator are analyzed. A gust propagation model of the gust field is investigated. An unsteady lift model is developed using a varying-fidelity approach which includes a 2D interference aerodynamic model of the combined gust generator-wing system. The computed integrated unsteady lift is compared to experimental data for validation of the unsteady lift model. Both the amplitude and transport delay are found to be accurately captured by the unsteady lift model.

Nhan, Nguyen↗

Scaled equation of state parameters for gases in the critical region

In the light of recent theoretical developments, the paper presents an accurate characterization of anomalous thermodynamic behavior of xenon, helium 4, helium 3, carbon dioxide, steam and oxygen in the critical region. This behavior is associated with long range fluctuations in the system and the physical properties depend primarily on a single variable, namely, the correlation length. A description of the thermodynamic behavior of fluids in terms of scaling laws is formulated, and the two successfully used scaled equations of state (NBS equation and Linear Model parametric equation) are compared. Methods for fitting both equations to experimental equation of state data are developed and formulated, and the optimum fit for each of the two scaled equations of the above gases are presented and the results are compared. By extending the experimental data for the above one-component fluids to partially miscible binary liquids, superfluid liquid helium, ferromagnets and solids exhibiting order-disorder transitions, the principle of universality is concluded. Finally by using this principle, the critical regions for nine additional fluids are described.

Sengers, J. M. H. L.↗

A solar-wind model including proton thermal anisotropy.

A model of the solar wind is presented which divides the interplanetary space into two regions. It assumes that the solar wind is one-fluid in an 'inner region' of radius less than 0.4 A, and two-fluid and collisionless in an 'outer region' of radius greater than 0.4 A. The second and third moment equations of the Vlasov equation together with conservation equations are used to govern the solar wind flow in the outer region. The model produces solutions for all macroscopic quantities as well as for the microscopic proton distribution function, both as functions of heliocentric distances up to 50 A. Calculated results at 1 A include the solar wind velocity, electron temperature, proton temperature, proton thermal anisotropy and magnetic field magnitude. They are in good agreement with observations.

Whang, Y. C.↗

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗

Rheological effects on friction in elastohydrodynamic lubrication

An analytical and experimental investigation is presented of the friction in a rolling and sliding elastohydrodynamic lubricated contact. The rheological behavior of the lubricant is described in terms of two viscoelastic models. These models represent the separate effects of non-Newtonian behavior and the transient response of the fluid. A unified description of the non-Newtonian shear rate dependence of the viscosity is presented as a new hyperbolic liquid model. The transient response of viscosity, following the rapid pressure rise encountered in the contact, is described by a compressional viscoelastic model of the volume response of a liquid to an applied pressure step. The resulting momentum and energy equations are solved by an iterative numerical technique, and a friction coefficient is calculated. The experimental study was performed, with two synthetic paraffinic lubricants, to verify the friction predictions of the analysis. The values of friction coefficient from theory and experiment are in close agreement.

Trachman, E. G.↗

A numerical study of the effects of wind tunnel wall proximity on an airfoil model

A procedure was developed for modeling wind tunnel flows using computational fluid dynamics. Using this method, a numerical study was undertaken to explore the effects of solid wind tunnel wall proximity and Reynolds number on a two-dimensional airfoil model at low speed. Wind tunnel walls are located at varying wind tunnel height to airfoil chord ratios and the results are compared with freestream flow in the absence of wind tunnel walls. Discrepancies between the constrained and unconstrained flows can be attributed to the presence of the walls. Results are for a Mach Number of 0.25 at angles of attack through stall. A typical wind tunnel Reynolds number of 1,200,000 and full-scale flight Reynolds number of 6,000,000 were investigated. At this low Mach number, wind tunnel wall corrections to Mach number and angle of attack are supported. Reynolds number effects are seen to be a consideration in wind tunnel testing and wall interference correction methods. An unstructured grid Navier-Stokes code is used with a Baldwin-Lomax turbulence model. The numerical method is described since unstructured flow solvers present several difficulties and fundamental differences from structured grid codes, especially in the area of turbulence modeling and grid generation.

Potsdam, Mark↗

Accretion of rotating fluids by barytropes - Numerical results for white-dwarf models

Numerical sequences of rotating axisymmetric nonmagnetic equilibrium models are constructed which represent the evolution of a barytropic star as it accretes material from a rotating medium. Two accretion geometries are considered - one approximating accretion from a rotating cloud and the other, accretion from a Keplerian disk. It is assumed that some process, such as Ekman spin-up or nonequilibrium oscillations, maintains nearly constant angular velocity along cylinders about the rotation axis. Transport of angular momentum in the cylindrically radial direction by viscosity is included. Fluid instabilities and other physical processes leading to enhancement of this transport are discussed. Particular application is made to zero-temperature white-dwarf models, using the degenerate electron equation of state. An initially nonrotating 0.566-solar-mass white dwarf is followed during the accretion of more than one solar mass of material. Applications to degenerate stellar cores, to mass-transfer binary systems containing white dwarfs, such as novae and dwarf novae, to Type I supernovae, and to galactic X-ray sources are considered.

Durisen, R. H.↗

Analysis and testing of propellant feed system priming process

This paper presents the analytical and the experimental results pertaining to the priming process of the propellant feed system with initial line pressures starting from 0 psia and greater. The analytical methods employ the method of characteristics to solve for the one-dimensional liquid transients in the liquid-full segments and the lumped inertia technique to model the dynamics of the partially filled segments or the two-phase segments. The advantages of these methods are (1) fluid compressibility and piping flexibility are accounted for in the solution; and (2) the characteristics method can be used in the solution of a complex system. The analytical results were compared with test results. Excellent correlation was obtained between predictions and test results verifying the essential aspects of the analytical modeling techniques and providing the foundation for analyzing a complicated network system.

Lin, T. Y.↗

CAS2D: FORTRAN program for nonrotating blade-to-blade, steady, potential transonic cascade flows

An exact, full-potential-equation (FPE) model for the steady, irrotational, homentropic and homoenergetic flow of a compressible, homocompositional, inviscid fluid through two dimensional planar cascades of airfoils was derived, together with its appropriate boundary conditions. A computer program, CAS2D, was developed that numerically solves an artificially time-dependent form of the actual FPE. The governing equation was discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field was discretized by providing a boundary-fitted, nonuniform computational mesh. The mesh was generated by using a sequence of conforming mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the FPE was solved iteratively by using successive line overrelaxation. The possible isentropic shocks were correctly captured by adding explicitly an artificial viscosity in a conservative form. In addition, a three-level consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two dimensional cascade flows.

Dulikravich, D. S.↗

Instability of baroclinic flows with horizontal shear along topography

The stability of baroclinic flows with horizontal shear over sloping topography is analyzed with special emphasis on the structure and energetics of the unstable perturbations. The study is conducted by using a linearized two-layer quasi-geostrophic channel model for different topography profiles and distributions of the basic velocity field. Interactions between the two fluid layers and the energy conversions by the unstable perturbations are described. It is found that topography sloping as (opposed to) the fluid interface contributes to enhance the perturbation amplitude in the upper (lower) layer relative to the lower (upper) layer. The results for bottom topography with differing characteristics across the flow indicate pronounced localized effects on the energy conversions over the slopes and the meridional scale of the perturbations in the lower layer.

Mechoso, C. R.↗

Lithospheric flexure at fracture zones

Studies attempting to demonstrate that lithospheric flexure occurs across the Pioneer and Mendocino fracture zones, and that the flexural topography is a topographic expression at these fracture zones, are presented. The flexure is modelled and compared with predicted depths with five bathymetric profiles which cross the two fracture zones at different ages. The model uses a thin elastic plate overlying an incompressible fluid half-space, and incorporates a temperature-dependent effective elastic thickness. Several conclusions were derived from this study. First, it is found that no significant slip on the fossil fault planes of the Mendocino and Pioneer fracture zones exists. In addition, the flexural amplitude is determined to increase with age. Finally, it is concluded that there is elastic coupling between the Mendocino and Pioneer fracture zones since the separation is less than a flexural wavelength.

Sandwell, D.↗

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗

A general method for unsteady stagnation region heat transfer and results for model turbine flows

Recent experiments suggest that the heat-transfer characteristics of stator blades are influenced by the frequency of passing of upstream rotor blades. The calculation of these effects requires that the movement of the stagnation point with variations in freestream velocity is properly represented together with the possible effects of turbulence characteristics on the thin leading-edge boundary layer. A procedure to permit the achievement of these purposes is described for laminar flows in this paper together with results of its application to two model problems which demonstrate its abilities and quantify the influence of wake characteristics on fluid-dynamic and heat-transfer properties of the flow and their effects on surface heat transfer.

Cebeci, Tuncer↗