Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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,801 records · Page 100

Simulation of Methane Oxycombustion in Supercritical Carbon Dioxide

The Allam-Fetvedt Cycle (AFC) offers the potential for electricity generation with minimal or even negative carbon emissions by using oxycombustion in supercritical carbon dioxide (sCO2). Because of the extreme high temperatures and pressures encountered in these systems, simulation of oxycombustion requires use of a complex equation of state such as the modified Soave-Redlich-Kwong (SRK) equation, coupled with equally complex correlations for transport properties. The current work simulates a model combustor using both an ideal gas equation of state and the SRK equation of state. The model combustor features a jet of CH4 mixing with and reacting with O2 in the sCO2 ambient. Two jet configurations are considered: one of pure CH4, and one of CH4 mixed with sCO2. Comparison of results shows that key output properties such as combustor temperature and flow distribution are impacted by the choice of equations of state. Inspection of simulation results shows that mixing and reaction rates, as well as the uniformity of the flow, are impacted by the choice of equation of state. These results show that even though use of a complex equation of state greatly increases the computational costs of these simulations, it is critical to accurately capturing the physics of these systems.

computational fluid dynamics↗

A perspective of theoretical and applied computational fluid dynamics

The present effectiveness and future promise of computational fluid dynamics (CFD) are discussed. The CFD ideal of computer-designed flight vehicles and the progress that has been made toward its achievement, expressed in terms of primary and secondary pacing items, are also addressed. The primary pacing items listed by Chapman - grid generation, turbulence modeling, solution methodology development, and mainframe computer methods and architecture - are modified and reviewed. Important secondary pacing items in CFD development are also examined: algorithm development, complex geometry definition, and predata and postata processing. In addition, user requirements placed on CFD computer codes are defined and some problems associated with the dissemination of CFD technology are identified.

Kutler, P.↗

Recent experience in seeding transonic/supersonic flows at AEDC

The laser velocimeter has been utilized for several years at the Arnold Engineering and Development Center (AEDC) as a flow diagnostics tool. Most applications, following the initial proof-of-concept experiments, have involved relatively complex unknown flow fields in which the more conventional, intrusive techniques had either not been attempted or had yielded unsatisfactory results. A blunt-base nozzle-afterbody base flow study is discussed as a respresentative example of such applications. A wide variety of problems have been encountered during these tests, many of which have proven to be closely related to the size and/or size distribution of the seeding material within the fluid. Resulting measurement uncertainties could often not be conclusively resolved because of the unknown nature of the flow field. The other experiments listed were conducted to provide known aerodynamic conditions for comparison with the velocimeter results.

Heltsley, F. L.↗

Fast molecular shocks. I - Reformation of molecules behind a dissociative shock

The physical and chemical processes that operate in the cooling gas behind a fast, dissociative, single-fluid shock propagating in a dense interstellar cloud are discussed. The treatment extends previous theoretical work on fast molecular shocks by including the effects of the conversion of Ly-alpha photons into radiation of the two-photon continuum and into H2 Lyman band emission lines, the effects of CO photodissociation following line absorption, and the formation and destruction of molecules containing the elements nitrogen, silicon, and sulphur, and of the complex hydrocarbons. Abundance profiles for the molecular species of interest are presented. After molecular hydrogen begins to reform, by means of gas phase and grain surface processes, the neutral species OH, H2O, O2, CO, CN, HCN, N2, NO, SO, and SiO reach substantial abundances. The molecular ions HeH(+), OH(+), SO(+), CH(+), H2(+), and H3(+), are produced while the gas is still hot and partially ionized. Emissions from them provide a possible diagnostic probe of fast molecular shocks.

Neufeld, David A.↗

Tensoral: A system for post-processing turbulence simulation data

Many computer simulations in engineering and science -- and especially in computational fluid dynamics (CFD) -- produce huge quantities of numerical data. These data are often so large as to make even relatively simple post-processing of this data unwieldy. The data, once computed and quality-assured, is most likely analyzed by only a few people. As a result, much useful numerical data is under-utilized. Since future state-of-the-art simulations will produce even larger datasets, will use more complex flow geometries, and will be performed on more complex supercomputers, data management issues will become increasingly cumbersome. My goal is to provide software which will automate the present and future task of managing and post-processing large turbulence datasets. My research has focused on the development of these software tools -- specifically, through the development of a very high-level language called 'Tensoral'. The ultimate goal of Tensoral is to convert high-level mathematical expressions (tensor algebra, calculus, and statistics) into efficient low-level programs which numerically calculate these expressions given simulation datasets. This approach to the database and post-processing problem has several advantages. Using Tensoral the numerical and data management details of a simulation are shielded from the concerns of the end user. This shielding is carried out without sacrificing post-processor efficiency and robustness. Another advantage of Tensoral is that its very high-level nature lends itself to portability across a wide variety of computing (and supercomputing) platforms. This is especially important considering the rapidity of changes in supercomputing hardware.

Dresselhaus, Eliot↗

Numerical Investigation of Slag Behavior for RSRM

It is known that the flow field of the redesigned solid rocket motor (RSRM) is very complicated due to the complex characteristics of turbulent multi-phase flow, chemical reaction, particle combustion, evaporation, breakup and agglomeration etc. It requires multi-phase calculations, chemical reaction simulation, and particle combustion, evaporation, and breakup models to obtain a better understanding of thermophysics for the RSRM design using numerical methods. Also, the slag buildup due to the molten particles is another factor affecting the performance of the RSRM. To achieve this goal, the volume of fluid (VOF) method is used to capture the free surface motion so as to simulate the accumulation of the molten particles (slag) of the RSRM. A finite rate chemistry model is used to simulate the chemical reaction effects. For multi-phase calculations, the Hermsen combustion model is used for the aluminum particle combustion analysis and the Taylor Analogy Breakup (TAB) model is used for the particle breakup analysis. An interphase mas-exchange model introduced by Spalding is used for the evaporation calculation. The particle trajectories are calculated using a one-step implicit method for several groups of particle sizes by which the drag forces and heat fluxes are then coupled with the gas phase equations. The preliminary results predicted a reasonable physical simulation of the particle effects using a simple two dimensional solid rocket motor configuration. It shows that the AL/AL2O3 particle sizes are reduced due to the combustion, evaporation, and breakup. The flow field is disturbed by the particles. Mach number distributions in the nozzle are deformed due to the effect of particle concentrations away from the center line.

Liaw, P.↗

Doppler Global Velocimeter Measurements of the Vortical Flow Above a Thin Delta Wing

A new measurement technique is being developed by NASA to measure off-surface flow fields. This method, Doppler global velocimetry, will allow quantification of complex three-dimensional flow fields at video camera rates. The entire flow field structure within a selected plane is measured simultaneously rather than by scanned, point-by-point measurements using conventional laser velocimetry. Data obtained using this technique will be used to correlate with other data sets for verification, and following verification, provide a quantified, highly detailed definition of the flow field. This will help to improve the understanding of fluid physics, supplement and broaden the database required to validate and refine computational fluid dynamics (CFD) models, and improve aircraft design methodology. To assess the capability of the technique, velocity measurements of the vortical flow field above a thin 75-degree delta wing were made in the NASA - Langley Basic Aerodynamics Research Tunnel. Preliminary comparisons of the results were made with similar measurements obtained using a three component laser velocimeter indicate that this technique is capable of describing the entire three - component velocity flow field simultaneously within a measurement plane in real time.

Usry, Jimmy W.↗

Utilizing Radioisotope Power System Waste Heat for Spacecraft Thermal Management

One of the advantages of using a Radioisotope Power System (RPS) for deep space or planetary surface missions is the readily available waste heat, which can be used to maintain electronic components within a controlled temperature range, to warm propulsion tanks and mobility actuators, and to gasify liquid propellants. Previous missions using Radioisotope Thermoelectric Generators (RTGs) dissipated a very large quantity of waste heat due to the relatively low efficiency of the thermoelectric conversion technology. The next generation RPSs, such as the 110-watt Stirling Radioisotope Generator (SRG110) will have much higher conversion efficiencies than their predecessors and therefore may require alternate approaches to transferring waste heat to the spacecraft. RTGs, with efficiencies of approx. 6 to 7% and 200 C housing surface temperatures, would need to use large and heavy radiator heat exchangers to transfer the waste heat to the internal spacecraft components. At the same time, sensitive spacecraft instruments must be shielded from the thermal radiation by using the heat exchangers or additional shields. The SRG110, with an efficiency around 22% and 50 C nominal housing surface temperature, can use the available waste heat more efficiently by more direct heat transfer methods such as heat pipes, thermal straps, or fluid loops. The lower temperatures allow the SRG110 much more flexibility to the spacecraft designers in configuring the generator without concern of overheating nearby scientific instruments, thereby eliminating the need for thermal shields. This paper will investigate using a high efficiency SRG110 for spacecraft thermal management and outline potential methods in several conceptual missions (Lunar Rover, Mars Rover, and Titan Lander) to illustrate the advantages with regard to ease of assembly, less complex interfaces, and overall mass savings.

Pantano, David R.↗

Methodology for Thermodynamic Analysis Coupled with Computational Fluid Dynamics Modeling for Casting a Novel Aluminum–Cerium Alloy

In this paper, we present the results of computational fluid dynamics (CFD) analysis to assess castability, porosity, molten metal fluidity and other technological properties of a new Al-Ce alloy. A thermodynamic analysis of a new Al-Ce alloy is done to obtain the physical properties as a function of temperature across both solid and liquid phases. These properties are then used to build the CFD of model the full casting process from initial pouring through to final solidification of the part, here an example of a heavy duty mor mount is used. While several excellent commercial CFD codes exists this process shows that CFD can be used to assess the capability of the alloy to properly fill the mold, as well as give predictions where scattered porosity or large-scale defects may occur in the casting. Further, unlike the commercially available software (e.g., ProCAST, SOLIDCast , MAGMASOFT®) the complex 3D-analysis of the stress / strain fields in cast parts is not performed at this time. However, the availability of free software for assessing the required thermodynamic / thermo-physical properties of new alloys (OpenCALPHAD) and CFD codes such as OpenFOAM® makes the developed option attractive and economical, especially for the analysis of new Al-Ce alloys, for which the available data does not exist.

36 MATERIALS SCIENCE↗

Mixing of cryogenic fluids and predicted detonation properties for multiphase liquid oxygen and liquid hydrogen

The mixing and explosion of LOX and LH sub 2 is a concern for Space Transportation System operations. To understand this problem, cryogenic mixing is experimentally studied by pouring 1,2,2 trichloro 2,1,1 trifluoroethane (Freon 113) into LN sub 2, LN sub 2 into LH sub 2, and LH sub 2 into into LN sub 2 in a 1 m by 15 cm cylindrical glass vessel. Data from these experiments is compared with previous studies and a hypothesis is advanced that LOX/LH sub 2 mixing will result in a complex, heterogenous, multiphase aggregation including LOX, SOX, LH sub 2, and VH sub 2 using the multiphase hypothesis. Visual and x ray observations of the process and mass measurements are reported.

Bishop, C. V.↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction of hemiwicking dynamics in micropillar arrays

Dynamic hemiwicking behavior is observable in both nature and a wide range of industrial applications ranging from biomedical devices to thermal management. We present a semi-analytical modeling framework (without empirical fitting coefficients) to predict transient capillary-driven hemiwicking behavior of a liquid through a nano/microstructured surface, specifically a micropillar array. In our model framework, the liquid domain is discretized into micropillar unit cells to enable the time marching of the hemiwicking front. A simplified linear pressure drop is assumed along the hemiwicking length such that the local meniscus curvature, contact angle, and effective liquid height are determined at each time step in our transient model. This semi-analytical model is validated with experimental data from our own experiments and from published literature for different fluids. Our model predicts hemiwicking dynamics with <20% error over a broad range of micropillar geometries with height-to-pitch ratio ranging between ≈0.34 and 6.7 and diameter-to-pitch ratio in the range of ≈0.25–0.7 and without any fitting parameters. For lower diameter-to-pitch ratio data points related to sparse micropillar array arrangements, we suggest modifications to the semi-analytical model. This work sheds light on complex and dynamic solid–liquid–vapor interfacial interactions which could serve as a guide for the design of textured surfaces for wicking enhancement in multi-phase thermal and mass transport technologies and applications.

Mechanics↗

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗

The annular flow electrothermal ramjet

The annular flow, electrothermal, plug ramjet is examined as a possible means of achieving rapid projectile acceleration to velocities for such applications as direct launch of spacebound payloads. The performance of this ramjet operating with hydrogen propellant is examined for cases where this working fluid is treated: (1) as a perfect gas, and (2) as a gas that is allowed to dissociate and ionize and then recombine with finite reaction rates in the nozzle. Performance results for these cases are compared to the performance of a conventional ramjet operating with perfect gas hydrogen propellant. The performance of the conventional ramjet is superior to that of the annular flow, electrothermal ramjet. However, it is argued that the mechanical complexities associated with conventional ramjet operation are difficult to attain, and for this reason the annular flow, electrothermal ramjet is more desirable as a launch system. Models are presented which describe both electrothermal plug ramjet and conventional ramjet operation, and it is shown that for a given flight velocity there is a rate of heat addition per unit propellant mass for which ramjet operation is optimized.

Shaw, B. D.↗

Geometric description and grid generation for space vehicles

The understanding and analysis of three dimensional fluid behavior around reentry vehicles, like the Space Shuttle Orbiter, is of considerable importance for flight applications. Vehicle designers need pressure, heat transfer, and shear distributions to support design activities. Traditionally these have been obtained by experimental data. Significant reduction in design costs can be obtained if the calculated results based on finite difference procedures are used to verify and supplement experimental results. Computing the flow field over Space Shuttle Orbiter, which has complex cross sectional geometries is a major problem. A first step in this direction is development of body definition and grid generation procedure. In this project a digitized body definition and grid generation procedure, developed earlier at JSC, is modified and applied to the Shuttle Orbiter X24C-10D. Two dimensional grid and body definition for various sections of the space vehicle is obtained using the differential equation method. This study demonstrates that this geometry and grid generation procedure can be applied to other aerospace vehicles such as aerospaceplanes.

Krishna, Madakasira G.↗

Wavenumber Vacillation in Weakly-Stratified Baroclinic Flows

The wavenumber selection in an unstable baroclinic flow has long been a subject of investigation in geophysical fluid dynamics. For a weakly unstable flow, the linearly most unstable wave is always the final realized wave. As the flow becomes more unstable, the dominant wave gradually shifts from the most unstable wave to a longer wave. However, when the supercriticality is reduced from the state with a dominant longwave, it may remain in that state even beyond the threshold where it first developed. Hysteresis in a quasigeostrophic system also displays a complex behavior: it occurs not only between the states of different dominant wavenumbers, but also between the states of identical dominant wavenumber but of different dynamic characteristics. In a recent numerical experiment using a high-resolution spectral Eady model with asymmetric Ekman dissipation, it is found that in a weakly stratified flow with strong meridional forcing the flow does not show a preference for any single wave. Instead, the dominant wave periodically varies among the linearly most unstable wave and the longer waves. The mechanism of wavenumber vacillation will be examined in terms of meridional heat flux, the wave PV gradient, and wave saturation. The importance of wave-wave and wave-mean flow interactions will also be investigated.

Chou, Shih-Hung↗

Refine your search to explore more results.