Search NASA⌕ Search

SEARCH · Search NASA

Results for “constitutive flow 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 73 records · Page 4

Continuum theories for fluid-particle flows: Some aspects of lift forces and turbulence

A general framework is outlined for the modeling of fluid particle flows. The momentum exchange between the constituents embodies both lift and drag forces, constitutive equations for which can be made explicit with reference to known single particle analysis. Relevant results for lift are reviewed, and invariant representations are posed. The fluid and particle velocities and the particle volume fraction are then decomposed into mean and fluctuating parts to characterize turbulent motions, and the equations of motion are averaged. In addition to the Reynolds stresses, further correlations between concentration and velocity fluctuations appear. These can be identified with turbulent transport processes such as eddy diffusion of the particles. When the drag force is dominant, the classical convection dispersion model for turbulent transport of particles is recovered. When other interaction forces enter, particle segregation effects can arise. This is illustrated qualitatively by consideration of turbulent channel flow with lift effects included.

Mctigue, David F.↗

Verification Test Suite for Spalart-Allmaras QCR2000 Turbulence Model

The paper presents three benchmark cases for verification of Reynolds-averaged Navier-Stokes solvers. The verification studies focus on a one-equation Spalart-Allmaras model, SA-[neg]-QCR2000, that uses a version of quadratic constitutive relations.The benchmark cases are a two-dimensional subsonic flow around a Joukowski airfoil, a three-dimensional supersonic flow through a square duct, and a three-dimensional flow over a wing-fuselage configuration. The turbulence-model formulation, geometry, flow conditions, grids, and the expected output are described in detail. Reference solutions computed by several established codes are shown

Computational Aerodynamics↗

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Biaxial experiments supporting the development of constitutive theories for advanced high-temperature materials

Complex states of stress and strain are introduced into components during service in engineering applications. It follows that analysis of such components requires material descriptions, or constitutive theories, which reflect the tensorial nature of stress and strain. For applications involving stress levels above yield, the situation is more complex in that material response is both nonlinear and history dependent. This has led to the development of viscoplastic constitutive theories which introduce time by expressing the flow and evolutionary equation in the form of time derivatives. Models were developed here which can be used to analyze high temperature components manufactured from advanced composite materials. In parallel with these studies, effort was directed at developing multiaxial testing techniques to verify the various theories. Recent progress in the development of constitutive theories from both the theoretical and experimental viewpoints are outlined. One important aspect is that material descriptions for advanced composite materials which can be implemented in general purpose finite element codes and used for practical design are verified.

Ellis, J. R.↗

Juncture Flow Computations using kL-Based Turbulence Models

The development and implementation of kL-based Reynolds-Average Navier-Stokes (RANS) two-equation turbulence models are reported in this paper. The kL model is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta showed that a reliable transport equation can be formed from the turbulence length scale L, and the turbulence kinetic energy k. Rotta's kL equation is well suited for a term-by-term modeling and displays useful features compared to other scale formulations. 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 k-kL scheme has been modified to include the Algebraic Reynolds Stress Model (ARSM) and the Quadratic Constitutive Relation (QCR) as nonlinear models. These models have gone through extensive validations using two-dimensional benchmark problems. The present study documents the application of the k-kL schemes to simulate flow around the Juncture Flow Model (JFM). The k-kL prediction results show generally good comparisons with measurements. The results from this formulation are similar to, or better than results using the SA-RC-QCR2000 two-equation turbulence model. The ARSM and the QCR formulations of k-kL show promise with a similar level of computational resources as basic two-equation turbulence models.

K S Abdol-Hamid↗

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗

SAM-ML: Integrating data-driven closure with nuclear system code SAM for improved modeling capability

Advanced reactors often involve complicated thermal-fluid (T-F) phenomena. Modeling such phenomena with the traditional one-dimensional (1-D) system code is a challenging task. The System Analysis Module (SAM), a modern nuclear system code, has developed a coarse mesh multi-dimensional (multi-D) flow model to capture the spatial effect of T-F phenomena in advanced reactors. As a coarse mesh solver, constitutive relations are required for SAM's multi-D model for unresolved fine-scale physics, such as turbulence. Here this work presents a novel approach that integrates neural networks as data-driven closure for SAM's multi-D flow model. The data-driven closure is trained with fine-resolution data to ensure its accuracy while maintaining a coarse mesh setup to ensure its efficiency and consistency with SAM. We demonstrate the applicability of this SAM-ML capability in an open volume thermal stratification problem, where a neural network model serves as the eddy viscosity closure. A customized interface between the neural network and SAM is developed to ensure flexible and efficient data exchange. The SAM-ML results demonstrate superior performance compared to SAM's built-in zero-equation eddy viscosity closure. The case study shows that although the generalization capability of the data-driven closure still needs to be improved for different transient case or different geometric setup, SAM -ML demonstrates good potential for challenging simulation problems with improved accuracy and computational efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗

Vadose Zone Flow and Transport Parameters Data Package for the Hanford Site Composite Analysis

This report provides a description of the basis for the development and implementation of a conceptual model for vadose zone flow and transport for the Hanford Site composite analysis (CA) groundwater pathway analysis. The parameterization for a numerical model is intimately linked to the conceptual model framework. This report describes the basis for the selection of hydraulic and transport parameters for the hydrostratigraphic units (HSUs) identified in the 200 East and 200 West Areas. The following three-step process was used to develop the hydraulic properties and transport parameters for the CA modeling: 1. Breakdown of existing database on hydraulic properties by HSUs 2. Constitutive model parameters for laboratory, core-scale hydraulic properties 3. Upscaling for macroscopic, field-scale flow, and transport parameters This document provides data and information on the following items: A brief description and technical basis for the conceptual model selected to represent vadose zone flow and transport for the CA modeling. Variation in sediment particle size distribution for representative sites in the 200 East and West Areas. A summary of effective (upscaled) moisture retention, and saturated and unsaturated hydraulic conductivity estimates for the HSUs in the 200 East and West Areas. A summary of effective transport parameters, including sediment bulk density and macrodispersivity estimates for various HSUs. The 200 East Area B Complex perched water aquifer hydraulic properties

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Predictive Flow Simulation with the P2R Model for the 200-IA-1 Preliminary Remediation Goal Saturated Zone Abstraction

In order to develop and evaluate alternatives for remediation of waste sites constituting the 200-IA-1 Operable Unit (OU), a focused feasibility study (FFS) is currently underway. The 200-IA-1 OU FFS requires preliminary remediation goals (PRGs) protective of groundwater be evaluated which includes estimates of fate and transport of contaminants in the groundwater from multiple sources within the Central Plateau Inner Area. This document details the application of the Plateau-to-River (P2R) Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the 200-IA-1 OU PRG calculation. The simulated flow field will support the simulation of fate and transport of contaminants for use in estimating peak unit-length concentration as part of the 200-IA-1 OU PRG calculation presented in a separate environmental calculation file (ECF) (ECF-200IA1-20-0108, Determination of Preliminary Remediation Goals for 200-IA-1 Source Operable Unit).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Kinetic Model Development of the Carbon Vapor Deposition Reactor to Predict Performance Versus Design

To enable complete closure of a Sabatier-based atmosphere revitalization system where oxygen recovery is predicated by carbon dioxide reduction to methane and water, there must be a means of recovering the hydrogen from the byproduct methane. Honeywell’s carbon vapor deposition (CVD) reactor is a potential technology that provides a thermochemical means by which methane decomposes at high temperatures (> 1,000 °C) into hydrogen, solid carbon, and a small percentage of hydrocarbons. This advanced methane pyrolysis reactor utilizes a carbon fiber substrate to capture and immobilize the solid pyrolytic carbon rather than allowing continuation toward undesirable soot. The CVD reactor is a highly dynamic system whose characteristic geometry changes while operating. The reactor model was developed to provide predictive performance assessment capabilities at different operating conditions, geometries, and configurations to aid in technology development. Modeling of the CVD reactor considered the multitude of simultaneous and interdependent physicochemical processes that are entailed, including: (1) the homogeneous gas-phase pyrolysis reactions, (2) the heterogeneous deposition reactions, and (3) the transport of gaseous species in the reactor domains. A 1-D CVD reactor model was developed based on a general packed-bed, axial-dispersion plug flow reactor model with a solid substrate balance to account for carbon deposition. The model includes 59 species and 243 reactions constituting the methane pyrolysis reaction network, carbon deposition, and interconversions of oxygenate species. The reactor model was compared to brassboard test data, where it was shown to be able to capture the pertinent physicochemical phenomena as evidenced by its agreement with short-duration run data (i.e., looking at different temperature and flow rate conditions), as well as trend agreement with the brassboard reactor’s long-duration test data. Recommendations are provided on further model development steps as well as its use cases.

Air Revitalization↗

Agricultural Impacts on Nitrogen Cycling: Climate and Air Pollution (Final Technical Report)

We have successfully added a parameterization of ammonia emissions from synthetic fertilizer and manure from agriculture to Earth System Models. We have evaluated how emissions change in the future and the processes driving these changes. We have examined the impact of agricultural ammonia emissions on present-day air quality and have evaluated the impact of resolving the emissions with high temporal resolution. We have also examined the impact of how Earth System Models treat the cycling of emissions in soils. Our code, the Flows of Agricultural Nitrogen – version 2 – model (FANv2) has been run in both the CLM framework and the ELM framework. To my knowledge the code has not yet been fully implemented into ELM. Volatilization of ammonia (NH3) from fertilizers and livestock wastes forms a significant pathway of nitrogen losses in agricultural ecosystems and constitutes the largest source of atmospheric emissions of NH3. We have provided a major update to the process model FAN (Flow of Agricultural Nitrogen), which evaluates NH3 emissions interactively within an Earth system model; in this work, the Community Earth System Model (CESM) is used. The updated version (FANv2) includes a more detailed treatment of both physical and agricultural processes, which allows the model to differentiate between the volatilization losses from animal housings, manure storage, grazed pastures, and the application of manure and different types of mineral fertilizers. The modeled ammonia emissions are first evaluated at a local scale against experimental data for various types of fertilizers and manure, and they are subsequently run globally to evaluate NH3 emissions for 2010–2015 based on gridded datasets of fertilizer use and livestock populations. Comparison of regional emissions shows that FANv2 agrees with previous inventories for North America and Europe and is within the range of previous inventories for China. However, due to higher NH3 emissions in Africa, India, and Latin America, the global emissions simulated by FANv2(48 Tg N) are 30 %–40 % higher than in the existing inventories.

54 ENVIRONMENTAL SCIENCES↗

Turbine Engine Hot Section Technology, 1984

Presentations were made concerning the hot section environment and behavior of combustion liners, turbine blades, and waves. The presentations were divided into six sessions: instrumentation, combustion, turbine heat transfer, structural analysis, fatigue and fracture, and surface properties. The principal objective of each session was to disseminate research results to date, along with future plans. Topics discussed included modeling of thermal and fluid flow phenomena, structural analysis, fatigue and fracture, surface protective coatings, constitutive behavior, stress-strain response, and life prediction methods.

Source record↗

Rheological properties of nondilute suspensions of deformable particles

The creeping flow equations have been solved to account for effects of shear on drop shape when the influence of neighboring drops is considered. A cell model is used to account for finite concentrations of the dispersed phase. From these results a constitutive equation is derived and implications for several flows are discussed. The shape and orientation of drops are computed and compared with existing data. Finally, dynamic viscosity and rigidity are obtained for nondilute suspensions of viscoelastic particles in an oscillatory shear flow.-

Choi, S. J.↗

Strain shift measured from stress-controlled oscillatory shear: Evidence for a continuous yielding transition and new techniques to determine recovery rheology measures

Understanding the yielding of complex fluids is an important rheological challenge that affects our ability to engineer and process materials for a wide variety of applications. Common theoretical understandings of yield stress fluids follow the Oldroyd–Prager formalism in which the material behavior below the yield stress is treated as solidlike, and above the yield stress as liquidlike, with an instantaneous transition between the two states. This formalism was built on a quasi-static approach to the yield stress, while most applications, ranging from material processing to end user applications, involve a transient approach to yielding over a finite timescale. Using stress-controlled oscillatory shear experiments, we show that yield stress fluids flow below their yield stresses. This is quantified through measuring the strain shift, which is the value about which the strain oscillates during a stress-controlled test and is a function of only the unrecoverable strain. Measurements of the strain shift are, therefore, measurements of flow having taken place. These experimental results are compared to the Herschel–Bulkley form of the Saramito model, which utilizes the Oldroyd–Prager formalism, and the recently published Kamani–Donley–Rogers (KDR) model, in which one constitutive equation represents the entire range of material responses. Scaling relationships are derived, which allow us to show why yield stress fluids will flow across all stresses, above and below their yield stress. Finally, derivations are presented that show strain shift can be used to determine average metrics previously attainable only through recovery rheology, and these are experimentally verified.

Mechanics↗

Anisotropic constitutive modeling for nickel-base single crystal superalloys

An anisotropic constitutive model was developed based on crystallographic slip theory for nickel base single crystal superalloys. The constitutive equations developed utilizes drag stress and back stress state variables to model the local inelastic flow. Specially designed experiments were conducted to evaluate the existence of back stress in single crystal superalloy Rene N4 at 982 C. The results suggest that: (1) the back stress is orientation dependent; and (2) the back stress state variable is required for the current model to predict material anelastic recovery behavior. The model was evaluated for its predictive capability on single crystal material behavior including orientation dependent stress-strain response, tension/compression asymmetry, strain rate sensitivity, anelastic recovery behavior, cyclic hardening and softening, stress relaxation, creep and associated crystal lattice rotation. Limitation and future development needs are discussed.

Sheh, Michael Y.↗

Turbine Engine Hot Section Technology 1986

On October 21 and 22, 1986, the Turbine Engine Hot Section Technology (HOST) Project of the NASA Lewis Research Center sponsored a workshop to discuss current research pertinent to turbine engine durability problems. Presentations were made concerning the hot-section environment and the behavior of combustion liners, turbine blades, and turbine vanes. The presentations were divided Into six sessions: Instrumentation, Combustion, Turbine Heat Transfer, Structural Analysis, Fatigue and Fracture, and Surface Protection. The principal objective of each session was to disseminate research results to date. Topics discussed included modeling of thermal and fluid-flow phenomena, structural analysis, fatigue and fracture, surface protective coatings, constitutive behavior of materials, stress-strain response, and life-prediction methods. Researchers from industry, academia, and government presented results of their work sponsored by the HOST Project. This publication contains extended abstracts and visual material presented during the workshop.

Turbine Engine Hot Section Technology (HOST)↗