Search NASA⌕ Search

SEARCH · Search NASA

Results for “derivative discontinuity”

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 181 records · Page 10

Strongly coupled stress waves in heterogeneous plates.

Consideration of coupled stress waves generated by an impulsive load applied at one end of a semiinfinite plate. For the field equations governing the one-dimensional coupled waves a hyperbolic system of equations is obtained in which a strong coupling in the second derivatives exists. The method of characteristics described by Chou and Mortimer (1967) is extended to cover the case of strong coupling, and a study is made of the transient stress waves in a semiinfinite plate subjected to an initial step input. Coupled discontinuity fronts are found to propagate at different velocities. The normal plate stress and the bending moment at different time regimes are illustrated by graphs.

Wang, A. S. D.↗

Context dependent anti-aliasing image reconstruction

Image Reconstruction has been mostly confined to context free linear processes; the traditional continuum interpretation of digital array data uses a linear interpolator with or without an enhancement filter. Here, anti-aliasing context dependent interpretation techniques are investigated for image reconstruction. Pattern classification is applied to each neighborhood to assign it a context class; a different interpolation/filter is applied to neighborhoods of differing context. It is shown how the context dependent interpolation is computed through ensemble average statistics using high resolution training imagery from which the lower resolution image array data is obtained (simulation). A quadratic least squares (LS) context-free image quality model is described from which the context dependent interpolation coefficients are derived. It is shown how ensembles of high-resolution images can be used to capture the a priori special character of different context classes. As a consequence, a priori information such as the translational invariance of edges along the edge direction, edge discontinuity, and the character of corners is captured and can be used to interpret image array data with greater spatial resolution than would be expected by the Nyquist limit. A Gibb-like artifact associated with this super-resolution is discussed. More realistic context dependent image quality models are needed and a suggestion is made for using a quality model which now is finding application in data compression.

Beaudet, Paul R.↗

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗

A new Lagrangian method for three-dimensional steady supersonic flows

In this report, the new Lagrangian method introduced by Loh and Hui is extended for three-dimensional, steady supersonic flow computation. The derivation of the conservation form and the solution of the local Riemann solver using the Godunov and the high-resolution TVD (total variation diminished) scheme is presented. This new approach is accurate and robust, capable of handling complicated geometry and interactions between discontinuous waves. Test problems show that the extended Lagrangian method retains all the advantages of the two-dimensional method (e.g., crisp resolution of a slip-surface (contact discontinuity) and automatic grid generation). In this report, we also suggest a novel three dimensional Riemann problem in which interesting and intricate flow features are present.

Loh, Ching-Yuen↗

Trapped particles at a magnetic discontinuity

At a tangential discontinuity between two constant magnetic fields a layer of trapped particles can exist, this work examines the conditions under which the current carried by such particles tends to maintain the discontinuity. Three cases are examined. If the discontinuity separates aligned vacuum fields, the only requirement is that they be antiparallel. With arbitrary relative orientations, the field must have equal intensities on both sides. Finally, with a guiding center plasma on both sides, the condition reduces to a relation which is also derivable from hydromagnetic theory. Arguments are presented for the occurrence of such trapped modes in the magnetopause and for the non-existence of specular particle reflection.

Stern, D. P.↗

Surface fitting three-dimensional bodies

The geometry of general three-dimensional bodies is generated from coordinates of points in several cross sections. Since these points may not be smooth, they are divided into segments and general conic sections are curve fit in a least-squares sense to each segment of a cross section. The conic sections are then blended in the longitudinal direction by fitting parametric cubic-spline curves through coordinate points which define the conic sections in the cross-sectional planes. Both the cross-sectional and longitudinal curves may be modified by specifying particular segments as straight lines and slopes at selected points. Slopes may be continuous or discontinuous and finite or infinite. After a satisfactory surface fit has been obtained, cards may be punched with the data necessary to form a geometry subroutine package for use in other computer programs. At any position on the body, coordinates, slopes and second partial derivatives are calculated. The method is applied to a blunted 70 deg delta wing, and it was found to generate the geometry very well.

Dejarnette, F. R.↗

An IUE survey of interstellar H I Ly alpha absorption. 1: Column densities

We measure Galactic interstellar neutral hydrogen column densities by analyzing archival interstellar Ly alpha absorption line data toward 554 B2 and hotter stars observed at high resolution with the IUE satellite. This study more than doubles the number of lines of sight with measures of N(H I) based on Ly alpha. We have included the scattered light background correction algorithm of Bianchi and Bohlin in our data reduction. We use the correlation between the Balmer discontinuity (c sub 1) index and the stellar Ly alpha absorption in order to assess the effects of stellar Ly alpha contamination. Approximately 40% of the B stars with measured (c sub 1) index, exhibit serious stellar Ly alpha contamination. One table contains the derived values of the interstellar N(H I) for 393 stars with at most small amounts of stellar contamination. Another lists the observed values of total N(H I) for 161 stars with suspected stellar Ly alpha contamination and/or uncertain stellar parameters.

Diplas, Athanassios↗

Practical formulation of a positively conservative scheme

Approximate Riemann solvers have been highly successful for computing the Euler/Navier-Stokes equations, but linearized Riemann solvers are known to fail occasionally by predicting non-physical states with negative density or internal energy. Positively conservative schemes, in contrast, guarantee physical solutions from realistic input. The Harten-Lax-van Leer-Einfeldt (HLLE) scheme is a typical example of a positively conservative scheme. However, the HLLE scheme is highly dissipative at contact discontinuities and shear layers and thus it is not applicable to practicle simulations. An existing modification to the HLLE scheme, known as HLLEM, enhances the resolution to that of the Roe scheme. However, this modification violates the positivity of density and internal energy. Precise derivation of the modification yields a quatratic inequality and thus requires a case-by-case treatment. This Note describes a new, modified HLLE scheme that satisfies the positively conservative condition approximately. Sample computationa are included to demonstrate the resolution and the robustness of the scheme.

Obayashi, Shigeru↗

Modified Geometric Truncation of the Scattering Phase Function

Phase function of light scattering on large atmospheric particles has very strong peak in forward direction constituting a challenge for accurate numerical calculations of radiance required in remote sensing problems. Scaling transformation replaces original phase function with a sum of the delta function and a new regular smooth phase function. Geometric truncation is one of the ways to construct such a smooth function. The replacement phase function coincides with the original one outside the forward cone and preserves the asymmetry parameter. It has discontinuity at the cone. Another simple functional form of the replacement phase function within the cone is suggested. It enables continuity and allows for a number of modifications. Three of them are considered in this study: preserving asymmetry parameter, providing continuity of the 1st derivative of the phase function, and preserving mean scattering angle. Yet another problem addressed in this study is objective selection of the width of the forward cone. That angle affects truncation fraction and values of the phase function within the cone. A heuristic approach providing unambiguous criterion of selection of the truncation angle is proposed. The approach has easy numerical implementation. Suggested modifications were tested on cloud phase function using discrete ordinates and Monte Carlo methods. It was shown that the modifications provide better accuracy of the radiance computation compare to the original geometric truncation with discrete ordinates while continuous derivative approach provides significant gain in computer time with Monte Carlo simulations.

Alexander Radkevich↗

Modeling of switching regulator power stages with and without zero-inductor-current dwell time

State space techniques are employed to derive accurate models for buck, boost, and buck/boost converter power stages operating with and without zero-inductor-current dwell time. A generalized procedure is developed which treats the continuous-inductor-current mode without the dwell time as a special case of the discontinuous-current mode, when the dwell time vanishes. An abrupt change of system behavior including a reduction of the system order when the dwell time appears is shown both analytically and experimentally.

Lee, F. C.↗

Luna 24 - Opaque mineral chemistry of gabbroic and basaltic fragments from Mare Crisium

Spinels and ilmenites are relatively sparse in the Luna 24 gabbro and basalts. Spinel compositions show some affinities to those of spinels in Apollo 12, Apollo 14 and Luna 16 basalts; a characteristic feature is high Al2O3, reaching a maximum of 19.8 wt%. A comparison of spinels in the Luna 24 gabbro with those in other deep-seated lunar intrusive rocks shows a characteristic trend for Fe/Mg. This trend is systematic from gabbro to anorthosite to troctolite and is interpreted to be P-T dependent. Compositions of spinels in the gabbro fall within the Cr/Al trend defined by the spinels of the basalts, but form a Fe/Mg trend parallel to that of the basalts; this relationship suggests that both the gabbro and the basalts are derived from a closely similar source region, with the basalts originating at a slightly greater depth than the gabbro. The spinels in both rock types are considered to have formed at high crustal levels, at low pressures. The Luna 24 data suggest that the compositional discontinuities which exist between chromian spinels and titanian spinels in a large proportion of mare basalts are the result of nucleation of chromian spinels at high crustal levels prior to eruption, and of titanian spinels during melt crystallization at the lunar surface.

Haggerty, S. E.↗

An extended classical solution of the droplet growth problem

Problems of applying the classical kinetic theory to the growth of small droplets from vapor are examined. A solution for the droplet growth equation is derived which is based on the assumption of a diffusive field extending to the drop surface. The method accounts for partial thermal and mass accommodation at the interface and the kinetic limit to the mass and heat fluxes, and it avoids introducing the artifact of a discontinuity in the thermal and vapor field near the droplet. Consideration of the environmental fields in spherical geometry utilizing directional fluxes yields boundary values in terms of known parameters and a new Laplace transform integral.

Anderson, B. J.↗

The Effects of a 1998 Observing System Change on MERRA-2-Based Ozone Profile Simulations

Model simulations of ozone (O3) driven by meteorological reanalyses are useful for filling observational gaps and interpreting observed O3 variability and trends. However, the transport circulation of reanalysis products is impacted by changes to the observing system (the data assimilated into the reanalyses). We examine the impacts of these changes on simulated O3 from two models, Global Modeling Initiative (GMI) Chemistry Transport Model (GMI CTM) and Modern‐Era Retrospective Analysis for Research Applications version 2 (MERRA‐2) GMI Replay (M2 GMI Replay) simulation, using observations from global ozonesondes (>50,000 profiles) and satellites from 1980 to 2016. Both models are constrained by meteorology from the NASA MERRA‐2 reanalysis, and both use versions of NASA's GMI chemical mechanism. We focus on an observing system change affecting simulated O3 after 1998, associated with the assimilation of temperature and humidity data from new microwave profiling satellites. A large post‐1998 O3 increase, mainly confined to 15–20‐km altitude, of ~10 Dobson units (DU) in midlatitudes occurs in the GMI CTM, worsening the bias compared to observations. In contrast, an increase in M2 GMI Replay simulation O3 of ~10 DU is observed only near −60° latitude, reducing the bias compared to observations. The GMI CTM O3 high biases display a Quasi‐Biennial Oscillation (QBO)‐like periodicity that result from excessive transport from the tropical stratosphere to the midlatitude lower stratosphere during the QBO westerly phase. We quantify O3 discontinuities caused by MERRA‐2 observing system changes and demonstrate how the MERRA‐2 Global Modeling Initiative Replay simulation dampens the effects of these changes and QBO‐driven artifacts on simulated lower stratospheric and total O3. We caution against using simulations driven by a reanalysis to derive multidecadal O3 trends, especially prior to 1998.

Replay↗

DORIS Satellite Phase Center Determination and Consequences on the Derived Scale of the Terrestrial Reference Frame

This slide presentation reports on the analysis of several years of Delft Object-oriented Radar Interferometric Software (Doris) data undertaken to estimate daily determination of satellite antenna phase center corrections. This was done on a satellite by satellite basis. For each DORIS satellite, we considered long-term time series of such individual estimations, looking for possible biases, discontinuities, trends or annual signals. The analysis compared DORIS to GPS estimates for common satellites (Jason and TOPEX/Poseidon).

terrestrial reference frame scale↗

Flux splitting algorithms for two-dimensional viscous flows with finite-rate chemistry

The Roe flux-difference splitting method has been extended to treat two-dimensional viscous flows with nonequilibrium chemistry. The derivations have avoided unnecessary assumptions or approximations. For spatial discretization, the second-order Roe upwind differencing is used for the convective terms and central differencing for the viscous terms. An upwind-based TVD scheme is applied to eliminate oscillations and obtain a sharp representation of discontinuities. A two-stage Runge-Kutta method is used to time integrate the discretized Navier-Stokes and species transport equations for the asymptotic steady solutions. The present method is then applied to two types of flows: the shock wave/boundary layer interaction problems and the jet in cross flows.

Shuen, Jian-Shun↗

Flux splitting algorithms for two-dimensional viscous flows with finite-rate chemistry

The Roe flux difference splitting method was extended to treat 2-D viscous flows with nonequilibrium chemistry. The derivations have avoided unnecessary assumptions or approximations. For spatial discretization, the second-order Roe upwind differencing is used for the convective terms and central differencing for the viscous terms. An upwind-based TVD scheme is applied to eliminate oscillations and obtain a sharp representation of discontinuities. A two-state Runge-Kutta method is used to time integrate the discretized Navier-Stokes and species transport equations for the asymptotic steady solutions. The present method is then applied to two types of flows: the shock wave/boundary layer interaction problems and the jet in cross flows.

Shuen, Jian-Shun↗

The Stress-strain Behavior of Polymer-Nanotube Composites from Molecular Dynamics Simulations

Stress-strain curves of polymer-carbon nanotube composites are derived from molecular dynamics simulations of a single-walled carbon nanotube embedded in polyethylene. A comparison is made between the response to mechanical loading of a composite with a long, continuous nanotube (replicated via periodic boundary conditions) and the response of a composite with a short, discontinuous nanotube. Both composites are mechanically loaded in the direction of and transverse to the NT axis. The long-nanotube composite shows an increase in the stiffness relative to the polymer and behaves anisotropically under the different loading conditions. The short-nanotube composite shows no enhancement relative to the polymer, most probably because of its low aspect ratio. The stress-strain curves are compared with rule-of-mixtures predictions.

Frankland, S. J. V.↗

Capillary surface discontinuities above reentrant corners

A particular configuration of a vertical capillary tube for which S is the equilibrium interface between two fluids in the presence of a downward pointing gravitational field was investigated. S is the graph a function u whose domain is the (horizontal) cross section gamma of the tube. The mean curvature of S is proportional to its height above a fixed reference plane and lambda is a prescribed constant and may be taken between zero and pi/2. Domains gamma for which us is a bounded function but does not extend continuously to d gamma are sought. Simple domains are found and the behavior of u in those domains is studied. An important comparison principle that has been used in the literature to derive many of the results in capillarity is reviewed. It allows one to deduce the approximate shape of a capillary surface by constructing comparison surfaces with mean curvature and contact angle close to those of the (unknown) solution surface. In the context of nonparametric problems the comparison principle leads to height estimates above and below for the function u. An example from the literature where these height estimates have been used successfully is described. The promised domains for which the bounded u does not extend continuously to the boundary are constructed. The point on the boundary at which u has a jump discontinuity will be the vertext of a re-entrant corner having any interior angle theta pi. Using the comparison principle the behavior of u near this point is studied.

Korevaar, H. J.↗