Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

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 631 records · Page 35

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

Pulsating Hydrodynamic Instability and Thermal Coupling in an Extended Landau/Levich Model of Liquid-Propellant Combustion: Inviscid Analysis - 1

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a non-zero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wavenumber, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a non-steady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, Stephen B.↗

Computer program for definition of transonic axial-flow compressor blade rows

A method is presented for designing axial-flow compressor blading from blade elements defined on cones which pass through the blade-edge streamline locations. Each blade-element centerline is composed of two segments which are tangent to each other. The centerline and surfaces of each segment have constant change of angle with path distance. The stacking line for the blade elements can be leaned in both the axial and tangential directions. The output of the computer program gives coordinates for fabrication and properties for aeroelastic analysis for planar blade sections. These coordinates and properties are obtained by interpolation across conical blade elements. The program is structured to be coupled with an aerodynamic design program.

Crouse, J. E.↗

A model for helicopter guidance on spiral trajectories

A point mass model is developed for helicopter guidance on spiral trajectories. A fully coupled set of state equations is developed and perturbation equations suitable for 3-D and 4-D guidance are derived and shown to be amenable to conventional state variable feedback methods. Control variables are chosen to be the magnitude and orientation of the net rotor thrust. Using these variables reference controls for nonlevel accelerating trajectories are easily determined. The effects of constant wind are shown to require significant feedforward correction to some of the reference controls and to the time. Although not easily measured themselves, the controls variables chosen are shown to be easily related to the physical variables available in the cockpit.

Mendenhall, S.↗

A time-dependent theoretical model of the polar wind Preliminary results

The coupled, time-dependent continuity, momentum, and energy equations of a two-ion O(+) and H(+) quasi-neutral plasma were solved in order to extend understanding of polar wind behavior. This numerical code allows studies of the time dependent behavior of polar wind-type flows into and out of the ionosphere. Initial studies indicate that the typical time constants for electron and ion temperature changes are of the order of minutes and tens of minutes, respectively. The response time of the minor high altitude ion O(+) is less than an hour, whereas that of the major ion, H(+), is many hours. The initial test runs also demonstrate the fact that temporary supersonic flows of both O(+) and H(+) are possible, especially in the presence of significant ion heating.

Gombosi, T. I.↗

A spectral view of nonlinear fluxes and stationary-transient interaction in the atmosphere

Nonlinear spectral transfers of kinetic energy and enstrophy, and stationary-transient interaction, are studied using global FGGE data for January 1979. The spectral transfers are found to be dominated by two dynamically distinct components, one due to nonlinear self-interactions within the transient flow and the other arising from coupled interactions between the stationary and transient parts of the flow. The former seems to be attributable to two-dimensional homogeneous turbulence, while the second cannot be understood within that theory. It is argued that the mixed stationary-transient interactions can be understood physically as a process of shear-induced spectral transfer of transient enstrophy along lines of constant zonal wavenumber. The transients act mainly to strengthen the stationary flow, and this 'forcing' operates on a timescale comparable to presumed viscous timescales.

Shepherd, Theodore G.↗

Time-dependent modulation of galactic cosmic rays by merged interaction regions

Models that solve the one-dimensional, solar modulation equation have reproduced the 11-year galactic cosmic ray using functional representations of global merged interaction regions (MIRs). This study extends those results to the solution of the modulation equation with explicit time dependence. The magnetometers on Voyagers 1 and 2 provide local magnetic field intensities at regular intervals, from which one calculates the ratio of the field intensity to the average local field. These ratios in turn are inverted to form diffusion coefficients. Strung together in radius and time, these coefficents then fall and rise with the strength of the interplanetary magnetic field, becoming representations of MIRs. These diffusion coefficients, calculated locally, propagate unchanged from approx. 10 AU to the outer boundary (120 AU). Inside 10 AU, all parameters, including the diffusion coefficient are assumed constant in time and space. The model reproduces the time-intensity profiles of Voyager 2 and Pioneer 10. Radial gradient data from 1982-1990 between Pioneer 10 and Voyager 2 are about the same magnitude as those calculated in the model. It is also shows agreement in rough magnitude with the radial gradient between Pioneer 10 and 1 AU. When coupled with enhanced, time-dependent solar wind speed at the probe's high latitude, as measured by independent observers, the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source for the observed negative latitudinal gradients. The model exhibits the 11-year cyclical cosmic ray intensity behavior at all radii, including 1 AU, not just at the location of the spacecraft where the magnetic fields are measured. In addition, the model's point of cosmic ray maximum correctly travels at the solar wind speed, illustrating the well-known propagation of modulation. Finally, at least in the inner heliosphere this model accounts for the delay experienced by lower-rigidity protons in reaching their time-intensity peak. The actual delays in this model, however, are somewhat smaller than the data. In the outer heliosphere the models sees no delays, and the data are ambiguous as to their existence. It appears that strong magnetic field compression regions (merged interaction regions) that are 3-4 times the average field strength can, at least in a helioequatorial band, disrupt effects, such as drifts, that could dominate in quieter magnetic fields. The question remains: Is the heliosphere ever quiet enough to allow such effects to be unambiguously measured, at least in the midlatitudes?

Perko, J. S.↗

Interaction of Ethyl Alcohol Vapor with Sulfuric Acid Solutions

We investigated the uptake of ethyl alcohol (ethanol) vapor by sulfuric acid solutions over the range approx.40 to approx.80 wt % H2SO4 and temperatures of 193-273 K. Laboratory studies used a fast flow-tube reactor coupled to an electron-impact ionization mass spectrometer for detection of ethanol and reaction products. The uptake coefficients ((gamma)) were measured and found to vary from 0.019 to 0.072, depending upon the acid composition and temperature. At concentrations greater than approx.70 wt % and in dilute solutions colder than 220 K, the values approached approx.0.07. We also determined the effective solubility constant of ethanol in approx.40 wt % H2SO4 in the temperature range 203-223 K. The potential implications to the budget of ethanol in the global troposphere are briefly discussed.

sulfate aerosols↗

The Attenuation of a Detonation Wave by an Aircraft Engine Axial Turbine Stage

A Constant Volume Combustion Cycle Engine concept consisting of a Pulse Detonation Combustor (PDC) followed by a conventional axial turbine was simulated numerically to determine the attenuation and reflection of a notional PDC pulse by the turbine. The multi-stage, time-accurate, turbomachinery solver TURBO was used to perform the calculation. The solution domain consisted of one notional detonation tube coupled to 5 vane passages and 8 rotor passages representing 1/8th of the annulus. The detonation tube was implemented as an initial value problem with the thermodynamic state of the tube contents, when the detonation wave is about to exit, provided by a 1D code. Pressure time history data from the numerical simulation was compared to experimental data from a similar configuration to verify that the simulation is giving reasonable results. Analysis of the pressure data showed a spectrally averaged attenuation of about 15 dB across the turbine stage. An evaluation of turbine performance is also presented.

VanZante, Dale↗

Solid-state current transformer

A signal transformation network which is uniquely characterized to exhibit a very low input impedance while maintaining a linear transfer characteristic when driven from a voltage source and when quiescently biased in the low microampere current range is described. In its simplest form, it consists of a tightly coupled two transistor network in which a common emitter input stage is interconnected directly with an emitter follower stage to provide virtually 100 percent negative feedback to the base input of the common emitter stage. Bias to the network is supplied via the common tie point of the common emitter stage collector terminal and the emitter follower base stage terminal by a regulated constant current source, and the output of the circuit is taken from the collector of the emitter follower stage.

Farnsworth, D. L.↗

Optimization of shallow arches against instability using sensitivity derivatives

The author discusses the problem of optimization of shallow frame structures which involve a coupling of axial and bending responses. A shallow arch of a given shape and of given weight is optimized such that its limit point load is maximized. The cross-sectional area, A(x) and the moment of inertia, I(x) of the arch obey the relationship I(x) = rho A(x) sup n, n = 1,2,3 and rho is a specified constant. Analysis of the arch for its limit point calculation involves a geometric nonlinear analysis which is performed using a corotational formulation. The optimization is carried out using a second-order projected Lagrangian algorithm and the sensitivity derivatives of the critical load parameter with respect to the areas of the finite elements of the arch are calculated using implicit differentation. Results are presented for an arch of a specified rise to span ratio under two different loadings and the limitations of the approach for the intermediate rise arches are addressed.

Kamat, Manohar P.↗

Reference level winds from balloon platforms

The superpressure balloon was developed to provide a method of obtaining global winds at all altitudes from 5 to 30 km. If a balloon could be made to fly for several weeks at a constant altitude, and if it could be tracked accurately on its global circuits, the balloon would provide a tag for the air parcel in which it was embedded. The Lagrangian data on the atmospheric circulation would provide a superior data input to the numerical model. The Global Atmospheric Research Program (GARP) was initiated in large part based on the promise of this technique coupled with free-floating ocean buoys and satellite radiometers. The initial name proposed by Charney for GARP was SABABURA 'SAtellite BAlloon BUoy RAdiometric system' (Charney, 1966). However, although the superpressure balloon exceeded its designers' expectations for flight duration in the stratosphere (longest flight duration of 744 days), flight duration below 10 km was limited by icing in super-cooled clouds to a few days. The balloon was relegated to a secondary role during the GARP Special Observing Periods. The several major superpressure balloon programs for global wind measurement are described as well as those new developments which make the balloon once again an attractive vehicle for measurement of global winds as a reference and bench-mark system for future satellite systems.

Lally, Vincent E.↗

Effect of Dynamic Rolling Oscillations on Twin Tail Buffet Response

The effect of dynamic rolling oscillations of delta-wing/twin-tail configuration on twin-tail buffet response is investigated. The computational model consists of a sharp-edged delta wing of aspect ratio one and swept-back flexible twin tail with taper ratio of 0.23. The configuration model is statically pitched at 30 deg. angle of attack and then forced to oscillate in roll around the symmetry axis at a constant amplitude of 4 deg. and reduced frequency of pi and 2(pi). The freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. This multidisciplinary problem is solved using three sets of equations on a dynamic multi-block grid structure. The first set is the unsteady, full Navier-Stokes equations, the second set is the aeroelastic equations for coupled bending and torsion vibrations of the tails, and the third set is the grid-displacement equations. The configuration is investigated for inboard position of the twin tails which corresponds to a separation distance between the twin tails of 33% wing span. The computed results are compared with the results of stationary configuration, which previously have been validated using experimental data. The results conclusively showed that the rolling oscillations of the configuration have led to higher loads, higher deflections, and higher excitation peaks than those of the stationary configuration. Moreover, increasing the reduced frequency has led to higher loads and excitation peaks and lower bending and torsion deflections and acceleration.

Sheta, Essam F.↗

Nystagmus induced by electrical stimulation of the vestibular and prepositus hypoglossi nuclei in the monkey: evidence for site of induction of velocity storage

Electrical stimulation of the vestibular nuclei (VN) and prepositus hypoglossi nuclei (PPH) of alert cynomolgus monkeys evoked nystagmus and eye deviation while they were in darkness. At some sites in VN, nystagmus and after-nystagmus were induced with characteristics suggesting that velocity storage had been excited. We analyzed these responses and compared them to the slow component of optokinetic nystagmus (OKN) and to optokinetic after-nystagmus (OKAN). We then recorded unit activity in VN and determined which types of nystagmus would be evoked from the sites of recording. Nystagmus and eye deviations were also elicited by electrical stimulation of PPH, and we characterized the responses where unit activity was recorded in PPH. Horizontal slow phase velocity of the VN "storage" responses was contralateral to the side of stimulation. The rising time constants and peak steady-state velocities were similar to those of OKN, and the falling time constants of the after-nystagmus and of OKAN were approximately equal. Both the induced after-nystagmus and OKAN were habituated by stimulation of the VN. When horizontal after-nystagmus was evoked with animals on their sides, it developed yaw and pitch components that tended to shift the vector of the slow phase velocity toward the spatial vertical. Similar "cross-coupling" occurs for horizontal OKAN or for vestibular post-rotatory nystagmus elicited in tilted positions. Thus, the storage component of nystagmus induced by VN stimulation had the same characteristics as the slow component of OKN and the VOR. Positive stimulus sites for inducing nystagmus with typical storage components were located in rostral portions of VN. They lay in caudal ventral superior vestibular nucleus (SVN), dorsal portions of central medial vestibular nucleus (MVN) caudal to the abducens nuclei and in adjacent lateral vestibular nucleus (LVN). More complex stimulus responses, but with contralateral after-nystagmus, were induced from surrounding regions of ventral MVN and LVN, rostral descending vestibular nucleus and the marginal zone between MVN and PPH. Vestibular-only (VO), vestibular plus saccade (VPS) and tonic vestibular pause (TVP) units were identified by extracellular recording. Stimulation near type I lateral and vertical canal-related VO units elicited typical "storage" responses with after-nystagmus in 23 of 29 tracks (79%). Stimulus responses were more complex from the region of neurons with oculomotor-related signals, i.e., TVP or VPS cells, although after-nystagmus was also elicited from these sites. Effects of vestibular nerve and nucleus stimulation were compared.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗

Doppler phase shifting using dual, switched phase shifting devices

A system of inducing a phase shift using moving reflector elements. The moving reflectors can be moving mirrors or an acousto-optical filter. The moving reflectors oscillate i.e. the move first in a first direction and then in a second direction. Two different reflectors are used so that the light can be switched between the reflectors. During a first portion of the cycle the light is coupled to the first modulator which moves the reflector in the first direction. The second modulator is out of phase with the first modulator, and the light is switched to that second modulator during a second portion of the cycle. The second modulator is also moving in the first direction when the light is applied thereto. In this way, the light obtains a constant direction Doppler shift.

Gutierrez, Roman C.↗