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 145 records · Page 8

Magnetopause shapes - General solutions

The problem of determining the shape of a magnetopause, namely, the tangential discontinuity separating (typically) the magnetic field of a source from external plasma, is reduced to solving an integral equation. No symmetry assumptions whatsoever are used in the derivation, consequently, realistic problems such as model shapes for the earth's magnetosphere (tilted dipole, flowing plasma, etc.) can in principle be treated. Moreover, the magnetic field source itself need not be a simple dipole. A variational method is suggested whereby the optimum parameters for any finite term trial shape are obtained.

Michel, F. C.↗

Correction Factor for Determining the London Penetration Depth from Strip Resonators

A significant disagreement is often seen between the theoretical temperature dependent magnetic penetration depth profile and experimentally derived calculations based on stripline type resonators. This short paper shows that the disagreement can be attributed to the susceptance coupled into the resonator from the gap discontinuity as well as the feed line. When the effect is taken into account, the natural resonant frequency of the resonator is increased, and the frequency shift due to kinetic inductance can be calculated much more accurately. While it is necessary to include this effect to determine the penetration depth, it is shown that the impact on unloaded quality factor is generally negligible. The situation when the strip characteristic impedance is not matched to the generator is included.

Romanofsky, Robert R.↗

A new family of stable elements for the Stokes problem based on a mixed Galerkin/least-squares finite element formulation

Adding to the classical Hellinger-Reissner formulation, a residual form of the equilibrium equation, a new Galerkin/least-squares finite element method is derived. It fits within the framework of a mixed finite element method and is stable for rather general combinations of stress and velocity interpolations, including equal-order discontinuous stress and continuous velocity interpolations which are unstable within the Galerkin approach. Error estimates are presented based on a generalization of the Babuska-Brezzi theory. Numerical results (not presented herein) have confirmed these estimates as well as the good accuracy and stability of the method.

Franca, Leopoldo P.↗

The Radiated Field Generated by a Monopole Source in a Short, Rigid, Rectangular Duct

This thesis develops a method to model the acoustic field generated by a monopole source placed in a moving rectangular duct. The walls of the duct are assumed to be infinitesimally thin and the source is placed at the center of the duct. The total acoustic pressure is written in terms of the free-space pressure, or incident pressure, and the scattered pressure. The scattered pressure is the augmentation to the incident pressure due to the presence of the duct. It satisfies a homogeneous wave equation and is discontinuous across the duct walls. Utilizing an integral representation of the scattered pressure, a set of singular boundary integral equations governing the unknown jump in scattered pressure is derived. This equation is solved by the method of collocation after representing the jump in pressure as a double series of shape functions. The solution obtained is then substituted back into the integral representation to determine the scattered pressure, and the total acoustic pressure at any point in the field. A few examples are included to illustrate the influence of various geometric and kinematic parameters on the radiated sound field.

Lakota, Barbara Anne↗

Congruence Approximations for Entrophy Endowed Hyperbolic Systems

Building upon the standard symmetrization theory for hyperbolic systems of conservation laws, congruence properties of the symmetrized system are explored. These congruence properties suggest variants of several stabilized numerical discretization procedures for hyperbolic equations (upwind finite-volume, Galerkin least-squares, discontinuous Galerkin) that benefit computationally from congruence approximation. Specifically, it becomes straightforward to construct the spatial discretization and Jacobian linearization for these schemes (given a small amount of derivative information) for possible use in Newton's method, discrete optimization, homotopy algorithms, etc. Some examples will be given for the compressible Euler equations and the nonrelativistic MHD equations using linear and quadratic spatial approximation.

Barth, Timothy J.↗

A Multi-scale Refined Zigzag Theory for Multilayered Composite and Sandwich Plates with Improved Transverse Shear Stresses

The Refined Zigzag Theory (RZT) enables accurate predictions of the in-plane displacements, strains, and stresses. The transverse shear stresses obtained from constitutive equations are layer-wise constant. Although these transverse shear stresses are generally accurate in the average, layer-wise sense, they are nevertheless discontinuous at layer interfaces, and thus they violate the requisite interlaminar continuity of transverse stresses. Recently, Tessler applied Reissner's mixed variational theorem and RZT kinematic assumptions to derive an accurate and efficient shear-deformation theory for homogeneous, laminated composite, and sandwich beams, called RZT(m), where "m" stands for "mixed". Herein, the RZT(m) for beams is extended to plate analysis, where two alternative assumptions for the transverse shear stresses field are examined: the first follows Tessler's formulation, whereas the second is based on Murakami's polynomial approach. Results for elasto-static simply supported and cantilever plates demonstrate that Tessler's formulation results in a powerful and efficient structural theory that is well-suited for the analysis of multilayered composite and sandwich panels.

Iurlaro, Luigi↗

Cloud Climate Data Records from Geostationary Satellites to Support CERES: Challenges and Recent Progress

Cloud properties are critical for understanding the Earth’s radiation budget and cloud feedbacks. At NASA Langley Research Center, the Satellite ClOud and Radiative Property retrieval System (SatCORPS) provides real-time and historical analyses of clouds derived from Geostationary satellite (GEOsat) data for weather and climate applications. For the Clouds and the Earth’s Radiant Energy System (CERES) program, the global constellation of GEOsats has been analyzed since 2000 to help characterize and account for the diurnal cycle of clouds and their radiative impacts in the CERES climate data record. Obtaining consistent cloud properties over the GEOsat data record at all times of day during the CERES era is a major objective but a significant challenge considering the diversity of imaging capabilities deployed during that time. A brief overview of the current CERES/GEO cloud property data record constructed from three generations of GEOsat imagers for CERES Edition-4 is presented. In this version, different algorithms were applied to different satellites to take advantage of as much spectral information as possible for each satellite. Under this approach, the derived cloud properties were found to become accurate (and more similar to MODIS) for the GEOsat imagers with more spectral information. However, this approach also led to some discontinuities in cloud properties across platforms. Significant differences were also found between daytime and nighttime cloud properties. This paper will highlight the major issues in the Edition-4 record and discuss new strategies and early results from the next version of the SatCORPS GEOsat cloud property time series in development for CERES Edition-5 with the goal of improving day-night and cross-platform consistency during the CERES record.

clouds↗

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

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗

Gas Density Discontinuities in Merging Clusters

In a number of galaxy clusters, Chandra has discovered a new phenomenon, the sharp gas density edges. Depending on the sign of the temperature jump across the edge, these features may either be bow shocks or cold fronts. While the merger origin of the bow shocks is unambiguous, what causes cold fronts is not entirely clear, as they are observed both in mergers and in relaxed clusters. The purpose of the study of A3376, an interesting cluster with density edges, is to understand the origin of cold fronts and to look for possible shocks. This cluster also has a radio halo, and comparison of an X-ray temperature map and radio image may shed light on the nature of the cluster radio halos, which is our secondary goal. The XMM data for A3376 have been obtained, processed and cleaned. Most of the preliminary analysis, including deriving a temperature map, is complete. We also have Chandra data for this cluster, and intend to compare and combine the two datasets. We are well underway toward determining the dynamical state of A3376 and the nature of its gas density discontinuities.

Mushotzky, Richard↗

On wind-type flows in astrophysical jets. III - Temporal evolution of perturbations and the formation of shocks

The temporal evolution of disturbances in a spherically symmetric polytropic wind from a central object is studied. Such disturbances may be due to localized momentum addition/subtraction, as, for example, by MHD waves, heating/cooling mechanisms in the outflow, or localized deviations from spherical symmetric expansion. The evolution of an initial perturbed state to a continuous or discontinuous final equilibrium state, as predicted by previous analytic calculations for stationary flows, is followed. It is shown that some of the predicted discontinuous equilibrium states are not physically accessible, while the attainment of the remaining equilibrium states depends on both the temporal and the spatial parameters characterizing the perturbation. The results are derived for solar conditions, but in fact can be applied to outflows in other astrophysical systems. In particular, applications to the solar wind and flows in astrophysical jets are discussed.

Trussoni, E.↗

Pressure regimes and core formation in the accreting earth

Recent work suggests that a large degree of melting is required to segregate metal from silicates, suggesting a connection with the formation of magma oceans. At low pressures metallic liquids do not wet silicate minerals, preventing the metal from aggregating into large masses that can sink. At high pressures, above 25 GPa, the dihedral angles of grains in contact with oxygen-rich metallic liquids may be reduced enough to allow percolation of metal, but this has not been confirmed. Physical models of core formation and accretion may therefore involve the formation of magma oceans and the segregation of metal at both high and low pressures. Models of core formation involving different pressure regimes are discussed as well as chemical evidence bearing on the models. Available geophysical data is ambiguous. The nature of the 670 km boundary (chemical difference or strictly phase change) between the upper and lower mantle is in doubt. There is some evidence that plumes are derived from the lower mantle, and seismic tomography strongly indicates that penetration of subducting oceanic crust into the lower mantle, but the tomography data also indicates that the 670 km discontinuity is a significant barrier to general mantle convection. The presence of the D' layer at the base of the lower mantle could be a reaction zone between the mantle and core indicating core-mantle disequilibrium, or D' layer could be subducted material. The abundance of the siderophile elements in the mantle could provide clues to the importance of high pressure processes in Earth, but partition coefficients at high pressures are only beginning to be measured.

Newsom, H. E.↗

[Cosmic Microwave Background (CMB) Anisotropies]

One of the main areas of research is the theory of cosmic microwave background (CMB) anisotropies and analysis of CMB data. Using the four year COBE data we were able to improve existing constraints on global shear and vorticity. We found that, in the flat case (which allows for greatest anisotropy), (omega/H)0 less than 10(exp -7), where omega is the vorticity and H is the Hubble constant. This is two orders of magnitude lower than the tightest, previous constraint. We have defined a new set of statistics which quantify the amount of non-Gaussianity in small field cosmic microwave background maps. By looking at the distribution of power around rings in Fourier space, and at the correlations between adjacent rings, one can identify non-Gaussian features which are masked by large scale Gaussian fluctuations. This may be particularly useful for identifying unresolved localized sources and line-like discontinuities. Levin and collaborators devised a method to determine the global geometry of the universe through observations of patterns in the hot and cold spots of the CMB. We have derived properties of the peaks (maxima) of the CMB anisotropies expected in flat and open CDM models. We represent results for angular resolutions ranging from 5 arcmin to 20 arcmin (antenna FWHM), scales that are relevant for the MAP and COBRA/SAMBA space missions and the ground-based interferometer. Results related to galaxy formation and evolution are also discussed.

Silk, Joseph↗

A Method of Efficient Inclination Changes for Low-thrust Spacecraft

The evolution of low-thrust propulsion technologies has reached a point where such systems have become an economical option for many space missions. The development of efficient, low trip time control laws has received an increasing amount of attention in recent years, though few studies have examined the subject of inclination changing maneuvers in detail. A method for performing economical inclination changes through the use of an efficiency factor is derived front Lagrange's planetary equations. The efficiency factor can be used to regulate propellant expenditure at the expense of trip time. Such a method can be used for discontinuous-thrust transfers that offer reduced propellant masses and trip-times in comparison to continuous thrust transfers, while utilizing thrusters that operate at a lower specific impulse. Performance comparisons of transfers utilizing this approach with continuous-thrust transfers are generated through trajectory simulation and are presented in this paper.

Falck, Robert↗

Analysis and closed-form solutions of circular and rectangular apertures in the transverse plane of a circular waveguide

An analysis is carried out to determine the discontinuity susceptance of a circular or rectangular aperture in the transverse plane of a circular waveguide. Closed-form solutions based on the scattered amplitude calculation and aperture coupling theory are derived. Example solutions for a small circular aperture and a resonant rectangular aperture are given. The calculated results agree very well with the experiments. The results have many applications in the design of circular waveguide components and circular waveguide-backed aperture antennas.

Eastham, Gary B.↗

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗

On empirical models of the upper atmosphere in the polar regions.

The expression for the exospheric temperature in Jacchia's static diffusion models of the upper atmosphere has a discontinuous gradient at the poles. Therefore it cannot describe the true state of the upper atmosphere in the polar regions. Furthermore, it cannot be used to calculate quantities that depend on the derivative of the exospheric temperature, or the density, like pressure gradients, horizontal forces or horizontal heat fluxes. A modified expression for the exospheric temperature is suggested. This modification yields variables of state of the upper atmosphere that deviate little from Jacchia's values, but it has continuous gradient at the poles and is therefore more suitable for treating dynamical problems like the global wind pattern.

Blum, P. W.↗

Shape Measurement of Large Aerospace Structures Using Sensitivity Electrical TDR Distributed Strain Sensor

Electrical time domain reflectometry (ETDR) sensing technique can be best described as "closed-loop radar," where the information is derived from the reflections of a voltage pulse sent through a transmission medium. The ETDR sensing technique is a well-developed method and has been widely used to locate and evaluate discontinuities in long coaxial power transmission cables. The ETDR technique provides a true distributed sensing capability which can not only sense the distributed loading condition of the structure but also can pin-point the location of disturbance, such as the locations of stress concentration and structural damages. Proof-of-concept experiments have been conducted using photoelastic specimens with embedded commercial coaxial cables, i.e., RG85/U and RG174, to demonstrate the stress/strain sensing capability of ETDR sensors for structural health monitoring application. Although the test results showed that the ETDR sensor signals capture specimen deformation pattern both in bending and tension and indicate the location and type of crack damages of the photoelastic specimen; yet, the low signal-to-noise ratio of the sensor signal smears the details of the strain measurement that the ETDR signals can convey. A high-sensitivity ETDR coaxial strain sensor prototype newly developed at Clark Atlanta University will be presented. The construction of the prototype sensing cable as well as its electrical properties relevant to distributed strain sensing application will be shown in details. Test results of the sensitivity and tension responses of the ETDR signal of the prototype sensor will be presented and compared with those of commercial coaxial cables. Promising potentials of the ETDR distributed strain sensing method for shape measurement application of large aerospace structures will also be demonstrated using long slender beam with surface-bonded ETDR distributed strain sensor.

Lin, Mark W.↗