Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 1,135 records · Page 63

On turbulent spots in plane Poiseuille flow

Turbulence characteristics inside a turbulent spot in plume Poiseuille flow are investigated by analyzing a data base obtained from a direct numerical simulation. The spot is found to consist of two distinct regions - a turbulent area and a wave area. The flow inside the turbulent area has a strong resemblance to that found in the fully developed turbulent channel. Suitably defined mean and rms fluctuations as well as the internal shear layer structures are found to be similar to the turbulent counterpart. In the wave area, the inflexional mean spanwise profiles cause a rapid growth of oblique waves, which break down to turbulence. The breakdown process of the oblique waves is reminiscent of the secondary instability observed during transition to turbulence in channel and boundary layer flows. Other detailed characteristics associated with the Poiseuille spot are presented and are compared with experimental results.

Henningson, Dan S.↗

Energy buildup in sheared force-free magnetic fields

Photospheric displacement of the footpoints of solar magnetic field lines results in shearing and twisting of the field, and consequently in the buildup of electric currents and magnetic free energy in the corona. The sudden release of this free energy may be the origin of eruptive events like coronal mass ejections, prominence eruptions, and flares. An important question is whether such an energy release may be accompanied by the opening of magnetic field lines that were previously closed, for such open field lines can provide a route for matter frozen into the field to escape the sun altogether. This paper presents the results of numerical calculations showing that opening of the magnetic field is permitted energetically, in that it is possible to build up more free energy in a sheared, closed, force-free magnetic field than is in a related magnetic configuration having both closed and open field lines. Whether or not the closed force-free field attains enough energy to become partially open depends on the form of the shear profile; the results presented compare the energy buildup for different shear profiles. Implications for solar activity are discussed briefly.

Wolfson, Richard↗

Metallographic cooling rates of L-group ordinary chondrites

Shock metamorphism appears to be a ubiquitous feature in L-group ordinary chondrites. Brecciation and heterogeneous melting obscure much of the early history of this meteorite group and have caused confusion as to whether L chondrites have undergone thermal metamorphism within onion-shell or rubble-pile parent bodies. Employing the most recent shock criteria, we have examined 55 Antarctic and 24 non-Antarctic L chondrites in order to identify those which have been least affected by post-accretional shock. Six low-shock samples (those with shock grade less than S4) of petrographic types L3-L5 were selected from both populations and metallographic cooling rates were obtained following the technique of Willis and Goldstein. All non-Antarctic L6 chondrites inspected were too heavily shocked to be included in this group. However, 4 shocked L6 chondrites were analyzed in order to determine what effects shock may impose on metallographic cooling rates. Metallographic cooling rates were derived by analyzing the cores of taenite grains and then measuring the distance to the nearest grain edge. Taenites were identified using backscatter imaging on a Cameca SX-50 electron microprobe. Using backscatter we were able to locate homogeneous, rust-free, nearly spherical grains. M-shaped profiles taken from grain traverses were also used to help locate the central portions of selected grains. All points which contained phosphorus above detection limits were discarded. Plots of cooling-rate data are summarized and data from the high-shock samples are presented. The lack of coherency of cooling rates for individual samples is indicative of heterogeneous cooling following shock. The data confirms the statement expressed by numerous workers that extreme care must be taken when selecting samples of L chondrites for cooling-rate studies. Data for the 6 non-Antarctic low-shock samples are also presented. The samples display a general trend in cooling rates. The lowest metamorphic grade yielded the slowest cooling rates and an increase in grade follows an increase in cooling rate. This is the opposite relationship to that predicted by the onion-shell model.

Bennett, Marvin E.↗

Hydrodynamic Analyses and Evaluation of New Fluid Film Bearing Concepts

Over the past several years, numerical and experimental investigations have been performed on a waved journal bearing. The research work was undertaken by Dr. Florin Dimofte, a Senior Research Associate in the Mechanical Engineering Department at the University of Toledo. Dr. Theo Keith, Distinguished University Professor in the Mechanical Engineering Department was the Technical Coordinator of the project. The wave journal bearing is a bearing with a slight but precise variation in its circular profile such that a waved profile is circumscribed on the inner bearing diameter. The profile has a wave amplitude that is equal to a fraction of the bearing clearance. Prior to this period of research on the wave bearing, computer codes were written and an experimental facility was established. During this period of research considerable effort was directed towards the study of the bearing's stability. The previously developed computer codes and the experimental facility were of critical importance in performing this stability research. A collection of papers and reports were written to describe the results of this work. The attached captures that effort and represents the research output during the grant period.

Keith, Theo G., Jr.↗

Terrestrial Planet Formation Constrained By Mars and the Structure of the Asteroid Belt

Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disc between 1 and 3 au. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 au and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas drag, as has been proposed to explain the formation of close-in super-Earths. Here we show that the asteroid belt’s orbital excitation provides a crucial constraint against this scenario for the Solar system. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from discs of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Discs with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogues significantly more massive than the real planet. In contrast, a disc with a surface density gradient proportional to r(exp −5.5) reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disc profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.

Numerical methods↗

Computer program for calculating laminar, transitional, and turbulent boundary layers for a compressible axisymmetric flow

A finite-difference program is described for calculating the viscous compressible boundary layer flow over either planar or axisymmetric surfaces. The flow may be initially laminar and progress through a transitional zone to fully turbulent flow, or it may remain laminar, depending on the imposed boundary conditions, laws of viscosity, and numerical solution of the momentum and energy equations. The flow may also be forced into a turbulent flow at a chosen spot by the data input. The input may contain the factors of arbitrary Reynolds number, free-stream Mach number, free-stream turbulence, wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and wall roughness. The solution may start from an initial Falkner-Skan similarity profile, an approximate equilibrium turbulent profile, or an initial arbitrary input profile.

Albers, J. A.↗

Simulator study of a pictorial display for general aviation instrument flight

A simulation study of a computer drawn pictorial display involved a flight task that included an en route segment, terminal area maneuvering, a final approach, a missed approach, and a hold. The pictorial display consists of the drawing of boxes which either move along the desired path or are fixed at designated way points. Two boxes may be shown at all times, one related to the active way point and the other related to the standby way point. Ground tracks and vertical profiles of the flights, time histories of the final approach, and comments were obtained from time pilots. The results demonstrate the accuracy and consistency with which the segments of the flight are executed. The pilots found that the display is easy to learn and to use; that it provides good situation awareness, and that it could improve the safety of flight. The small size of the display, the lack of numerical information on pitch, roll, and heading angles, and the lack of definition of the boundaries of the conventional glide slope and localizer areas were criticized.

Adams, J. J.↗

Modelling of TGS growth in space

Attention is called to the necessity of programming the temperature downward in order to maintain a constant growth rate in the absence of convection. Trial and error numerical computations are performed to find the first approximation for the required isothermal dissolution period, linear ramp, and polynomial temperature variation period. The linear ramp rate is limited by the specified maximum temperature gradient. The isothermal dissolution turns out to be unnecessary. The linear ramp must be stopped before the polynomial period is begun in order to avoid overshooting the specified maximum growth rate. After a few hours, the temperature profile approaches steady-state behavior.

Liu, L. C.↗

Surface cracks in a plate of finite width under extension or bending

In this paper the problem of a finite plate containing collinear surface cracks is considered. The problem is solved by using the line spring model with plane elasticity and Reissner's plate theory. The main purpose of the study is to investigate the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors and to provide extensive numerical results which may be useful in applications. First, some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks and two corner cracks for wide range of relative dimensions. Particularly in corner cracks the agreement with the finite element solution is surprisingly very good. The results are obtained for semielliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.

Erdogan, F.↗

A new convective parameterization scheme for global models

A convective adjustment scheme is presented for parameterization of the vertical temperature and moisture structure for large-scale numerical weather models. The temperature and moisture structure is adjusted toward observed quasi-equilibrium structures to maintain a realistic vertical temperature and moisture structure in grid-scale global models during convection. Two reference thermodynamic structures are modeled analytically to account for shallow and deep convection. The thermodynamic structure obtained is permitted to change as a function of large-scale advective terms, radiation and surface fluxes. The reference profiles are applied when an energy integral provides a description of the shallow or deep convective processes occurring. Comparisons with GATE, BOMEX and ATEX data and a 10 day forecast for weather over the tropical ocean illustrate the climate and precipitation forecast accuracy improvements available with the adjustment scheme relative to standard parameterization schemes.

Miller, M. J.↗

On turbulent spots in plane Poiseuille flow

Turbulence characteristics inside a turbulent spot in plane Poiseuille flow are investigated by analyzing a database obtained from a direct numerical simulation. The spot is found to consist of two distinct regions - a turbulent area and a wave area. The flow inside the turbulent area has a strong resemblance to that found in the fully developed turbulent channel. Suitably defined mean and r.m.s. fluctuations as well as the internal shear-layer structures are found to be similar to the turbulent counterpart. In the wave area the inflexional mean spanwise profiles cause a rapid growth of oblique waves, which break down to turbulence. The breakdown process of the oblique waves is reminiscent of the secondary instability observed during transition to turbulence in channel and boundary-layer flows. Other detailed characteristics associated with the Poiseuille spot are presented and are compared with experimental results.

Henningson, Dan S.↗

On the three-dimensional instability of laminar boundary layers on concave walls

A study is made of the stability of laminar boundary-layer profiles on slightly curved walls relative to small disturbances that result from vortices whose axes are parallel to the principal direction of flow. The result is an eigenvalue problem by which, for a given undisturbed flow at a prescribed wall, the amplification or decay is computed for each Reynolds number and each vortex thickness. For neutral disturbances (zero amplification) a critical Reynolds number is determined for each vortex distribution. The numerical calculation produces amplified disturbances on concave walls only.

FLOW, LAMINAR↗

Propagation model for the Land Mobile Satellite channel in urban environments

This paper presents the major characteristics of a simulation package capable of performing a complete narrow and wideband analysis of the mobile satellite communication channel in urban environments for any given orbital configuration. The wavelength-to-average urban geometrical dimension ratio has required the use of the Geometrical Theory of Diffraction (GTD). For the RF frequency range, the model has been designed to be (1 up to 60 GHz) extended to include effects of non-perfect conductivity and surface roughness. Taking advantage of the inherent capabilities of such a high frequency method, we are able to provide a complete description of the electromagnetic field at the mobile terminal. Using the information made available at the ray-tracer and GTD solver outputs, the Land Mobile Satellite (LMS) urban model can also give a detailed description of the communication channel in terms of power delay profiles, Doppler spectra, channel scattering functions, and so forth. Statistical data, e.g. cumulative distribution functions, level crossing rates or distributions of fades are also provided. The user can access the simulation tool through a Design-CAD user-friendly interface by means of which she can effectively design her own urban layout and run consequently all the envisaged routines. The software is optimized in its execution time so that numerous runs can be achieved in a considerably short time.

Sforza, M.↗

Superfluid helium sloshing dynamics induced oscillations and fluctuations of angular momentum, force and moment actuated on spacecraft driven by gravity gradient or jitter acceleration associated with slew motion

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by the gravity gradient and jitter accelerations associated with slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with the Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) for slew motion which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid-vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers have also been derived. Examples are also given for cases applicable to the AXAF-S spacecraft sloshing dynamics associated with slew motion.

Hung, R. J.↗

Numerical studies of the surface tension effect of cryogenic liquid helium

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose of performing scientific observation during the normal spacecraft operation is investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics has been based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers, have been derived.

Hung, R. J.↗

The Response of Disks to Oscillatory Modes in Galaxies

Recent studies suggest that galaxies can oscillate in normal modes with essentially no damping over a Hubble time. These modes may play an important role in the structure and evolution of disk/halo systems. Motivated by the possibility that normal mode oscillations exist in real galaxies, we are investigating the response of galactic disks to halo oscillations. The goal of these investigations is to ascertain whether or not observational signatures exist for such oscillations. Our approach is to perform numerical experiments on the response of a self-gravitating disk to a time-varying halo potential. We assume that a significant fraction of the mass in a galaxy is in a dark halo. The halo oscillates and the luminous disk material responds to these oscillations. Preliminary results are reported for disks embedded in a radially oscillating gravitational potential. The equilibrium initial disk is represented by an exponential density profile. Considerable care was taken to build an initial disk model that was "stable" over long time scales. A control experiment was run with the disk in a static halo potential. The disk responds to the time-varying potential by developing a ring structure, which forms and disappears during each halo oscillation cycle. The density of stars becomes depressed in an annular region at the radius where the disk epicycle frequency is equal to the halo oscillation period. This pattern of response persists over time periods approaching a Hubble time. In the oscillating potential, a bar develops in the inner disk. This bar is absent when the halo remains static. Specific targets of this study include the implications for large-scale disk structure, the gas dynamical response of the interstellar medium in such systems, and the inflow of material into the central regions of the galaxy.

Smith, Bruce F.↗

A numerical investigation of thermal convection in a heat-generating fluid layer

Finite difference solutions of the equations governing thermal convection driven by uniform volumetric energy sources are presented for two-dimensional flows in a rectangular domain. The boundary conditions are a rigid (i.e, zero slip), zero heat-flux lower surface, rigid adiabatic sides, and either a rigid or free (i.e., zero shear) isothermal upper surface. Computations are carried out for Prandtl numbers from 0.05 to 20 and Rayleigh numbers from 5 x 10 to the 4th to 5 x 10 to the 8th. Nusselt numbers and average temperature profiles within the layer are in good agreement with experimental data for rigid-rigid boundaries. For rigid-free boundaries, Nusselt numbers are larger than in the former case. The structure of the flow and temperature fields in both cases is dominated by rolls, except at larger Rayleigh numbers where large-scale eddy transport occurs. Generally, low velocity upflows over broad regions of the layer are balanced by higher velocity downflows when the flow exhibits a cellular structure. The hydrodynamic constraint at the upper surface and the Prandtl number are found to influence only the detailed nature of flow and temperature fields. No truly steady velocity and temperature fields are found despite the fact that average Nusselt numbers reach steady values.

Emara, A. A.↗

Numerical study of low pressure nuclear thermal rockets

The flowfields and performance of low pressure nuclear thermal rockets, which use hydrogen as a propellant, are studied by solving the Navier-Stokes equations and the species equations. A finite-rate chemistry model is used in the species equations, and the turbulence is simulated by the Baldwin-Lomax turbulence model with a modified van Driest's damping constant. The calculated results for the chamber temperatures of 3200 K and 4000 K with a chamber pressure range of 0.1 atm to 6 atm are presented as contours, centerline variations, and exit profiles. The performance values from the present calculations, such as the vacuum specific impulse and thrust, are compared with those from the 1D, inviscid equilibrium and frozen flow code.

Kim, Suk C.↗