Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow simulation”

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 685 records · Page 38

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗

Laboratory formation of a simulated comet

A laboratory simulation experiment to study the interaction between a cometary plasma and the solar wind has been performed using the UCR-T1 space simulation facility (dia. 1.3 m, length 11 m). Intense plasma flow simulating the solar wind interacts with another light-emitting plasma composed of Ba, Sr, and/or C by a plasma emitter which simulates a cometary coma. The purposes of this experiment are to investigate how the solar wind parameters contribute to the formation of the cometary ion tail and to determine the magnetic field structure of a comet. In order to estimate the solar wind parameters by ground-based observations of actual comets, knowledges of such relationships are essential. The experimental results show that the interplanetary magnetic field of the solar wind is very important in forming the cometary tail.

Minami, S.↗

Localized Ignition And Subsequent Flame Spread Over Solid Fuels In Microgravity

Localized ignition is initiated by an external radiant source at the middle of a thin solid sheet under external slow flow, simulating fire initiation in a spacecraft with a slow ventilation flow. Ignition behavior, subsequent transition simultaneously to upstream and downstream flame spread, and flame growth behavior are studied theoretically and experimentally. There are two transition stages in this study; one is the first transition from the onset of the ignition to form an initial anchored flame close to the sample surface, near the ignited area. The second transition is the flame growth stage from the anchored flame to a steady fire spread state (i.e. no change in flame size or in heat release rate) or a quasi-steady state, if either exists. Observations of experimental spot ignition characteristics and of the second transition over a thermally thin paper were made to determine the effects of external flow velocity. Both transitions have been studied theoretically to determine the effects of the confinement by a relatively small test chamber, of the ignition configuration (ignition across the sample width vs spot ignition), and of the external flow velocity on the two transitions over a thermally thin paper. This study is currently extending to two new areas; one is to include a thermoplastic sample such poly(methymethacrylate), PMMA, and the other is to determine the effects of sample thickness on the transitions. The recent results of these new studies on the first transition are briefly reported.

Kashiwagi, T.↗

Movement and collision of Lagrangian particles in hydro-turbine intakes: a case study

Studies of the stress/survival of migratory fish during downstream passage through operating hydro-turbines are normally conducted to determine the fish-friendliness of units. One field approach consisting of recording extreme hydraulics with autonomous sensors is largely sensitive to the conditions of sensor release and the initial trajectories at the turbine intake. This study applies a modelling strategy based on flow simulations using computational fluid dynamics and Lagrangian particle tracking to represent the travel of live fish and autonomous sensor devices through hydro-turbine intakes. For the flow field calculation, the simulations were conducted with both a time-averaging turbulence model and an eddy-resolving technique. For the particle tracking calculation, different modelling assumptions for turbulence forcing, mass formulation, buoyancy, and release condition were tested. The modelling assumptions are evaluated with respect to data sets collected using a laboratory physical model and an autonomous sensor device deployed at Ice Harbor Dam (Snake River, State of Washington, U.S.A.) at the same discharge and release point as in the present computer simulations. We found an acceptable agreement between the simulated results and observed data and discuss relevant features of Lagrangian particle movement that are critical in turbine design and in the experimental design of field studies.

CFD↗

Methods Of Simulation Of Incompressible Flow

Report discusses methods for numerical simulation of Navier-Stokes equations of viscous, incompressible flow, with emphasis on pseudocompressibility method. Advanced computational techniques reviewed contribute to improved designs by rapidly providing detailed information about conditions at all points in flow fields of proposed designs.

Dwak, D.↗

The Dynamics of Miscible Interfaces: Simulations

The goal of this experimental/computational investigation (joint with Prof Maxworthy at USC) has been to study the dynamics of miscible interfaces, both from a scientific and a practical point of view, and to prepare a related experiment to be flown on the International Space Station. In order to address these effects, we have focused experimental and computational investigations on miscible displacements in cylindrical capillary tubes, as well as in Hele-Shaw cells. Regarding the flow in a capillary tube, the question was addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations and experiments for miscible flows in capillary tubes. An estimate of the vorticity and streamfunction fields induced by the Kortewegs stresses was derived, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, were seen to be very small, and they have a negligible influence on the tip velocity. As a result, it can be concluded that they are not responsible for the discrepancies between experiments and conventional Stokes simulations. A further part of our investigation focussed on the development of high-accuracy three-dimensional spectral element simulation techniques for miscible flows in capillary tubes, including the effects of variable density and viscosity. Towards this end, the conservation equations are treated in cylindrical coordinates.

Meiburg, Eckart↗

Evaluating geothermal and hydrogeologic controls on regional groundwater temperature distribution

Described is a one-dimensional (1-D) analytic solution is developed for heat transport through an aquifer system where the vertical temperature profile in the aquifer is nearly uniform. The general anisotropic form of the viscous heat generation term is developed for use in groundwater flow simulations. The 1-D solution is extended to more complex geometries by solving the equation for piece-wise linear or uniform properties and boundary conditions. A moderately complex example, the Eastern Snake River Plain (ESRP), is analyzed to demonstrate the use of the analytic solution for identifying important physical processes. For example, it is shown that viscous heating is variably important and that heat conduction to the land surface is a primary control on the distribution of aquifer and spring temperatures. Use of published values for all aquifer and thermal properties results in a reasonable match between simulated and measured groundwater temperatures over most of the 300 km length of the ESRP, except for geothermal heat flow into the base of the aquifer within 20 km of the Yellowstone hotspot. Previous basal heat flow measurements (110 mW/m 2 ) made beneath the ESRP aquifer were collected at distances of > 50 km from the Yellowstone Plateau, but a higher basal heat flow of 150 mW/m 2 is required to match groundwater temperatures near the Plateau. As a result, the ESRP example demonstrates how the new tool can be used during preliminary analysis of a groundwater system, allowing efficient identification of the important physical processes that must be represented during more-complex 2-D and 3-D simulations of combined groundwater and heat flow.

58 GEOSCIENCES↗

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↗

Implementation of a Transition Model in a NASA Code and Validation Using Heat Transfer Data on a Turbine Blade

The purpose of this report is to summarize and document the work done to enable a NASA CFD code to model laminar-turbulent transition process on an isolated turbine blade. The ultimate purpose of the present work is to down-select a transition model that would allow the flow simulation of a variable speed power turbine to be accurately performed. The flow modeling in its final form will account for the blade row interactions and their effects on transition which would lead to accurate accounting for losses. The present work only concerns itself with steady flows of variable inlet turbulence. The low Reynolds number k- model of Wilcox and a modified version of the same model will be used for modeling of transition on experimentally measured blade pressure and heat transfer. It will be shown that the k- model and its modified variant fail to simulate the transition with any degree of accuracy. A case is thus made for the adoption of more accurate transition models. Three-equation models based on the work of Mayle on Laminar Kinetic Energy were explored. The three-equation model of Walters and Leylek was thought to be in a relatively mature state of development and was implemented in the Glenn-HT code. Two-dimensional heat transfer predictions of flat plate flow and two-dimensional and three-dimensional heat transfer predictions on a turbine blade were performed and reported herein. Surface heat transfer rate serves as sensitive indicator of transition. With the newly implemented model, it was shown that the simulation of transition process is much improved over the baseline k- model for the single Reynolds number and pressure ratio attempted; while agreement with heat transfer data became more satisfactory. Armed with the new transition model, total-pressure losses of computed three-dimensional flow of E3 tip section cascade were compared to the experimental data for a range of incidence angles. The results obtained, form a partial loss bucket for the chosen blade. In time the loss bucket will be populated with losses at additional incidences. Results obtained thus far will be discussed herein.

Ameri, Ali A.↗

An experimental study of jet exhaust simulation

Afterbody drag predictions for jet aircraft are usually made experimentally with the jet exhaust flow simulated. The physical gas properties of the fluid used for the model jet exhaust can affect the accuracy of simulation of the airplane's jet exhaust plume. The effect of the accuracy of this simulation on afterbody drag was investigated by wind-tunnel tests with single engine model. In addition to unheated air as the exhaust gas, the decomposition products of three different concentrations of hydrogen peroxide were utilized. The air jet simulation consistently resulted in higher boattail drag than hydrogen peroxide simulation. The differences in drag for the various exhaust gases are attributed to different plume shapes and entrainment properties of the gases. The largest differences in drag due to exhaust gas properties were obtained for the combination of high transonic Mach numbers and high boattail angles. For these conditions, the current data indicate that the use of air to simulate a nonafterburning turbojet exhaust can result in an increase in afterbody amounting to 20 percent of the nonafterburning turbojet value.

Compton, W. B., III↗

Developments and applications of dynamic models for large eddy simulation of complex flows

The dynamic modeling procedure for large eddy simulation of turbulent flows is reviewed and recent developments in the theoretical aspects and applications are described. Methods for inclusion of backscatter of energy from small to large scale motions are presented. New formulations of the dynamic procedure are proposed which are optimized based on the subgrid scale flux vector or the energy dissipation rate instead of the subgrid scale stress tensor. Recent results from application of the model to forced isotropic turbulence with an inertial subrange, flow over a backward facing step at Reynolds number of 28000, and flow over a concave curved surface are presented.

Moin, P.↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1996 Summer Program - Part 6

Topics considered include: New approach to turbulence modeling; Second moment closure analysis of the backstep flow database; Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers; Turbulent flame propagation in partially premixed flames; Ensemble averaged dynamic modeling. Also included a study of the turbulence structures of wall-bounded shear flows; Simulation and modeling of the elliptic streamline flow.

Source record↗

Assessment of three numerical methods for the computation of a low-density plume flow

Results from three numerical methods including one based on the Navier-Stokes equations, one based on kinetic theory using the DSMC method, and one based on the Boltzmann equation with a Krook-type collision term are compared to each other and to experimental data for a model problem of heated nitrogen flow in a conical nozzle expanding into a vacuum. The problem simulates flow in a resistojet, a low-thrust, electrothermal rocket. The continuum method is applied to both the internal flow and near-field plume. The DSMC and Boltzmann methods are applied primarily to the plume. Experimental measurements of Pitot pressure and flow angle, taken with an apparatus that duplicates the model nozzle flow, are used in the comparisons.

Penko, Paul F.↗

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three dimensional unsteady results were produced and then time averaged for the collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Collection efficiencies were generated for three spherical particle sizes, 100, 20, and 5 micron in diameter, using the codes LEWICE3D and LEWICE2D. The free stream Mach number was 0.27, representing a velocity of approximately 86 ms. It was observed that a reduction in velocity of about 15-20 occurred as the flow entered the shroud of the probe.Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the effect of the protective shroud.

Aircraft Icing↗

Asymmetry Effects in Numerical Simulation of Supersonic Flows with Upstream Separated Regions

The present paper studies the numerical simulation of flows with shock/boundary-layer upstream interaction, under conditions of symmetry in geometry, boundary conditions, and grid. For this purpose, a series of two- and three-dimensional numerical test-cases were carried out. The tests showed that standard numerical schemes, which appear to be symmetry-preserving under most flow configurations, produce nonsymmetric perturbations when large separated regions are present. These perturbations are amplified when the core flow is under compression. If the flow-blockage due to separation is sufficiently large, the symmetry of the flow may collapse altogether. Experimental evidence of this numerical behavior is also considered.

Rodriguez, C. G.↗

Asymmetry Effects in Numerical Simulation of Supersonic Flows with Upstream Separated Regions

The present paper studies the numerical simulation of flows with shock/boundary-layer upstream interaction, under conditions of symmetry in geometry, boundary conditions, and grid. For this purpose, a series of two- and three-dimensional numerical test-cases were carried out. The tests showed that standard numerical schemes, which appear to be symmetry preserving under most flow configurations, produce nonsymmetric perturbations when large separated regions are present. These perturbations are amplified when the core flow is under compression. If the flow-blockage due to separation is sufficiently large, the symmetry of the flow may collapse altogether. Experimental evidence of this numerical behavior is also considered.

Riodriguez, C. G.↗

Simulation of blood flow through an artificial heart

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady state and unsteady flow calculations are done by solving the incompressible Navier-Stokes equations in 3-D generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Kiris, Cetin↗