SEARCH · Search NASA
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.
Experimental data on the interaction mean free path of hadrons with energy E 20 TeV in lead obtained with multilayer X-ray emulsion chambers(XEC) are compared with results of simulation of nuclear-electromagnetic cascades in lead chamber. It is shown that, to explain experimental data, the value of the inelasticity coefficients, K .8 should be assumed.
The current computational methods for solving transport equations of turbulent reacting single-phase flows are critically reviewed, with primary attention given to those methods that lead to model-free simulations. In particular, consideration is given to direct numerical simulations using spectral (Galerkin) and pseudospectral (collocation) methods, spectral element methods, and Lagrangian methods. The discussion also covers large eddy simulations and turbulence modeling.
Work is devoted to study of free-stream vorticity normal to leading edge interaction with boundary layer over plate and resulting flow distortion influence on laminar-turbulent transition. In experiments made the wake behind the vertically stretched wire was used as a source of vortical disturbances and its effect on the boundary layer over the horizontally mounted plate with various leading edge shapes was investigated. The purpose of experiments was to check the predictions of theoretical works of M.E. Goldstein, et. al. This theory shows that small free-stream inhomogeneity interacting with leading edge produces considerable distortion of boundary layer flow. In general, results obtained confirms predictions of Goldstein's theory, i.e., the amplification of steady vortical disturbances in boundary layer caused by vortex lines stretching was observed. Experimental results fully coincide with predictions of theory for large Reynolds number, relatively sharp leading edge and small disturbances. For large enough disturbances the flow distortion caused by symmetric wake unexpectedly becomes antisymmetric in spanwise direction. If the leading edge is too blunt the maximal distortion takes place immediately at the nose and no further amplification was observed. All these conditions and results are beyond the scope of Goldstein's theory.
Heating requirements for satisfactory cyclic de-icing over a wide range of icing and operating conditions have been determined for a gas-heated, 36deg swept airfoil of 6.9-foot chord with a partial-span leading-edge slat. Comparisons of heating requirements and effectiveness were made between the slatted and unslatted portions of the airfoil. Studies were also made comparing cyclic de-icing with continuous anti-icing, and cycll.cde-icing systems with and without leading-edge ice-free parting strips. De-icing heat requirements were approximately the same with either heated or unheated parting strips because of the aerodynamic effects of the 36deg sweep angle and the spanwise saw-tooth profile of leading-edge glaze-ice deposits. Cyclic de-icing heat-source requirements were found to be one-fourth or less of the heat requirements for complete anti-icing. The primary factors that affected the performance of the cyclic de-icing heating system were ambient air temperature, heat distribution, and thermal lag.
In this study, we have assessed the effectiveness of the use of existing observing systems in the lower troposphere in the GEOS hybrid–4DEnVar data assimilation system through a set of observing system experiments. The results show that microwave radiances have a large impact in the Southern Hemisphere and Tropical ocean, but the large influence is mostly observed above 925 hPa and dissipates relatively quickly with longer forecast lead times. Conventional data information holds better in the forecast ranging from the surface to 100 hPa, depending on the field evaluated, in the Northern Hemisphere and lowest model levels in the Tropics. Infrared radiances collectively have much less impact in the lower troposphere. Removing surface observations has small but persistent impact on specific humidity in the upper atmosphere, but small or negligible impact on planetary boundary layer (PBL) height and temperature. The model responses to the incremental analysis update (IAU) forcing are also analyzed. In the IAU assimilation window, the physics responds strongly to the IAU forcing in the lower troposphere, and the changes of physics tendency in the lower troposphere and hydrodynamics tendency in the mid- and upper troposphere are viewed as beneficial to the reduction of state error covariance. In the subsequent forecast, the model tendencies continue to deviate further from the original free forecast with forecast lead times around 300–400 hPa, but physics tendency has showed signs of returning to its original free forecast mechanisms at 1-day forecast in the lower troposphere.
By rolling up the surfaces of discontinuity originating from the leading edge of delta wings, free vortex cores are formed above the wing. In case of greater angles of incidence, the flow in these vortex cores shows an instability which abruptly produces strong turbulence. In the present paper an explanation is given of this instability being a "frictionless instability" of the vortex core flow by increasing helical interference vortices. The occurring vortex core flows are calculated and investigated for stability by means of a stability criterion concerning flows with helical streamlines given by H. Ludwieg.
Experimental and theoretical studies of low speed leading edge boundary layer receptivity to free-stream vorticity produced by upstream wires normal to the leading edge are discussed. Data include parametric variations in leading edge configuration and details of the incident disturbance field including single and multiple wakes. The induced disturbance amplitude increases with increases in the leading edge diameter and wake interactions. Measurements agree with the theory of M. E. Goldstein.
Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].
The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.
An analysis is presented of the effects of airfoil thickness and mean aerodynamic loading on boundary-layer receptivity in the leading-edge region. The case of acoustic free-stream disturbances, incident on a thin cambered airfoil with a parabolic leading edge in a low Mach number flow, is considered. An asymptotic analysis based on large Reynolds number is developed, supplemented by numerical results. The airfoil thickness distribution enters the theory through a Strouhal number based on the nose radius of the airfoil, S = (omega)tau(sub n)/U, where omega is the frequency of the acoustic wave and U is the mean flow speed. The influence of mean aerodynamic loading enters through an effective angle-of-attack parameter ti, related to flow around the leading edge from the lower surface to the upper. The variation of the receptivity level is analyzed as a function of S, mu, and characteristics of the free-stream acoustic wave. For an unloaded leading edge, a finite nose radius dramatically reduces the receptivity level compared to that for a flat plate, the amplitude of the instability waves in the boundary layer being decreased by an order of magnitude when S = 0.3. Modest levels of aerodynamic loading are found to further decrease the receptivity level for the upper surface of the airfoil, while an increase in receptivity level occurs for the lower surface. For larger angles of attack close to the critical angle for boundary layer separation, a local rise in the receptivity level occurs for the upper surface, while for the lower surface the receptivity decreases. The effects of aerodynamic loading are more pronounced at larger values of S. Oblique acoustic waves produce much higher receptivity levels than acoustic waves propagating downstream parallel to the airfoil chord.
The problem concerned in this work is that of calculating magnetic field configurations in which the Lorentz force (vector)j x (vector)B is everywhere zero, subject to specified boundary conditions. The magnetic field is represented in terms of Clebsch variables in the form (vector)B = del alpha x del beta. These variables are constant on any field line. The most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. It is proposed that the field lines move in the direction of local Lorentz force and relax towards a force-free field configuration. This concept leads to an iteration procedure for modifying the variables alpha and beta that tends asymptotically towards the force-free state. This method is first applied to a simple problem in two rectangular dimensions; the calculation shows that the convergence of magnetic field energy to a minimum state (force-free) is close to exponential. This method is then applied to study some astrophysical force-free magnetic fields, such as the structures and evolution of magnetic fields of rotating sunspots and accretion disks. The implication of the results, as related to the mechanisms of solar flares, extragalactic radio sources and radio jets, are discussed.
A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, finite-volume, spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free to-roll case by simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.
In this brief paper we examine a Riccati equation decomposition due to Reid and Lainiotis and apply the result to the continuous time-invariant linear filtering problem. Exploitation of the time-invariant structure leads to integration-free covariance recursions which are of use in covariance analyses and in filter implementations. A super-linearly convergent iterative solution to the algebraic Riccati equation (ARE) is developed. The resulting algorithm, arranged in a square-root form, is thought to be numerically stable and competitive with other ARE solution methods. Certain covariance relations that are relevant to the fixed-point and fixed-lag smoothing problems are also discussed.
Explore the source record for details and available documents.
Electrophoresis has long been recognized as an effective analytic technique for the separation of proteins and other charged species, however attempts at scaling up to accommodate commercial volumes have met with limited success. In this report we describe a novel electrophoretic separation technique - Binary Oscillatory Crossflow Electrophoresis (BOCE). Numerical simulations indicate that the technique has the potential for preparative scale throughputs with high resolution, while simultaneously avoiding many problems common to conventional electrophoresis. The technique utilizes the interaction of an oscillatory electric field and a transverse oscillatory shear flow to create an active binary filter for the separation of charged protein species. An oscillatory electric field is applied across the narrow gap of a rectangular channel inducing a periodic motion of charged protein species. The amplitude of this motion depends on the dimensionless electrophoretic mobility, alpha = E(sub o)mu/(omega)d, where E(sub o) is the amplitude of the electric field oscillations, mu is the dimensional mobility, omega is the angular frequency of oscillation and d is the channel gap width. An oscillatory shear flow is induced along the length of the channel resulting in the separation of species with different mobilities. We present a model that predicts the oscillatory behavior of charged species and allows estimation of both the magnitude of the induced convective velocity and the effective diffusivity as a function of a in infinitely long channels. Numerical results indicate that in addition to the mobility dependence, the steady state behavior of solute species may be strongly affected by oscillating fluid into and out of the active electric field region at the ends of the cell. The effect is most pronounced using time dependent shear flows of the same frequency (cos((omega)t)) flow mode) as the electric field oscillations. Under such conditions, experiments indicate that solute is drawn into the cell from reservoirs at both ends of the cell leading to a large mass build up. As a consequence, any initially induced mass flux will vanish after short times. This effect was not captured by the infinite channel model and hence numerical and experimental results deviated significantly. The revised model including finite cell lengths and reservoir volumes allowed quantitative predictions of the time history of the concentration profile throughout the system. This latter model accurately describes the fluxes observed for both oscillatory flow modes in experiments using single protein species. Based on the results obtained from research funded under NASA grant NAG-8-1080.S, we conclude that binary separations are not possible using purely oscillatory flow modes because of end effects associated with the cos((omega)t) mode. Our research shows, however, that a combination of cos(2(omega)t) and steady flow should lead to efficient separation free of end effects. This possibility is currently under investigation.
A small total-pressure tube resting against a flat-plate surface was used as a Stanton tube and calibrated as a skin-friction meter at various subsonic and supersonic speeds. Laminar flow was maintained for the supersonic runs at a Mach number M(sub infinity) of 2. At speeds between M(sub infinity) = 1.33 and M(sub infinity) = 1.87, the calibrations were carried-out in a turbulent boundary layer. The subsonic flows were found to be in transition. The skin-friction readings of a floating-element type of balance served as the reference values against which the Stanton tube was calibrated. A theoretical model was developed which, for moderate values of the shear parameter tau, accurately predicts the performance of the Stanton tube in subsonic and supersonic flows. A "shear correction factor" was found to explain the deviations from the basic model when T became too large. Compressibility effects were important only in the case of turbulent supersonic flows, and they did not alter the form of the calibration curve. The test Reynolds numbers, based on the distance from the leading edge and free-stream conditions, ranged from 70,000 to 875,000. The turbulent-boundary-layer Reynolds numbers, based on momentum thickness, varied between 650 and 2,300. Both laminar and turbulent velocity profiles were taken and the effect of pressure gradient on the calibration was investigated.
Tropospheric aerosols are of considerable environmental importance. They modify the radiative budget of Earth by scattering and absorbing radiation, and by providing nuclei for cloud formation. Additionally, they provide surfaces for heterogeneous and multiphase reactions that affect tropospheric chemistry. For example, Dentener and Crutzen (1993) showed that reactions of N2O5 and NO3 with sulfate aerosols may significantly alter the tropospheric concentrations of NO(x), O3, and OH by converting NOx to HNO3 which is rapidly removed by precipitation. Zhang et al. (1994) assumed these same reactions would occur on dust aerosols and showed that dust outbreaks may reduce NO(x) levels by up to 50%. Dentener et al. (1996) studied the possible effect of reactions on dust on sulfate, nitrate, and O3 concentration. Heterogeneous and multiphase reactions on aerosols may also perturb the sulfur cycle the chlorine cycle and the bromine cycle. Because these reactions can release free chlorine and free bromine they might lead to the destruction of ozone in the marine boundary layer that may be important to include in models of tropospheric chemistry. The goal of our proposed work is to examine the role of heterogeneous and multiphase reactions in the tropospheric cycles of reactive nitrogen and sulfur.
How simple membrane peptides performed such essential proto-cellular functions as transport of ions and organic matter across membranes separating the interior of the cell from the environment, capture and utilization of energy, and transduction of environmental signals, is a key question in protobiological evolution. On the basis of detailed, molecular-level computer simulations we investigate how these peptides insert into membranes, self-assemble into higher-order structures and acquire functions. We have studied the insertion of an a-helical peptide containing leucine (L) and serine (S) of the form (LSLLLSL)S into a model membrane. The transmembrane state is metastable, and approximately 15 kcal/mol is required to insert the peptide into the membrane. Investigations of dimers formed by (LSLLLSL)S and glycophorin A demonstrate how the favorable free energy of helix association can offset the unfavorable free energy of insertion, leading to self- assembly of peptide helices in the membrane. An example of a self-assembled structure is the tetrameric transmembrane pore of the influenza virus M2 protein, which is an efficient and selective voltage-gated proton channel. Our simulations explain the gating mechanism and provide guidelines how to reengineering the channel to act as a simple proton pump. In general, emergence of integral membrane proteins appears to be quite feasible and may be easier to envision than the emergence of water-soluble proteins.