Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 325 records · Page 18

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.↗

A search for periodicity in the x ray spectrum of black hole candidate A0620-00

The archived data from the SAS-3 observations of the X-ray nova A0620-00, the best of the stellar blackhole candidates, were exhaustively examined for evidence of variable phenomena correlated with the orbital motion of the binary system of which it is a member. The original analysis of these data was completed before discovery of the binary companion and determination of the orbital period of the system. New interest was drawn to the task of a reexamination of the archive data by the recent discovery of the massive nature of the X-ray source through analysis of the Doppler variations and ellipsoidal light variations of the faint K-star companion by McClintock and Remillard. The archive research, carried out under the supervision of the principal investigator, was the topic of the thesis submitted to the MIT Department of Physics by Kenneth Plaks in partial fulfillment of the requirements for the degree of Master of Science. Plaks' effort was focused on the elimination of fluctuations in the data due to errors in attitude solutions and other extraneous causes. The first products of his work were long-term light curves of the X-ray intensities in the various energy channels as functions of time during the time from outbursts in August 1975 to quiescence approximately 6 months later. These curves, are refined versions of the preliminary results published in 1976 (Matilsky et al. 1976). Smooth exponentials were fitted to these long term light curves to provide the basis for detrending the data, thereby permitting a calculation of residuals derived by subtracting the fitted curve from the data. The residuals were then analyzed by Fourier analysis to search for variations with the period of the binary orbit, namely 7.75 hours. No evidence of an orbital periodicity was found. However, the refined light curve provides a much clearer picture of the outburst and subsequent decay of the X-ray luminosity. In fact, there were two outbursts, each followed by an exponential decay with similar time constants of about 25 days. Previous evidence of a three-oscillation variation with a 7.8 day period were confirmed. Substantial theoretical effort has been devoted to attempts to account for the decay characteristics as the result of the gradual eating up of an accretion disk by a stellar-mass blackhole (e.g., Huang and Wheeler 1989). The improved decay curves will provide significant new constraints on the theoretical analyses.

Clark, George W.↗

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.↗

A multigrid algorithm for the cell-centered finite difference scheme

In this article, we discuss a non-variational V-cycle multigrid algorithm based on the cell-centered finite difference scheme for solving a second-order elliptic problem with discontinuous coefficients. Due to the poor approximation property of piecewise constant spaces and the non-variational nature of our scheme, one step of symmetric linear smoothing in our V-cycle multigrid scheme may fail to be a contraction. Again, because of the simple structure of the piecewise constant spaces, prolongation and restriction are trivial; we save significant computation time with very promising computational results.

Ewing, Richard E.↗

Extreme ultraviolet spectroscopy of PKS 2155-304

We present the extreme ultraviolet (75-110 A) spectrum of the BL Lac object PKS 2155-304, the first spectrum of an extragalactic source obtained with the Extreme Ultraviolet Explorer (EUVE). The spectrum shows a generally smooth continuum, which can be modeled by a single power law plus interstellar absorption, and possibly an absorption feature at approximately 80 A. The best fit to the data suggests that the EUV spectrum can be interpreted as a simple extrapolation of the X-ray continuum, with an energy index alpha is approximately = 1.6; however, shallower or steeper power laws with indices between -0.4 and 2.7 cannot be ruled out by the existing EUV data alone. The data provide strong constraints on the interstellar neutral H and He along the line of sight. Using a column density of 1.36 x 10(exp 20)/sq cm for the Galactic neutral hydrogen along the PKS 2155-304 line of sight, the neutral helium column density is constrained to be 9%-10% of the hydrogen amount.

Fruscione, Antonella↗

An experimental study of changes in the structure of a turbulent boundary layer due to surface geometry changes

This investigation attempts to observe changes in the flow structure of a nominally zero pressure gradient turbulent boundary layer by altering the surface of a smooth plate with small triangular riblets oriented in the direction of flow and protruding just above the viscous sublayer. Hot-wire measurements show a reduction in the local skin friction over the riblet valleys by about 40% compared to the smooth surface near the end of the plate. Local skin friction increases by about 10% over the riblet peaks. Results showing the effects of the riblets on the mean, root-mean-square and skewness factor distributions of the streamwise velocity are presented. The mean profile displays characteristics similar to those of drag reducing additives; the streamwise turbulence intensity is reduced by approximately 10-15% near the surface by the riblets; the skewness factor distribution is virtually unchanged. A comparison of results from an application of the VITA 'burst' detection technique to smooth and riblet surface data is also presented although the results are rather inconclusive. Flow visualization of 'bursting' over both surfaces, while dramatically revealing the structure, reveals little difference resulting from the presence of the riblet.

Source record↗

The far-infrared spectrum of the core of Sagittarius B2

A Michelson interferometer aboard NASA's Kuiper Airborne Observatory has been used to measure the spectrum of Sgr B2 from 40 to 200 kaysers with 5-kayser resolution in a 1.4-arcmin beam. The measured spectrum is smooth and featureless with a broad maximum at about 85 kaysers. The data can be fitted analytically with a model corresponding to thermal emission by a uniform sla of dust filling the beam, with an average temperature of approximately 32 K, an optical depth at 100 microns of about 1.6, and a spectral index of the dust emissivity about 1.5. The absence of features implies either that the source is optically thick or that the emission spectrum of the individual grains is smooth in the passband. The possible physical significance of this model is discussed.

Erickson, E. F.↗

A Colebrook equation for impinging radial wall jets

Here, in this study, we evaluate the skin coefficient of friction for steady turbulent radial wall jets across smooth and rough surfaces. Although the Colebrook equation has been used successfully for many decades to evaluate friction factors for flows through smooth and rough pipes, how roughness affects the skin friction coefficient for steady turbulent radial wall jets remains unclear. Here we explore a Colebrook-type equation for skin friction coefficients associated with single-phase turbulent radial wall jets arising from orthogonally impinging circular jets. The fully iterative solution, based on well-established concepts of turbulent wall-bounded flow, is presented along with a power-law approximation and a non-iterative approximation for the friction coefficient derived therefrom. We find the skin coefficient of friction defined on the peak radial velocity to be a function of position over rough but not smooth surfaces in contrast to pipe friction factors that remain independent of axial position. These results follow expected trends, explain prior heterogeneity in power-law expressions for the skin friction coefficient and have significant implications for the industrial use of jets in mixing vessels.

friction losses↗

The early light curve of A0620-00

The Ariel-5 All-Sky Monitor measured the 3-6 keV X-ray intensity of 0620-00 for two days shortly after peak emission in August 1975, and continuously throughout September 1975. The effective exposure each day for this source (and every other source in the 80% of the celestial sphere covered by the monitor) is approximately 250 sq cm sec. The light curve obtained through October 1 (when the spin axis was again pointed to A0620-00) is shown. The intensity difference at maximum is obviously a manifestation of the very soft spectrum of the source. Clearly, A0620-00 was approximately four times as bright as Sco X-1 at maximum in the band 3-6 keV. The decay is quite smooth, but cannot be fit with a single e-folding time. The interval between the early All-Sky Monitor points and the onset of continuous coverage has an inferred e-folding time of 22 days, but it is continually increasing throughout September. A0620-00 is apparently similar in its X-ray character to the very strong, long-lasting transient X-ray sources which presently number six. In contrast, only three can be sensibly reconciled with the lower-intensity, shorter-duration hard-spectrum transients which have been found to pulse on a time scale of minutes.

Kaluzienski, L. J.↗

Cyclic plasticity and failure of structural components

An analytical approach for low-cycle fatigue prediction is presented. The approach combines a cyclic plasticity model with the finite element method and a damage accumulation criterion for ductile metals. The cyclic plasticity model is based on the concept of the combination of several yield surfaces. The surfaces are related to the material uniaxial stress-strain curve idealized by piecewise linear segments. The damage criterion is based on the Coffin-Manson formulae modified for the mean stress variation effect. It is extended to the multiaxial varying stress-strain field and applied for both the crack initiation and the crack growth processes. The stable slow crack growth rate is approximated by the damage accumulation gradient computed from the cracked finite element models. This procedure requires fatigue testing data of only smooth specimens under constant strain amplitudes. The present approach is illustrated by numerical examples of an aircraft wing stiffened panel subjected to compression, which causes material yielding and residual tension.

Kalev, I.↗

Ray tracing of broadband bursty radio emissions from Uranus

To determine the source position of the broadband bursty emission, rays of X-mode emissin were traced from source positions along magnetic field lines with footprints that form a large grid centered approximately on the south magnetic pole of Uranus. For large wave normal angles, source regions different from those producing b-smooth emission were found. The emission observed prior to closest approach has a source along field lines that are distinct from those which generate emissions observed after closest approach.

Curran, D. B.↗

A Hubble Space Telescope ultraviolet spectrum of SN 1993J

We obtained a Hubble Space Telescope (HST) ultraviolet spectrum of the Type II supernova SN 1993J in M81 on 1993 April 15. The approximately 1650-2900 A region is smoother than observed for SN 1987A and SN 1992A and lacks strong P Cygni absorptions caused by iron peak element lines. Synthetic spectra calculated using a parameterized Local Thermodynamic Equilibrium (LTE) procedure and a simple model atmosphere do not fit the UV region. Radio observations suggest that SN 1993J is embedded in a thick circumstellar envelope. The UV spectra of other supernovae that are believed to have thick circumstellar envelopes also have approximately 1650-2900 A regions lacking in strong P Cygni absorptions. Interaction of supernova ejecta and circumstellar matter may cause the smooth UV spectrum. If so, UV observations of supernovae will provide insight into the circumstellar environment of the supernova progenitors.

Jeffery, David J.↗

Hydrodynamic and Aerodynamic Characteristics of a Model of a Supersonic Multijet Water-Based Aircraft Equipped with Supercavitating Hydrofoils

The hydrodynamic and aerodynamic characteristics of a model of a multijet water-based Mach 2.0 aircraft equipped with hydrofoils have been determined. Takeoff stability and spray characteristics were very good, and sufficient excess thrust was available for takeoff in approximately 32 seconds and 4,700 feet at a gross weight of 225,000 pounds. Longitudinal and lateral stability during smooth-water landings were good. Lateral stability was good during rough-water landings, but forward location of the hydrofoils or added pitch damping was required to prevent diving. Hydrofoils were found to increase the aerodynamic lift-curve slope and to increase the aerodynamic drag coefficient in the transonic speed range, and the maximum lift-drag ratio decreased from 7.6 to 7.2 at the cruise Mach number of 0.9. The hydrofoils provided an increment of positive pitching moment over the Mach number range of the tests (0.6 to 1.42) and reduced the effective dihedral and directional stability.

McKann, Robert E.↗

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme.

Solomonoff, A.↗

Thermally generated magnetic fields in laser-driven compressions and explosions

The evolution of thermally generated magnetic fields in a plasma undergoing a nearly spherically symmetric adiabatic compression or expansion is calculated. The analysis is applied to obtain approximate results for the development of magnetic fields in laser-driven compression and explosion of a pellet of nuclear fuel. Localized sources, such as those occurring at composition boundaries in structured pellets or at shock fronts, give stronger fields than those deriving from smoothly distributed asymmetries. Although these fields may approach 10 million G in the late stages of compression, this is not expected to present difficulties for the compression process. Assuming ignition of a nuclear explosion occurs, the sources become much stronger, and values of approximately 10 billion G are obtained at tamper boundaries assuming a 20% departure from spherical symmetry during the explosion.

Tidman, D. A.↗

Nuclear cusps and singularities in the nonadditive kinetic potential bifunctional from analytical inversion

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗

Physically interpretable approximations of many-body spectral functions

The rational function approximation provides a natural and interpretable representation of response functions such as the many-body spectral functions. We apply the vector fitting (VFIT) algorithm to fit a variety of spectral functions calculated from the Holstein model of electron-phonon interactions. We show that the resulting rational functions are highly efficient in their fitting of sharp features in the spectral functions, and could provide a means to infer physically relevant information from a spectral data set. The position of the peaks in the approximated spectral function are determined by the location of poles in the complex plane. Additionally, we developed a variant of VFIT that incorporates regularization to improve the quality of fits. With this procedure, we demonstrate it is possible to achieve accurate spectral function fits that vary smoothly as a function of physical conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GradeIT (Road Grade Inference Tool)

The National Renewable Energy Laboratory (NREL) has developed the Road Grade Inference Tool, "GradeIT". This tool is a modular Python package the leverages the United States Geological Survey’s (USGS) 1/3 arc-second Digital Elevation Model (DEM), an open and free public dataset. Roadway gradient (slope) has a strong influence on vehicle energy consumption behavior and longitudinal dynamics. Therefore, any tools, models, or analyses that consider vehicle energy consumption or driving profiles is incomplete without accurate elevation and grade information. Despite the acute need for this data, it is rarely collected alongside GPS (latitude and longitude) data, and tools to post-process GPS data to append elevation and grade are difficult to come by. GradeIT satisfies this significant need in the research community. Utilization of the USGS DEM makes the underlying data available to any potential user and ensures continuous data anywhere in the continental United States. The USGS DEM is a great resource; however, in order to be useful for automotive and traffic engineering applications, significant interpolation and filtering is required to achieve an accurate representation of gradient on the roadway. GradeIT employs a two-dimensional interpolation to approximate elevation at a point on the road from the nearest nodes in the DEM. Then, GradeIT applies multiple filtering techniques to smooth the elevation profile to a realistic road topology. Finally, road grade is derived from the elevation profile and returned as the primary output.

Holden, Jacob↗