Results for “constitutive flow model” — Page 6

NASA NTRSDate not supplied

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in the NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency. In contrast to the legacy FUN3D iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear Courant-Friedrichs-Lewy number (CFL) ramping, HANIM is based upon a hierarchy of modules including preconditioner, generalized conjugate residual, realizability check, nonlinear control, and CFL adaption modules. HANIM performance is systematically compared with the performance of the legacy solver of FUN3D and a baseline solver based on a preconditioner alone. Iterative solutions are compared for three benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-fuselage juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmaras model and its nonlinear extension based on quadratic constitutive relations.

CFD↗

NASA NTRS1991

On the thermodynamic framework of generalized coupled thermoelastic-viscoplastic-damage modeling

A complete potential based framework using internal state variables is put forth for the derivation of reversible and irreversible constitutive equations. In this framework, the existence of the total (integrated) form of either the (Helmholtz) free energy or the (Gibbs) complementary free energy are assumed a priori. Two options for describing the flow and evolutionary equations are described, wherein option one (the fully coupled form) is shown to be over restrictive while the second option (the decoupled form) provides significant flexibility. As a consequence of the decoupled form, a new operator, i.e., the Compliance operator, is defined which provides a link between the assumed Gibb's and complementary dissipation potential and ensures a number of desirable numerical features, for example the symmetry of the resulting consistent tangent stiffness matrix. An important conclusion reached, is that although many theories in the literature do not conform to the general potential framework outlined, it is still possible in some cases, by slight modifications of the used forms, to restore the complete potential structure.

Arnold, S. M.↗

NASA NTRS1977

COMOC 2: Two-dimensional aerodynamics sequence, computer program user's guide

The COMOC finite element fluid mechanics computer program system is applicable to diverse problem classes. The two dimensional aerodynamics sequence was established for solution of the potential and/or viscous and turbulent flowfields associated with subsonic flight of elementary two dimensional isolated airfoils. The sequence is constituted of three specific flowfield options in COMOC for two dimensional flows. These include the potential flow option, the boundary layer option, and the parabolic Navier-Stokes option. By sequencing through these options, it is possible to computationally construct a weak-interaction model of the aerodynamic flowfield. This report is the user's guide to operation of COMOC for the aerodynamics sequence.

Manhardt, P. D.↗

NASA NTRSDate not supplied

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency of FUN3Dsolutions. In contrast to the baseline iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear CFL ramping, HANIM is based upon a hierarchy of modules including pre conditioner, generalized conjugate residual, realizability check, nonlinear control,and CFL adaption modules. HANIM performance is systematically compared with the performance of the baseline solver. The iterative solutions are compared for three aerodynamic benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-body juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmar as model and its nonlinear extension based on quadratic constitutive relations.

Li Wang↗

NASA NTRS1984

The rheology of icy satellites

High-temperature creep in orthoenstatite under conditions of controlled oxygen fugacity was studied. It was found that creep was conttrolled by the extremely thin layer of SiO2 which wetted the grain boundaries. Slight reduction of the (Mg, Fe)SiO3 enstatite during hot pressing produced microscopic particles of Fe and the thin film of intergranular SiO2. This result highlights another complication in determining the flow properties of iron bearing silicates which constitute the bulk of terrestrial planets and moons. The Phenomenon may be important in the ductile formation of any extraterrestrial body which is formed in a reducing environment. The rheology of dirty ice was studied. This involves micromechanical modeling of hardening phenomena due to contamination by a cosmic distribution of silicate particles. The larger particles are modeled by suspension theory. In order to handle the distribution of particles sizes, the hardening is readed as a critical phenomenon, and real space renormalization group techniques are used. Smaller particles interact directly with the dislocations. The particulate hardening effect was studied in metals. The magnitude of such hardening in ice and the defect chemistry of ice are studied to assess the effects of chemical contamination by methane, ammonia, or other likely contaminants.

Sammis, C. G.↗

NASA NTRS1999

Development and Validation of a Supersonic Helium-Air Coannular Jet Facility

Data are acquired in a simple coannular He/air supersonic jet suitable for validation of CFD (Computational Fluid Dynamics) codes for high speed propulsion. Helium is employed as a non-reacting hydrogen fuel simulant, constituting the core of the coannular flow while the coflow is composed of air. The mixing layer interface between the two flows in the near field and the plume region which develops further downstream constitute the primary regions of interest, similar to those present in all hypersonic air breathing propulsion systems. A computational code has been implemented from the experiment's inception, serving as a tool for model design during the development phase.

Carty, Atherton A.↗

NASA NTRS2016

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed.

Impact↗

NASA NTRS2016

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed

Impact↗

NASA NTRS2013

Material Ignition and Suppression Test (MIST) in Space Exploration Atmospheres, Summary of Research

The Material Ignition and Suppression Test (MIST) project has had the objective of evaluating the ease of ignition and the fire suppression of materials used in spacecraft under environmental condition expected in a spacecraft. For this purpose, an experimental and theoretical research program is being conducted on the effect of space exploration atmospheres (SEA) on the piloted ignition of representative combustible materials, and on their fire suppression characteristics. The experimental apparatus and test methodology is derived from the Forced Ignition and Flame Spread Test (FIST), a well-developed bench scale test designed to extract material properties relevant to prediction of material flammability. In the FIST test, materials are exposed to an external radiant flux and the ignition delay and critical mass flux at ignition are determined as a function of the type of material and environmental conditions. In the original MIST design, a small-scale cylindrical flow duct with fuel samples attached to its inside wall was heated by a cylindrical heater located at the central axis of the cylinder. However, as the project evolved it was decided by NASA that it would be better to produce an experimental design that could accommodate other experiments with different experimental concepts. Based on those instructions and input from the requirements of other researchers that may share the hardware in an ISS/CIR experiment, a cylindrical design based on placing the sample at the center of an optically transparent tube with heaters equally spaced along the exterior of the cylinder was developed. Piloted ignition is attained by a hot wire igniter downstream of the fuel sample. Environment variables that can be studied via this experimental apparatus include: external radiant flux, oxidizer oxygen concentration, flow velocity, ambient pressure, and gravity level (if flown in the ISS/CIR). This constitutes the current experimental design, which maintains fairly good consistency with Dr Tien's and Dr Olson's project approaches. A further goal of the project has been to develop a combined solid/gas phase numerical model based on the MIST test methodology to predict the flammability behavior of practical materials in spacecraft.

flammability↗

NASA NTRS2003

Implementation of an Associative Flow Rule Including Hydrostatic Stress Effects Into the High Strain Rate Deformation Analysis of Polymer Matrix Composites

A previously developed analytical formulation has been modified in order to more accurately account for the effects of hydrostatic stresses on the nonlinear, strain rate dependent deformation of polymer matrix composites. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical J2 plasticity theory definitions of effective stress and effective inelastic strain, along with the equations used to compute the components of the inelastic strain rate tensor, are appropriately modified. To verify the revised formulation, the shear and tensile deformation of two representative polymers are computed across a wide range of strain rates. Results computed using the developed constitutive equations correlate well with experimental data. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite for several fiber orientation angles across a variety of strain rates. The computed values compare well to experimentally obtained results.

Goldberg, Robert K.↗

NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗

NASA NTRS1994

Particle Distribution in Suspension Shear Flow

An earlier proposed constitutive relation for normal stresses originated by random particle fluctuations is used to describe a joint effect of thermal and shear-induced fluctuations on concentrational distributions in suspension flow. Averaged products of components of the fluctuation velocity are evaluated on a basis of the rational mechanics approach combined with a simple kinematic consideration. The equation of momentum conservation of the dispersed phase of a suspension closed with this constitutive relation is applied to unidirectional shear flow in the gravity field and to rotational Couette flow. Coupling of the thermal and shear-induced fluctuations results in that the ability of shear flow to suspend particles has a minimum at a certain particle size, all other things being equal. The developed model provides also for a reasonable explanation of particle distributions observed in Couette flow. The approach based on the consideration of momentum balance for the dispersed phase is proved to lead to an effective equation of convective diffusion of the suspended particles. Coefficients of mutual diffusion due to both thermal and shear-induced fluctuations are drastically different from corresponding self-diffusivities as regards both their scaling and their concentrational dependence.

Buyevich, Yu A.↗

NASA NTRS2024

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

NASA NTRS2023

Occurrence of Large Geomagnetically Induced Currents Within the EPRI SUNBURST Monitoring Network

Space weather, a natural hazard, can adversely impact human technological assets. High-voltage electric power transmission grids constitute one of the most critical technological systems vulnerable to space weather driven geomagnetically induced currents (GICs). One of the major challenges pertaining to the study of GICs over the continental United States has been the availability of GIC measurements, which are critical for validation of geoelectric field and power flow models, for example. In this study, we analyze GIC measurements collected at 17 Electrical Power Research Institute (EPRI) SUNBURST transformer locations across the United States for which a GIC value of 10 A or greater was recorded. This data set includes 52 individual geomagnetic storms with Kp index 6 and above during the period from 2010 to 2021. The analysis confirms that there is a good correlation between the number of geomagnetic storms per year and the number of recorded GIC events. Our results also show that about 76% of the top 17 GIC events are associated with the storm main phase, while only 24% are attributed to storm sudden commencements. In addition, it is shown, for the first time, that mid-latitude positive bays can cause large GICs over the continental United States. Finally, this study shows that the largest measured GIC event in the data set was associated with a localized intense dB/dt structure, which could be attributed to substorm activity.

Chigomezyo M. Ngwira↗

NASA NTRS1988

Hardware development for the surface tension driven convection experiment aboard the USML-1 spacelab mission

The Surface Tension Driven Convection Experiment is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the USML-1 Spacelab mission planned for March 1992. Hardware is under development to establish the experimental conditions and perform the specified measurements, for both ground based research and the flight experiment in a Spacelab single rack. Major development areas include an infrared thermal imaging system for surface temperature measurement, a CO2 laser and control system for surface heating, and for flow visualization, a He-Ne laser and optical system in conjunction with an intensified video camera. For ground based work the components of each system were purchased or designed, and tested individually. The three systems will be interfaced with the balance of the experimental hardware and will constitute a working engineering model. A description of the three systems and examples of the component performance is given along with the plans for the development of flight hardware.

Pline, A. D.↗

NASA NTRS1989

Hardware development for the Surface Tension Driven Convection Experiment aboard the USML-1 Spacelab mission

The Surface Tension Driven Convection Experiment is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the USML-1 Spacelab mission planned for March 1992. Hardware is under development to establish the experimental conditions and perform the specified measurements, for both ground based research and the flight experiment in a Spacelab single rack. Major development areas include an infrared thermal imaging system for surface temperature measurement, a CO2 laser and control system for surface heating, and for flow visualization, a He-Ne laser and optical system in conjunction with an intensified video camera. For ground based work the components of each system were purchased or designed, and tested individually. The three systems will be interfaced with the balance of the experimental hardware and will constitute a working engineering model. A description of the three systems and examples of the component performance is given along with the plans for the development of flight hardware.

Pline, A. D.↗

NASA NTRS1992

An annual model of SSM/I radiobrightness for dry soil

An annual model is presented of the temperature structure within a homogeneous, dry soil halfspace that is subject to both diurnal and annual insolation, radiant heating from the atmosphere, sensible heat exchange with the atmosphere, and radiant cooling. The thermal constitutive properties of the soil are assumed to be constant so that the heat flow equation can be solved analytically. For computational economy, a variable time interval Laplace transform method is developed to predict the temperature.

Liou, Yuei-An↗

NASA NTRS1994

A New Reynolds Stress Algebraic Equation Model

A general turbulent constitutive relation is directly applied to propose a new Reynolds stress algebraic equation model. In the development of this model, the constraints based on rapid distortion theory and realizability (i.e. the positivity of the normal Reynolds stresses and the Schwarz' inequality between turbulent velocity correlations) are imposed. Model coefficients are calibrated using well-studied basic flows such as homogeneous shear flow and the surface flow in the inertial sublayer. The performance of this model is then tested in complex turbulent flows including the separated flow over a backward-facing step and the flow in a confined jet. The calculation results are encouraging and point to the success of the present model in modeling turbulent flows with complex geometries.

Shih, Tsan-Hsing↗