SEARCH · Search NASA
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.
Elastic wave propagation, deriving characteristic equations in generalized curvilinear coordinates using linear elastic, isotropic and homogeneous constitutive equations
A steady-state analysis is conducted to examine the basic flow structure of a non-Newtonian fluid in a domain including an inflow region, a contraction region, and an outflow region. A Cartesian grid system is used throughout the entire flow domain, including the contraction region, thus creating an irregular grid cell structure adjacent to the curved boundary. At node points adjacent to the curved boundary symmetry conditions are derived for the different flow variables in order to solve the governing difference equations. Attention is given to the motion and non-Newtonian constitutive equations, the boundary conditions, the numerical modeling of the non-Newtonian equations, the stream function contour lines for the non-Newtonian fluid, the vorticity contour lines for the non-Newtonian fluid, the velocity profile across the contraction, and the shear stress contour lines for the non-Newtonian fluid.
Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.
A new hybrid stress finite element algorithm suitable for analyses of large quasistatic deformation of inelastic solids is presented. Principal variables in the formulation are the nominal stress rate and spin. The finite element equations which result are discrete versions of the equations of compatibility and angular momentum balance. Consistent reformulation of the constitutive equation and accurate and stable time integration of the stress are discussed at length. Examples which bring out the feasibility and performance of the algorithm conclude the work.
The advent of large deformation finite element methodologies is beginning to permit the numerical simulation of hot working processes whose design until recently has been based on prior industrial experience. Proper application of such finite element techniques requires realistic constitutive equations which more accurately model material behavior during hot working. A simple constitutive model for hot working is the single scalar internal variable model for isotropic thermal elastoplasticity proposed by Anand. The model is recalled and the specific scalar functions, for the equivalent plastic strain rate and the evolution equation for the internal variable, presented are slight modifications of those proposed by Anand. The modified functions are better able to represent high temperature material behavior. The monotonic constant true strain rate and strain rate jump compression experiments on a 2 percent silicon iron is briefly described. The model is implemented in the general purpose finite element program ABAQUS.
Characteristic equations are derived in generalized curvilinear coordinates. Linear elastic, isotropic, and homogeneous constitutive equations have been used in the derivation. The generalized characteristic equations readily lend themselves to any requirements of space dimension and geometry. A simple boundary value problem is solved to indicate the applicability of these equations.
Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).
We present a theoretical stability analysis for an expanding accretion shock that does not involve a rarefaction wave behind it. The dispersion equation that determines the eigenvalues of the problem and the explicit formulae for the corresponding eigenfunction profiles are presented for an arbitrary equation of state and finite-strength shocks. For spherically and cylindrically expanding steady shock waves, we demonstrate the possibility of instability in a literal sense, a power-law growth of shock-front perturbations with time, in the range of $h_c< h<1+2 {\mathcal {M}}_2$ , where $h$ is the D'yakov-Kontorovich parameter, $h_c$ is its critical value corresponding to the onset of the instability and ${\mathcal {M}}_2$ is the downstream Mach number. Shock divergence is a stabilizing factor and, therefore, instability is found for high angular mode numbers. As the parameter $h$ increases from $h_c$ to $1+2 {\mathcal {M}}_2$ , the instability power index grows from zero to infinity. This result contrasts with the classic theory applicable to planar isolated shocks, which predicts spontaneous acoustic emission associated with constant-amplitude oscillations of the perturbed shock in the range $h_c< h<1+2 {\mathcal {M}}_2$ . Examples are given for three different equations of state: ideal gas, van der Waals gas and three-terms constitutive equation for simple metals.
A state-of-the-art review of applicable constitutive models with selection of two for detailed comparison with a wide range of experimental tests was conducted. The experimental matrix contained uniaxial and biaxial tensile, creep, stress relaxation, and cyclic fatigue tests at temperatures to 1093 C and strain rates from .0000001 to .001/sec. Some nonisothermal cycles will also be run. The constitutive models will be incorporated into the MARC finite element structural analysis program with a demonstration computation made for advanced turbine blade configuration. In the code development work, particular emphasis is being placed on developing efficient integration algorithms for the highly nonlinear and stiff constitutive equations. Another area of emphasis is the appropriate and efficient methodology for determing constitutive constants from a minimum extent of experimental data.
We are concerned with the design of inhomogeneous, all dielectric (lossless) periodic structures which act as filters. Dielectric filters made as stacks of inhomogeneous gratings and layers of materials are being used in optical technology, but are not common at microwave frequencies. The problem is then finding the periodic cell's geometric configuration and permittivity values which correspond to a specified reflectivity/transmittivity response as a function of frequency/illumination angle. This type of design can be thought of as an inverse-source problem, since it entails finding a distribution of sources which produce fields (or quantities derived from them) of given characteristics. Electromagnetic sources (electric and magnetic current densities) in a volume are related to the outside fields by a well known linear integral equation. Additionally, the sources are related to the fields inside the volume by a constitutive equation, involving the material properties. Then, the relationship linking the fields outside the source region to those inside is non-linear, in terms of material properties such as permittivity, permeability and conductivity. The solution of the non-linear inverse problem is cast here as a combination of two linear steps, by explicitly introducing the electromagnetic sources in the computational volume as a set of unknowns in addition to the material unknowns. This allows to solve for material parameters and related electric fields in the source volume which are consistent with Maxwell's equations. Solutions are obtained iteratively by decoupling the two steps. First, we invert for the permittivity only in the minimization of a cost function and second, given the materials, we find the corresponding electric fields through direct solution of the integral equation in the source volume. The sources thus computed are used to generate the far fields and the synthesized triter response. The cost function is obtained by calculating the deviation between the synthesized value of reflectivity/transmittivity and the desired one. Solution geometries for the periodic cell are sought as gratings (ensembles of columns of different heights and widths), or combinations of homogeneous layers of different dielectric materials and gratings. Hence the explicit unknowns of the inversion step are the material permittivities and the relative boundaries separating homogeneous parcels of the periodic cell.
Single crystal materials are being used in gas turbine airfoils and are candidates for other hot section components because of their increased temperature capabilities and resistance to thermal fatigue. Development of a constitutive model which assesses the inelastic behavior of these materials has been studied in 2 NASA programs: Life Prediction and Constitutive Models for Engine Hot Section Anisotropic Materials and Biaxial Constitutive Equation Development for Single Crystals. The model has been fit to a large body of constitutive data for single crystal PWA 1480 material. The model uses a unified approach for computing total inelastic strains (creep plus plasticity) on crystallographic slip systems reproducing observed directional and strain rate effects as a natural consequence of the summed slip system quantities. The model includes several of the effects that have been reported to influence deformation in single crystal materials, such as shear stress, latent hardening, and cross slip. The model is operational in a commercial Finite Element code and is being installed in a Boundary Element Method code.
A new high-fidelity micromechanics model has been developed under funding from the NASA Glenn Research Center for predicting the response of multiphase materials with arbitrary periodic microstructures. The model's analytical framework is based on the homogenization technique, but the method of solution for the local displacement and stress fields borrows concepts previously employed in constructing the higher order theory for functionally graded materials. The resulting closed-form macroscopic and microscopic constitutive equations, valid for both uniaxial and multiaxial loading of periodic materials with elastic and inelastic constitutive phases, can be incorporated into a structural analysis computer code. Consequently, this model now provides an alternative, accurate method.
The objectives of the research were as follows: (1) Extend the Representative Elementary Area (RE) concept, first proposed and developed in Wood et al, (1988), to the water balance fluxes of the interstorm period (redistribution, evapotranspiration and baseflow) necessary for the analysis of long-term water balance processes. (2) Derive spatially averaged water balance model equations for spatially variable soil, topography and vegetation, over A RANGE OF CLIMATES. This is a necessary step in our goal to derive consistent hydrologic results up to GCM grid scales necessary for global climate modeling. (3) Apply the above macroscale water balance equations with remotely sensed data and begin to explore the feasibility of parameterizing the water balance constitutive equations at GCM grid scale.
The objective was to develop unified constitutive equations which can model a variety of nonlinear material phenomena observed in Rene 80 at elevated temperatures. A constitutive model was developed based on back stress and drag stress. The tensorial back stress was used to model directional effects; whereas, the scalar drag stress was used to model isotropic effects and cyclic hardening or softening. A flow equation and evolution equations for the state variables were developed in multiaxial form. Procedures were developed to generate the material parameters. The model predicted very well the monotonic tensile, cyclic, creep, and stress relaxation behavior of Rene 80 at 982 C. The model was then extended to 871, 760, and 538 C. It was shown that strain rate dependent behavior at high temperatures and strain rate independent behavior at the lower temperatures could be predicted very well. A large number of monotonic tensile, creep, stress relation, and cyclic experiments were predicted. The multiaxial capabilities of the model were verified extensively for combined tension/torsion experiments. The prediction of the model agreed very well for proportional, nonproportional, and pure shear cyclic loading conditions at 982 and 871 C.
In the present treatment of a viscoplastic composite wafer formed from a viscoplastic matrix that is reinforced by elastic or viscoplastic fibers attached to its plane surfaces, one constitutive equation is established by considering the viscoplastic behavior of the matrix as determined by an integral-type constitutive law whose intrinsic time-measure is pertinent to endochronic viscoplasticity. Attention is given to asymptotic cases where fiber stiffnesses and the hydrostatic modulus of the wafer are much larger than the latter's shear modulus. An explicit calculation is used when the stress field is uniaxial.
A tensor formulation of the equation of radiative transfer is derived in a seven-dimensional Riemannian space such that the resulting equation constitutes a divergence in any coordinate system. After being transformed to a spherically symmetric comoving coordinate system, the transfer equation contains partial derivatives in angle and frequency, as well as optical depth due to the effects of aberration and the Doppler shift. However, by virtue of the divergence form of this equation, the divergence theorem may be applied to yield a numerical differencing scheme which is expected to be stable and to conserve luminosity. It is shown that the equation of transfer derived by this method in a Lagrangian coordinate system may be reduced to that given by Castor (1972), although it is, of course, desirable to leave the equation in divergence form.
A viscous-inviscid interaction algorithm is developed for prediction of two-dimensional mean and fluctuating velocity distributions in the wake immediately downstream of an airfoil trailing edge. A composite pressure field is defined, and a Poisson equation solved for transverse pressure variations. A parabolized form of the time-averaged steady Navier-Stokes equations are solved in conjunction with a viscous-augmented two-dimensional inviscid potential flow analysis. A tensor constitutive equation is employed to predict Reynolds stress distributions from solutions of a turbulence kinetic energy two equation closure model. Numerical predictions compared favorably with detailed experimental data for mean and fluctuating velocities, and Reynolds shear stress distributions, in the trailing edge region of a NACA 63-012 airfoil.
Direct Boundary Element formulations and their numerical implementation for periodic and transient elastic as well as inelastic transient dynamic analyses of two-dimensional, axisymmetric and three-dimensional solids are presented. The inelastic formulation is based on an initial stress approach and is the first of its kind in the field of Boundary Element Methods. This formulation employs the Navier-Cauchy equation of motion, Graffi's dynamic reciprocal theorem, Stokes' fundamental solution, and the divergence theorem, together with kinematical and constitutive equations to obtain the pertinent integral equations of the problem in the time domain within the context of the small displacement theory of elastoplasticity. The dynamic (periodic, transient as well as nonlinear transient) formulations have been applied to a range of problems. The numerical formulations presented here are included in the BEST3D and GPBEST systems.