Search NASA⌕ Search

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 1,225 records · Page 68

ALFVEN WAVES IN SHEAR FLOWS REVISITED

We revisit our earlier study of the evolution of an initial propagating Alfven wave in a magnetic-field-aligned flow with a cross-field velocity shear. Our goal is to show how the Alfven wave drives up plasma density fluctuations which might be observed and serve as a signature of the presence of Alfven waves in regions such as the solar corona which are inaccessible to direct observations. Here, we introduce a new initial condition which takes into account the initial distortion of the streamlines by the Alfven wave, and we present new analytical results for the driven waves. We find that the density fluctuations of a properly placed linearly polarized Alfven wave in a shear flow are much smaller than we originally estimated.

79 ASTRONOMY AND ASTROPHYSICS↗

AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.

79 ASTRONOMY AND ASTROPHYSICS↗

Influences of δB contribution and parallel inertial term of energetic particles on MHD-kinetic hybrid simulations: a case study of the 1/1 internal kink mode

The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model (Park et al 1992 Phys. Fluids B 4 2033–7) has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling scheme is adopted in this model. However, two noteworthy issues arise in the model application: firstly, the coupled term introduced in the pressure-coupling scheme, (∇•P h ) ⟂ , is often simplified by ∇•P h , which is equivalent to neglecting the parallel inertial term of EPs; secondly, besides the $δf$ contribution caused by changing in the EP distribution function, the magnetic field perturbation (the $δB$ contribution) generated during development of the instabilities should also be considered, but it is often ignored in existing hybrid simulations. In this paper, we derive the analytical formulations under these two coupling schemes and then numerically study the representative case of the linear stability of the $m/n$ = $1/1$ internal kink mode (IKM) (Fu et al 2006 Phys. Plasmas 13 052517) by using the CLT-K code. Further, it is found that the approximated models can still yield reasonable results when EPs are isotopically distributed. But it fails completely in cases with anisotropic EP distributions. In addition, we further investigate the influence of EP's orbit width on the stability of IKM and verify the equivalence between pressure-coupling scheme and the current-coupling scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CMB and BAO constraints for an induced gravity dark energy model with a quartic potential

We study the predictions for structure formation in an induced gravity dark energy model with a quartic potential. By developing a dedicated Einstein-Boltzmann code, we study self-consistently the dynamics of homogeneous cosmology and of linear perturbations without using any parametrization. By evolving linear perturbations with initial conditions in the radiation era, we accurately recover the quasi-static analytic approximation in the matter dominated era. We use PLANCK 2013 data and a compilation of baryonic acoustic oscillation (BAO) data to constrain the coupling γ to the Ricci curvature and the other cosmological parameters. By connecting the gravitational constant in the Einstein equation to the one measured in a Cavendish-like experiment, we find γ < 0.0012 at 95% CL with PLANCK 2013 and BAO data. This is the tightest cosmological constraint on γ and on the corresponding derived post-Newtonian parameters. Because of a degeneracy between γ and the Hubble constant H{sub 0}, we show how larger values for γ are allowed, but not preferred at a significant statistical level, when local measurements of H{sub 0} are combined in the analysis with PLANCK 2013 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Hall viscosity and nonlocal conductivity of gapped graphene

We calculate the Hall viscosity and the nonlocal (i.e., dependent on wave vector q) Hall conductivity of “gapped graphene” (a nontopological insulator with two valleys) in the presence of a strong perpendicular magnetic field. Using the linear-response theory at zero temperature within the Dirac approximation for the Landau levels, we present analytical expressions for both valley and total Hall viscosities and conductivity up to q 2 at all frequencies. Although the final formulas for total Hall viscosity and conductivity are similar to the ones previously obtained for gapless graphene, the derivation reveals a significant difference between the two systems. First of all, both the Hall viscosity and the Hall conductivity vanish when the Fermi level lies in the gap that separates the lowest Landau level in the conduction band from the highest Landau level in the valence band. It is only when the Fermi level is not in the gap that the familiar formulas of gapless graphene are recovered. Second, in the case of gapped graphene, it is not possible (at least, within our present approach) to define a single-valley Hall viscosity: this quantity diverges with a strength proportional to the magnitude of the gap. It is only when both valleys are included that the diverging terms, having opposite signs in the two valleys, cancel out, and the familiar result is recovered. In contrast to this, the nonlocal Hall conductivity is finite in each valley. These results indicate that the Hoyos-Son formula [Phys. Rev. Lett. 108, 066805 (2012)] connecting the Hall viscosity to the coefficient of q 2 in the small-q expansion of the q-dependent Hall conductivity cannot be applied to each valley, but only to the system as a whole. As a result, the problem of defining a “valley Hall viscosity” remains open.

36 MATERIALS SCIENCE↗

Role of Sink Density in Nonequilibrium Chemical Redistribution in Alloys

Nonequilibrium chemical redistribution in open systems submitted to external forces, such as particle irradiation, leads to changes in the structural properties of the material, potentially driving the system to failure. Such redistribution is controlled by the complex interplay between the production of point defects, atomic transport rates, and the sink character of the microstructure. In this work, we analyze this interplay by means of a kinetic Monte Carlo (KMC) framework with an underlying atomistic model for the Fe-Cr model alloy to study the effect of ideal defect sinks on Cr concentration profiles, with a particular focus on the role of interface density. We observe that the amount of segregation decreases linearly with decreasing interface spacing. Within the framework of the thermodynamics of irreversible processes, a general analytical model is derived and assessed against the KMC simulations to elucidate the structure-property relationship of this system. Interestingly, in the kinetic regime where elimination of point defects at sinks is dominant over bulk recombination, the solute segregation does not directly depend on the dose rate but only on the density of sinks. Furthermore, this model provides new insight into the design of microstructures that mitigate chemical redistribution and improve radiation tolerance.

36 MATERIALS SCIENCE↗

Recovery of the characteristics of plasma turbulence from the radial correlation backscattering diagnostics

Signals of the backscattering radial correlation Doppler diagnostics of plasma density fluctuations in the presence of the cutoff of the probing wave are analyzed theoretically with allowance for the curvature of magnetic surfaces. The scattering of the probing electromagnetic wave is considered in the linear (Born) approximation with respect to the amplitude of fluctuations. Using the Wentzel−Kramers−Brillouin approach, analytical expressions for the scattered signal and the correlation function of two scattered signals corresponding to oblique probing at different frequencies are derived. A criterion is obtained for the tilt angle of the antenna pattern at which the two-point turbulence correlation function can be measured directly. A method is proposed to recover the spectrum of plasma density fluctuations from the data on the radial wavenumbers even if this criterion is violated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SWinzip v. 2.0

SWinzip is a Matlab and C++ library for scientific lossy data compression and reconstruction using compressed sensing and tree-wavelets transforms. These methods are known for their large compression and usefulness in data analytics such as features extraction. Compressed sensing and wavelets methods rely heavily on sparse and dense linear algebra operations implemented through the Boost codes in our library. SWinzip accommodates data represented on both regular grids (e.g. image data) and point-clouds (e.g. unstructured meshes). In SWinzip v2.0, there is more focus on wavelet compression and the following features are added: - The matrix-free (direct) wavelet transform is implemented in both C++ and Matlab. A threshold-based compression can now be used. Examples for distributed compression are also added. - Python implementation of Alpert wavelet transforms - Field data comparison using wavelets in Matlab with global and local metrics

Solloum, Maher↗

Lunar electromagnetic scattering. I - Propagation parallel to the diamagnetic cavity axis.

A general analytic solution is obtained for the interaction of the moon and its downstream cavity with a linearly polarized plane electromagnetic wave propagating parallel to the cavity axis. The solution is formulated in terms of a spherical moon model with arbitrary radially dependent electromagnetic parameters and a nonconducting cylindrical downstream cavity. Use is made of a number of approximations that are consistent with the physical nature of the interaction between the moon and the solar wind.

Schwartz, K.↗

High supersonic aerodynamic characteristics of five irregular planform wings with systematically varying wing fillet geometry tested in the NASA/LaRC 4-foot UPWT (LEG 2) (LA45A/B)

An experimental and analytical aerodynamic program to develop predesign guides for irregular planform wings is reported. The benefits are linearization of subsonic lift curve slope to high angles of attack and avoidance of subsonic pitch instabilities at high lift by proper tailoring of the planform fillet wing combination while providing the desired hypersonic trim angle and stability. The two prime areas of concern are to optimize shuttle orbiter landing and entry characteristics. Basic longitudinal aerodynamic characteristics at high supersonic speeds are developed.

Source record↗

Numerical solutions and laser-Doppler measurements of spin-up

The spin-up flow in a cylinder of homogeneous fluid has been examined both experimentally and numerically. A series of laser-Doppler measurements was made of the zonal flow over a range of Ekman numbers and Rossby numbers at various locations in the interior of the flow. These measurements exceed previous ones in accuracy. The weak inertial modes excited by the impulsive start are detectable. The numerical simulations used the primitive equations in axisymmetric form and employed finite-difference techniques on both constant and variable grids. The number of grid points necessary to resolve the Ekman layers was determined. A thorough comparison of the simulations and the experimental measurements is made which includes the details of the amplitude and frequency of the inertial modes. Agreement to within the experimental tolerance is achieved. Analytical results for conditions identical to those in the experiments are not available but some similar linear and nonlinear theories are also compared with the experiments.

Warn-Varnas, A.↗

Detailed gravity anomalies from GEOS-3 satellite altimetry data

A technique for deriving mean gravity anomalies from dense altimetry data was developed. A combination of both deterministic and statistical techniques was used. The basic mathematical model was based on the Stokes' equation which describes the analytical relationship between mean gravity anomalies and geoid undulations at a point; this undulation is a linear function of the altimetry data at that point. The overdetermined problem resulting from the excessive altimetry data available was solved using Least-Squares principles. These principles enable the simultaneous estimation of the associated standard deviations reflecting the internal consistency based on the accuracy estimates provided for the altimetry data as well as for the terrestrial anomaly data. Several test computations were made of the anomalies and their accuracy estimates using GOES-3 data.

Gopalapillai, G. S.↗

Design curves for non-linear analysis of simply-supported, uniformly-loaded rectangular plates

Design curves for the non-linear analysis of simply-supported rectangular plates subjected to uniform normal pressure loads have been developed. These curves yield the center deflection, center stress and corner stress in non-dimensionalized form plotted against a dimensionless parameter describing the load intensity. The results presented are based on extensive non-linear finite element analysis employing the ARGUS structural analysis program. Plates with length to width ratios of 1, 1.5, 2, 3 and 4 are included. The load range considered extends to 1000 times the load at which the behavior of the plate becomes significantly non-linear. Over the load range considered, the analysis shows that the ratio of center deflection to plate thickness for a square plate is less than 16 to 1, whereas linear theory would predict a center deflection 400 times the plate thickness. Likewise, the stress is markedly lower than would be predicted by linear theory. The present results are shown to be in excellent agreement with the classical linear theory up to a central deflection to plate thickness ratio of about one-half. In the non-linear regime the present results for deflection and stress are in very good agreement with the analytical and experimental work of other investigators.

Moore, D.↗

Effect of streamwise vortices on Tollmien-Schlichting waves

The method of multiple scales is used to determine a first-order uniform expansion for the effect of counter-rotating steady streamwise vortices in growing boundary layers on oblique Tollmien-Schlichting waves. The results show that such vortices have a strong tendency to amplify oblique Tollmien-Schlichting waves having a spanwise wavelength that is twice the wavelength of the vortices. An analytical expression is derived for the growth rates of these waves. These exponential growth rates increase linearly with increasing amplitudes of the vortices. Numerical results are presented. They suggest that this mechanism may dominate the instability.

Nayfeh, A. H.↗

On the design derivatives of eigenvalues and eigenvectors for distributed parameter systems

In this paper, analytic expressions are obtained for the design derivatives of eigenvalues and eigenfunctions of self-adjoint linear distributed parameter systems. Explicit treatment of boundary conditions is avoided by casting the eigenvalue equation into integral form. Results are expressed in terms of the linear operators defining the eigenvalue problem, and are therefore quite general. Sufficiency conditions appropriate to structural optimization of eigenvalues are obtained.

Reiss, R.↗

Identifying approximate linear models for simple nonlinear systems

This paper addresses the identification (realization) of approximate linear models from response data for certain nonlinear dynamic systems. Response characteristics for several typical nonlinear joints are analyzed mathematically and represented by series expansions. The parameters of the series expansion are then compared with the modal parameters of a linear model identified by the Eigensystem Realization Algorithm. The agreement of the identified model and the analytically derived representation is excellent for the cases studied. Also laboratory data from a model which exhibited stiffening behavior was analyzed using the Eigensystem Realization algorithm and Fast Fourier Transform. The laboratory experiment demonstrated the ability of the technique to recover the model characteristics using real data.

Horta, L. G.↗

Recursive mass matrix factorization and inversion: An operator approach to open- and closed-chain multibody dynamics

This report advances a linear operator approach for analyzing the dynamics of systems of joint-connected rigid bodies.It is established that the mass matrix M for such a system can be factored as M=(I+H phi L)D(I+H phi L) sup T. This yields an immediate inversion M sup -1=(I-H psi L) sup T D sup -1 (I-H psi L), where H and phi are given by known link geometric parameters, and L, psi and D are obtained recursively by a spatial discrete-step Kalman filter and by the corresponding Riccati equation associated with this filter. The factors (I+H phi L) and (I-H psi L) are lower triangular matrices which are inverses of each other, and D is a diagonal matrix. This factorization and inversion of the mass matrix leads to recursive algortihms for forward dynamics based on spatially recursive filtering and smoothing. The primary motivation for advancing the operator approach is to provide a better means to formulate, analyze and understand spatial recursions in multibody dynamics. This is achieved because the linear operator notation allows manipulation of the equations of motion using a very high-level analytical framework (a spatial operator algebra) that is easy to understand and use. Detailed lower-level recursive algorithms can readily be obtained for inspection from the expressions involving spatial operators. The report consists of two main sections. In Part 1, the problem of serial chain manipulators is analyzed and solved. Extensions to a closed-chain system formed by multiple manipulators moving a common task object are contained in Part 2. To retain ease of exposition in the report, only these two types of multibody systems are considered. However, the same methods can be easily applied to arbitrary multibody systems formed by a collection of joint-connected regid bodies.

Rodriguez, G.↗

Theoretical analysis and meteorological interpretation of the role of raindrop shape on microwave attenuation and propagation phase shifts - Implications for the radar measurement of rain

The effects of raindrop shape on the attenuation and propagation phase shift at vertical and horizontal linear polarizations is examined, focusing on the implications of these effects on rain measurements by radars. Analytic expressions for the imaginary and real components of the forward scattering amplitude function are dervied for a wavelength of 2.2 cm. It is shown that the summation attenuation rate and the rate of summation propagation phase shift are both insensitive to the shapes of raindrops.

Jameson, A. R.↗