Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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,405 records · Page 78

Second-Order Moller-Plesset Perturbation Theory for Molecular Dirac-Hartree-Fock Wave Functions

Moller-Plesset perturbation theory is developed to second order for a selection of Kramers restricted Dirac-Hartree-Fock closed and open-shell reference wave functions. The open-shell wave functions considered are limited to those with no more than two electrons in open shells, but include the case of a two-configuration SCF reference. Denominator shifts are included in the style of Davidson's OPT2 method. An implementation which uses unordered integrals with labels is presented, and results are given for a few test cases.

Dyall, Kenneth G.↗

Assessment of autonomic response by broad-band respiration

We present a technique for introducing broad-band respiratory perturbations so that the response characteristics of the autonomic nervous system can be determined noninvasively over a wide range of physiologically relevant frequencies. A subject's respiratory bandwidth was broadened by breathing on cue to a sequence of audible tones spaced by Poisson intervals. The transfer function between the respiratory input and the resulting instantaneous heart rate was then computed using spectral analysis techniques. Results using this method are comparable to those found using traditional techniques, but are obtained with an economy of data collection.

NASA Discipline Cardiopulmonary↗

Methods studies toward simplified rotor-body dynamics, part 1

This report is directed to the problem of developing an adequate but not overly complex linear flight dynamics analytical model of a rotorcraft to study stability, control, gust and random turbulence responses. Since the conventional flight dynamics analysis using quasisteady rotor derivatives is adequate for the long period modes like the phugoid mode, only short time responses are considered here, where rotor-body coupling is of importance. Thus the body motion consists of pitch, roll and vertical motion, omitting linear longitudinal and lateral and yaw perturbations. Five analytical models of varying degree of sophistication are applied to a hypothetical hingeless compound helicopter operating up to .8 rotor advance ratio. Stability and response data are obtained for the basic helicopter and for the vehicle with two simple control feedback systems.

Hohenemser, K. H.↗

Uncertainty in Heart Rate Complexity Metrics Caused by R-Peak Perturbations

Heart rate complexity (HRC) is a proven metric for gaining insight into human stress and physiological deterioration. To calculate HRC, the detection of the exact instance of when the heart beats, the R-peak, is necessary. Electrocardiogram (ECG) signals can often be corrupted by environmental noise (e.g., from electromagnetic interference, movement artifacts), which can potentially alter the HRC measurement, producing erroneous inputs which feed into decision support models. Current literature has only investigated how HRC is affected by noise when R-peak detection errors occur (false positives and false negatives). However, the numerical methods used to calculate HRC are also sensitive to the specific location of the fiducial point of the R-peak. This raises many questions regarding how this fiducial point is altered by noise, the resulting impact on the measured HRC, and how we can account for noisy HRC measures as inputs into our decision models. This work uses Monte Carlo simulations to systematically add white and pink noise at different permutations of signal-to-noise ratios (SNRs), time segments, sampling rates, and HRC measurements to characterize the influence of noise on the HRC measure by altering the fiducial point of the R-peak. Using the generated information from these simulations provides improved decision processes for system design which address key concerns such as permutation entropy being a more precise, reliable, less biased, and more sensitive measurement for HRC than sample and approximate entropy.

Napoli, Nicholas J.↗

Optical Modeling Activities for the James Webb Space Telescope (JWST) Project: Determining Image Motion and Wavefront Error Over an Extended Field of View with a Segmented Optical System - II

This is part two of a series on the optical modeling activities for JWST. Starting with the linear optical model discussed in part one, we develop centroid and wavefront error sensitivities for the special case of a segmented optical system such as JWST, where the primary mirror consists of 18 individual segments. Our approach extends standard sensitivity matrix methods used for systems consisting of monolithic optics, where the image motion is approximated by averaging ray coordinates at the image and residual wavefront error is determined with global tip/tilt removed. We develop an exact formulation using the linear optical model, and extend it to cover multiple field points for performance prediction at each instrument aboard JWST. This optical model is then driven by thermal and dynamic structural perturbations in an integrated modeling environment. Results are presented.

Howard, Joseph M.↗

The ground signatures of the expansion phase during multiple onset substorms

The paper examines the signatures of multiple-expansion-phase substorm onsets in night-time magnetograms as well as in recordings of auroral activity and energetic electron precipitation. Individual onsets of Pi2 magnetic pulsations observed at three widely spaced stations are used to define and time each onset accurately in order to distinguish between local variations in bay activity and fully developed substorm onsets; this method is equivalent to defining each new onset in terms of the brightening of an auroral arc and the formation of a westward-travelling surge. It is found that the formation of multiple auroral surges appears to be a fundamental feature of multiple-onset substorms, that each surge seems to be associated with a localized field-aligned current system which moves westward with the surge and perturbs the preexisting current wedge, and that this gives rise to the multiple-onset signatures observed on subauroral and low-latitude magnetograms. Since certain findings contradict existing models of multiple-onset events, an alternative model is proposed which is based on the fundamental role of surge activity and localized current wedges.

Pytte, T.↗

Optical phase distortion due to compressible flow over laser turrets

Analytical models for optical phase distortion due to compressible flow over a laser turret are developed. Phase distortion is calculated for both blunt and small perturbation turrets. For the blunt turret, the Janzen-Rayleigh technique is used to determine the flow field. Phase distortions of 2.2 wavelengths at 3.8 microns are calculated for the blunt turret. For small perturbation turrets, a versatile analytical model is developed for a turret on a fuselage with circular cross section. With a two dimensional Fourier series representation of the turret, any shape can be considered. Both subsonic and supersonic flows can be calculated. Phase distortions of 1.2 wavelengths at 3.8 microns are calculated for one turret at high subsonic Mach number. In addition to being of value for laser turrets, the methods are applicable to reconnaissance aircraft using photographic equipment and cruise missiles using celestial navigation.

Fuhs, A. E.↗

Learned pattern recognition using synthetic-discriminant-functions

A method of using synthetic-discriminant-functions to facilitate learning in a pattern recognition system is discussed. Learning is accomplished by continually adding images to the training set used for synthetic discriminant functions (SDF) construction. Object identification is performed by efficiently searching a library of SDF filters for the maximum optical correlation. Two library structures are discussed - binary tree and multilinked graph - along with maximum ascent, back-tracking, perturbation, and simulated annealing searching techniques. By incorporating the distortion invariant properties of SDFs within a library structure, a robust pattern recognition system can be produced.

Jared, David A.↗

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 1: Code validation and stationary accretors

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing stationary spheres of varying sizes (from 10.0 down to 0.02 Bondi radii) accrete matter from a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3 or 1.2). To accommodate the long-range gravitational forces, the extent of the computational volume is typically a factor of 100 larger than the radius of the accretor. We compare the numerical mass accretion rates with the theoretical predictions of Bondi, to assess the validity of the code. The hydrodynamics is modeled by the piecewise parabolic method. No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (6-8) grids around the stationary sphere, each finer grid being a factor of 2 smaller spatially than the next coarser grid. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid while at the same time evolving the gas on the coarser grids. The accretion rates derived numerically are in in very good agreement (to about 10% over several orders of magnitude) with the values given by Bondi for a stationary accretor within a hydrodynamic medium. However, the equations have to be changed in order to include the finite size of the accretor (in some cases very large compared to the sonic point or even to the Bondi radius).

Ruffert, Maximilian↗

Time-Separating Heating and Sensor Functions of Thermistors in Precision Thermal Control Applications

A method allows combining the functions of a heater and a thermometer in a single device, a thermistor, with minimal temperature read errors. Because thermistors typically have a much smaller thermal mass than the objects they monitor, the thermal time to equilibrate the thermometer to the temperature of the object is typically much shorter than the thermal time of the object to change its temperature in response to an external perturbation.

Cho, Hyung J.↗

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

The perfectly matched layer (PML) technique is developed in the context of a high- order spectral-element Discontinuous-Galerkin (DG) method. The technique is applied to a range of test cases and is shown to be superior compared to other approaches, such as those based on using characteristic boundary conditions and sponge layers, for treating the inflow and outflow boundaries of computational domains. In general, the PML technique improves the quality of the numerical results for simulations of practical flow configurations, but it also exhibits some instabilities for large perturbations. A preliminary analysis that attempts to understand the source of these instabilities is discussed.

Spectral-Element↗

Study of modal coupling procedures for the shuttle: A matrix method for damping synthesis

The damping method was applied successfully to real structures as well as analytical models. It depends on the ability to determine an appropriate modal damping matrix for each substructure. In the past, modal damping matrices were assumed diagonal for lack of being able to determine the coupling terms which are significant in the general case of nonproportional damping. This problem was overcome by formulating the damped equations of motion as a linear perturbation of the undamped equations for light structural damping. Damped modes are defined as complex vectors derived from the complex frequency response vectors of each substructure and are obtained directly from sinusoidal vibration tests. The damped modes are used to compute first order approximations to the modal damping matrices. The perturbation approach avoids ever having to solve a complex eigenvalue problem.

Hasselman, T. K.↗

Evaluation of turbulence induced noise in coherent anti-Stokes Raman scattering

The effect of turbulence in a transonic wind tunnel on coherent anti-Stokes Raman scattering is considered. The driving pump and Stokes waves are taken to be coaxially propagating Gaussian beam waves which are focused on the Raman active medium through the turbulent boundary layer of the flow tube. The random index of refraction variations in the layer are modeled as phase perturbations of the driving waves which cause a reduction of the mean on-axis field and an increase in the mean diameter of the beams. Effective Gaussian beam parameters are developed and the radiated anti-Stokes power calculated as a function of the phase screen parameters. A significant reduction in signal strength occurs for realistic estimates of the phase screen parameter appropriate to a confined transonic flow. A method for estimating the signal degradation which could be applied to other experimental situations is presented.

Elliott, R. A.↗

The Normalization of Surface Anisotropy Effects Present in SEVIRI Reflectances by Using the MODIS BRDF Method

A modified version of the MODerate resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) algorithm is presented for use in the angular normalization of surface reflectance data gathered by the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) aboard the geostationary Meteosat Second Generation (MSG) satellites. We present early and provisional daily nadir BRDFadjusted reflectance (NBAR) data in the visible and near-infrared MSG channels. These utilize the high temporal resolution of MSG to produce BRDF retrievals with a greatly reduced acquisition period than the comparable MODIS products while, at the same time, removing many of the angular perturbations present within the original MSG data. The NBAR data are validated against reflectance data from the MODIS instrument and in situ data gathered at a field location in Africa throughout 2008. It is found that the MSG retrievals are stable and are of high-quality across much of the SEVIRI disk while maintaining a higher temporal resolution than the MODIS BRDF products. However, a number of circumstances are discovered whereby the BRDF model is unable to function correctly with the SEVIRI observations-primarily because of an insufficient spread of angular data due to the fixed sensor location or localized cloud contamination.

SEVIRI↗

Time transformations and Cowell's method

The precise numerical integration of Cowell's equations of satellite motion is frequently performed with an independent variable s defined by an equation of the form dt = cr to the n-th power ds, where t represents time, r the radial distance from the center of attraction, c is a constant, and n is a parameter. This has been primarily motivated by the 'uniformizing' effects of such a transformation resulting in desirable 'analytic' stepsize control for elliptical orbits. This report discusses the 'proper' choice of the parameter n defining the independent variable s for various types of orbits and perturbation models, and develops a criterion for its selection.

Velez, C. E.↗

Mechanical forces and their second messengers in stimulating cell growth in vitro

Mechanical forces play an important role in modulating the growth of a number of different tissues including skeletal muscle, smooth muscle, cardiac muscle, bone, endothelium, epithelium, and lung. As interest increases in the molecular mechanisms by which mechanical forces are transduced into growth alterations, model systems are being developed to study these processes in tissue culture. This paper reviews the current methods available for mechanically stimulating tissue cultured cells. It then outlines some of the putative 'mechanogenic' second messengers involved in altering cell growth. Not surprisingly, many mechanogenic second messengers are the same as those involved in growth factor-induced cell growth. It is hypothesized that from an evolutionary standpoint, some second messenger systems may have initially evolved for unicellular organisms to respond to physical forces such as gravity and mechanical perturbation in their environment. As multicellular organisms came into existence, they appropriated these mechanogenic second messenger cascades for cellular regulation by growth factors.

Vandenburgh, Herman H.↗

Classifying Multi-Model Wheat Yield Impact Response Surfaces Showing Sensitivity to Temperature and Precipitation Change

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (minus 2 to plus 9 degrees Centigrade) and precipitation (minus 50 to plus 50 percent). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Classification; Climate change; Crop model; Ensemb↗

Aura/MLS Observes and SD-WACCM-X Simulates the Seasonality, Quasi-Biennial Oscillation and El Niño-Southern Oscillation of the Migrating Diurnal Tide Driving Upper Mesospheric CO Primarily Through Vertical Advection

This work uses 17 years of upper mesospheric carbon monoxide (CO) and temperature observations by the Microwave Limb Sounder (MLS) on-board the Aura satellite to present and explain the seasonal and interannual variability of the migrating diurnal tide (DW1) component of upper mesospheric CO. This work then compares these observations to simulations by the Specified Dynamics – Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X). Results show that, for all seasons, MLS CO local-time perturbations peaks above 85 km and has a latitude structure resembling the (1,1) mode in temperature. On the other hand, SD-WACCM-X DW1 also peaks above 85 km and has a latitude structure resembling the (1,1) mode but it simulates 2 local maximum of the (1,1) mode between 85 km and 92 km. Despite the differences in altitude structure, a tendency analysis and the adiabatic displacement method revealed that, on seasonal and interannual timescales, observed and modelled CO’s (1,1) component can be reproduced solely using vertical advection. It was also found that both observed and modelled CO’s (1,1) component contains interannual oscillations with periodicities close to that of the Quasi-biennial Oscillation and the El Niño–Southern Oscillation. From these results, this work concludes that on seasonal and interannual timescales, the observed and modelled (1,1) mode affects the global structure of upper mesospheric CO primarily through vertical advection.

Cornelius Csar Jude H Salinas↗