Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical solution of the incompressible Navier-Stokes equations

The current work is initiated in an effort to obtain an efficient, accurate, and robust algorithm for the numerical solution of the incompressible Navier-Stokes equations in two- and three-dimensional generalized curvilinear coordinates for both steady-state and time-dependent flow problems. This is accomplished with the use of the method of artificial compressibility and a high-order flux-difference splitting technique for the differencing of the convective terms. Time accuracy is obtained in the numerical solutions by subiterating the equations in psuedo-time for each physical time step. The system of equations is solved with a line-relaxation scheme which allows the use of very large pseudo-time steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. Numerous laminar test flow problems are computed and presented with a comparison against analytically known solutions or experimental results. These include the flow in a driven cavity, the flow over a backward-facing step, the steady and unsteady flow over a circular cylinder, flow over an oscillating plate, flow through a one-dimensional inviscid channel with oscillating back pressure, the steady-state flow through a square duct with a 90 degree bend, and the flow through an artificial heart configuration with moving boundaries. An adequate comparison with the analytical or experimental results is obtained in all cases. Numerical comparisons of the upwind differencing with central differencing plus artificial dissipation indicates that the upwind differencing provides a much more robust algorithm, which requires significantly less computing time. The time-dependent problems require on the order of 10 to 20 subiterations, indicating that the elliptical nature of the problem does require a substantial amount of computing effort.

Rogers, Stuart E.↗
NASA NTRS1986

Viscous vortex flows

Several computational studies are currently being pursued that focus on various aspects of representing the entire lifetime of the viscous trailing vortex wakes generated by an aircraft. The formulation and subsequent near-wing development of the leading-edge vortices formed by a delta wing are being calculated at modest Reynolds numbers using a three-dimensional, time-dependent Navier-Stokes code. Another computational code was developed to focus on the roll-up, trajectory, and mutual interaction of trailing vortices further downstream from the wing using a two-dimensional, time-dependent, Navier-Stokes algorithm. To investigate the effect of a cross-wind ground shear flow on the drift and decay of the far-field trailing vortices, a code was developed that employs Euler equations along with matched asymptotic solutions for the decaying vortex filaments. And finally, to simulate the conditions far down stream after the onset of the Crow instability in the vortex wake, a full three-dimensional, time-dependent Navier-Stokes code was developed to study the behavior of interacting vortex rings.

Weston, R. P.↗
NASA NTRS1993

Alternative descriptions of wave and particle aspects of the harmonic oscillator

The dynamical properties of the wave and particle aspects of the harmonic oscillator can be studied with the help of the time-dependent Schroedinger equation (SE). Especially the time-dependence of maximum and width of Gaussian wave packet solutions allow to show the evolution and connections of those two complementary aspects. The investigation of the relations between the equations describing wave and particle aspects leads to an alternative description of the considered systems. This can be achieved by means of a Newtonian equation for a complex variable in connection with a conservation law for a nonclassical angular momentum-type quantity. With the help of this complex variable, it is also possible to develop a Hamiltonian formalism for the wave aspect contained in the SE, which allows to describe the dynamics of the position and momentum uncertainties. In this case the Hamiltonian function is equivalent to the difference between the mean value of the Hamiltonian operator and the classical Hamiltonian function.

Schuch, Dieter↗
NASA NTRS2000

Visualization of In-Flight Flow Phenomena Using Infrared Thermography

Infrared thermography was used to obtain data on the state of the boundary layer of a natural laminar flow airfoil in supersonic flight. In addition to the laminar-to-turbulent transition boundary, the infrared camera was able to detect shock waves and present a time dependent view of the flow field. A time dependent heat transfer code was developed to predict temperature distributions on the test subject and any necessary surface treatment. A commercially available infrared camera was adapted for airborne use in this application. Readily available infrared technology has the capability to provide detailed visualization of various flow phenomena in subsonic to hypersonic flight regimes.

Banks, D. W.↗
NASA NTRS1993

Theoretical and Experimental Investigations on Droplet Evaporation and Droplet Ignition at High Pressures

The evaporation of fuel droplets under high ambient pressure and temperature in normal gravity and microgravity has been investigated experimentally. For subcritical ambient conditions, droplet evaporation after a heat-up period follows the d(exp 2)-law. For all data the evaporation constant increases as the ambient temperature increases. At identical ambient conditions the evaporation constant under microgravity is smaller compared to normal gravity. This effect can first be observed at 1 bar and increases with ambient pressure. Preliminary experiments on ignition delay for self-igniting fuel droplets have been performed. Above a 1 s delay time, at identical ambient conditions, significant differences in the results of the normal and microgravity data are observed. Self-ignition occurs within different temperature ranges due to the influence of gravity. The time dependent behavior of the droplet is examined theoretically. In the calculations two different approaches for the gas phase are applied. In the first approach the conditions at the interface are given using a quasi steady theory approximation. The second approach uses a set of time dependent governing equations for the gas phase which are then evaluated. In comparison, the second model shows a better agreement with the drop tower experiments. In both cases a time dependent gasification rate is observed.

Ristau, R.↗
NASA NTRS1989

Comparison of parameterized nitric acid rainout rates using a coupled stochastic-photochemical tropospheric model

A major tropospheric loss of soluble species such as nitric acid results from scavenging by water droplets. Several theoretical formulations have been advanced which relate an effective time-independent loss rate for soluble species to statistical properties of precipitation such as the wet fraction and length of a precipitation cycle. In this paper, various 'effective' loss rates that have been proposed are compared with the results of detailed time-dependent model calculations carried out over a seasonal time scale. The model is a stochastic precipitation model coupled to a tropospheric photochemical model. The results of numerous time-dependent seasonal model runs are used to derive numerical values for the nitric acid residence time for several assumed sets of preciptation statistics. These values are then compared with the results obtained by utilizing theoretical 'effective' loss rates in time-independent models.

Stewart, Richard W.↗
NASA NTRS1996

Stress and Reliability Analysis of a Metal-Ceramic Dental Crown

Interaction of mechanical and thermal stresses with the flaws and microcracks within the ceramic region of metal-ceramic dental crowns can result in catastrophic or delayed failure of these restorations. The objective of this study was to determine the combined influence of induced functional stresses and pre-existing flaws and microcracks on the time-dependent probability of failure of a metal-ceramic molar crown. A three-dimensional finite element model of a porcelain fused-to-metal (PFM) molar crown was developed using the ANSYS finite element program. The crown consisted of a body porcelain, opaque porcelain, and a metal substrate. The model had a 300 Newton load applied perpendicular to one cusp, a load of 30ON applied at 30 degrees from the perpendicular load case, directed toward the center, and a 600 Newton vertical load. Ceramic specimens were subjected to a biaxial flexure test and the load-to-failure of each specimen was measured. The results of the finite element stress analysis and the flexure tests were incorporated in the NASA developed CARES/LIFE program to determine the Weibull and fatigue parameters and time-dependent fracture reliability of the PFM crown. CARES/LIFE calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/Or proof test loading. This program is an extension of the CARES (Ceramics Analysis and Reliability Evaluation of Structures) computer program.

Anusavice, Kenneth J↗
NASA NTRS1990

The behavior of elastic anisotropic laminated composite flat structures subjected to deterministic and random loadings

Within this research project, the following topics were studied: (1) foundation of the refined theory of flat cross-ply laminated composite flat and curved panels as well as their static and dynamic response analysis; (2) foundation of a geometrically-nonlinear shear-deformable theory of composite laminated flat panels including the effect of initial geometric imperfections and its application in the postbuckling analysis; (3) the study of the dynamic response of shear deformable elastic laminated composite panels to deterministic time-dependent external excitations as the sonic boom and explosive blast type-loadings; (4) the study of the dynamic response of shear deformable elastic laminated composite panels to random excitation as e.g. the one produced by a jet noise or by any time-dependent external excitation whose characteristics are expressed in a statistical sense; and (5) the dynamic stability of fiber-reinforced composite flat panels whose materials (due to e.g. an ambient high temperature field) exhibit a time-dependent physical behavior.

Librescu, Liviu↗
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗
NASA NTRS1993

Numerical Simulation of a High Mach Number Jet Flow

The recent efforts to develop accurate numerical schemes for transition and turbulent flows are motivated, among other factors, by the need for accurate prediction of flow noise. The success of developing high speed civil transport plane (HSCT) is contingent upon our understanding and suppression of the jet exhaust noise. The radiated sound can be directly obtained by solving the full (time-dependent) compressible Navier-Stokes equations. However, this requires computational storage that is beyond currently available machines. This difficulty can be overcome by limiting the solution domain to the near field where the jet is nonlinear and then use acoustic analogy (e.g., Lighthill) to relate the far-field noise to the near-field sources. The later requires obtaining the time-dependent flow field. The other difficulty in aeroacoustics computations is that at high Reynolds numbers the turbulent flow has a large range of scales. Direct numerical simulations (DNS) cannot obtain all the scales of motion at high Reynolds number of technological interest. However, it is believed that the large scale structure is more efficient than the small-scale structure in radiating noise. Thus, one can model the small scales and calculate the acoustically active scales. The large scale structure in the noise-producing initial region of the jet can be viewed as a wavelike nature, the net radiated sound is the net cancellation after integration over space. As such, aeroacoustics computations are highly sensitive to errors in computing the sound sources. It is therefore essential to use a high-order numerical scheme to predict the flow field. The present paper presents the first step in a ongoing effort to predict jet noise. The emphasis here is in accurate prediction of the unsteady flow field. We solve the full time-dependent Navier-Stokes equations by a high order finite difference method. Time accurate spatial simulations of both plane and axisymmetric jet are presented. Jet Mach numbers of 1.5 and 2.1 are considered. Reynolds number in the simulations was about a million. Our numerical model is based on the 2-4 scheme by Gottlieb & Turkel. Bayliss et al. applied the 2-4 scheme in boundary layer computations. This scheme was also used by Ragab and Sheen to study the nonlinear development of supersonic instability waves in a mixing layer. In this study, we present two dimensional direct simulation results for both plane and axisymmetric jets. These results are compared with linear theory predictions. These computations were made for near nozzle exit region and velocity in spanwise/azimuthal direction was assumed to be zero.

Hayder, M. Ehtesham↗
NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗
NASA NTRS1995

Neutral solar wind evolution during solar cycle

The time dependent model of the expected fluxes of the neutral H and He components of the solar wind in the inner heliosphere is discussed. The model takes into account typical temporal evolution of the distribution of neutral interstellar gases (hydrogen and helium) in the interplanetary space due to solar cycle effects and the long term variability of the solar wind. The contribution of different charge exchange processes to the production of particular NSW element is presented. The distribution of the NSW flux is analysed with respect to the heliocentric distance and azimuthal angle from the Interstellar Wind apex. It demonstrates significant, time-dependent upwind/downwind H and He flux asymmentries. It is shown that the most pronounced modulation of the NSW flux is expected around the solar maximum epoch, when a strong decrease of the energetic H flux by two three orders of magnitude at 1 AU is predicted. The computations show that in the inner solar system (approx. 1 AU) energetic helium atoms production in the downwind region usually dominates the production of the hydrogen component This leads to the conclusion that the NSW composition at the Earth orbit strongly depends on time and the position of the observation point in reference to the apex direction.

Bzowski, M.↗
NASA NTRS1985

Drag characteristics of unsteady, perturbed boundary layer flows

A series of time-dependent numerical computations have been performed for the flow of an unsteady boundary layer over an embedded cavity, which is aligned normal to the flow direction. The unsteady flow is the result of the imposition of a Stuart vortex, at the inflow boundary of the computational domain, onto an otherwise unperturbed two-dimensional laminar boundary layer flow. This produces a strongly vortical, but spatially monochromatic, motion at inflow which is allowed to evolve downstream through the time-dependent numerical solution of appropriate conservation equations. As the imposed vortex evolves over the solid boundary, a weaker induced vortical motion is produced and together these perturbing motions interact with the fluid motion in the embedded cavity. An analysis of the pertinent dynamic variables as well as the time-dependent drag characteristics for this type of non-planar geometry is performed; with particular emphasis on the relative contribution of both pressure drag and frictional drag at various times in the evolution process.

Gatski, T. B.↗
NASA NTRS1987

Linear baroclinic instability in the presence of large scale topography

The effect of a planetary-scale, wavenumber 2 topography on baroclinically active disturbances is investigated for a channel domain in a two-layer, quasi-geostrophic context. When the lower-layer zonal velocity is nonzero, the topography influences the disturbances by forcing a stationary wave, and the topography and the forced wave influence the growth rates and the spatial structures of the time-dependent solutions. The case of zero zonal velocity in the lower layer was also investigated, for which no forced wave exists. Asymptotic forms of the equations, valid when the topographic effect (governed by the ratio of the nondimensional topographic height to the rotational Froude number) is small, are used to obtain both the stationary and time-dependent solutions. The latter are also obtained using a numerical approach, in which is determined the eigenvalues and eigenfunctions of a matrix representing the dynamical equations. Agreement is good between the two approaches. Recent laboratory experiments with a baroclinic annulus in which there is a false bottom with wavenumber 2 topography, are used to select governing parameters. The simultaneous presence of a stationary forced wave of wavenumber 2 and a time-dependent baroclinic wave of wavenumber 4, which has wavenumber 2 and 6 sidebands due to the topography, yields a flow field that exhibits some principal features of the laboratory experiments. The position of the forced wave and the location of an excursion in latitude of the storm track show qualitative resemblance to those features observed in the atmosphere.

Reynolds, Nathaniel Dunton↗
NASA NTRS1979

An energy-balance model with multiply-periodic and quasi-chaotic free oscillations

A slightly modified version of the one-dimensional time-dependent energy-balance climate model of Ghil and Bhattacharya (1978) is presented. The albedo-temperature parameterization has been reformulated and the smoothing of the temperature distribution in the tropics has been eliminated. The model albedo depends on time-lagged temperature in order to account for finite growth and decay time of continental ice sheets. Two distinct regimes of oscillatory behavior which depend on the value of the albedo-temperature time lag are considered.

Bhattacharya, K.↗
NASA NTRS1997

Science Support for Space-Based Droplet Combustion: Drop Tower Experiments and Detailed Numerical Modeling

This program supports the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies. Experimental emphasis is on the study of simple alcohols (methanol, ethanol) and alkanes (n-heptane, n-decane) as fuels with time dependent measurements of drop size, flame-stand-off, liquid-phase composition, and finally, extinction. Experiments have included bench-scale studies at Princeton, studies in the 2.2 and 5.18 drop towers at NASA-LeRC, and both the Fiber Supported Droplet Combustion (FSDC-1, FSDC-2) and the free Droplet Combustion Experiment (DCE) studies aboard the shuttle. Test matrix and data interpretation are performed through spherically-symmetric, time-dependent numerical computations which embody detailed sub-models for physical and chemical processes. The computed burning rate, flame stand-off, and extinction diameter are compared with the respective measurements for each individual experiment. In particular, the data from FSDC-1 and subsequent space-based experiments provide the opportunity to compare all three types of data simultaneously with the computed parameters. Recent numerical efforts are extending the computational tools to consider time dependent, axisymmetric 2-dimensional reactive flow situations.

Marchese, Anthony J.↗
NASA NTRS1995

The simulation of complete 11 and 12 year modulation cycles for cosmic rays in the heliosphere using a drift model with global merged interaction regions

Two-dimensional, time-dependent drift models have done exceptionally well in explaining major modulation features, especially during the A less than 0 magnetic polarity cycle of the heliospheric magnetic field when positively charge particles are drifting in along the heliospheric neutral sheet (HNS). These models were found to do well when the heliospheric 'tilt angles' alpha less than approximately 30 deg (le Roux & Potgeiter). However, they seem to do less well when alpha greater than approximately 30 deg during A less than 0 cycles seem to fail when this happens in A greater than 0 cycles. Progress was made in understanding these phases of the modulation cycle when merged interaction regions (MIRs) were incorporated in time-dependent drift models (Potgieter et al.). It was also explicitly shown that in obtaining large step decreases in cosmic rays, the MIRs had to be global, i.e., having a latitudinal extent of more than approximately 60 deg. Other classes of MIRs, such as local MIRs and co-rotating MIRs were found to be of secondary importance for establishing long-term modulation. In a previous paper we studied the effects of two consecutive, identical global MIRs, together with a changing wavy HNS, on long-term modulation (Potegieter & le Roux). This approach gave a very natural and convincing explanation for the observed step decreases in cosmic-ray modulation. Emphasis was placed in the declining and recovery phases of the 11 yr modulation cycle. In this paper, four consecutive, nonidentical global MIRs, in combination with a varying wavy HNS, were included in our time-dependent drift model in order to do simulations closer to what was observed between 1977 and 1987. By doing this we could model, for the first time, complete 11 and 22 yr cycles in the heliospheric modulation of galactic cosmic-rays, including the solar magnetic polarity reversals.

Le Roux, J. A.↗
NASA NTRS2001

Measurement of Microscopic Growth Rates in Float-Zone Silicon Crystals

Time dependent convective flows during crystal growth of doped semiconductors lead to fluctuations of the composition, so called dopant striations. In general, it is difficult to decide which is the main mechanism for the generation of these striations, it might be either the fluctuation of the concentration field in the melt and the extent of the solute boundary layer ahead of the solid-liquid interface or a variation of the growth velocity. Direct access to the concentration field is rather complicated to achieve, especially considering the high process temperature and the chemical activity of liquid silicon. The contribution of growth rate fluctuations to the formation of compositional fluctuations can be determined by measuring microscopic growth rates. The classical method of current pulses requires electrical feed-throughs and good electrical contacts, both are critical issues for the growth of high purity silicon crystals. Using a radiation based heating system, the heating power can be modulated very fast and effectively. We added to the normal heater power a sinusoidal off-set in the frequency range of 1 to 10 Hz, generating a narrow spaced weak rippling in the grown crystals which are superposed to the dopant striations caused by natural and by thermocapillary convection. The pulling speed was varied between 1 and 4mm/min. The microscope images of etched crystals slices have been analyzed by peak-search algorithms (measuring the spacing between each artificially induced marker) and by FFT. Performing growth experiments under a time-dependent flow regime, fluctuations of the microscopic growth velocity of Delta(v)/v(sub average) up to 50% have been measured. Damping the time-dependent convection by the use of an axial, static magnetic field of 500mT, the microscopic growth rate became constant within the resolution limit of this method. The results will be discussed using analytical methods for the calculation of microscopic growth velocity and by comparing them with measurements of temperature fluctuations in the melt during growth experiments itself.

Dold, P.↗