Engineering PapersSearch

SEARCH · Engineering Papers

Results for “consistent discretization”

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 55 records · Page 3

Contouring Trivariate Data

In many applications, the data consist of discrete 3-D points at which one or more parameters are given. To display contours, the data are represented by a continuous function which is evaluated at any point as needed for contouring. Contouring results are presented which are applicable both to artitrarily spaced data and to data which lie on a topologically rectangular three dimensional grid. the contours are assumed to be described by 3-D display lists for viewing on a dynamic color graphics device; that is, they will not simply be projected into 2-D and viewed as a static image on a frame buffer. Dynamic viewing of color contours may be essential to the proper interpretation of results. The gridded data lie on a topologically rectangular grid although two or more nodes may be the same point. Parametric tensor product methods may be used to fit the gridded data and to generate the contours. Rectangular elements are convenient but not necessary. For example, there are other methods which are effective for contouring over tetrahedral elements.

Stead, S. E.

The detection of high-velocity outflows from M8E-IR

A high-resolution (0.059/cm) M band (4.6 micron) spectrum of the embedded young stellar object M8E-IR is presented and discussed. The spectrum shows strong absorption to large blueshifts in the rotational lines of the fundamental vibrational band, v = 1-0, of CO. The absorption is interpreted as being due to gas near to, and flowing from, the central object. The outflowing gas is warm (95-330 K) and consists of discrete velocity components with the very high velocities of 90, 130, 150, and 160 km/s. On the basis of a simple model, it is estimated that the observed outflows are less than 100 yr old.

Mitchell, George F.

A search for water on the moon at the Reiner Gamma Formation - A possible site of cometary coma impact

Earth-based telescopic measurements of the Reiner Gamma Formation in the 3-micron region were made to search for evidence of water at this proposed site of cometary interaction with the lunar surface. Comparison of the spectra of Reiner Gamma and laboratory measurements of a variety of minerals with adsorbed water or structural OH(-) show that there is no evidence for water at the locations measured. While these measurements rule out the presence of abundant bound H2O at these locations, they do not exclude the presence of OH(-) bearing minerals that have absorptions not detectable through the wet terrestrial atmosphere. The possible effects of high lunar temperatures on the reflectance spectra were examined, and it was concluded that these effects would not obscure a bound water absorption if it were present. Upper limits on the amount of water that could be present at the locations observed and remain undetected are 2.25 average wt pct for a mixture consisting of discrete patches of water-bearing and water-free minerals and 0.01 wt pct for an intimate mixture of these two materials.

Roush, Ted L.

Enhanced element-specific modal formulations for flexible multibody dynamics

The accuracy of current flexible multibody formalisms based on assumed modes is examined in the context of standard spacecraft motions involving structural components undergoing both slow and fast overall translational and rotational motions as well as small deformations. Limitations of current techniques in treating: (1) element-specific coupling behavior of large motion and small deformation, and (2) motion-induced structural stiffness variations, are noted. The roles of nonlinear and linear elastic structural theories in accurately predicting transient large-displacement dynamic behavior of flexible multibody systems are examined in detail. Coupling effects between deformation and overall motion are carefully scrutinized in the context of assumed-mode discretization techniques. Consistently linearized beam, plate, and shell formulations involving in-plane stretch variables are proposed and shown to yield very accurate simulation results and extremely fast modal convergence for most motions involving small strains. In some particular cases, however, in which membrane stiffness dominates bending stiffness, a nonlinear strain formulation is required in order to capture proper coupling between deformation and overall motion.

Ryan, Robert R.

Scattering from randomly oriented scatterers with strong permittivity fluctuations

Strong permittivity fluctuation theory is used to solve the problem of scattering from a medium composed of completely randomly oriented scatterers under a low frequency limit. Gaussian statistics are not assumed for the renormalized scattering sources. Numerical results on effective permittivity are illustrated for oblate and prolate spheroidal scatterers and compared with the results for spherical scatterers. The results are consistent with discrete scatterer theory. The effective permittivity of a random medium embedded with nonspherical scatterers shows a higher imaginary part than the spherical scatterer case with equal correlation volume. Under the distorted Born approximation, the polarimetric covariance matrix for the backscattered electric field is calculated for half-space randomly oriented scatterers. The nonspherical geometry of the scatterers shows significant effects on the cross-polarized backscattering returns, and the correlation coefficient between HH and VV returns. The polarimetric backscattering coefficients can provide useful information in distinguishing the geometry of scatterers.

Yueh, S. H.

An Algorithm for Integrated Subsystem Embodiment and System Synthesis

Consider the statement,'A system has two coupled subsystems, one of which dominates the design process. Each subsystem consists of discrete and continuous variables, and is solved using sequential analysis and solution.' To address this type of statement in the design of complex systems, three steps are required, namely, the embodiment of the statement in terms of entities on a computer, the mathematical formulation of subsystem models, and the resulting solution and system synthesis. In complex system decomposition, the subsystems are not isolated, self-supporting entities. Information such as constraints, goals, and design variables may be shared between entities. But many times in engineering problems, full communication and cooperation does not exist, information is incomplete, or one subsystem may dominate the design. Additionally, these engineering problems give rise to mathematical models involving nonlinear functions of both discrete and continuous design variables. In this dissertation an algorithm is developed to handle these types of scenarios for the domain-independent integration of subsystem embodiment, coordination, and system synthesis using constructs from Decision-Based Design, Game Theory, and Multidisciplinary Design Optimization. Implementation of the concept in this dissertation involves testing of the hypotheses using example problems and a motivating case study involving the design of a subsonic passenger aircraft.

Lewis, Kemper

Delaying Trains of Short Light Pulses in WGM Resonators

Suitably configured whispering-gallery-mode (WGM) optical resonators have been proposed as delay lines for trains of short light pulses. Until now, it has been common practice to implement an optical delay line as a coiled long optical fiber, which is bulky and tends to be noisy. An alternative has been to implement an optical delay line as a coupled-resonator optical waveguide (a chain of coupled optical resonators), which is compact but limits the width of the pulse spectrum to the width of an optical resonance and thereby places a lower limit on the duration of a pulse. In contrast, a delay line according to the proposal could be implemented as a single WGM resonator, and the pulses delayed by the resonator could be so short that their spectral widths could greatly exceed the spectral width of any single resonance. The proposal emerged from theoretical and experimental studies of the propagation of a pulse train in a WGM resonator. An important element of the theoretical study was recognition that the traditional definition of group velocity in effect, the velocity of a single pulse comprising a packet of waves propagating in a medium, the responsivity of which is a monotonous function of frequency does not necessarily apply in the case of a WGM resonator or other medium having a spectrum consisting of discrete resonance peaks at different frequencies. A new definition of group velocity, applicable to a train of pulses propagating in such a medium, was introduced and found to lead to the discovery of previously unknown features of propagation. Notably, it was found that in a micro - sphere optical resonator that supports a suitable combination of WGM modes, the group velocity for a train of light pulses could be positive, zero, or negative. A positive group velocity could be so small that the delay could be much longer than the ring-down time of the resonator; a delay of such great length is impossible for a single pulse interacting with either a linearly responding lossless resonator or a coupled- resonator optical waveguide.

Matsko, Andrey

Smooth Phase Interpolated Keying

Smooth phase interpolated keying (SPIK) is an improved method of computing smooth phase-modulation waveforms for radio communication systems that convey digital information. SPIK is applicable to a variety of phase-shift-keying (PSK) modulation schemes, including quaternary PSK (QPSK), octonary PSK (8PSK), and 16PSK. In comparison with a related prior method, SPIK offers advantages of better performance and less complexity of implementation. In a PSK scheme, the underlying information waveform that one seeks to convey consists of discrete rectangular steps, but the spectral width of such a waveform is excessive for practical radio communication. Therefore, the problem is to smooth the step phase waveform in such a manner as to maintain power and bandwidth efficiency without incurring an unacceptably large error rate and without introducing undesired variations in the amplitude of the affected radio signal. Although the ideal constellation of PSK phasor points does not cause amplitude variations, filtering of the modulation waveform (in which, typically, a rectangular pulse is converted to a square-root raised cosine pulse) causes amplitude fluctuations. If a power-efficient nonlinear amplifier is used in the radio communication system, the fluctuating-amplitude signal can undergo significant spectral regrowth, thus compromising the bandwidth efficiency of the system. In the related prior method, one seeks to solve the problem in a procedure that comprises two major steps: phase-value generation and phase interpolation. SPIK follows the two-step approach of the related prior method, but the details of the steps are different. In the phase-value-generation step, the phase values of symbols in the PSK constellation are determined by a phase function that is said to be maximally smooth and that is chosen to minimize the spectral spread of the modulated signal. In this step, the constellation is divided into two groups by assigning, to information symbols, phase values that result in equal numbers of clockwise and counter-clockwise phase rotations for equally likely symbols. The purpose served by assigning phase values in this way is to prevent unnecessary generation of spectral lines and prevent net shifts of the carrier signal. In the phase-interpolation step, the smooth phase values are interpolated over a number, n, of consecutive symbols (including the present symbol) by means of an unconventional spline curve fit.

Borah, Deva K.

A Structural Model Decomposition Framework for Hybrid Systems Diagnosis

Nowadays, a large number of practical systems in aerospace and industrial environments are best represented as hybrid systems that consist of discrete modes of behavior, each defined by a set of continuous dynamics. These hybrid dynamics make the on-line fault diagnosis task very challenging. In this work, we present a new modeling and diagnosis framework for hybrid systems. Models are composed from sets of user-defined components using a compositional modeling approach. Submodels for residual generation are then generated for a given mode, and reconfigured efficiently when the mode changes. Efficient reconfiguration is established by exploiting causality information within the hybrid system models. The submodels can then be used for fault diagnosis based on residual generation and analysis. We demonstrate the efficient causality reassignment, submodel reconfiguration, and residual generation for fault diagnosis using an electrical circuit case study.

structural model decomposition

Introduction to COFFE: The Next-Generation HPCMP CREATE-AV CFD Solver

HPCMP CREATE-AV Conservative Field Finite Element (COFFE) is a modular, extensible, robust numerical solver for the Navier-Stokes equations that invokes modularity and extensibility from its first principles. COFFE implores a flexible, class-based hierarchy that provides a modular approach consisting of discretization, physics, parallelization, and linear algebra components. These components are developed with modern software engineering principles to ensure ease of uptake from a user's or developer's perspective. The Streamwise Upwind/Petrov-Galerkin (SU/PG) method is utilized to discretize the compressible Reynolds-Averaged Navier-Stokes (RANS) equations tightly coupled with a variety of turbulence models. The mathematics and the philosophy of the methodology that makes up COFFE are presented.

Glasby, Ryan S.

Qualitative Fault Isolation of Hybrid Systems: A Structural Model Decomposition-Based Approach

Quick and robust fault diagnosis is critical to ensuring safe operation of complex engineering systems. A large number of techniques are available to provide fault diagnosis in systems with continuous dynamics. However, many systems in aerospace and industrial environments are best represented as hybrid systems that consist of discrete behavioral modes, each with its own continuous dynamics. These hybrid dynamics make the on-line fault diagnosis task computationally more complex due to the large number of possible system modes and the existence of autonomous mode transitions. This paper presents a qualitative fault isolation framework for hybrid systems based on structural model decomposition. The fault isolation is performed by analyzing the qualitative information of the residual deviations. However, in hybrid systems this process becomes complex due to possible existence of observation delays, which can cause observed deviations to be inconsistent with the expected deviations for the current mode in the system. The great advantage of structural model decomposition is that (i) it allows to design residuals that respond to only a subset of the faults, and (ii) every time a mode change occurs, only a subset of the residuals will need to be reconfigured, thus reducing the complexity of the reasoning process for isolation purposes. To demonstrate and test the validity of our approach, we use an electric circuit simulation as the case study.

Hybrid Systems Diagnosis

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan

A prefiltering version of the Kalman filter with new numerical integration formulas for Riccati equations

A prefiltering version of the Kalman filter is derived for both discrete and continuous measurements. The derivation consists of determining a single discrete measurement that is equivalent to either a time segment of continuous measurements or a set of discrete measurements. This prefiltering version of the Kalman filter easily handles numerical problems associated with rapid transients and ill-conditioned Riccati matrices. Therefore, the derived technique for extrapolating the Riccati matrix from one time to the next constitutes a new set of integration formulas which alleviate ill-conditioning problems associated with continuous Riccati equations. Furthermore, since a time segment of continuous measurements is converted into a single discrete measurement, Potter's square root formulas can be used to update the state estimate and its error covariance matrix. Therefore, if having the state estimate and its error covariance matrix at discrete times is acceptable, the prefilter extends square root filtering with all its advantages, to continuous measurement problems.

Womble, M. E.

Finite difference solutions of the Euler equations in the vicinity of sharp edges

Attempts have been made to explain why finite difference solutions of the Euler equations can describe flows with large vortical structures around sharp-edged bodies. The present paper is concerned with the influence of a singular sharp edge on the truncation error for a set of discretized Euler equations. An analysis is conducted of the distribution of the truncation error of one finite difference approximation of the Euler equations near a sharp edge of a thin plate. The analysis leads to a determination of the size of the region of the neighborhood of such a singularity. Attention is given to the consistency of a discretization of the Euler equations, and numerical experiments.

Hartwich, P.-M.

Electromagnetic Scattering by Fully Ordered and Quasi-Random Rigid Particulate Samples

In this paper we have analyzed circumstances under which a rigid particulate sample can behave optically as a true discrete random medium consisting of particles randomly moving relative to each other during measurement. To this end, we applied the numerically exact superposition T-matrix method to model far-field scattering characteristics of fully ordered and quasi-randomly arranged rigid multiparticle groups in fixed and random orientations. We have shown that, in and of itself, averaging optical observables over movements of a rigid sample as a whole is insufficient unless it is combined with a quasi-random arrangement of the constituent particles in the sample. Otherwise, certain scattering effects typical of discrete random media (including some manifestations of coherent backscattering) may not be accurately replicated.

Radiative transfer

The Suzaku View of Highly Ionized Outflows in AGN: Statistical Detection and Global Absorber Properties - 1

We present the results of a new spectroscopic study of Fe K-band absorption in active galactic nuclei (AGN). Using data obtained from the Suzaku public archive we have performed a statistically driven blind search for Fe XXV Healpha and/or Fe XXVI Lyalpha absorption lines in a large sample of 51 Type 1.0−1.9 AGN. Through extensive Monte Carlo simulations we find that statistically significant absorption is detected at E greater than or approximately equal to 6.7 keV in 20/51 sources at the P(sub MC) greater than or equal tov 95 per cent level, which corresponds to approximately 40 per cent of the total sample. In all cases, individual absorption lines are detected independently and simultaneously amongst the two (or three) available X-ray imaging spectrometer detectors, which confirms the robustness of the line detections. The most frequently observed outflow phenomenology consists of two discrete absorption troughs corresponding to Fe XXV Healpha and Fe XXVI Lyalpha at a common velocity shift. From xstar fitting the mean column density and ionization parameter for the Fe K absorption components are log (N(sub H) per square centimeter)) is approximately equal to 23 and log (Xi/erg centimeter per second) is approximately equal to 4.5, respectively. Measured outflow velocities span a continuous range from less than1500 kilometers per second up to approximately100 000 kilometers per second, with mean and median values of approximately 0.1 c and approximately 0.056 c, respectively. The results of this work are consistent with those recently obtained using XMM-Newton and independently provides strong evidence for the existence of very highly ionized circumnuclear material in a significant fraction of both radio-quiet and radio-loud AGN in the local universe.

Gofford, Jason

Learning in stochastic neural networks for constraint satisfaction problems

Researchers describe a newly-developed artificial neural network algorithm for solving constraint satisfaction problems (CSPs) which includes a learning component that can significantly improve the performance of the network from run to run. The network, referred to as the Guarded Discrete Stochastic (GDS) network, is based on the discrete Hopfield network but differs from it primarily in that auxiliary networks (guards) are asymmetrically coupled to the main network to enforce certain types of constraints. Although the presence of asymmetric connections implies that the network may not converge, it was found that, for certain classes of problems, the network often quickly converges to find satisfactory solutions when they exist. The network can run efficiently on serial machines and can find solutions to very large problems (e.g., N-queens for N as large as 1024). One advantage of the network architecture is that network connection strengths need not be instantiated when the network is established: they are needed only when a participating neural element transitions from off to on. They have exploited this feature to devise a learning algorithm, based on consistency techniques for discrete CSPs, that updates the network biases and connection strengths and thus improves the network performance.

Johnston, Mark D.

Development of Position-Sensitive Magnetic Calorimeters for X-Ray Astronomy

Metallic magnetic calorimeters (MMC) are one of the most promising devices to provide very high energy resolution needed for future astronomical x-ray spectroscopy. MMC detectors can be built to large detector arrays having thousands of pixels. Position-sensitive magnetic (PoSM) microcalorimeters consist of multiple absorbers thermally coupled to one magnetic micro calorimeter. Each absorber element has a different thermal coupling to the MMC, resulting in a distribution of different pulse shapes and enabling position discrimination between the absorber elements. PoSMs therefore achieve the large focal plane area with fewer number of readout channels without compromising spatial sampling. Excellent performance of PoSMs was achieved by optimizing the designs of key parameters such as the thermal conductance among the absorbers, magnetic sensor, and heat sink, as well as the absorber heat capacities. Micro fab ri - cation techniques were developed to construct four-absorber PoSMs, in which each absorber consists of a two-layer composite of bismuth and gold. The energy resolution (FWHM full width at half maximum) was measured to be better than 5 eV at 6 keV x-rays for all four absorbers. Position determination was demonstrated with pulse-shape discrimination, as well as with pulse rise time. X-ray microcalorimeters are usually designed to thermalize as quickly as possible to avoid degradation in energy resolution from position dependence to the pulse shapes. Each pixel consists of an absorber and a temperature sensor, both decoupled from the cold bath through a weak thermal link. Each pixel requires a separate readout channel; for instance, with a SQUID (superconducting quantum interference device). For future astronomy missions where thousands to millions of resolution elements are required, having an individual SQUID readout channel for each pixel becomes difficult. One route to attaining these goals is a position-sensitive detector in which a large continuous or pixilated array of x-ray absorbers shares fewer numbers of temperature sensors. A means of discriminating the signals from different absorber positions, however, needs to be built into the device for each sensor. The design concept for the device is such that the shape of the temperature pulse with time depends on the location of the absorber. This inherent position sensitivity of the signal is then analyzed to determine the location of the event precisely, effectively yielding one device with many sub-pixels. With such devices, the total number of electronic channels required to read out a given number of pixels is significantly reduced. PoSMs were developed that consist of four discrete absorbers connected to a single magnetic sensor. The design concept can be extended to more than four absorbers per sensor. The thermal conductance between the sensor and each absorber is different by design and consequently, the pulse shapes are different depending upon which absorber the xrays are received, allowing position discrimination. A magnetic sensor was used in which a paramagnetic Au:Er temperature-sensitive material is located in a weak magnetic field. Deposition of energy from an x-ray photon causes an increase in temperature, which leads to a change of magnetization of the paramagnetic sensor, which is subsequently read out using a low noise dc-SQUID. The PoSM microcalorimeters are fully microfabricated: the Au:Er sensor is located above the meander, with a thin insulation gap in between. For this position-sensitive device, four electroplated absorbers are thermally linked to the sensor via heat links of different thermal conductance. One pixel is identical to that of a single-pixel design, consisting of an overhanging absorber fabricated directly on top of the sensor. It is therefore very strongly thermally coupled to it. The three other absorbers are supported directly on a silicon-nitride membrane. These absorbers are thermally coupled to the sensor via Ti (5 nm)/Au250 nm) metal links. The strength of the links is parameterized by the number of gold squares making up the link. For detector performance, experimentally different pulse-shapes were demonstrated with 6 keV x-rays, which clearly show different rise times for different absorber positions. For energy resolution measurement, the PoSM was operated at 32 mK with an applied field that was generated using a persistent current of 50 mA. Over the four pixels, energy resolution ranges from 4.4 to 4.7 eV were demonstrated.

Bandler, SImon