Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient problem”

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 235 records · Page 13

A Taylor-Galerkin finite element algorithm for transient nonlinear thermal-structural analysis

A Taylor-Galerkin finite element method for solving large, nonlinear thermal-structural problems is presented. The algorithm is formulated for coupled transient and uncoupled quasistatic thermal-structural problems. Vectorizing strategies ensure computational efficiency. Two applications demonstrate the validity of the approach for analyzing transient and quasistatic thermal-structural problems.

Thornton, E. A.↗

Combining DPG in space with DPG time-marching scheme for the transient advection-reaction equation

In this article, we present a general methodology to combine the Discontinuous PetrovGalerkin (DPG) method in space and time in the context of methods of lines for transient advection-reaction problems. We first introduce a semidiscretization in space with a DPG method redefining the ideas of optimal testing and practicality of the method in this context. Then, we apply the recently developed DPG-based time-marching scheme, which is of exponential-type, to the resulting system of Ordinary Differential Equations (ODEs). We also discuss how to efficiently compute the action of the exponential of the matrix coming from the space semidiscretization without assembling the full matrix. Finally, we verify the proposed method for 1D+time advection-reaction problems showing optimal convergence rates for smooth solutions and more stable results for linear conservation laws comparing to the classical exponential integrators.

97 MATHEMATICS AND COMPUTING↗

The leaking mode problem in atmospheric acoustic-gravity wave propagation

The problem of predicting the transient acoustic pressure pulse at long horizontal distances from large explosions in the atmosphere is examined. Account is taken of poles off the real axis and of branch line integrals in the general integral governing the transient waveform. Perturbation techniques are described for the computation of the imaginary ordinate of the poles and numerical studies are described for a model atmosphere terminated by a halfspace with c = 478 m/sec above 125 km. For frequencies less than 0.0125 rad/sec, the GR sub 1 mode, for example, is found to have a frequency dependent amplitude decay of the order of 0.0001 nepers/km. Examples of numerically synthesized transient waveforms are exhibited with and without the inclusion of leaking modes. The inclusion of leaking modes results in waveforms with a more marked beginning rather than a low frequency oscillating precursor of gradually increasing amplitude. Also, the revised computations indicate that waveforms invariably begin with a pressure rise, a result supported by other theoretical considerations and by experimental data.

Kinney, W. A.↗

A simulation model of time-dependent plasma-spacecraft interactions

A plasma simulation code is presented that models the time-dependent plasma properties in the vicinity of a spherical, charged spacecraft. After showing agreement with analytic, steady-state theories and ATS-6 satellite data, the following three problems are treated: (1) transient pulses from photoemission at various emission temperatures and ambient plasma conditions, (2) spacecharge limited emission, and (3) simulated plasma oscillations in the long wavelength limit.

Rothwell, P. L.↗

Interpretation methodology and analysis of in-flight lightning data

A methodology is presented whereby electromagnetic measurements of inflight lightning stroke data can be understood and extended to other aircraft. Recent measurements made on the NASA F106B aircraft indicate that sophisticated numerical techniques and new developments in corona modeling are required to fully understand the data. Thus the problem is nontrivial and successful interpretation can lead to a significant understanding of the lightning/aircraft interaction event. This is of particular importance because of the problem of lightning induced transient upset of new technology low level microcircuitry which is being used in increasing quantities in modern and future avionics. Inflight lightning data is analyzed and lightning environments incident upon the F106B are determined.

Rudolph, T.↗

Thermal-Analysis Program

Temperature changes with time modeled with traditional electrical analog. General Thermal Analyzer program solves transient and steadystate thermal problems using desk-top computers. Program written in BASIC.

Levine, D. I.↗

Heat flux calibration facility capable of SSME conditions

There is a need to more thoroughly characterize the hostile space shuttle main engine (SSME) turbopump environment. It has been estimated that component surface heat flux in the hot-gas environment is about 10 MW/square meter, and this is about 50 times that encountered in aircraft engines. Also, material temperature transients can be as high as 1000 K in about 1 second. These transients can cause durability problems such as material cracking. Heat flux sensors placed in the turbopump components can partially characterize this environment by measuring surface heat flux. These heat flux data can be used to verify analytical-stress, boundary-layer, and heat-transfer design models. Preliminary plans were discussed at the first SSME durablity conference for designing and fabricating a new facility for the calibration and durability testing of prototype heat flux sensors for the SSME. This facility, which is necessary for assessment of new heat flux gauge concepts needed in the hostile SSME turbopump environment, is described.

Liebert, Curt H.↗

Development of advanced modal methods for calculating transient thermal and structural response

Higher-order modal methods for predicting thermal and structural response are evaluated. More accurate methods or ones which can significantly reduce the size of complex, transient thermal and structural problems are desirable for analysis and are required for synthesis of real structures subjected to thermal and mechanical loading. A unified method is presented for deriving successively higher-order modal solutions related to previously-developed, lower-order methods such as the mode displacement and mode-acceleration methods. A new method, called the force-derivative method, is used to obtain higher-order modal solutions for both uncoupled (proportionally-damped) structural problems as well as thermal problems and coupled (non-proportionally damped) structural problems. The new method is called the force-derivative method because, analogous to the mode-acceleration method, it produces a term that depends on the forcing function and additional terms that depend on the time derivatives of the forcing function.

Camarda, Charles J.↗

A novel pulse technique for thermal diffusivity determination of high-temperature levitated materials - Inverse analysis and experimental set-up

In order to determine the thermal diffusivity of materials at high temperatures, a two-step extended flash technique, which is applicable to levitated spherical samples, is proposed. The containerless flash method is modeled as an axisymmetric transient conduction heat transfer problem within the sphere subjected to a radiative boundary condition on its surface. The problem of nonlinearity which arises from the radiative heat transfer boundary condition is solved by replacing it with the measured time-dependent surface temperature data, thus giving rise to an 'Inverse Analysis'. Upon obtaining the analytic solution for the temperature field, the determination of the thermal diffusivity turns into a minimization problem, whereby the difference between the actual surface boundary condition and the analytic solution is minimized. In performing the proposed experiments, there is a need to undertake a cool-down experiment which should be governed by the lumped parameter analysis.

Shen, F.↗

Observations of Methane and Ethylene Diffusion Flames Stabilized Around a Blowing Porous Sphere Under Microgravity Conditions

This paper presents the experimental and theoretical results for expanding methane and ethylene diffusion flames in microgravity. A small porous sphere made from a low-density and low-heat-capacity insulating material was used to uniformly supply fuel at a constant rate to the expanding diffusion flame. A theoretical model which includes soot and gas radiation is formulated but only the problem pertaining to the transient expansion of the flame is solved by assuming constant pressure infinitely fast one-step ideal gas reaction and unity Lewis number. This is a first step toward quantifying the effect of soot and gas radiation on these flames. The theoretically calculated expansion rate is in good agreement with the experimental results. Both experimental and theoretical results show that as the flame radius increases, the flame expansion process becomes diffusion controlled and the flame radius grows as gamma t. Theoretical calculations also show that for a constant fuel mass injection rate a quasi-steady state is developed in the region surrounded by the flame and the mass flow rate at any location inside this region equals the mass injection rate.

Atreya, Arvind↗

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗

Theoretical Studies of Accreting Neutron Stars

Among the newly discovered classes of X-ray sources which have attracted wide attention are close binary systems in which mass is transferred via Roche lobe overflow from a low mass donor star to its neutron star companion. Many of these sources exhibit intense bursts of X-ray radiation as well as periodic and quasi-periodic phenomena. Intensive analysis of these sources as a class has provided insight into the accretion process in binary star systems and into the magnetic field, rotational, and nuclear evolution of the underlying neutron star. In this proposal we have focused on theoretical studies of the hydrodynamical and nuclear processes that take place on the surface of accreting neutron stars in these systems. The investigation of these processes is critical for providing an understanding of a number of outstanding problems related to their transient behavior and evolution.

Taam, Ronald E.↗

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

Transmission line inversion and synthesis from the point of view of transient response

The transmission line inversion problem is applied to the speech communication process. The problem involves acquiring information necessary for synthesis of speech and for making models of the speech production mechanism. Formulation for the problem is derived from an analysis of the human vocal tract functions. It is known that speech sounds are produced by acoustic excitation of the tract by means of quasi-periodic pulses from the vocal cords during vowel sounds and turbulent air flow at various points along the tract during fricative sounds. By assuming plane wave propagation, differential equations are obtained relating the cross sectional area of the tract, the pressure, and the volume velocity. By replacing the pressure factor with a voltage factor and the volume velocity by current, the problem becomes a transmission line problem.

Sondhi, M. M.↗

Optimum data utilization for parameter identification with application to lifting rotors

The work is concerned with determining the minimum quantity of data needed for achieving best possible accuracy of identified parameters in transient testing of a lifting rotor blade. Specifically, the problem is that of determining from blade flapping transients caused by blade pitch inputs the equivalent Lock number and the equivalent collective pitch setting. For a given time-dependent blade pitch input function, the running time for the test less than which insufficient parameter accuracy is obtained and more than which little accuracy improvement is achieved was calculated with the aid of the Cramer-Rao lower bound for the parameter covariance matrix.

Banerjee, D.↗

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact

A numerical transient solution of the hydrodynamically lubricated point contact problem is obtained using the ball-on-plane model. Results, which include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency, are compared with the analytic solution of the transient step bearing problem with the same dynamic loading function.

Lim, Sang G.↗

Deriving Transmissibility Functions from Finite Elements for Specifications

This work explores deriving transmissibility functions for a missile from a measured location at the base of the fairing to a desired location within the payload. A pressure on the outside of the fairing and the rocket motor’s excitation creates an acceleration at a measured location and a desired location. Typically, the desired location is not measured. In fact, it is typical that the payload may change, but measured acceleration at the base of the fairing is generally similar to previous test flights. Given this knowledge, it is desired to use a finite-element model to create a transmissibility function which relates acceleration from the previous test flight’s measured location at the base of the fairing to acceleration at a location in the new payload. Four methods are explored for deriving this transmissibility, with the goal of finding an appropriate transmissibility when both the pressure and rocket motor excitation are equally present. In conclusion, these methods are assessed using transient results from a simple example problem, and it is found that one of the methods gives good agreement with the transient results for the full range of loads considered.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Finite element analysis of large transient elastic-plastic deformations of simple structures, with application to the engine rotor fragment containment/deflection problem

Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.

Wu, R. W.↗