Search NASASearch

SEARCH · Search NASA

Results for “Analytic linearization”

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 37 records · Page 2

Design sensitivity analysis of nonlinear structural response

A unified theory is described of design sensitivity analysis of linear and nonlinear structures for shape, nonshape and material selection problems. The concepts of reference volume and adjoint structure are used to develop the unified viewpoint. A general formula for design sensitivity analysis is derived. Simple analytical linear and nonlinear examples are used to interpret various terms of the formula and demonstrate its use.

Cardoso, J. B.

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang

Impact of Cross-Axis Structural Dynamics on Validation of Linear Models for Space Launch System

A feasibility study was performed to examine the advisability of incorporating a set of Programmed Test Inputs (PTIs) during the Space Launch System (SLS) vehicle flight. The intent of these inputs is to provide validation to the preflight models for control system stability margins, aerodynamics, and structural dynamics. During October 2009, Ares I-X program was successful in carrying out a series of PTI maneuvers which provided a significant amount of valuable data for post-flight analysis. The resulting data comparisons showed excellent agreement with the preflight linear models across the frequency spectrum of interest. However unlike Ares I-X, the structural dynamics associated with the SLS boost phase configuration are far more complex and highly coupled in all three axes. This presents a challenge when implementing this similar system identification technique to SLS. Preliminary simulation results show noticeable mismatches between PTI validation and analytical linear models in the frequency range of the structural dynamics. An alternate approach was examined which demonstrates the potential for better overall characterization of the system frequency response as well as robustness of the control design.

Pei, Jing

Development of a computational testbed for numerical simulation of combustion instability

A synergistic hierarchy of analytical and computational fluid dynamic techniques is used to analyze three-dimensional combustion instabilities in liquid rocket engines. A mixed finite difference/spectral procedure is employed to study the effects of a distributed vaporization zone on standing and spinning instability modes within the chamber. Droplet atomization and vaporization are treated by a variety of classical models found in the literature. A multi-zone, linearized analytical solution is used to validate the accuracy of the numerical simulations at small amplitudes for a distributed vaporization region. This comparison indicates excellent amplitude and phase agreement under both stable and unstable operating conditions when amplitudes are small and proper grid resolution is used. As amplitudes get larger, expected nonlinearities are observed. The effect of liquid droplet temperature fluctuations was found to be of critical importance in driving the instabilities of the combustion chamber.

Grenda, Jeffrey

Coronal heating by the resonant absorption of Alfven waves: Wavenumber scaling laws.

The importance of global modes in coronal loop heating is well established. In the present work the scaling of the global-mode resonant heating rate with the perturbation wavenumbers is studied with the numerical solution of the linearized time-dependent MHD equations for a full compressible, low-beta, resistive plasma using an implicit integration scheme. The numerical simulations demonstrate that the dissipation on inhomogeneties in the background Alfven speed occurs in narrow resonant layer with the highest heating rate at the global-mode frequency. The global-mode heating rate H (sub r) was found to scale as H (Sub r) approximately k (sub y) (exp 1.03) when k (sub z) = 0.1, and as H (sub r) approximately k (sub y) (exp -1.93) when k (sub z) = 0.75, where k (sub y) and k (sub z) are the wavenumbers in the perpendicular and parallel to the magnetic field directions, respectively, while the dependence of H (sub r) on k (sub z) is more complex. The quality factor Q of the MHD resonance cavity scales as Q approximately k (sub y) (exp -1.8) for k (sub z) = 0.75 and as Q approximately k (sub y) (exp -1.46) for k (sub z) = 0.1. The numerically determined heating rate scaling, the global-mode fequency, and the quality factor are in good agreement with the analytical linear theory. The magnitude of the perturbed velocities was found to decrease with k (sub y). Assuming typical coronal loop parameters (B (sub 0) = 100-200 G, upsilon (sub A) = 2000-4000 km/s), the Alfven waves can supply the required heating to a low-Q loops.

Ofman, L.

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics

Numerical solution of laminar jet mixing with and without free stream.

Systematic numerical solutions of two-dimensional and axisymmetrical laminar jet of an incompressible fluid with and without free stream have been obtained. For the case without free stream, the exact numerical solutions have been obtained for different initial velocity profiles and compared with experimental results and similarity solutions. For the case with free stream, the numerical solutions have been compared with the linearized analytical solutions. An approximate numerical solution for the axial velocity distribution of the three-dimensional laminar jet is proposed. The accuracy of the method has been determined from the exact solutions of the two limiting cases of three-dimensional jets, i.e., the two-dimensional and axisymmetric cases.

Pai, S. I.

Theoretical performance of cross-wind axis turbines with results for a catenary vertical axis configuration

A general analysis capable of predicting performance characteristics of cross-wind axis turbines was developed, including the effects of airfoil geometry, support struts, blade aspect ratio, windmill solidity, blade interference and curved flow. The results were compared with available wind tunnel results for a catenary blade shape. A theoretical performance curve for an aerodynamically efficient straight blade configuration was also presented. In addition, a linearized analytical solution applicable for straight configurations was developed. A listing of the computer program developed for numerical solutions of the general performance equations is included in the appendix.

Muraca, R. J.

Analysis and testing of large space structures

Considerations and approaches to the ground testing of large space structures are discussed. The large size combined with the loading due to gravity makes testing of the complete structure difficult. Gravitational stiffening, suspension effects, virtual air mass, preloads, and air damping alter the dynamic characteristics. Low resonant frequencies and high modal densities within the frequency range of interest combine with small motions and accelerations to make testing difficult. Mechanism complexities and nonlinearities associated with space-erected/assembled structures cause structural complexity regardless of other considerations. Ground test approaches include scale models, element and substructure tests, and structural linearization. Analytical approaches are also discussed.

Stahle, C. V.

Application of unsteady laminar triple-deck theory to viscous-inviscid interactions from an oscillating flap in supersonic and subsonic flow

Unsteady triple-deck theory is applied to analyze the local viscous-inviscid interaction of an idealized oscillating flap with a laminar boundary layer in either supersonic or subsonic external flow. For small flap amplitudes and small-to-moderate nondimensional frequencies, linearized analytical solutions by means of Fourier transformation are given for the pressure and shear distributions ahead of and behind the flap hinge. In the supersonic case, the predicted unsteady viscous effects reduce the pressure amplitude and spread it out upstream while causing it to lag the quasi-steady inviscid prediction (flap motion); this also results in an unsteady delay of incipient separation at the hinge to a higher flap amplitude. In contrast, the lowered unsteady viscous pressure distribution in purely subsonic flow slightly leads the flap motion, resulting in incipient separation at a smaller flap amplitude than the quasi-steady value. These trends are shown to be in qualitative agreement with available experimental data.

Inger, G. R.

LEO high voltage solar array arcing response model

A series of mathematical models were developed that describe the electrical behavior of a large solar cell array floating electrically in the low Earth orbit (LEO) space plasma and struck by an arc at a point of negative bias. There are now three models in this series: ARCII, which is a fully analytical, linearized model; ARCIII, which is an extension of ARCIII that includes solar cell inductance as well as load reactance; Nonlinear ARC, which is a numerical model able to treat effects such as non-linearized, i.e., logarithmic solar cell I/V characteristics, conductance switching as a solar cell crosses plasma ground on a voltage excursion and non-ohmic plasma leakage current collection.

Metz, Roger N.

On the origin of the tropical intraseasonal oscillation

This study provides an explanation for the origin of the tropical intraseasonal (40-50 day) oscillation (TIO) based on a simple generalization of Gill's linear analytic model for tropical large-scale heat-induced circulation. The solution contains a convective region that excites an eastward-moving Kelvin wave and a westward-moving Rossby wave. The entire system moves eastward as a response to the circulation it excites at a speed at which the latent heat energy in the tropics is best extracted. The TIO speed is a weighted mean of the speed of the Kelvin wave and that of the Rossby wave.

Chao, Winston C.

A generic multi-flex-body dynamics, controls simulation tool for space station

An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.

London, Ken W.

A nested grid, nonhydrostatic, elastic model using a terrain-following coordinate transformation - The radiative-nesting boundary conditions

This study presents a nested-grid nonhydrostatic and elastic model using a terrain-following coordinate transformation as well as a unique application of grid-nesting techniques to the time-splitting elastic model. A simulation of the 10-m-high Witch of Agnesi Mountain provides the control to test this new model. The results show that the model produces the same solution as that derived from a simple linear analytic model. It is demonstrated that the new nested-grid model improves model resolution without resorting to the costly method of placing a fine-resolution grid over the entire domain. Since the wave reflection from the boundaries of the fine-grid model is well controlled, the boundary of the nested fine-grid model can be located even at the wave-active region. The model can be used to simulate various weather systems in which scale interactions are important.

Chen, Chaing

Effects of Meteoroid Erosion in Planetary Rings

This grant supported continuing studies of the effects of ballistic transport on the evolution of Saturn's rings. Ballistic transport, as used in this context, refers to the net transport of mass and angular momentum caused by the exchange of meteoroid impact ejecta between neighboring ring regions (Ip 1983, 1984, Morfill et al. 1983, Lissauer 1984, Durisen 1984a,b). The characteristic time scale associated with this process is the gross erosion time t(sub g) the time it would take a ring region to be completed eroded if all impact ejecta were lost. This time scale is estimated to be about 10(exp 5) to 10(exp 6) years for a ring region with normal optical depth tau approximately 1. Earlier work by myself and collaborators developed the physical theory and simulation techniques to model this process (Durisen et al. 1989, Cuzzi and Durisen 1990). Detailed simulations, supported in part by this grant, have demonstrated that ballistic transport can produce observed structures in Saturn's rings, especially at and near the inner edges of the A and B Rings (Durisen et al. 1992, 1996). The structures of interest in the real rings are illustrated in Figures 1 and 2. Most of these structures were previously unexplained. The computational results plus analytic treatments place useful constraints on fundamental ring properties, including the indication of a relatively young ring age less than or equal 10(exp 8) years (see reviews by Nicholson and Dones 1991, Esposito 1993, Cuzzi 1995, and Porco 1995). This grant also supported development of the faster computational algorithms necessary to permit longer evolutions. Resulting simplications in the ballistic transport equations permitted an analytic linear stability analysis (Durisen 1995), which has provided considerable insight into ballistic transport processes and applications to Saturn's rings. All these accomplishments are described in more detail below.

Durisen, Richard H.

Latitudinal Dependence of the Radial IMF Component - Interplanetary Imprint

Ulysses measurements have confirmed that there is no significant gradient with respect to heliomagnetic latitude in the radial component, B(sub r,), of the interplanetary magnetic field. There are two processes responsible for this observation. In the corona, the plasma beta is much less than 1, except directly above streamers, so both longitudinal and latitudinal (meridional) gradients in field strength will relax, due to the transverse magnetic pressure gradient force, as the solar wind carries magnetic flux away from the Sun. This happens so quickly that the field is essentially uniform by 5 solar radius. Beyond 10 solar radius, beta is greater than 1 and it is possible for a meridional thermal pressure gradient to redistribute magnetic flux - an effect apparently absent in Ulysses and earlier ICE and Interplanetary Magnetic Physics (IMP) data. We discuss this second effect here, showing that its absence is mainly due to the perpendicular part of the anisotropic thermal pressure gradient in the interplanetary medium being too small to drive significant meridional transport between the Sun and approx. 4 AU. This is done using a linear analytic estimate of meridional transport. The first effect was discussed in an earlier paper.

Suess, S. T.

Cellular Spacing Selection During the Directional Solidification of Binary Alloys. A Numerical Approach

The evolution of cellular solid/liquid interfaces from an initially unstable planar front was studied by means of a two-dimensional computer simulation. The developed numerical model makes use of an interface tracking procedure and has the capability to describe the dynamics of the interface morphology based on local changes of the thermodynamic conditions. The fundamental physics of this formulation was validated against experimental microgravity results and the predictions of the analytical linear stability theory. The performed simulations revealed that in certain conditions, based on a competitive growth mechanism, an interface could become unstable to random perturbations of infinitesimal amplitude even at wavelengths smaller than the neutral wavelength, lambda(sub c), predicted by the linear stability theory. Furthermore, two main stages of spacing selection have been identified. In the first stage, at low perturbations amplitude, the selection mechanism is driven by the maximum growth rate of instabilities while in the second stage the selection is influenced by nonlinear phenomena caused by the interactions between the neighboring cells. Comparison of these predictions with other existing theories of pattern formation and experimental results will be discussed.

Catalina, Adrian V.