Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relaxation”

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,441 records · Page 80

Matrix dominated stress/strain behavior in polymeric composites: Effects of hold time, nonlinearity and rate dependency

In order to understand matrix dominated behavior in laminated polymer matrix composites, an elastic/viscoplastic constitutive model was developed and used to predict stress strain behavior of off-axis and angle-ply symmetric laminates under in-plane, tensile axial loading. The model was validated for short duration tests at elevated temperatures. Short term stress relaxation and short term creep, strain rate sensitivity, and material nonlinearity were accounted for. The testing times were extended for longer durations, and periods of creep and stress relaxation were used to investigate the ability of the model to account for long term behavior. The model generally underestimated the total change in strain and stress for both long term creep and long term relaxation respectively.

Gates, Thomas S.↗

Resolution of differences between collision number definitions in particle and continuum simulations

The use of the same symbol, Z, representing a 'collision number' for thermal relaxation, has led to confusion regarding its definition in the context of both continuum and particle simulations. Examination of the relaxation mechanics employed in particle simulations demonstrates that these definitions differ by a numerical factor that depends upon the intermolecular potential. Particle and continuum simulations employing appropriate definitions of Z lead to identical results during isothermal and adiabatic stationary relaxation.

Lumpkin, Forrest E., III↗

A nonlinear relaxation/quasi-Newton algorithm for the compressible Navier-Stokes equations

A highly efficient implicit method for the computation of steady, two-dimensional compressible Navier-Stokes flowfields is presented. The discretization of the governing equations is hybrid in nature, with flux-vector splitting utilized in the streamwise direction and central differences with flux-limited artificial dissipation used for the transverse fluxes. Line Jacobi relaxation is used to provide a suitable initial guess for a new nonlinear iteration strategy based on line Gauss-Seidel sweeps. The applicability of quasi-Newton methods as convergence accelerators for this and other line relaxation algorithms is discussed, and efficient implementations of such techniques are presented. Convergence histories and comparisons with experimental data are presented for supersonic flow over a flat plate and for several high-speed compression corner interactions. Results indicate a marked improvement in computational efficiency over more conventional upwind relaxation strategies, particularly for flowfields containing large pockets of streamwise subsonic flow.

Edwards, Jack R.↗

Development of a new model for vibration-dissociation coupling in nitrogen

A new model for vibration-dissociation coupling is derived using knowledge obtained from a simulation of the vibrational relaxation of heated nitrogen by collisional processes. This bilevel model is based on the assumption that nitrogen can be represented by two molecular species. The vibrational energy relaxation is represented by a modified form of the Landau-Teller equation. The effects of the bilevel model on a typical hypervelocity flowfield computation are examined and compared to the Park TTv model. Specifically, it is found that the bilevel model exhibits a decreased dissociation rate. This is believed to be due to the inhibition of dissociation by nonequilibrium vibrational relaxation.

Landrum, D. B.↗

Constraints on crustal rheology and age of deformation from models of gravitational spreading in Ishtar Terra, Venus

Gravitational spreading is expected to lead to rapid relaxation of high relief due to the high surface temperature and associated weak crust on Venus. In this study, we use new Magellan radar and altimetry data to determine the extent of gravitational relaxation in Ishtar Terra, which contains the highest relief on Venus as well as areas of extremely high topographic slope. Within Ishtar Terra the only mountain belts found on Venus, Akna, Danu, Freyja, and Maxwell Montes, nearly encircle the smooth, high (3-4 km) plateau of Lakshmi Planum. Finite-element models of this process give expected timescales for relaxation of relief and failure at the surface. From these modeling results we attempt to constrain the strength of the crust and timescales of deformation in Ishtar Terra. Below we discuss observational evidence for gravitational spreading in Ishtar Terra, results from the finite-element modeling, independent age constraints, and implications for the rheology and timing of deformation.

Smrekar, Suzanne E.↗

Gravitational spreading of high terrain in Ishtar Terra, Venus

A plane-strain viscoelastic finite element algorithm derived by Melosh and Rafesky (1980, 1983) is used to model the gravitational spreading of the high topography of Ishtar Terra. Brittle failure and relaxation of topographic relief are predicted as functions of time for ranges of crustal thickness and thermal gradient and for observed ranges in topographic elevation and slope. For a crustal thickness greater than 10 km, a thermal gradient of 15 K/km or more, and average values of relief (3 km) and slope (3 deg), the topography relaxes to 25 percent of its original height in less than 10 m.y. Although initial failure, in the form of normal faulting on the highlands and margin and occasional shallow thrusting in the lowlands, is predicted to occur much earlier, values of horizontal surface strain that are likely to be observable (about 1 percent) do not accumulate until significant relaxation of relief begins.

Smrekar, Suzanne E.↗

The H + OCS hot atom reaction - CO state distributions and translational energy from time-resolved infrared absorption spectroscopy

Time-resolved infrared diode laser spectroscopy has been used to probe CO internal and translational excitation from the reaction of hot H atoms with OCS. Product distributions should be strongly biased toward the maximum 1.4 eV collision energy obtained from 278 nm pulsed photolysis of HI. Rotations and vibrations are both colder than predicted by statistical density of states theory, as evidenced by large positive surprisal parameters. The bias against rotation is stronger than that against vibration, with measurable population as high as v = 4. The average CO internal excitation is 1920/cm, accounting for only 13 percent of the available energy. Of the energy balance, time-resolved sub-Doppler line shape measurements show that more than 38 percent appears as relative translation of the separating CO and SH fragments. Studies of the relaxation kinetics indicate that some rotational energy transfer occurs on the time scale of our measurements, but the distributions do not relax sufficiently to alter our conclusions. Vibrational distributions are nascent, though vibrational relaxation of excited CO is unusually fast in the OCS bath, with rates approaching 3 percent of gas kinetic for v = 1.

Nickolaisen, Scott L.↗

Simple Creep Test For Ceramic Fibers

Simple bend-stress-relaxation test yields information on creep-related properties of polycrystalline ceramic fibers. Determination of these properties important part of efforts to develop ceramic composite materials that retain mechanical strength and resistance to creep at high temperatures. Present test measures effects of time, temperature, and applied strain on creep-related relaxation of bend stress in ceramic fiber of almost any diameter in almost any environment, without need for contact sensors. Degree of relaxation of bend stress determined from radii of curvature.

Dicarlo, James A.↗

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on enhancement and application of the direct simulation Monte Carlo (DSMC) particle method for computing hypersonic flows of re-entry vehicles. Enhancement efforts dealt with modeling gas-gas interactions for thermal non-equilibrium relaxation processes and gas-surface interactions for prediction of vehicle surface temperatures. Both are important for application to problems of engineering interest. The code was employed in a parametric study to improve future applications, and in simulations of aeropass maneuvers in support of the Magellan mission. Detailed comparisons between continuum models for internal energy relaxation and DSMC models reveals that several discrepancies exist. These include definitions of relaxation parameters and the methodologies for implementing them in DSMC codes. These issues were clarified and all differences were rectified in a paper (Appendix A) submitted to Physics of Fluids A, featuring several key figures in the DSMC community as co-authors and B. Haas as first author. This material will be presented at the Fluid Dynamics meeting of the American Physical Society on November 21, 1993. The aerodynamics of space vehicles in highly rarefied flows are very sensitive to the vehicle surface temperatures. Rather than require prescribed temperature estimates for spacecraft as is typically done in DSMC methods, a new technique was developed which couples the dynamic surface heat transfer characteristics into the DSMC flow simulation code to compute surface temperatures directly. This model, when applied to thin planar bodies such as solar panels, was described in AIAA Paper No. 93-2765 (Appendix B) and was presented at the Thermophysics Conference in July 1993. The paper has been submitted to the Journal of Thermophysics and Heat Transfer. Application of the DSMC method to problems of practical interest requires a trade off between solution accuracy and computational expense and limitations. A parametric study was performed and reported in AIAA Paper No. 93-2806 (Appendix C) which assessed the accuracy penalties associated with simulations of varying grid resolution and flow domain size. The paper was also presented at the Thermophysics Conference and will be submitted to the journal shortly. Finally, the DSMC code was employed to assess the pitch, yaw, and roll aerodynamics of the Magellan spacecraft during entry into the Venus atmosphere at off-design attitudes. This work was in support of the Magellan aerobraking maneuver of May 25-Aug. 3, 1993. Furthermore, analysis of the roll characteristics of the configuration with canted solar panels was performed in support of the proposed 'Windmill' experiment. Results were reported in AIAA Paper No. 93-3676 (Appendix D) presented at the Atmospheric Flight Mechanics Conference in August 1993, and were submitted to Journal of Spacecraft and Rockets.

Haas, Brian L.↗

Modelling of Surfaces. Part 2: Metallic Alloy Surfaces Using the BFS Method

Using BFS, a new semiempirical method for alloys, we study the surface structure of fcc ordered binary alloys. We concentrate on the calculation of surface energies and surface relaxations for the L1(sub 0) and L1(sub 2) ordered structures. Different terminations of the low-index faces are studied. Also, we present results for the interlayer relaxations for planes close to the surface, revealing different relaxations for atoms of different species producing a rippled surface layer.

Bozzolo, Guillermo↗

Time-domain approach for the transient responses in stratified viscoelastic Earth models

We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Hanyk, L.↗

Research activity at the shock tube facility at NASA Ames

The real gas phenomena dominate the relaxation process occurring in the flow around hypersonic vehicles. The air flow around these vehicles undergoes vibrational excitation, chemical dissociation, and ionization. These chemical and kinetic phenomena absorb energy, change compressibility, cause temperature to fall, and density to rise. In high-altitude, low density environments, the characteristic thicknesses of the shock layers can be smaller than the relaxation distances required for the gas to attain chemical and thermodynamic equilibrium. To determine the effects of chemical nonequilibrium over a realistic hypersonic vehicle, it would be desirable to conduct an experiment in which all aspects of fluid flow are simulated. Such an experiment is extremely difficult to setup. The only practical alternative is to develop a theoretical model of the phenomena and to compute the flow around the vehicle including the chemical nonequilibrium, and compare the results with the experiments conducted in the facilities under conditions where only a portion of the flow phenomena is simulated. Three types of experimental data are needed to assist the aerospace community in this model development process: (1) data which will enhance our phenomenological understanding of the relaxation process, (2) data on rate reactions for the relevant reactions, and (3) data on bulk properties, such as spectral radiation emitted by the gas, for a given set of aerodynamic conditions. NASA Ames is in a process of collecting such data by simulating the required aerothermochemical conditions in an electric arc driven shock tube.

Sharma, Surendra P.↗

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

The Extensional Rheology of Non-Newtonian Materials

The evolution of the transient extensional stresses in dilute and semi-dilute viscoelastic polymer solutions are measured with a filament stretching rheometer of a design similar to that first introduced by Sridhar, et al. The solutions are polystyrene-based (PS) Boger fluids that are stretched at constant strain rates ranging from 0.6 less than or equal to epsilon(0) less than or equal to 4s(exp -1) and to Hencky strains of epsilon greater than 4. The test fluids all strain harden and Trouton ratios exceeding 1000 are obtained at high strains. The experimental data strain hardens at lower strain levels than predicted by bead-spring FENE models. In addition to measuring the transient tensile stress growth, we also monitor the decay of the tensile viscoelastic stress difference in the fluid column following cessation of uniaxial elongation as a function of the total imposed Hencky strain and the strain rate. The extensional stresses initially decay very rapidly upon cessation of uniaxial elongation followed by a slower viscoelastic relaxation, and deviate significantly from FENE relaxation predictions. The relaxation at long times t is greater than or equal to 5 s, is compromised by gravitational draining leading to non-uniform filament profiles. For the most elastic fluids, partial decohension of the fluid filament from the endplates of the rheometer is observed in tests conducted at high strain rates. This elastic instability is initiated near the rigid endplate fixtures of the device and it results in the progressive breakup of the fluid column into individual threads or 'fibrils' with a regular azimuthal spacing. These fibrils elongate and bifurcate as the fluid sample is elongated further. Flow visualization experiments using a modified stretching device show that the instability develops as a consequence of an axisymmetry-breaking meniscus instability in the nonhomogeneous region of highly deformed fluid near the rigid endplate.

Spiegelberg, Stephen H.↗

Comparison of Nonequilibrium Solution Algorithms Applied to Chemically Stiff Hypersonic Flows

Three solution algorithms, explicit under-relaxation, point implicit, and lower-upper symmetric Gauss-Seidel, are used to compute nonequilibrium flow around the Apollo 4 return capsule at the 62-km altitude point in its descent trajectory. By varying the Mach number, the efficiency and robustness of the solution algorithms were tested for different levels of chemical stiffness.The performance of the solution algorithms degraded as the Mach number and stiffness of the flow increased. At Mach 15 and 30, the lower-upper symmetric Gauss-Seidel method produces an eight order of magnitude drop in the energy residual in one-third to one-half the Cray C-90 computer time as compared to the point implicit and explicit under-relaxation methods. The explicit under-relaxation algorithm experienced convergence difficulties at Mach 30 and above. At Mach 40 the performance of the lower-upper symmetric Gauss-Seidel algorithm deteriorates to the point that it is out performed by the point implicit method. The effects of the viscous terms are investigated. Grid dependency questions are explored.

Palmer, Grant↗

Statistics of Experiments on Cluster Formation and Transport in a Gravitational Field

Metastable state relaxation in a gravitational field is investigated in the case of non-critical binary solutions. A relaxation description is presented in terms of the time-dependent Ginzburg-Landau formalism for a non-conserved order parameter. A new ansatz for solution of the corresponding partial nonlinear stochastic differential equation is discussed. It is proved that, for the supersaturated solution under consideration, the metastable state relaxation in a gravitational field leads to formation of solute concentration gradients due to the sedimentation of subcritical solute clusters. The pure discussion of the possible methods to compare theoretical results and experimental data related to solute sedimentation in a gravitational field is presented. It is shown that in order to describe these experiments it is necessary to deal both with the value of the solute concentration gradient and with its formation rate. The stochastic nature of the sedimentation process is shown.

Izmailov, Alexander F.↗

Laboratory Kinetic Studies of OH and CO2 Relevant to Upper Atmospheric Radiation Balance

The purpose of this project was to quantify the rates of two processes which are crucial to our understanding of radiative energy balance in the upper atmosphere. The first process is radiative emission from vibrationally hot OH radicals following the H + O3 reaction in the upper mesosphere. The importance of this process depends strongly on the OH radiative emission coefficients. Our goal was to measure the OH permanent dipole moment in excited vibrational states and to use these measurements to construct an improved OH dipole moment function and improved radiative emission coefficients. Significant progress was made on these experiments including the construction of a supersonic jet source for vibrationally excited OH radicals. Unfortunately, our efforts to transport the OH radicals into a second lower pressure vacuum chamber were not successful, and we were unable to make improved dipole moment measurements for OH. The second key kinetic process which we attempted to quantify during this project is the rate of relaxation of bend-excited CO2 by oxygen atoms. Since excitation of the bending vibrational mode of CO2 is the major cooling mechanism in the upper mesosphere/lower thermosphere, the cooling rate of this region depends crucially on the rate of energy transfer out of this state. It is believed that the most efficient transfer mechanism is via atomic oxygen but the rate for this process has not been directly measured in the laboratory at appropriate temperatures and even the room temperature rate remains controversial. We attempted to directly measure the relaxation rate Of CO2 (010) by oxygen atoms using the discharge flow technique. This experiment was set up at Aerodyne Research. Again, significant progress was achieved in this experiment. A hot CO2 source was set up, bend excited CO2 was detected and the rate of relaxation of bend excited CO2 by He atoms was measured. Unfortunately, the project ran out of time before the oxygen atom kinetic studies could be implemented.

Nelson, David D.↗

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗