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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Non-newtonian Effects in Viscous Flows

Revision of the mathematical formalism of fluid dynamics suggests that some physical inconsistencies (infinite time of approaching equilibrium and fully deterministic solutions to the Navier-Stokes equations) can be removed by relaxing the Lipschitz conditions, i.e., the boundedness of the derivatives, in the constitutive equations. Physically such a modification can be interpreted as an incorporation of an infinitesimal static friction in the constitutive law. A modified version of the Navier-Stokes equations is introduced, discussed, and illustrated by examples. It is demonstrated that all the new effects in the modified model emerge within vanishingly small neighborhoods of equilibrium states which are the only domains where the governing equations are different from classical.

Zak, Michail↗

NASA NTRS1989

Computer simulation of the mathematical modeling involved in constitutive equation development: Via symbolic computations

Development of new material models for describing the high temperature constitutive behavior of real materials represents an important area of research in engineering disciplines. Derivation of mathematical expressions (constitutive equations) which describe this high temperature material behavior can be quite time consuming, involved and error prone; thus intelligent application of symbolic systems to facilitate this tedious process can be of significant benefit. A computerized procedure (SDICE) capable of efficiently deriving potential based constitutive models, in analytical form is presented. This package, running under MACSYMA, has the following features: partial differentiation, tensor computations, automatic grouping and labeling of common factors, expression substitution and simplification, back substitution of invariant and tensorial relations and a relational data base. Also limited aspects of invariant theory were incorporated into SDICE due to the utilization of potentials as a starting point and the desire for these potentials to be frame invariant (objective). Finally not only calculation of flow and/or evolutionary laws were accomplished but also the determination of history independent nonphysical coefficients in terms of physically measurable parameters, e.g., Young's modulus, was achieved. The uniqueness of SDICE resides in its ability to manipulate expressions in a general yet predefined order and simplify expressions so as to limit expression growth. Results are displayed when applicable utilizing index notation.

Arnold, S. M.↗

NASA NTRS2000

Time-Dependent Reversible-Irreversible Deformation Threshold Determined Explicitly by Experimental Technique

Structural materials for the design of advanced aeropropulsion components are usually subject to loading under elevated temperatures, where a material's viscosity (resistance to flow) is greatly reduced in comparison to its viscosity under low-temperature conditions. As a result, the propensity for the material to exhibit time-dependent deformation is significantly enhanced, even when loading is limited to a quasi-linear stress-strain regime as an effort to avoid permanent (irreversible) nonlinear deformation. An understanding and assessment of such time-dependent effects in the context of combined reversible and irreversible deformation is critical to the development of constitutive models that can accurately predict the general hereditary behavior of material deformation. To this end, researchers at the NASA Glenn Research Center at Lewis Field developed a unique experimental technique that identifies the existence of and explicitly determines a threshold stress k, below which the time-dependent material deformation is wholly reversible, and above which irreversible deformation is incurred. This technique is unique in the sense that it allows, for the first time, an objective, explicit, experimental measurement of k. The underlying concept for the experiment is based on the assumption that the material s time-dependent reversible response is invariable, even in the presence of irreversible deformation.

Castelli, Michael G.↗

NASA NTRS1995

Sensitivity analysis for the dynamic response of thermoviscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric, fully-coupled, thermoviscoplastic response of shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of large-strain, transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be viscoplasticity with strain hardening and thermal hardening, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used for the shell equations with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. The energy-balance equation is solved using a Galerkin procedure, with the temperature as the fundamental unknown. Spatial discretization is performed in one dimension (meridional direction) for the momentum and constitutive equations of the shell, and in two dimensions (meridional and thickness directions) for the energy-balance equation. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) for the momentum equation; a predictor-corrector version of the trapezoidal rule is used for the energy-balance equation; and an explicit scheme consistent with the central difference method is used to integrate the constitutive equations. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading.

Kulkarni, Makarand↗

NASA NTRS2002

Constitutive Relation in Transitional Granular Flows

To study the constitutive behavior of granular materials, the presence of gravity is detrimental. Although empirical relations have been obtained for engineering designs to control granular flows on Earth, it is not known how well these Earthbound relations can be used in another gravity field. Fundamental understanding must be derived to reliably design for granular flows in space exploration. There are two extremes of granular flows of which significant amount of knowledge is available. One deals with a dense and quasi-static situation where the deformation rate nearly vanishes. The other deals with dilute and rapidly fluctuating grain velocities where particle inertia dominates. This project, funded by the NASA Microgravity Fluid Physics Program, aims to study this transitional regime via physical experiments and computer simulations. A conceptual model has been established as described below. There are two natural time scales in a granular flow. One is the travel time between two consecutive collisions and the other is the duration of a collision contact. At a very low shear-rate, the shear-induced particle velocity is low. Hence the travel time between collisions is longer than the contact time between colliding particles. Binary collisions prevail. As the shear-rate increases, the traveling time between collisions reduces and the probability of multiple collisions goes up. These particle groups disperse shortly after and new groups form. When shear-rate is further increased, clusters grow in size due to an increasing chance for free particles to join before groups have the time to disperse. The maximum cluster size may depend on the global concentration and material properties. As the solid concentration approaches zero, the cluster size goes to one particle diameter. The maximum possible cluster size under any condition is the container size, provided that the shear flow is inside a container. The critical shear-rates that dictate the initiation of the multiple contacts, and the size and lifetime of the collision clusters, are functions of the concentration also. A 'regime' theory has been proposed by Babic et al. This theory suggested that both the solid concentration, C, and the non-dimensional shear-rate, B, are important in determining the regimes of the granular constitutive law. Additional information is included in the original extended abstract.

Shen, Hayley H.↗

NASA NTRS1990

The application of CFD to rotary wing flow problems

Rotorcraft aerodynamics is especially rich in unsolved problems, and for this reason the need for independent computational and experimental studies is great. Three-dimensional unsteady, nonlinear potential methods are becoming fast enough to enable their use in parametric design studies. At present, combined CAMRAD/FPR analyses for a complete trimmed rotor soltution can be performed in about an hour on a CRAY Y-MP (or ten minutes, with multiple processors). These computational speeds indicate that in the near future many of the large CFD problems will no longer require a supercomputer. The ability to convect circulation is routine for integral methods, but only recently was it discovered how to do the same with differential methods. It is clear that the differential CFD rotor analyses are poised to enter the engineering workplace. Integral methods already constitute a mainstay. Ultimately, it is the users who will integrate CFD into the entire engineering process and provide a new measure of confidence in design and analysis. It should be recognized that the above classes of analyses do not include several major limiting phenomena which will continue to require empirical treatment because of computational time constraints and limited physical understanding. Such empirical treatment should be included, however, into the developing CFD, engineering level analyses. It is likely that properly constructed flow models containing corrections from physical testing will be able to fill in unavoidable gaps in the experimental data base, both for basic studies and for specific configuration testing. For these kinds of applications, computational cost is not an issue. Finally, it should be recognized that although rotorcraft are probably the most complex of aircraft, the rotorcraft engineering community is very small compared to the fixed-wing community. Likewise, rotorcraft CFD resources can never achieve fixed-wing proportions and must be used wisely. Therefore the fixed-wing work must be gleaned for many of the basic methods.

Caradonna, F. X.↗

NASA NTRS1991

Circulation of the gyres of the world ocean: Observation modeling using TOPEX/Poseidon altimeter data

The overall objectives of the proposed investigation are to study the dynamics of the large-scale recirculating cells of water in the ocean, which are loosely defined as 'gyres' in this study. A gyre is normally composed of a swift western boundary current (e.g., the Gulf Stream and the Kuroshio), a tight recirculating cell attached to the current, and a large-scale sluggish return flow. The water, of course, is not entirely recirculating within a gyre. The exchange of water among gyres is an important process in maintaining the meridional heat transport of the ocean. The gyres constitute a major mode of water movement in the ocean and play significant roles in the global climate system.

Fu, Lee-Lueng↗

NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗

NASA NTRS2008

Transient Effects in Planar Solidification of Dilute Binary Alloys

The initial transient during planar solidification of dilute binary alloys is studied in the framework of the boundary integral method that leads to the non-linear Volterra integral governing equation. An analytical solution of this equation is obtained for the case of a constant growth rate which constitutes the well-known Tiller's formula for the solute transient. The more physically relevant, constant ramping down temperature case has been studied both numerically and analytically. In particular, an asymptotic analytical solution is obtained for the initial transient behavior. A numerical technique to solve the non-linear Volterra equation is developed and the solution is obtained for a family of the governing parameters. For the rapid solidification condition, growth rate spikes have been observed even for the infinite kinetics model. When recirculating fluid flow is included into the analysis, the spike feature is dramatically diminished. Finally, we have investigated planar solidification with a fluctuating temperature field as a possible mechanism for frequently observed solute trapping bands.

Mazuruk, K.↗

NASA NTRS2008

Transient Effects in Planar Solidification of Dilute Binary Alloys

Mazuruk, Konstantin↗

NASA NTRS2024

Computational Analysis of the X-57 Maxwell Airplane, the Landing Configuration with High-Lift Blowing and Aileron Deflections (Preliminary Fuselage)

The X-57 Maxwell is an all-electric airplane with a distributed electric propulsion system used for a high-lift system at takeoff and landing conditions. The Kestrel and USM3D flow solvers were used at NASA Langley to investigate the performance of the X-57 Maxwell in the development of an aerodynamic database. The configuration investigated in this paper had a 30° flap deflection, and in addition the pilot’s right aileron were deflected and undeflected in different cases. The solutions were computed at an airspeed of 58 KEAS, for an altitude of 6000 feet, and a flight Reynolds number of 0.588e+06 per foot. To evaluate the high-lift distributed electric propulsion system without aileron deflections, the solutions were computed for an angle-of-attack sweep from −2° to 20°, with the high-lift propellers blowing and the cruise propellers excluded from the simulation. To investigate the aileron effectiveness, the solutions were computed at angles of attack of −2° and 14°, for aileron deflections from −25° to 18° with the high-lift propellers blowing and the cruise propellers windmilling at idle-power. The high-lift propellers and the idle-power cruise propellers were modeled with an actuator disk. Results show negligible differences in lift, drag and pitching moment whether the idle-powered cruise propellers were included or excluded from the simulation. In general, the Kestrel and USM3D codes compared well for lift, drag and pitching moment for the landing configuration with no aileron control. The Kestrel code, using the Spalart-Allmaras turbulence model with the rotation corrections terms, predicted an increasing lift coefficient up to maximum lift coefficient of 4.65 at a 15° angle of attack. The USM3D code, using the Spalart-Allmaras turbulence model turbulence model with the Quadratic Constitutive Relation predicted an increasing lift coefficient up to maximum lift coefficient of 4.5 at a 12° angle of attack, with the lift remaining constant through a 14◦ angle of attack. A possible difference in lift coefficient between the codes for high angles of attack may result from the different available options used with the standard Spalart-Allmaras turbulence model. The Kestrel code predicted more drag across the range of angle of attack than the USM3D code. The codes compared well for pitching moment coefficient across the range of angle of attack. The Kestrel and USM3D codes compared well for aileron effectiveness at a 2° angle of attack. The Kestrel code predicts better aileron effectiveness than USM3D at a 14° angle of attack. There is more flow separation in the region outboard of the last high-lift nacelle for the USM3D solutions at a 14° angle of attack, which diminishes the ability of the aileron to be effective.

X-57 Maxwell Airplane↗

NASA NTRS2001

Coupled Growth in Hypermonotectics

The overall objective of this project is to obtain a fundamental understanding of the physics controlling solidification processes in immiscible alloy systems. The investigation involves both experimentation and the development of a model describing solidification in monotectic systems. The experimental segment was designed to first demonstrate that it is possible to obtain interface stability and steady state coupled growth in hypermonotectic alloys through microgravity processing. Microgravity results obtained to date have verified this possibility. Future flights will permit experimental determination of the limits of interface stability and the influence of alloy composition and growth rate on microstructure. The objectives of the modeling segment of the investigation include prediction of the limits of interface stability, modeling of convective flow due to residual acceleration, and the influence of surface tension driven flows at the solidification interface. The study of solidification processes in immiscible alloy systems is hindered by the inherent convective flow that occurs on Earth and by the possibility of sedimentation of the higher density immiscible liquid phase. It has been shown that processing using a high thermal gradient and a low growth rate can lead to a stable macroscopically planar growth front even in hypermonotectic alloys. Processing under these growth conditions can avoid constitutional supercooling and prevent the formation of the minor immiscible liquid phase in advance of the solidification front. However, the solute depleted boundary layer that forms in advance of the solidification front is almost always less dense than the liquid away from the solidification front. As a result, convective instability is expected. Ground based testing has indicated that convection is a major problem in these alloy systems and leads to gross compositional variations along the sample and difficulties maintaining interface stability. Sustained low gravity processing conditions are necessary in order to minimize these problems and obtain solidification conditions which approach steady state.

Andrews, J. Barry↗

NASA NTRS2004

Associative Flow Rule Used to Include Hydrostatic Stress Effects in Analysis of Strain-Rate-Dependent Deformation of Polymer Matrix Composites

designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and the strain-rate dependence of the composite response are due primarily to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. By applying micromechanics techniques along with given fiber properties, one can also determine the effects of the hydrostatic stresses in the polymer on the overall composite deformation response. First efforts to account for the hydrostatic stress effects in the composite deformation applied purely empirical methods that relied on composite-level data. In later efforts, to allow polymer properties to be characterized solely on the basis of polymer data, researchers at the NASA Glenn Research Center developed equations to model the polymers that were based on a non-associative flow rule, and efforts to use these equations to simulate the deformation of representative polymer materials were reasonably successful. However, these equations were found to have difficulty in correctly analyzing the multiaxial stress states found in the polymer matrix constituent of a composite material. To correct these difficulties, and to allow for the accurate simulation of the nonlinear strain-rate-dependent deformation analysis of polymer matrix composites, in the efforts reported here Glenn researchers reformulated the polymer constitutive equations from basic principles using the concept of an associative flow rule. These revised equations were characterized and validated in an experimental program carried out through a university grant with the Ohio State University, wherein tensile and shear deformation data were obtained for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. Tensile deformation data also were obtained over a variety of strain rates and fiber orientation angles for a representative polymer matrix composite composed using the polymer.

Goldberg, Robert K.↗

NASA NTRS1991

Stirling engine: Available tools for long-life assessment

A review is presented for the durability approaches applicable to long-time life assessment of Stirling engine hot-section components. The crucial elements are experimental techniques for generating long-time materials property data (both monotonic and cyclic flow and failure properties); analytic representations of slow strain rate material stress-strain response characteristics (monotonic and cyclic constitutive relations) at high temperatures and low stresses and strains; analytic creep-fatigue-environmental interaction life prediction methods applicable to long lifetimes at high temperatures and small stresses and strains; and experimental verification of life predictions. Long-lifetime design criteria for materials of interest are woefully lacking. Designing against failures due to creep, creep-rupture, fatigue, environmental attack, and creep-fatigue-environmental interaction will require considerable extrapolation. Viscoplastic constitutive models and time-temperature parameters will have to be calibrated for the hot-section materials of interest. Analysis combined with limited verification testing in a short-time regime will be required to build confidence in long-lifetime durability models.

Halford, Gary R.↗

NASA NTRS1991

Stirling engine - Available tools for long-life assessment

A review is presented for the durability approaches applicable to long-time life assessment of Stirling engine hot-section components. The crucial elements are experimental techniques for generating long-time materials property data (both monotonic and cyclic flow and failure properties); analytic representations of slow strain rate material stress-strain response characteristics (monotonic and cyclic constitutive relations) at high temperatures and low stresses and strains; analytic creep-fatigue-environmental interaction life prediction methods applicable to long lifetimes at high temperatures and small stresses and strains; and experimental verification of life predictions. Long-lifetime design criteria for materials of interest are woefully lacking. Designing against failures due to creep, creep-rupture, fatigue, environmental attack, and creep-failure-environmental interaction will require considerable extrapolation. Viscoplastic constitutive models and time-temperature parameters will have to be calibrated for the hot-section materials of interest. Analysis combined with limited verification testing in a short-time regime will be required to build confidence in long-lifetime durability models.

Halford, Gary R.↗

NASA NTRS1993

Physiological water model development

The water of the human body can be categorized as existing in two main compartments: intracellular water and extracellular water. The intracellular water consists of all the water within the cells and constitutes over half of the total body water. Since red blood cells are surrounded by plasma, and all other cells are surrounded by interstitial fluid, the intracellular compartment has been subdivided to represent these two cell types. The extracellular water, which includes all of the fluid outside of the cells, can be further subdivided into compartments which represent the interstitial fluid, circulating blood plasma, lymph, and transcellular water. The interstitial fluid surrounds cells outside of the vascular system whereas plasma is contained within the blood vessels. Avascular tissues such as dense connective tissue and cartilage contain interstitial water which slowly equilibrates with tracers used to determine extracellular fluid volume. For this reason, additional compartments are sometimes used to represent these avascular tissues. The average size of each compartment, in terms of percent body weight, has been determined for adult males and females. These compartments and the forces which cause flow between them are presented. The kidneys, a main compartment, receive about 25 percent of the cardiac output and filters out a fluid similar to plasma. The composition of this filtered fluid changes as it flows through the kidney tubules since compounds are continually being secreted and reabsorbed. Through this mechanism, the kidneys eliminate wastes while conserving body water, electrolytes, and metabolites. Since sodium accounts for over 90 percent of the cations in the extracellular fluid, and the number of cations is balanced by the number of anions, considering the renal handling sodium and water only should sufficiently describe the relationship between the plasma compartment and kidneys. A kidney function model is presented which has been adapted from a previous model of normal renal function in man. To test the validity of the proposed kidney model, results predicted by the model will be compared to actual data involving injected or ingested fluids and subsequent urine flow rates. Comparison of the model simulation to actual data following the ingestion of 1 liter of water is shown. The model simulation is also shown with actual data following the intravenous infusion of hypertonic saline.

Doty, Susan↗

NASA NTRS2018

Leverage Points for System Health Management of Autonomous Systems

Systems Health Management (SHM) is one of three basic functionalities that constitute an autonomous capability of a system. The other two functionalities are Planning & Scheduling, and Task Execution. In an autonomous system, variable autonomy is often distinct from variable authority to sense, decide, and act. There are quantifiable Levels of Autonomy that can be achieved by tuning different portions of the Observe-Orient-Decide-Act loop to provide flexibility and control. This approach is tabulated for multiple domains such as spacecraft and aerial vehicles. Examining SHM through a Systems Thinking lens helps us understand its stocks and flows, loops, and delays. Systems thinking, and modeling, is a useful way to understand change and complexity of systems of many types. There are certain archetypes that underlie well-known autonomy architectures. And there often are leverage points - best places to intervene in a system - that can resolve or mitigate some fundamental challenges in the design and deployment of autonomous systems. I identify these levers and present the ones that have been successfully used in NASA missions.

Systems Thinking↗

NASA NTRS1979

Visualization of turbulence in flow about a bluff body

A visual diagnostic investigation of turbulence in flow about a circular cylinder was conducted in order to gain a physical insight into the model for the structure of turbulence in flow around a bluff body advanced by the vorticity-amplification theory. A frame-by-frame qualitative and quantitative examination of selected movie strips over time intervals up to more than 800 ms was performed. The selective stretching of cross-vortex tubes, their streamwise tilting, the emergence of an organized turbulent flow pattern near the stagnation zone and the growth of a turbulent boundary layer were clearly revealed by the visualization. In particular, the visualization indicated that the cross-vortex tubes conveyed by the diverging stagnation flow constitute a coherent substructure within the overall turbulent flow that is triggered to its fullest manifestation by the stretching mechanism.

Sadeh, W. Z.↗