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 379 records · Page 21

The Aqua-planet Experiment (APE): Response to Changed Meridional SST Profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea- ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

AGCM↗

Challenges in Modelling of Lightning-Induced Delamination; Effect of Temperature-Dependent Interfacial Properties

Lightning is a major cause of damage in laminated composite aerospace structures during flight. Due to the dielectric nature of Carbon fiber reinforced polymers (CFRPs), the high energy induced by lightning strike transforms into extreme, localized surface temperature accompanied with a high-pressure shockwave resulting in extensive damage. It is crucial to develop a numerical tool capable of predicting the damage induced from a lightning strike to supplement extremely expensive lightning experiments. Delamination is one of the most significant failure modes resulting from a lightning strike. It can be extended well beyond the visible damage zone, and requires sophisticated techniques and equipment to detect. A popular technique used to model delamination is the cohesive zone approach. Since the loading induced from a lightning strike event is assumed to consist of extreme localized heating, the cohesive zone formulation should additionally account for temperature effects. However, the sensitivity to this dependency remains unknown. Therefore, the major focus point of this work is to investigate the importance of this dependency via defining various temperature dependency profiles for the cohesive zone properties, and analyzing the corresponding delamination area. Thus, a detailed numerical model consisting of multidirectional composite plies with temperature-dependent cohesive elements in between is subjected to lightning (excessive amount of heat and pressure) and delamination/damage expansion is studied under specified conditions.

interlaminar damage↗

Calculation of turbulent shear stress in supersonic boundary layer flows

An analysis of turbulent boundary layer flow characteristics and the computational procedure used are discussed. The integrated mass and momentum flux profiles and differentials of the integral quantities are used in the computations so that local evaluation of the streamwise velocity gradient is not necessary. The computed results are compared with measured shear stress data obtained by using hot wire anemometer and laser velocimeter techniques. The flow measurements were made upstream and downstream of an adiabatic unseparated interaction of an oblique shock wave with the turbulent boundary layer on the flat wall of a two dimensional wind tunnel. A comparison of the numerical analysis and actual measurements is made and the effects of small differences in mean flow profiles on the computed shear stress distributions are discussed.

Sun, C. C.↗

Creating Synthetic Coronal Observational Data From MHD Models: The Forward Technique

We present a generalized forward code for creating simulated corona) observables off the limb from numerical and analytical MHD models. This generalized forward model is capable of creating emission maps in various wavelengths for instruments such as SXT, EIT, EIS, and coronagraphs, as well as spectropolari metric images and line profiles. The inputs to our code can be analytic models (of which four come with the code) or 2.5D and 3D numerical datacubes. We present some examples of the observable data created with our code as well as its functional capabilities. This code is currently available for beta-testing (contact authors), with the ultimate goal of release as a SolarSoft package

Rachmeler, Laurel A.↗

Numerical studies of motion and decay of vortex filaments

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of large axial velocity and of nonsimilar initial profiles in vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of asymptotic analysis, which is the effective vortical core size is much smaller than all other length scales in the flowfield, e.g., the radius of curvature and interfilament distance. The computations are continued until the basic assumption is no longer valid, that is, when the merging or intersection of filaments have begun. Various types of local or global merging or intersection of filaments are classified and demonstrated by numerical examples.

Liu, C. H.↗

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

Tropospheric Wind Profile Measurements with a Direct Detection Doppler Lidar

Research has established the importance of global tropospheric wind measurements for large scale improvements in numerical weather prediction. In addition, global wind measurements provide data that are fundamental to the understanding and prediction of global climate change. These tasks are closely linked with the goals of the NASA Earth Science Enterprise and Global Climate Change programs. NASA Goddard has been actively involved in the development of direct detection Doppler lidar methods and technologies to meet the wind observing needs of the atmospheric science community. In this paper we describe a recently developed prototype wind lidar system using a direct detection Doppler technique for measuring wind profiles from the surface through the troposphere. This system uses a pulsed ND:YAG laser operating at 1064 nm as the transmitter. The laser pulse is directed to the atmosphere using a 40 cm diameter scan mirror. The portion of the laser energy backscattered from aerosols and molecules is collected by a 40 cm diameter telescope and coupled via fiber optics into the Doppler receiver. Single photon counting APD's are used to detect the atmospheric backscattered signal. The principle element of the receiver is a dual bandpass tunable Fabry Perot etalon which analyzes the Doppler shift of the incoming laser signal using the double edge technique. The double edge technique uses two high resolution optical filters having bandpasses offset relative to one another such that the 'edge' of the first filter's transmission function crosses that of the second at the half power point. The outgoing laser frequency is located approximately at the crossover point. Due to the opposite going slopes of the edges, a Doppler shift in the atmospheric backscattered laser frequency produces a positive change in signal for one filter and a negative change in the second filter. Taking the ratio of the two edge channel signals yields a result which is directly proportional to the component of the wind along the line-of-sight of the laser. Measuring the radial wind in several directions provides sufficient information to determine the true wind speed and direction. The lidar has operated from our laboratory at Goddard since June, 1997. Wind profiles have been obtained to altitudes of 12 km with a vertical resolution of 330 in. Vector wind data are obtained by rotating the scan mirror to measure line-of-sight wind profiles for at least two azimuth angles at an elevation angle of 45 degrees. The precision of the data as determined from the standard deviation of multiple independent lidar profiles is in the range of 1 to 3 m/sec up to 10 km. Good agreement is obtained when the lidar data are compared with the upper air rawinsonde soundings taken at Dulles airport. Examples of the wind lidar data will be presented along with a description of the instrument and future developments.

Gentry, Bruce M.↗

Remote temperature sounding in the presence of cloud by zenith scanning

The principles of remote sounding of vertical temperature profiles in planetary atmospheres containing cloud layers are presented for an approach that uses measurements at different angles to the local vertical to eliminate the unknown cloud effects. A numerical example for a model of the earth's atmosphere shows that soundings separated by 40 deg in zenith angle permit the recovery of the temperature profile with an error that is not significantly greater than that obtained in the clear atmosphere case, and without serious horizontal smear. It is proposed that the method will be valuable for the interpretation of measurements made by instruments on spin-stabilized spacecraft on planetary missions, since in these cases the required angular coverage is obtained automatically for most trajectories.

Taylor, F. W.↗

An experimental and computational study of turbulent boundary layers in oscillating flows

An experimental-computational study of the behavior of turbulent boundary layers for oscillating air flows over a plane surface with a small favorable mean pressure gradient is described. Experimental data were obtained for boundary layers generated on the test section wall of a new type of facility that produces a mean flow velocity with a superposed sinusoidal component over a wide frequency range. Experimental results for unsteady velocity profiles and turbulence intensity obtained from hot wire measurements at a fixed stream-wise position on the surface are presented for frequencies ranging from 5 to 21 Hz. Numerical predictions were obtained using an unsteady boundary layer computational code and the Cebeci-Smith and Glushko turbulence models. Predicted quantities related to unsteady velocity profiles exhibit fair agreement with experiment when the Cebeci-Smith model is used.

Cook, W. J.↗

Optimal Area Profiles for Ideal Single Nozzle Air-Breathing Pulse Detonation Engines

The effects of cross-sectional area variation on idealized Pulse Detonation Engine performance are examined numerically. A quasi-one-dimensional, reacting, numerical code is used as the kernel of an algorithm that iteratively determines the correct sequencing of inlet air, inlet fuel, detonation initiation, and cycle time to achieve a limit cycle with specified fuel fraction, and volumetric purge fraction. The algorithm is exercised on a tube with a cross sectional area profile containing two degrees of freedom: overall exit-to-inlet area ratio, and the distance along the tube at which continuous transition from inlet to exit area begins. These two parameters are varied over three flight conditions (defined by inlet total temperature, inlet total pressure and ambient static pressure) and the performance is compared to a straight tube. It is shown that compared to straight tubes, increases of 20 to 35 percent in specific impulse and specific thrust are obtained with tubes of relatively modest area change. The iterative algorithm is described, and its limitations are noted and discussed. Optimized results are presented showing performance measurements, wave diagrams, and area profiles. Suggestions for future investigation are also discussed.

Paxson, Daniel E.↗

Comparing Experiment and Computation of Hypersonic Laminar Boundary Layers with Isolated Roughness

Streamwise velocity profile behavior in a hypersonic laminar boundary layer in the presence of an isolated roughness element is presented for an edge Mach number of 8.2. Two different roughness element types are considered: a 2-mm tall, 4-mm diameter cylinder, and a 2-mm radius hemisphere. Measurements of the streamwise velocity behavior using nitric oxide (NO) planar laser-induced fluorescence (PLIF) molecular tagging velocimetry (MTV) have been performed on a 20-degree wedge model. The top surface of this model acts as a flat-plate and is oriented at 5 degrees with respect to the freestream flow. Computations using direct numerical simulation (DNS) of these flows have been performed and are compared to the measured velocity profiles. Particular attention is given to the characteristics of velocity profiles immediately upstream and downstream of the roughness elements. In these regions, the streamwise flow can experience strong deceleration or acceleration. An analysis in which experimentally measured MTV profile displacements are compared with DNS particle displacements is performed to determine if the assumption of constant velocity over the duration of the MTV measurement is valid. This assumption is typically made when reporting MTV-measured velocity profiles, and may result in significant errors when comparing MTV measurements to computations in regions with strong deceleration or acceleration. The DNS computations with the cylindrical roughness element presented in this paper were performed with and without air injection from a rectangular slot upstream of the cylinder. This was done to determine the extent to which gas seeding in the MTV measurements perturbs the boundary layer flowfield.

Bathel, Brett F.↗

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the 1-D Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation for EC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with non-fusing hydrogen-1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗

Titan atmosphere models, 1973

The composition and structure of the atmosphere of Titan, based on theory and on spectroscopic and infrared data, is reviewed for the development of numerical engineering models. Light, nominal, and heavy atmospheres are described and tabulated, and their profiles of radius, temperature, pressure, and density are illustrated. Corresponding descriptions of atmospheric dynamics, condensates and surfaces are outlined.

Divine, N.↗

Foam-machining tool with eddy-current transducer

Three-cutter machining system for foam-covered tanks incorporates eddy-current sensor. Sensor feeds signal to numerical controller which programs rotational and vertical axes of sensor travel, enabling cutterhead to profile around tank protrusions.

Copper, W. P.↗

Numerical simulation of tornado wind loading on structures

A numerical simulation of a tornado interacting with a building was undertaken in order to compare the pressures due to a rotational unsteady wind with that due to steady straight winds used in design of nuclear facilities. The numerical simulations were performed on a two-dimensional compressible hydrodynamics code. Calculated pressure profiles for a typical building were then subjected to a tornado wind field and the results were compared with current quasisteady design calculations. The analysis indicates that current design practices are conservative.

Maiden, D. E.↗

Remote sounding of cloudy atmospheres. II - Multiple cloud formations

The dual frequency-range principle developed in Part I (Chahine, 1974) for infrared remote sounding of atmospheric temperature profiles in the presence of a single cloud layer is extended here to the case of multiple cloud formations. The approach requires no a priori knowledge of the spectral properties of the clouds or the number of cloud layers in the fields of view. The method of solution requires measurements over adjacent fields of view and leads to the determination of the clear-column atmospheric temperature profiles with the same degree of accuracy and vertical resolution permitted under cloudless conditions. Numerical verifications are carried out to illustrate the stability and accuracy of the method using simulated radiance data, from the 4.3 and 15 micron CO2 bands in the terrestrial atmosphere, in the presence of up to three cloud layers in the fields of view.

Chahine, M. T.↗

Acoustic resonances and sound scattering by a shear layer

The energy reflection coefficient is evaluated numerically for plane waves incident on a plane shear layer having a linear velocity profile. The shear layer is found to exhibit no resonances and no Brewster angles. The behavior of the reflection coefficient depends crucially on the parameter tau, a nondimensional measure of the disturbance Strouhal number with respect to the disturbance Mach number in the mean flow direction. For moderate values of tau, the amplified reflection regime degenerates into the total reflection one, whereas in the ordinary reflection regime the variation of the reflection coefficient with tau depends on whether or not the corresponding vortex sheet has a Brewster angle. The results indicate that caution should be exercised in uncritically modeling a finite thickness shear layer by a corresponding vortex sheet.

Koutsoyannis, S. P.↗

Ultraviolet observations of 31 and 32 Cygni

Some observations of 31 and 32 Cygni, started in September 1978, at phases 0.62 and 0.17, respectively, are presented. The UV spectrum of 32 Cyg is described both qualitatively and quantitatively. Numerous lines of Si II, O I, C II, Al II and III and Fe II appear with P Cygni characteristics. The LWR-2275 image is dominated by the Mg II 2800 A resonance doublet and numerous Fe II lines in the 2200 - 2700 A region, all showing P Cygni profiles. More quantitatively, ionization modeling for the B star environs was made by comparing observed and computed equivalent widths to deduce density, temperature and turbulence values. At the projected separation between the stars, the log particle density obtained was 7.2, neglecting dust and molecule formation.

Stencel, R. E.↗