Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational mode decomposition”

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.

32 records · Page 2

East Asian Summer Monsoon Rainfall: A Historical Perspective of the 1998 Flood over Yangtze River

One of the main factors that might have caused the disastrous flood in China during 1998 summer is long-term variations that include a trend indicating increasing monsoon rainfall over the Yangtze River Valley. China's 160-station monthly rainfall anomaly for the summers of 1955-98 is analyzed for exploring such long-term variations. Singular value decomposition (SVD) between the summer rainfall and the global sea surface temperature (SST) anomalies reveals that the rainfall over Yangtze River Valley is closely related to global and regional SST variabilities at both interannual and interdecadal timescales. SVD1 mode links the above normal rainfall condition in central China to an El Nino-like SSTA distribution, varying on interannual timescale modified by a trend during the period. SVD3 mode links positive rainfall anomaly in Yangtze River Valley to the warm SST anomaly in the subtropical western Pacific, varying on interannual timescales modified by interdecadal timescales. This link tends to be stronger when the Nino3 area becomes colder and the western subtropical Pacific becomes warmer. The 1998 summer is a transition season when the 1997/98 El Nino event was in its decaying phase, and the SST in the Nino3 area emerged below normal anomaly while the subtropical western Pacific SST above normal. Thus, the first and third SVD modes become dominant in 1998 summer, favoring more Asian summer monsoon rainfall over the Yangtze River Valley.

Weng, H.-Y.↗

On Holo-Hilbert Spectral Analysis: A Full Informational Spectral Representation for Nonlinear and Non-Stationary Data

The Holo-Hilbert spectral analysis (HHSA) method is introduced to cure the deficiencies of traditional spectral analysis and to give a full informational representation of nonlinear and non-stationary data. It uses a nested empirical mode decomposition and Hilbert-Huang transform (HHT) approach to identify intrinsic amplitude and frequency modulations often present in nonlinear systems. Comparisons are first made with traditional spectrum analysis, which usually achieved its results through convolutional integral transforms based on additive expansions of an a priori determined basis, mostly under linear and stationary assumptions. Thus, for non-stationary processes, the best one could do historically was to use the time- frequency representations, in which the amplitude (or energy density) variation is still represented in terms of time. For nonlinear processes, the data can have both amplitude and frequency modulations (intra-mode and inter-mode) generated by two different mechanisms: linear additive or nonlinear multiplicative processes. As all existing spectral analysis methods are based on additive expansions, either a priori or adaptive, none of them could possibly represent the multiplicative processes. While the earlier adaptive HHT spectral analysis approach could accommodate the intra-wave nonlinearity quite remarkably, it remained that any inter-wave nonlinear multiplicative mechanisms that include cross-scale coupling and phase-lock modulations were left untreated. To resolve the multiplicative processes issue, additional dimensions in the spectrum result are needed to account for the variations in both the amplitude and frequency modulations simultaneously. HHSA accommodates all the processes: additive and multiplicative, intra-mode and inter-mode, stationary and nonstationary, linear and nonlinear interactions. The Holo prefix in HHSA denotes a multiple dimensional representation with both additive and multiplicative capabilities.

Huang, Norden E.↗

Does The Nile Reflect Solar Variability?

Historical records of the Nile water level provide a unique opportunity to investigate the possibility that solar variability influences the Earth's climate. Particularly important are the annual records of the water level, which are uninterrupted for the years 622-1470 A.D. These records are non-stationary, so that standard spectral analyses cannot adequately characterize them. Here the Empirical Mode Decomposition technique, which is designed to deal with nonstationary, nonlinear time series, becomes useful. It allows the identification of two characteristic time scales in the water level data that can be linked to solar variability: the 88 year period and a time scale of about 200 years. These time scales are also present in the concurrent aurora data. Auroras are driven by coronal mass ejections and the rate of auroras is an excellent proxy for solar variability. Analysis of auroral data contemporaneous with the Nile data shows peaks at 88 years and about 200 years. This suggests a physical link between solar variability and the low frequency variations of the Nile water level. The link involves the influence of solar variability on the North Annual Mode of atmospheric variability and its North Atlantic and Indian Oceans patterns that affect rainfall over Eastern Equatorial Africa where the Nile originates.

statistical methods↗

Seasonal Cycles of Meridional Overturning and Heat Transport of the Indian Ocean

A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (x 10(exp 6) cubic m/s) for the overturning and 1.8 and 1.4 (x 10(exp 15) W) for heat transport near 10 S and 10 N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman How plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10 deg N and 10 deg S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution front vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15 S. Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean. indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively. Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads 10 unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.

Lee, Tong↗

Non-Linear Structural Dynamics Characterization using a Scanning Laser Vibrometer

This paper presents the use of a scanning laser vibrometer and a signal decomposition method to characterize non-linear dynamics of highly flexible structures. A Polytec PI PSV-200 scanning laser vibrometer is used to measure transverse velocities of points on a structure subjected to a harmonic excitation. Velocity profiles at different times are constructed using the measured velocities, and then each velocity profile is decomposed using the first four linear mode shapes and a least-squares curve-fitting method. From the variations of the obtained modal \ielocities with time we search for possible non-linear phenomena. A cantilevered titanium alloy beam subjected to harmonic base-excitations around the second. third, and fourth natural frequencies are examined in detail. Influences of the fixture mass. gravity. mass centers of mode shapes. and non-linearities are evaluated. Geometrically exact equations governing the planar, harmonic large-amplitude vibrations of beams are solved for operational deflection shapes using the multiple shooting method. Experimental results show the existence of 1:3 and 1:2:3 external and internal resonances. energy transfer from high-frequency modes to the first mode. and amplitude- and phase- modulation among several modes. Moreover, the existence of non-linear normal modes is found to be questionable.

Pai, P. F.↗

An efficient design sensitivity analysis of eigenvectors

Subspace iteration has been a major advance in solving large eigen problems when only a subset of eigen-pairs is required. The essence of this method is a transformation from displacement coordinates of an n-th order eigensystem to generalized coordinates of a smaller q-th order. The eigenvalue problem is then solved in the reduced space. The method was first developed by Clint and Jennings for real symmetric systems and was then called 'simultaneous iteration'. The success of the method prompted further research along this line and there have been many improved algorithms developed. This approach has been widely used by structural engineers for extracting the most useful natural frequencies and mode shapes of large-scale dynamic systems. This paper exploits into a new direction which is in the form of iterative process for simultaneously calculating eigenvector derivatives of many eigenvectors with respect to multi-variables. The method fully uses all the available information from preceding eigenvalue solution and, thus, effectively economizes computational efforts. It iterates through two equations derived from the first variation of the two fundamental equations used in subspace iteration method. There is no expensive large matrix decomposition required and the process converges to acceptable solution in a finite number of iterations. Therefore, the procedure increases its efficiency superiority over the others as the system size or the number of interested eigenvectors become larger and larger.

Ting, T.↗

Coherent fluctuations of extratropical geopotential height and tropical convection in intraseasonal time scales

Northern-Hemisphere twice-daily 500-mb-geopotential-height data for November-March 1975-1980 and November-December 1981 and NOAA satellite measurements of outgoing longwave radiation for the same period are processed to remove seasonal cycles and Fourier-decomposition bandpassed to study variations with 20-70-d periods. The results of analysis using correlation, complex EOF, and composite techniques are presented in extensive maps and graphs and characterized. Extratropical wavetrains are found to evolve systematically from Eurasia eastward to North America and the North Atlantic on 5-6-d time scales, while the intraseasonal variation in tropical convection is dominated by a dipolelike east-west feature propagating from the western Indian Ocean to the dateline with a quasi-period of 40-50 d. The possibility that normal modes coupled between the tropics and midlatitudes may be responsible for these phenomena is considered.

Lau, K.-M.↗

Principal modes of atmospheric circulation anomalies associated with global angular momentum fluctuations

This paper provides a description of the variability of global atmospheric angular momentum (GAM) and its relationship with principal modes of three-dimensional atmospheric circulation anomalies. The data used are 5-day mean global wind fields from the European Centre for Medium-Range Weather Forecasts initialized dataset for 1980-1989. Significant seasonal variation of GAM is observed with maxima in April and November and a minimum during late July. The amplitude of the annual cycle is largest in the upper troposphere and decreases toward the surface. Although the lower tropospheric contribution to the total angular momentum is relatively small, its annual cycle is out of phase with those of the upper atmosphere and GAM. Also identified is a distinct semiannual component, with double peaks appearing in April and November. This signal is most noticeable in the upper troposphere above the 300-mb level. The principal modes of zonal-mean angular momentum and meridional circulation anomalies and their coupled modes are obtained by using empirical orthogonal function analysis and singular value decomposition. It is shown that the leading modes of the angular momentum and meridional circulation are coupled with each other and are responsible for much of the variability in GAM. The coupled modes represent fluctuations of upper-level subtropical zonal flow, which are linked to the modulation of Hadley circulation intensity in both hemispheres. It is found that GAM is highly correlated with the first eigenvector of upper-level streamfunction anomalies, which consists of a superrotational flow in the tropics and subtropics, except over the central Pacific where a 'blocked' flow with two subtropical anticyclonic circulation cells straddling the equator is found. Much of the blocked flow is due to the establishment of dipole anomalies in the velocity potential with centers over the central Pacific and the Maritime Continent on the interannual time scale. On the intraseasonal time scale, GAM fluctuation is dominated by superrotational flow in the tropics, with the blocked flow present to a much lesser extent. The associated velocity potential anomaly has a weak dipole structure with centers over the Indian Ocean and the eastern Pacific. The implications of the above results on the total angular momentum balance of the earth-atmosphere system are also discussed.

Kang, In-Sik↗

Tidal decomposition of zonal neutral and ion flows in the earth's upper equatorial thermosphere

Evidence is presented for strong coupling between the diurnal components of zonal neutral winds and ion drifts, suggesting that the relative importance of the E- and F-region dynamos be reevaluated. Measurements of zonal neutral winds in the equatorial region of the earth's thermosphere at an average altitude of about 350 km show that the nighttime zonal winds are very similar to the zonal ion-drifts. That similarity is examined, comparing the corresponding tidal components of the 24 hr variations of these two parameters. The amplitude spectrum of the neutral winds exhibits primary and secondary maxima at the diurnal and ter-diurnal frequencies respectively, while the ion-drift spectrum shows only the diurnal maximum. It is found that the simularity between neutral winds and ion-drifts is strongest in the diurnal mode where the phases differ by less than one half hour, the amplitude of the ion-drift being between 70 percent and 80 percent that of the neutral wind, suggesting a first-order relation between the two quantities. The largest difference is found in the steady component representing superrotation; under similar conditions of solar activity, the ions superrotate with a velocity of about 30 m/s and the neutrals with 10 m/s. For the ions, the steady component, the phase of the semi-diurnal component and the amplitude of the ter-diurnal component appear to be sensitive to solar activity and are responsible for the observed solar cycle variations in the times of eastward-to-westward reversals between 0400 and 0700 LT. The ion-drift diurnal amplitude and phase are relatively insensitive to changes in solar activity.

Herrero, F. A.↗

An Alternate Method for Estimating Dynamic Height from XBT Profiles Using Empirical Vertical Modes

A technique is presented that applies modal decomposition to estimate dynamic height (0-450 db) from Expendable BathyThermograph (XBT) temperature profiles. Salinity-Temperature-Depth (STD) data are used to establish empirical relationships between vertically integrated temperature profiles and empirical dynamic height modes. These are then applied to XBT data to estimate dynamic height. A standard error of 0.028 dynamic meters is obtained for the waters of the Gulf of Alaska- an ocean region subject to substantial freshwater buoyancy forcing and with a T-S relationship that has considerable scatter. The residual error is a substantial improvement relative to the conventional T-S correlation technique when applied to this region. Systematic errors between estimated and true dynamic height were evaluated. The 20-year-long time series at Ocean Station P (50 deg N, 145 deg W) indicated weak variations in the error interannually, but not seasonally. There were no evident systematic alongshore variations in the error in the ocean boundary current regime near the perimeter of the Alaska gyre. The results prove satisfactory for the purpose of this work, which is to generate dynamic height from XBT data for coanalysis with satellite altimeter data, given that the altimeter height precision is likewise on the order of 2-3 cm. While the technique has not been applied to other ocean regions where the T-S relation has less scatter, it is suggested that it could provide some improvement over previously applied methods, as well.

Lagerloef, Gary S. E.↗

MGS Radio Science Electron Density Profiles: Interannual Variability and Implications for the Martian Neutral Atmosphere

Martian electron density profiles provided by the Mars Global Surveyor (MGS) Radio Science (RS) experiment over the 95-200 km altitude range indicate what the height of the electron peak and the longitudinal structure of the peak height are sensitive indicators of the physical state of the Mars lower and upper atmospheres. The present analysis is carried out on five sets of occultation profiles, all at high solar zenith angles (SZA). Variations spanning 2 Martian years are investigated near aphelion conditions at high northern latitudes (64.7 - 77.6 N) making use of four of these data sets. A mean ionospheric peak height of 133.5 - 135 km is obtained near SZA = 78 - 82 deg.; a corresponding mean peak density of 7.3 - 8.5 x l0(exp 4)/ qu cm is also measured during solar moderate conditions at Mars. Strong wave number 2 - 3 oscillations in peak heights are consistently observed as a function of longitude over the 2 Martian years. These observed ionospheric features are remarkably similar during aphelion conditions 1 Martian year apart. This year-to-year repeatability in the thermosphere-ionosphere structure is consistent with that observed in multiyear aphelion temperature data of the Mars lower atmosphere. Coupled Mars general circulation model (MGCM) and Mars thermospheric general circulation model (MTGCM) codes are run for Mars aphelion conditions, yielding mean and longitude variable ionospheric peak heights that reasonably match RS observations. A tidal decomposition of MTGCM thermospheric densities shows that observed ionospheric wave number 3 features are linked to a non-migrating tidal mode with semidiurnal period (sigma = 2) and zonal wave number 1 (s = -1) characteristics. The height of this photochemically determined ionospheric peak should be monitored regularly.

Bougher, S. W.↗

Parameterizing Coefficients of a POD-Based Dynamical System

A method of parameterizing the coefficients of a dynamical system based of a proper orthogonal decomposition (POD) representing the flow dynamics of a viscous fluid has been introduced. (A brief description of POD is presented in the immediately preceding article.) The present parameterization method is intended to enable construction of the dynamical system to accurately represent the temporal evolution of the flow dynamics over a range of Reynolds numbers. The need for this or a similar method arises as follows: A procedure that includes direct numerical simulation followed by POD, followed by Galerkin projection to a dynamical system has been proven to enable representation of flow dynamics by a low-dimensional model at the Reynolds number of the simulation. However, a more difficult task is to obtain models that are valid over a range of Reynolds numbers. Extrapolation of low-dimensional models by use of straightforward Reynolds-number-based parameter continuation has proven to be inadequate for successful prediction of flows. A key part of the problem of constructing a dynamical system to accurately represent the temporal evolution of the flow dynamics over a range of Reynolds numbers is the problem of understanding and providing for the variation of the coefficients of the dynamical system with the Reynolds number. Prior methods do not enable capture of temporal dynamics over ranges of Reynolds numbers in low-dimensional models, and are not even satisfactory when large numbers of modes are used. The basic idea of the present method is to solve the problem through a suitable parameterization of the coefficients of the dynamical system. The parameterization computations involve utilization of the transfer of kinetic energy between modes as a function of Reynolds number. The thus-parameterized dynamical system accurately predicts the flow dynamics and is applicable to a range of flow problems in the dynamical regime around the Hopf bifurcation. Parameter-continuation software can be used on the parameterized dynamical system to derive a bifurcation diagram that accurately predicts the temporal flow behavior.

Kalb, Virginia L.↗

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗

Impacts of ENSO on the South American Summer Monsoon During 1997-1999

Using the National Center for Environmental Prediction (NCEP) Reanalysis, and CPC Merged Analysis Product (CMAP) rainfall, we have compared and contrasted the anomalies of the South American Summer Monsoon (SASM) during two extreme years of 1997/98 (EI Nino) and 1998/99 (La Nina). The results are assessed against a "canonical" ENSO response (CER) pattern for the SASM obtained from empirical mode decomposition based on a previous period (1980-1995). Overall, the SASM anomalies compare well with CER, but with some important differences. Anomalies occurring in the warm phase of the 1997-98 El Nino are very significant and robust, while those occurring in 1998/99 La Nina, appear to be reversed from 1997/98, but are relatively weak and less well-defined. The most pronounced signal in DJF 1997/98 is the development of drought conditions in northern Brazil, excessive rainfall over northern Peru and Ecuador, and over Uruguay and southern Brazil. The tropical rainfall anomalies are associated with the eastward shift of the Walker circulation, which is represented by pronounced low-level anomalous westerlies over the equatorial eastern Pacific and easterlies over northern Brazil. The easterlies are deflected sharply southeastward by the steep topography of the Andes, enhancing the low-level jet (LLJ) along the eastern foothills of the Andes near 15-20 S. The LLJ penetrates deep into the extratropics, yielding rainfall anomalies further poleward compared to CER. During DJF 1997/98, the eastward expansion of the warm tropospheric temperature over the Nino-3 region causes anomalous geopotential height to develop in the upper troposphere above the Altiplano, leading to a strengthened Bolivian High. An upper-tropospheric jet anomaly maximum is found over the subtropical continent near 30 S, due to increasing meridional gradient of tropospheric temperature, as well as teleconnection patterns linking the South Pacific and the South Atlantic. Consistent with the CER, the South Pacific High is weakened, and the South Atlantic High is strengthened in DJF 1997/98. However, rather than appearing as a coherent large scale signal over the entire Atlantic as in CER, the South Atlantic High anomalies in DJF 1997/98 split into two separate high pressure centers; one located south of the Amazon basin, and another over the southeastern South Atlantic. The former is induced as a Rossby-wave response by ENSO-Induced sinking motion over nor-them and eastern Brazil, and the latter may be associated with extratropical teleconnection signals. Anomalous convection observed over equatorial western Africa may also have contributed to the different circulation response over the Atlantic in DFJ 1997/99 compared to the CER. Comparison of the evolution of the SASM anomalies in 1997-99 with their corresponding annual cycles, suggest that the anomalies are phase-locked to and represent either enhancement or reduction of the annual variations. Results suggest that, from a system perspective the SASM was strengthened during the 1997/98 El Nino and weakened, albeit to a lesser extent, during the 1998/99 La Nina.

Lau, K.-M.↗