Search NASA⌕ Search

SEARCH · Search NASA

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

A model for the noise radiated by turbulent boundary layers and its interaction with laminar layers in supersonic flow

Motivated by the noise-boundary layer transition problem in supersonic flow, the interaction between incoming aerodynamic noise and the laminar boundary layer has been studied. The noise field has been modeled as a Mach wave radiation field consisting of discrete waves emanating from coherent turbulent entities moving downstream within the supersonic turbulent boundary layer. The individual disturbances are likened to miniature sonic booms and the laminar boundary layer is strafed by the waves as the sources move downstream. The mean, autocorrelation, and power spectral density of the field are expressed in terms of the wave shapes and their average arrival rates. The emphasis in the interaction analysis is on the behavior of the shocklets in the noise field. The shocklets are shown to be focused by the laminar boundary layer in its outer region. Borrowing wave-propagation terminology, this region is termed the caustic region. Using scaling laws from sonic boom work, focus factors at the caustic are estimated to vary from 2 to 6 for incoming shocklet strengths of 1 to .01 percent of the free stream pressure level.

Schopper, M. R.↗
NASA NTRS1984

Interaction of aerodynamic noise with laminar boundary layers in supersonic wind tunnels

The interaction between incoming aerodynamic noise and the supersonic laminar boundary layer is studied. The noise field is modeled as a Mach wave radiation field consisting of discrete waves emanating from coherent turbulent entities moving downstream within the supersonic turbulent boundary layer. The individual disturbances are likened to miniature sonic booms and the laminar boundary layer is staffed by the waves as the sources move downstream. The mean, autocorrelation, and power spectral density of the field are expressed in terms of the wave shapes and their average arrival rates. Some consideration is given to the possible appreciable thickness of the weak shock fronts. The emphasis in the interaction analysis is on the behavior of the shocklets in the noise field. The shocklets are shown to be focused by the laminar boundary layer in its outer region. Borrowing wave propagation terminology, this region is termed the caustic region. Using scaling laws from sonic boom work, focus factors at the caustic are estimated to vary from 2 to 6 for incoming shocklet strengths of 1 to .01 percent of the free stream pressure level. The situation regarding experimental evidence of the caustic region is reviewed.

Schopper, M. R.↗
NASA NTRS1984

Effect of alloy chemistry and exposure conditions on the oxidation of titanium

Multiwall is a new thermal protection system concept for advanced space transportation vehicles. The system consists of discrete panels made up of multiple layers of foil gage metal. Titanium is the proposed candidate metal for multiwall panels in the reentry temperature range up to 675 C. Oxidation and embrittlement are the principal concerns related to the use of Ti in heat shield applications. The results of a broad study on the oxidation kinetics of several titanium alloys subjected to different exposure conditions are described. The alloys include commercially pure titanium, Ti-6Al-4V, and Ti-6Al-2Sn-4Zr-2Mo. Oxidation studies were performed on these alloys exposed at 704 C in 5-760 torr air pressure and 0 to 50% relative humidity. The resulting weight gains were correlated with oxide thickness and substrate contamination. The contamination depth and weight gains due to solid solutioning were obtained from microhardness depth profiles and hardness versus weight percent oxygen calibration data.

Unnam, J.↗
NASA NTRS1984

Coma morphology and dust-emission pattern of periodic Comet Halley. I - High-resolution images taken at Mount Wilson in 1910

A radial/rotational shift-difference algorithm has been developed to improve visibility of the May-June 1910 high-resolution images taken at Mount Wilson of Comet Halley. Dust features for a period of up to three days were identified, which consisted of discrete active areas emitting dust continuously from the sunlit hemisphere of the rotating nucleus. A lower limit of the comet's rotation period was approximated at one day, and expansion velocities of the dust features were determined to be in the range of 0.2-0.3 km/s. Based on relative photometry of a bright jet, it was proposed that the column density of dust ejecta in the jet exceeded the density in the coma by a factor greater than or approximately equal to two, and that, if narrow along the line of light, the jet might have had a particle number concentration much higher than the coma background. Also noted are numerous ion features.

Larson, S. M.↗
NASA NTRS1985

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.↗
NASA NTRS1988

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.↗
NASA NTRS1988

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.↗
NASA NTRS1989

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.↗
NASA NTRS1990

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.↗
NASA NTRS1997

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↗
NASA NTRS2008

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↗
NASA NTRS2007

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.↗
NASA NTRS2015

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↗
NASA NTRS2016

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.↗
NASA NTRS2016

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↗
DOE OSTI2026

Endogenous Interface Pricing for Consistent Transmission–Distribution Co-Optimization With Discrete Distribution Controls

This paper proposes an endogenous interface pricing model for day-ahead transmission–distribution co-optimization that co-determines the interface locational marginal price (LMP) and the transmission–distribution exchange, ensuring price–dispatch consistency while optimally scheduling discrete distribution controls. The formulation couples a DC optimal power flow (OPF) with a branch-flow AC OPF that schedules distributed energy resources (DERs), tap-changer settings, capacitor banks (CBs), and multi-period energy storage systems (ESSs) under feeder voltage and current limits, and is solved as a mixed-integer second-order cone program (MISOCP). In a T14–D33 system, coordinated device scheduling recovers about 90% of the distribution-to-transmission export achievable in a reference case that ignores distribution network (DN) limits, while satisfying a 1.05 p.u. voltage upper bound. In a T39–D34/D37/D123 system, a sequential decoupled benchmark produces interface LMP distortions up to 12.5% and a 7.28% mismatch in net export energy, whereas the proposed model removes these distortions and the associated settlement mismatches. Second-order cone (SOC) relaxation gaps remain below $10^{-3}$ in all cases.

Noh, Seung-Gil↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
NASA NTRS1988

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↗