Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Time Series”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗
NASA NTRS1987

Plume formation in the D-double prime-layer and the roughness of the core-mantle boundary

A series of very-high-resolution finite element calculations of plume formation in the D-double prime-layer has been performed for several plausible rheologies and boundary conditions in order to study both the early and later stages of boundary layer development. The results show that plumes are initiated by coalescence of small-scale convective instabilities within the low-viscosity region immediately above the core-mantle boundary (CMB). These instabilities support topographic roughness on the CMB having horizontal scales of 20-50 km and provide a source for scattered P-waves seen as precursors to the phases PKIKP and PKKP. The calculated structure of fully developed plumes emerging from the D-double prime-layer consists of 5-50 cm/yr flow confined to 50-100 km thick vertical conduits. With strongly temperature-dependent viscosity, plumes exhibit time-dependent behavior, including upward propagating solitary conduit waves, which may contribute to episodicity in hotspot volcanism.

Olson, Peter↗
NASA NTRS1995

Application of Space Shuttle photography to studies of upper ocean dynamics

Three studies have been conducted using space shuttle imagery to explain the dynamics behavior of internal waves in the Atlantic and Indian Oceans and to derive tide-related parameters for Delaware Bay. By interpreting space shuttle photographs taken during mission STS-40, a total of 34 internal wave packets on the continental shelf of the Middle Atlantic Bight have been recognized. Using the finite-depth theory we derived that the maximum amplitude of solitons is 5.6 m, the phase speed 0.42 m/s, and the period 23.8 min. Deep-ocean internal waves in the western equatorial Indian Ocean on photographs taken during mission STS-44 were also interpreted and analyzed. The internal waves occurred in the form of a multisoliton packet in which there are about a dozen solitons. The average wavelength of the solitons is 1.8 +/- 0.5 km. The crest lines are mostly straight and reach as long as 100 km. The distance between two adjacent packets is about 66 km. Using the deepwater soliton theory, we derived that the mean amplitude of the solitons is 25 m, the nonlinear phase speed 1.7 m/s, and the average period 18 min. For both cases, the semidiural tides are the principal generating mechanism. The tide-related parameters of Delaware Bay were derived from space shuttle time-series photographs taken during mission STS-40. The water area in the bay were measured from interpretation maps of the photographs. The corresponding tidal levels were calculated using the exposure time. From these data, an approximate function relating the water area to the tidal level at a reference point was determined. Then, the water areas of the Delaware Bay at mean high water (MHW) and mean low water (MLW), below 0 m, for the tidal zone, and the tidal flux were inferred. All parameters derived were reasonable and compared well with results of previous investigations.

Zheng, Quanan↗
NASA NTRS2025

Single Integrated Flux Files to Account for Spacecraft Attitude Motion in Meteoroid Risk Assessments Using the NASA MEM3 and Bumper Codes

Spacecraft meteoroid risk assessments can be performed using NASA’s Meteoroid Engineering Model 3 (MEM3) and the Bumper risk assessment code. The MEM3 meteoroid flux is highly directional in heliocentric space. To represent the directionality, MEM3 outputs the flux in commonly used coordinate frames. For example, if a spacecraft attitude is fixed in the orbital tangential frame, a single Bumper run using the MEM3 flux file in the body-fixed frame can be performed. Cases exist in which the spacecraft does not remain fixed in one of these frames during a mission. This could be handled by performing a series of Bumper runs at points along the trajectory. For each trajectory point, the spacecraft has to be rotated into the frame in which the MEM3 flux file is defined. Another option presented here is to generate a single integrated flux (SIF) file. A SIF file is an average over a trajectory of the time varying MEM3 flux files rotated from their native frame into the spacecraft bus frame. A SIF file enables accounting for spacecraft rotation (coupled with flux variation over the trajectory) in a single Bumper run. An example of an application of SIF files is the meteoroid risk assessment for the Mars Sample Return (MSR) campaign Capture Containment and Return System (CCRS). This report discusses how SIF files were generated for the CCRS risk assessment and presents metrics on the accuracy of the impinging fluence that is generated when using the files. The results show a maximum difference per element in a finite element model within 1 percent for seven of the eight examined cases, with the outlier being less than 2 percent. The results also show that the maximum standard deviation of differences across all elements is within 0.275 percent. The maximum difference in the total fluence across all the elements is 0.0463 percent. The maximum difference in the mean fluence across all the elements is 0.0407 percent.

Meteoroids↗
NASA NTRS1998

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗
NASA NTRS1996

Cascade Error Projection: A Learning Algorithm for Hardware Implementation

In this paper, we workout a detailed mathematical analysis for a new learning algorithm termed Cascade Error Projection (CEP) and a general learning frame work. This frame work can be used to obtain the cascade correlation learning algorithm by choosing a particular set of parameters. Furthermore, CEP learning algorithm is operated only on one layer, whereas the other set of weights can be calculated deterministically. In association with the dynamical stepsize change concept to convert the weight update from infinite space into a finite space, the relation between the current stepsize and the previous energy level is also given and the estimation procedure for optimal stepsize is used for validation of our proposed technique. The weight values of zero are used for starting the learning for every layer, and a single hidden unit is applied instead of using a pool of candidate hidden units similar to cascade correlation scheme. Therefore, simplicity in hardware implementation is also obtained. Furthermore, this analysis allows us to select from other methods (such as the conjugate gradient descent or the Newton's second order) one of which will be a good candidate for the learning technique. The choice of learning technique depends on the constraints of the problem (e.g., speed, performance, and hardware implementation); one technique may be more suitable than others. Moreover, for a discrete weight space, the theoretical analysis presents the capability of learning with limited weight quantization. Finally, 5- to 8-bit parity and chaotic time series prediction problems are investigated; the simulation results demonstrate that 4-bit or more weight quantization is sufficient for learning neural network using CEP. In addition, it is demonstrated that this technique is able to compensate for less bit weight resolution by incorporating additional hidden units. However, generation result may suffer somewhat with lower bit weight quantization.

Duong, Tuan A.↗
NASA NTRS2001

Rapid Prototyping Integrated With Nondestructive Evaluation and Finite Element Analysis

Most reverse engineering approaches involve imaging or digitizing an object then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. Rapid prototyping (RP) refers to the practical ability to build high-quality physical prototypes directly from computer aided design (CAD) files. Using rapid prototyping, full-scale models or patterns can be built using a variety of materials in a fraction of the time required by more traditional prototyping techniques (refs. 1 and 2). Many software packages have been developed and are being designed to tackle the reverse engineering and rapid prototyping issues just mentioned. For example, image processing and three-dimensional reconstruction visualization software such as Velocity2 (ref. 3) are being used to carry out the construction process of three-dimensional volume models and the subsequent generation of a stereolithography file that is suitable for CAD applications. Producing three-dimensional models of objects from computed tomography (CT) scans is becoming a valuable nondestructive evaluation methodology (ref. 4). Real components can be rendered and subjected to temperature and stress tests using structural engineering software codes. For this to be achieved, accurate high-resolution images have to be obtained via CT scans and then processed, converted into a traditional file format, and translated into finite element models. Prototyping a three-dimensional volume of a composite structure by reading in a series of two-dimensional images generated via CT and by using and integrating commercial software (e.g. Velocity2, MSC/PATRAN (ref. 5), and Hypermesh (ref. 6)) is being applied successfully at the NASA Glenn Research Center. The building process from structural modeling to the analysis level is outlined in reference 7. Subsequently, a stress analysis of a composite cooling panel under combined thermomechanical loading conditions was performed to validate this process.

Abdul-Aziz, Ali↗
NASA NTRS1994

An Adaptively-Refined, Cartesian, Cell-Based Scheme for the Euler and Navier-Stokes Equations

A Cartesian, cell-based scheme for solving the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal 'cut' cells are created. The geometry of the cut cells is computed using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded, with a limited linear reconstruction of the primitive variables used to provide input states to an approximate Riemann solver for computing the fluxes between neighboring cells. A multi-stage time-stepping scheme is used to reach a steady-state solution. Validation of the Euler solver with benchmark numerical and exact solutions is presented. An assessment of the accuracy of the approach is made by uniform and adaptive grid refinements for a steady, transonic, exact solution to the Euler equations. The error of the approach is directly compared to a structured solver formulation. A non smooth flow is also assessed for grid convergence, comparing uniform and adaptively refined results. Several formulations of the viscous terms are assessed analytically, both for accuracy and positivity. The two best formulations are used to compute adaptively refined solutions of the Navier-Stokes equations. These solutions are compared to each other, to experimental results and/or theory for a series of low and moderate Reynolds numbers flow fields. The most suitable viscous discretization is demonstrated for geometrically-complicated internal flows. For flows at high Reynolds numbers, both an altered grid-generation procedure and a different formulation of the viscous terms are shown to be necessary. A hybrid Cartesian/body-fitted grid generation approach is demonstrated. In addition, a grid-generation procedure based on body-aligned cell cutting coupled with a viscous stensil-construction procedure based on quadratic programming is presented.

Coirier, William John↗
NASA NTRS1950

Use of source distributions for evaluating theoretical aerodynamics of thin finite wings at supersonic speeds

A series of publications on the source-distribution methods for evaluating the aerodynamics of thin wings at supersonic speeds is summarized, extended, and unified. Included in the first part are the deviations of: (a) the linearized partial-differential equation for unsteady flow at a substantially constant Mach number. b) The source-distribution solution for the perturbation-velocity potential that satisfies the boundary conditions of tangential flow at the surface and in the plane of the wing; and (c) the integral equation for determining the strength and the location of sources to describe the interaction effects (as represented by upwash) of the bottom and top wing surfaces through the region between the finite wing boundary and the foremost Mach wave. The second part deals with steady-state thin-wing problems. The third part of the report approximates the integral equation for unsteady upwash and includes a solution of approximate equation. Expressions are then derived to evaluate the load distributions for time-dependent finite-wing motions.

Evvard, John C↗
NASA NTRS2004

Integrated Modeling Activities for the James Webb Space Telescope (JWST): Structural-Thermal-Optical Analysis

This is a continuation of a series of papers on modeling activities for JWST. The structural-thermal- optical, often referred to as "STOP", analysis process is used to predict the effect of thermal distortion on optical performance. The benchmark STOP analysis for JWST assesses the effect of an observatory slew on wavefront error. The paper begins an overview of multi-disciplinary engineering analysis, or integrated modeling, which is a critical element of the JWST mission. The STOP analysis process is then described. This process consists of the following steps: thermal analysis, structural analysis, and optical analysis. Temperatures predicted using geometric and thermal math models are mapped to the structural finite element model in order to predict thermally-induced deformations. Motions and deformations at optical surfaces are input to optical models and optical performance is predicted using either an optical ray trace or WFE estimation techniques based on prior ray traces or first order optics. Following the discussion of the analysis process, results based on models representing the design at the time of the System Requirements Review. In addition to baseline performance predictions, sensitivity studies are performed to assess modeling uncertainties. Of particular interest is the sensitivity of optical performance to uncertainties in temperature predictions and variations in metal properties. The paper concludes with a discussion of modeling uncertainty as it pertains to STOP analysis.

Johnston, John D.↗
NASA NTRS1994

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗
NASA NTRS1988

A methodology for the synthesis of robust feedback systems

A new methodology is developed for the synthesis of linear, time-variant (LTI) controllers for multivariable LTI systems. The resulting closed-loop system is nominally stable and exhibits a known level of performance. In addition, robustness of the feedback system is guaranteed, i.e., stability and performance are retained in the presence of multiple unstructured uncertainty blocks located at various points in the feedback loop. The design technique is referred to as the Causality Recovery Methodology (CRM). The CRM relies on the Youla parameterization of all stabilizing compensators to ensure nominal stability of the feedback system. A frequency-domain inequality in terms of the structured singular value mu defines the robustness specification. The optimal compensator, with respect to the mu condition, is shown to be noncausal in general. The aim of the CRM is to find a stable, causal transfer function matrix that approximates the robustness characteristics of the optimal solution. The CRM, via a series of infinite-dimensional convex programs, produces a closed-loop system whose performance robustness is at least as good as that of any initial design. The algorithm is approximated by a finite dimensional process for the purposes of implementation. Two numerical examples confirm the potential viability of the CRM concept; however, the robustness improvement comes at the expense of increased computational burden and compensator complexity.

Milich, David Albert↗
NASA NTRS1993

Three-dimensional numerical simulation of thermal tides in the Martian atmosphere

Analysis of the Mariner 9 IRIS and Viking Orbiter IRTM temperature data and of the time series of surface pressure from the Viking landers has indicated that atmospheric thermal tides are a significant component of the martian general circulation. Classical tidal theory, considering only the Sun synchronous components, has been shown to at least roughly account for the amplitudes of the observed diurnal and semidirunal surface pressure oscillations. In particular, Zurek and Leovy have demonstrated that classical tidal theory can reproduce the observed dependence of the strength of the semidiurnal tide on aerosol optical depth. The high-amplitude topograph of Mars, however, can cause longitudinal distortions in the thermotidal forcing, leading to additional tidal modes that are not Sun synchronous. In particular, observation and theory have suggested the likely presence of a diurnal, zonal wavenumber-1 K wave. Hamilton and Garcia noted that the wavenumber-1 K normal mode has a period of close to 24 hr and may be detectable in surface pressure observations. Zurek has discussed the possibility of this mode being resonantly excited for various atmospheric temperature studies. The semidiurnal wavenumber-2 K mode, with roughly a 12-hr period, could also be resonantly excited. A numerical model with finite amplitude topography and a self-consistent thermotidal forcing resulting from the daily heat flux at the surface and atmospheric absorption of solar radiation by aerosols is used here to explore the influence of topography on martian thermal tides and to examine the possibility of Kelvin wave resonance.

Wilson, R. John↗
NASA NTRS1987

Plasma depletions in the Jovian magnetosphere - Evidence of transport and solar wind interaction

A series of plasma voids ('dropouts') was observed by the Plasma Science (PLS) experiment in Jupiter's magnetosphere during the Voyager 2 encounter with that planet. A reexamination of Voyager 2 data has led to the conclusion that the dropout phenomenon cannot be a manifestation of a plasma wake produced by Ganymede. Rather, the appearance of the dropouts is attributed to changes in the upstream solar wind conditions and the global state of the magnetosphere; the proximity of Voyager 2 to Ganymede at the time is considered to be coincidental. It is suggested that these dropouts are evidence of a state of 'bubbling' of the magnetosphere that alternates with 'laminar' states in which, as in the case of the Voyager 1 encounter with Jupiter, voids are not present and that these states correspond to different processes by which plasma is transported out of the system. The nature of these states is related to changes in the magnitude of the upstream solar wind ram pressure. In the bubbling state, this pressure is higher than in the laminar state and drives an intermittent instability. The analysis presented is one of the first attempts to introduce, in space physics, recently acquired theoretical notions of the physics of the finite-beta plasmas of which the Jovian magnetospheric plasma is an important example.

Mcnutt, Ralph L., Jr.↗
NASA NTRS2020

Bi-Axial Load Testing of a Woven-Webbing Inflatable Space Habitat Restraint Layer Component

This report discusses testing of a bi-axial woven panel consisting of high-strength webbings that represent a portion of the structural restraint layer of an inflatable space structure. Inflatable softgoods vessels are being researched for human space missions as possible habitats, airlocks and tunnel elements. Understanding the complex behavior of the softgoods restraint layer and maturing finite element analysis capabilities to model these structures is critical to their successful implementation. The primary goal of this research is to study the load-up and load distribution in a weave of webbings, before and after the loss of tension in one of the webbings. In addition, a key objective is to evaluate the ability to convert strains measured via photogrammetry using digital image correlation to loads in the weave. The report gives an introduction and objectives for the test program and a description of the test fixture, setup and procedure. This is followed by a section focusing on the test data and associated discussion. To streamline presentation of the data, one focus case is detailed in the main text. (Comprehensive data set can be found in Appendix A for the three remaining test cases.) In addition to the load and time histories for one case, summary charts and tables that incorporate data from all four test cases are provided. The concluding remarks include major findings, lessons learned and recommendations for future work. Appendix B provides complementary full-field strain results when considering the existing woven webbing data as a fabric. These types of results could inform models intended for global representation. Appendix C contains a description from a series of exploratory instrumented hammer tap tests and results on the weave. These tests describe an area of research for the biaxial test fixture beyond the current testing studying the propagation and detection of vibrations in a tensioned weave that has applications in impact and damage detection.

Karen H Lyle↗
NASA NTRS1994

Particle-in-cell simulations of the critical ionization velocity effect in finite size clouds

The critical ionization velocity (CIV) mechanism in a finite size cloud is studied with a series of electrostatic particle-in-cell simulations. It is observed that an initial seed ionization, produced by non-CIV mechanisms, generates a cross-field ion beam which excites a modified beam-plasma instability (MBPI) with frequency in the range of the lower hybrid frequency. The excited waves accelerate electrons along the magnetic field up to the ion drift energy that exceeds the ionization energy of the neutral atoms. The heated electrons in turn enhance the ion beam by electron-neutral impact ionization, which establishes a positive feedback loop in maintaining the CIV process. It is also found that the efficiency of the CIV mechanism depends on the finite size of the gas cloud in the following ways: (1) Along the ambient magnetic field the finite size of the cloud, L (sub parallel), restricts the growth of the fastest growing mode, with a wavelength lambda (sub m parallel), of the MBPI. The parallel electron heating at wave saturation scales approximately as (L (sub parallel)/lambda (sub m parallel)) (exp 1/2); (2) Momentum coupling between the cloud and the ambient plasma via the Alfven waves occurs as a result of the finite size of the cloud in the direction perpendicular to both the ambient magnetic field and the neutral drift. This reduces exponentially with time the relative drift between the ambient plasma and the neutrals. The timescale is inversely proportional to the Alfven velocity. (3) The transvers e charge separation field across the cloud was found to result in the modulation of the beam velocity which reduces the parallel heating of electrons and increases the transverse acceleration of electrons. (4) Some energetic electrons are lost from the cloud along the magnetic field at a rate characterized by the acoustic velocity, instead of the electron thermal velocity. The loss of energetic electrons from the cloud seems to be larger in the direction of plasma drift relative to the neutrals, where the loss rate is characterized by the neutral drift velocity. It is also shown that a factor of 4 increase in the ambient plasma density, increases the CIV ionization yield by almost 2 orders of magnitude at the end of a typical run. It is concluded that a larger ambient plasma density can result in a larger CIV yield because of (1) larger seed ion production by non-CIV mechanisms, (2) smaller Alfven velocity and hence weak momentum coupling, and (3) smaller ratio of the ion beam density to the ambient ion density, and therefore a weaker modulation of the beam velocity. The simulation results are used to interpret various chemical release experiments in space.

Moghaddam-Taaheri, E.↗
NASA NTRS1986

Correlation of data on strain accumulation adjacent to the San Andreas Fault with available models

Theoretical and numerical studies of deformation on strike slip faults were performed and the results applied to geodetic observations performed in the vicinity of the San Andreas Fault in California. The initial efforts were devoted to an extensive series of finite element calculations of the deformation associated with cyclic displacements on a strike-slip fault. Measurements of strain accumulation adjacent to the San Andreas Fault indicate that the zone of strain accumulation extends only a few tens of kilometers away from the fault. There is a concern about the tendency to make geodetic observations along the line to the source. This technique has serious problems for strike slip faults since the vector velocity is also along the fault. Use of a series of stations lying perpendicular to the fault whose positions are measured relative to a reference station are suggested to correct the problem. The complexity of faulting adjacent to the San Andreas Fault indicated that the homogeneous elastic and viscoelastic approach to deformation had serious limitations. These limitation led to the proposal of an approach that assumes a fault is composed of a distribution of asperities and barriers on all scales. Thus, an earthquake on a fault is treated as a failure of a fractal tree. Work continued on the development of a fractal based model for deformation in the western United States. In order to better understand the distribution of seismicity on the San Andreas Fault system a fractal analog was developed. The fractal concept also provides a means of testing whether clustering in time or space is a scale-invariant process.

Turcotte, Donald L.↗
NASA NTRS2002

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗