Search NASA⌕ Search

SEARCH · Search NASA

Results for “Jump conditions”

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 199 records · Page 11

A vector-continuous loading concept for aerodynamic panel methods

An approach to the reduction of discretization errors in aerodynamic panel methods is presented. The approach is based on preventing the occurence of induced velocity singularities at panel slope discontinuities by maintaining continuity of the velocity jump vector across the panels. The approach was implemented in a two-dimensional incompressible panel method formulation and evaluated by application to several external and internal flow problems. The method is shown to exhibit a second order accuracy trend and to produce smaller errors with velocity component boundary conditions imposed on the real flow than with equipotential boundary conditions imposed on the imaginary flow behind the panels. For flows around airfoil sections with either sharp or blunt trailing edges, the method gives excellent agreement with results from a well developed finite difference method. The method is well behaved and is insensitive to irregularities in panel size distribution.

Kemp, W. B., Jr.↗

Multivalued velocity field model of turbulence

The considered study is based on new theoretical concepts regarding a post-instability model of a fluid discussed by Zak (1980). The model permits the completion of the governing equations of turbulence by introducing multivalued fields of velocities. Attention is given to the mechanism of energy dissipation, the characteristic wave propagation, a simplified model, the formation of turbulence around stagnation points, the formulation of boundary conditions, and the mechanism of turbulence formation. The mechanism of turbulence formation can be understood as propagation of initial discontinuities from the boundaries into a flow with the characteristic velocity which is defined by the normal (to the boundary) velocity components. These components emerge at the boundary as a result of jumps in the tangential components due to the continuity equation.

Zak, M.↗

Oxidation Behavior of GRCop-84 (Cu-8Cr-4Nb) at Intermediate and High Temperatures

The oxidation behavior of GRCop-84 (Cu-8 at %Cr-4 at %Nb) has been investigated in air and in oxygen, for durations of 0.5 to 50 hours and temperatures ranging from 500 to 900 C. For comparison, data was also obtained for the oxidation of Cu and NARloy-Z (Cu-3 wt% Ag-0.5 wt% Zr) under the same conditions. Arrhenius plots of those data showed that all three materials had similar oxidation rates at high temperatures (> 750 C). However, at intermediate temperatures (500 to 750 C) GRCop exhibited significantly higher oxidation resistance than Cu and NARloy-Z. The oxidation kinetics of GRCop-84 exhibited a sharp and discontinuous jump between the two regimes. Also, in the high temperature regime GRCop-84 oxidation rate was found to change from a high initial value to a significantly smaller terminal value at each temperature, with progress of oxidation; the two different oxidation rates were found to correlate with a porous intial oxide and a dense final oxide, respectively.

Thomas-Ogbuji, Linus U.↗

Secular rotational motions and the mechanical structure of a dynamical viscoelastic earth

A survey is presented of analytical methods for computing the linear responses of the rotational axis of a layered viscoelastic earth to surface loading. Theoretical research in this area is first summarized, and the differences between the mechanical boundary conditions to be applied at the interface separating the upper and lower mantles for an adiabatically and chemically stratified mantle are discussed. Some examples of polar wander and secular variation of the spin rate from glacial excitation are presented for various types of chemical and viscosity stratifications. The effects of an artificial density jump at the base of the lithosphere in models are examined, and certain issues concerning the fluid tidal Love number for different types of density stratification are addressed. The meaning of effective plate thickness over geological time scales for rotational dynamics is discussed.

Yuen, D. A.↗

Direct simulation of hypersonic flows over blunt slender bodies

Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.

Moss, J. N.↗

First-Order System Least-Squares for Second-Order Elliptic Problems with Discontinuous Coefficients

The first-order system least-squares methodology represents an alternative to standard mixed finite element methods. Among its advantages is the fact that the finite element spaces approximating the pressure and flux variables are not restricted by the inf-sup condition and that the least-squares functional itself serves as an appropriate error measure. This paper studies the first-order system least-squares approach for scalar second-order elliptic boundary value problems with discontinuous coefficients. Ellipticity of an appropriately scaled least-squares bilinear form of the size of the jumps in the coefficients leading to adequate finite element approximation results. The occurrence of singularities at interface corners and cross-points is discussed. and a weighted least-squares functional is introduced to handle such cases. Numerical experiments are presented for two test problems to illustrate the performance of this approach.

Manteuffel, Thomas A.↗

Intermittent Flow Regimes in a Transonic Fan Airfoil Cascade

A study was conducted in the NASA Glenn Research Center linear cascade on the intermittent flow on the suction surface of an airfoil section from the tip region of a modern low aspect ratio fan blade. Experimental results revealed that, at a large incidence angle, a range of transonic inlet Mach numbers exist where the leading-edge shock-wave pattern was unstable. Flush mounted high frequency response pressure transducers indicated large local jumps in the pressure in the leading edge area, which generates large intermittent loading on the blade leading edge. These measurements suggest that for an inlet Mach number between 0.9 and 1.0 the flow is bi-stable, randomly switching between subsonic and supersonic flows. Hence, it appears that the change in overall flow conditions in the transonic region is based on the frequency of switching between two stable flow states rather than on the continuous increase of the flow velocity. To date, this flow behavior has only been observed in a linear transonic cascade. Further research is necessary to confirm this phenomenon occurs in actual transonic fans and is not the byproduct of an endwall restricted linear cascade.

Lepicovsky, J.↗

Steady-state dynamic behavior of an auxiliary bearing supported rotor system

This paper investigates the steady-state responses of a rotor system supported by auxiliary bearings in which there is a clearance between the rotor and the inner race of the bearing. A simulation model based upon the rotor of a production jet engine is developed and its steady-state behavior is explored over a wide range of operating conditions for various parametric configurations. Specifically, the influence of rotor imbalance, support stiffness, and damping is studied. It is found that imbalance may change the rotor responses dramatically in terms of frequency contents at certain operating speeds. Subharmonic responses of 2nd order through 10th order are all observed except the 9th order. Chaotic phenomenon is also observed. Jump phenomena (or double-valued responses) of both hard-spring type and soft-spring type are shown to occur at low operating speeds for systems with low auxiliary bearing damping or large clearance even with relatively small imbalance. The effect of friction between the shaft and the inner race of the bearing is also discussed.

Xie, Huajun↗

Jump resonant frequency islands in nonlinear feedback control systems

A new type of jump resonance is predicted and observed in certain nonlinear feedback control systems. The new jump resonance characteristic is described as a 'frequency island' due to the fact that a portion of the input-output transfer characteristic is disjoint from the main body. The presence of such frequency islands was predicted by using a sinusoidal describing function characterization of the dynamics of an inertial gyro employing nonlinear ternary rebalance logic. While the general conditions under which such islands are possible has not been examined, a numerical approach is presented which can aid in establishing their presence. The existence of the frequency islands predicted for the ternary rebalanced gyro was confirmed by simulating the nonlinear system and measuring the transfer function.

Koenigsberg, W. D.↗

The numerical simulation of steady transonic rotational flow using a dual potential formulation

A finite-difference method is presented that simulates steady transonic rotational flow of an inviscid fluid by representing the velocity field as the sum of scalar and vector potentials. This dual potential velocity decomposition extends the validity of the scalar (full) velocity potential to include vorticity. The inclusion of a vector potential also permits an alternate treatment of lift that does not require a circulation wake cut. This is accomplished by specifying the vector potential as a constant on the airfoil surface in order to satisfy a Kutta condition. The governing equations are solved as iteratively decoupled scalar equations using approximate factorization techniques, and the overall efficiency approaches that of the full potential equation. The governing equations are able to convect entropy and vorticity throughout the flow field and are equivalent to the Euler equations in continuous flow domains, however at shocks the Rankine-Hugoniot entropy jump must be supplied. An entropy correction method is presented and verified with transonic airfoil solutions of the Euler equations.

Chaderjian, N. M.↗

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗

The Dryline on 22 May 2002 during IHOP_2002: Convective-Scale Measurements at the Profiling Site

A detailed analysis of the structure of a double dryline observed over the Oklahoma panhandle during the first International H2O Project (IHOP_2002) convective initiation (CI) mission on 22 May 2002 is presented. A unique and unprecedented set of high temporal and spatial resolution measurements of water vapor mixing ratio, wind, and boundary layer structure parameters were acquired using the National Aeronautics and Space Administration (NASA) scanning Raman lidar (SRL), the Goddard Lidar Observatory for Winds (GLOW), and the Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE), respectively. These measurements are combined with the vertical velocity measurements derived from the National Center for Atmospheric Research (NCAR) Multiple Antenna Profiler Radar (MAPR) and radar structure function from the high-resolution University of Massachusetts frequency-modulated continuous-wave (FMCW) radar to reveal the evolution and structure of the late afternoon double-dryline boundary layer. The eastern dryline advanced and then retreated over the Homestead profiling site in the Oklahoma panhandle, providing conditions ripe for a detailed observation of the small-scale variability within the boundary layer and the dryline. In situ aircraft data, dropsonde and radiosonde data, along with NCAR S-band dual-polarization Doppler radar (S-Pol) measurements, are also used to provide the larger-scale picture of the double-dryline environment. Moisture and temperature jumps of about 3 g kg(sup -1) and 1 -2 K, respectively, were observed across the eastern radar fine line (dryline), more than the moisture jumps (1-2 g kg(sup -1)) observed across the western radar fine line (secondary dryline). Most updraft plumes observed were located on the moist side of the eastern dryline with vertical velocities exceeding 3 m s(sup -1) and variable horizontal widths of 2-5 km, although some were as wide as 7-8 km. These updrafts were up to 1.5 g kg(sup -1) moister than the surrounding environment. Although models suggested deep convection over the Oklahoma panhandle and several cloud lines were observed near the dryline, the dryline itself did not initiate any storms over the intensive observation region (IOR). Possible reasons for this lack of convection are discussed. Strong capping inversion and moisture detrainment between the lifting condensation level and the level of free convection related to an overriding drier air, together with the relatively small near-surface moisture values (less than 10 g kg(sup -1)), were detrimental to CI in this case.

Demoz, Belay↗

The origin of the Kirkwood gaps - A mapping for asteroidal motion near the 3/1 commensurability

A mapping of the phase space onto itself with the same low-order resonance structure as the 3/1 commensurability in the planar-elliptic restricted three-body problem is obtained. This mapping is about 1,000 times faster than the usual method of numerically integrating the averaged equations of motion. It exhibits some surprising behavior that might provide a key to the origin of the Kirkwood gaps. It is noted that a test asteroid placed in the gap may evolve for a million years with low eccentricity (less than 0.05) and then suddenly jump to large eccentricity (greater than 0.3), becoming a Mars crosser. The removal of the asteroid by a close encounter with Mars would then be possible. As a first test of this hypothesis, a distribution of 300 test asteroids in the area of the 3/1 commensurability was evolved for two million years. When the Mars crossers are removed, the distribution of initial conditions reveals a gap at the location of the 3/1 Kirkwood gap.

Wisdom, J.↗

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations That Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + definite integral of g(x, t)F(t,y(t))dt with limits between 0 and 1,0 less than or equal to x les than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integral equations arise, e.g., when one applied Green's function techniques to nonlinear two-point boundary value problems of the form y "(x) =f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and y(l) = y(sub l), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trepezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal rule, thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations that Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗

Towards a Theory of Sampled-Data Piecewise-Deterministic Markov Processes

The analysis and design of practical control systems requires that stochastic models be employed. Analysis and design tools have been developed, for example, for Markovian jump linear continuous and discrete-time systems, piecewise-deterministic processes (PDP's), and general stochastic hybrid systems (GSHS's). These model classes have been used in many applications, including fault tolerant control and networked control systems. This paper presents initial results on the analysis of a sampled-data PDP representation of a nonlinear sampled-data system with a jump linear controller. In particular, it is shown that the state of the sampled-data PDP satisfies the strong Markov property. In addition, a relation between the invariant measures of a sampled-data system driven by a stochastic process and its associated discrete-time representation are presented. As an application, when the plant is linear with no external input, a sufficient testable condition for the convergence in distribution to the invariant delta Dirac measure is given.

Herencia-Zapana, Heber↗

Multi-Pass Sequential Mini-Batch Stochastic Gradient Descent Algorithms for Noise Covariance Estimation in Adaptive Kalman Filtering

Estimation of unknown noise covariances in a Kalman filter is a problem of significant practical interest in a wide array of applications. Although this problem has a long history, reliable algorithms for their estimation were scant, and necessary and sufficient conditions for identifiability of the covariances were in dispute until recently. Necessary and sufficient conditions for covariance estimation and a batch estimation algorithm were presented in our previous study. This paper presents stochastic gradient descent algorithms for noise covariance estimation in adaptive Kalman filters that are an order of magnitude faster than the batch method for similar or better root mean square error. More significantly, these algorithms are applicable to non-stationary systems where the noise covariances can occasionally jump up or down by an unknown magnitude. The computational efficiency of the new algorithms stems from adaptive thresholds for convergence, recursive fading memory estimation of the sample cross-correlations of the innovations, and accelerated stochastic gradient descent algorithms. The comparative evaluation of the proposed methods on a number of test cases demonstrates their computational efficiency and accuracy.

Adaptive Kalman filtering↗

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗