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 37 records · Page 2

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗

An eigenvalue method for solving transient heat conduction problems

The eigenvalue method, which has been used by researchers in structure mechanics, is applied to problems in heat conduction. Its formulation is decribed in terms of an examination of transient heat conduction in a square slab. Taking advantage of the availability of the exact solution, we compare the accuracy and other numerical properties of the eigenvalue method with those of existing numerical schemes. The comparsion shows that, overall, the eigenvalue method appears to be fairly attractive. Furthermore, only a few dominant eigenvalues and their corresponding eigenvectors need to be computed and retained to yield reasonably high accuracy. Greater savings are attained in the computation time for a transient problem with long time duration and a large computational domain.

Shih, T. M.↗

New computer program solves wide variety of heat flow problems

Boeing Engineering Thermal Analyzer /BETA/ computer program uses numerical methods to provide accurate heat transfer solutions to a wide variety of heat flow problems. The program solves steady-state and transient problems in almost any situation that can be represented by a resistance-capacitance network.

Almond, J. C.↗

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

Transient thermal stress problem for a circumferentially cracked hollow cylinder

The transient thermal stress problem for a hollow elasticity cylinder containing an internal circumferential edge crack is considered. It is assumed that the problem is axisymmetric with regard to the crack geometry and the loading, and that the inertia effects are negligible. The problem is solved for a cylinder which is suddenly cooled from inside. First the transient temperature and stress distributions in an uncracked cylinder are calculated. By using the equal and opposite of this thermal stress as the crack surface traction in the isothermal cylinder the crack problem is then solved and the stress intensity factor is calculated. The numerical results are obtained as a function of the Fourier number tD/b(2) representing the time for various inner-to-outer radius ratios and relative crack depths, where D and b are respectively the coefficient of diffusivity and the outer radius of the cylinder.

Nied, H. F.↗

Engineering thermal analyzer /BETA 2/

Computer program uses numerical methods to provide accurate heat transfer solutions to a wide variety of heat flow problems. This highly versatile program will solve steady state and transient problems in almost any situation that can be presented by a resistance-capacitance network.

Scates, J. H.↗

A theoretical analysis of steady-state photocurrents in simple silicon diodes

A theoretical analysis solves for the steady-state photocurrents produced by a given photo-generation rate function with negligible recombination in simple silicon diodes, consisting of a uniformly doped quasi-neutral region (called 'substrate' below) adjacent to a p-n junction depletion region (DR). Special attention is given to conditions that produce 'funneling' (a term used by the single-eventeffects community) under steady-state conditions. Funneling occurs when carriers are generated so fast that the DR becomes flooded and partially or completely collapses. Some or nearly all of the applied voltage, plus built-in potential normally across the DR, is now across the substrate. This substrate voltage drop affects substrate currents. The steady-state problem can provide some qualitative insights into the more difficult transient problem. First, it was found that funneling can be induced from a distance, i.e., from carriers generated at locations outside of the DR. Secondly, it was found that the substrate can divide into two subregions, with one controlling substrate resistance and the other characterized by ambipolar diffusion. Finally, funneling was found to be more difficult to induce in the p(sup +)/n diode than in the n(sup +)/p diode. The carrier density exceeding the doping density in the substrate and at the DR boundary is not a sufficient condition to collapse a DR.

Edmonds, L.↗

Solutions to divergence form equations using the method of partial implicitization

Partial implicitization is well suited for use on vector-processing computers since the method is an explicit unconditionally stable numerical technique. So far equations cast in nondivergence form have been used, though equations cast in the divergence form would have inherently greater accuracy. Since the method of partial implicitization is applicable only to relaxing problems to their steady state and is not applicable to true transient problems, a simple mathematical manipulation can be carried out which will produce the desired partial implicitization solution to divergence form equations.

Graves, R. A., Jr.↗

Sensitivity Equation Derivation for Transient Heat Transfer Problems

The focus of the paper is on the derivation of sensitivity equations for transient heat transfer problems modeled by different discretization processes. Two examples will be used in this study to facilitate the discussion. The first example is a coupled, transient heat transfer problem that simulates the press molding process in fabrication of composite laminates. These state equations are discretized into standard h-version finite elements and solved by a multiple step, predictor-corrector scheme. The sensitivity analysis results based upon the direct and adjoint variable approaches will be presented. The second example is a nonlinear transient heat transfer problem solved by a p-version time-discontinuous Galerkin's Method. The resulting matrix equation of the state equation is simply in the form of Ax = b, representing a single step, time marching scheme. A direct differentiation approach will be used to compute the thermal sensitivities of a sample 2D problem.

Hou, Gene↗

Transient thermal stress problem for a circumferentially cracked hollow cylinder

The paper is concerned with the transient thermal stress problem for a long hollow circular cylinder containing an internal axisymmetric circumferential edge crack that is suddenly cooled from inside. It is assumed that the transient thermal stress problem is quasi-static, i.e., the inertial effects are negligible. Also, all thermoelastic coupling effects and the possible temperature dependence of the thermoelastic constants are neglected. The problem is considered in two parts. The first part is the evaluation of transient thermal stresses in an uncracked cylinder; the second part is the isothermal perturbation problem for the cracked cylinder in which the crack surface tractions, equal and opposite to the thermal stresses obtained from the first problem, are the only external loads. The superposition of the two solutions gives results for the cracked cylinder.

Nied, H. F.↗

NASTRAN solutions of problems described by simultaneous parabolic differential equations

NASTRAN solution techniques are shown for a numerical analysis of a class of coupled vector flow processes described by simultaneous parabolic differential equations. To define one physical problem type where equations of this form arise, the differential equations describing the coupled transfers of heat and mass in mechanical equilibrium with negligible mass average velocity are presented and discussed. Also shown are the equations describing seepage when both electrokinetic and hydrodynamic forces occur. Based on a variational statement of the general problem type, the concepts of scalar transfer elements and parallel element systems are introduced. It is shown that adoptation of these concepts allows the direct use of NASTRAN's existing Laplace type elements for uncoupled flow (the heat transfer elements) for treating multicomponent coupled transfer. Sample problems are included which demonstrate the application of these techniques for both steady-state and transient problems.

Mason, J. B.↗

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

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗

Application of the Hughes-LIU algorithm to the 2-dimensional heat equation

An implicit explicit algorithm for the solution of transient problems in structural dynamics is described. The method involved dividing the finite elements into implicit and explicit groups while automatically satisfying the conditions. This algorithm is applied to the solution of the linear, transient, two dimensional heat equation subject to an initial condition derived from the soluton of a steady state problem over an L-shaped region made up of a good conductor and an insulating material. Using the IIT/PRIME computer with virtual memory, a FORTRAN computer program code was developed to make accuracy, stability, and cost comparisons among the fully explicit Euler, the Hughes-Liu, and the fully implicit Crank-Nicholson algorithms. The Hughes-Liu claim that the explicit group governs the stability of the entire region while maintaining the unconditional stability of the implicit group is illustrated.

Malkus, D. S.↗

A p-version finite element method for steady incompressible fluid flow and convective heat transfer

A new p-version finite element formulation for steady, incompressible fluid flow and convective heat transfer problems is presented. The steady-state residual equations are obtained by considering a limiting case of the least-squares formulation for the transient problem. The method circumvents the Babuska-Brezzi condition, permitting the use of equal-order interpolation for velocity and pressure, without requiring the use of arbitrary parameters. Numerical results are presented to demonstrate the accuracy and generality of the method.

Winterscheidt, Daniel L.↗

Finite element methods for integrated aerodynamic heating analysis

This report gives a description of the work which has been undertaken during the second year of a three year research program. The objectives of the program are to produce finite element based procedures for the solution of the large scale practical problems which are of interest to the Aerothermal Loads Branch (ALB) at NASA Langley Research Establishment. The problems of interest range from Euler simulations of full three dimensional vehicle configurations to local analyses of three dimensional viscous laminar flow. Adaptive meshes produced for both steady state and transient problems are to be considered. An important feature of the work is the provision of specialized techniques which can be used at ALB for the development of an integrated fluid/thermal/structural modeling capability.

Morgan, K.↗

Modeling of transient two-component flow using a four-point implicit method

The four-point, centered implicit scheme that is extensively used in open channel flow simulation is shown to be applicable to rapid and slow pressure transient problems in conduits with nearly single phase and two-phase flows. It is only necessary to choose the proper weighting factor value, theta, of the Courant number. For rapid pressure transients such as waterhammer, the implicit method can yield reasonable results with limited numerical dispersion and attenuation if theta is only slightly greater than the critical value of 0.5. For slower pressure gradients in single and two-phase flows, reasonable numerical solutions may be achieved for Courant number values as high as 20.

Wiggert, D. C.↗

Probabilistic finite elements for fatigue and fracture analysis

An overview of the probabilistic finite element method (PFEM) developed by the authors and their colleagues in recent years is presented. The primary focus is placed on the development of PFEM for both structural mechanics problems and fracture mechanics problems. The perturbation techniques are used as major tools for the analytical derivation. The following topics are covered: (1) representation and discretization of random fields; (2) development of PFEM for the general linear transient problem and nonlinear elasticity using Hu-Washizu variational principle; (3) computational aspects; (4) discussions of the application of PFEM to the reliability analysis of both brittle fracture and fatigue; and (5) a stochastic computational tool based on stochastic boundary element (SBEM). Results are obtained for the reliability index and corresponding probability of failure for: (1) fatigue crack growth; (2) defect geometry; (3) fatigue parameters; and (4) applied loads. These results show that initial defect is a critical parameter.

Belytschko, Ted↗