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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

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

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

Learning protocols for the fast and efficient control of active matter

Exact analytic calculation shows that optimal control protocols for passive molecular systems often involve rapid variations and discontinuities. However, similar analytic baselines are not generally available for active-matter systems, because it is more difficult to treat active systems exactly. Here we use machine learning to derive efficient control protocols for active-matter systems, and find that they are characterized by sharp features similar to those seen in passive systems. We show that it is possible to learn protocols that effect fast and efficient state-to-state transformations in simulation models of active particles by encoding the protocol in the form of a neural network. We use evolutionary methods to identify protocols that take active particles from one steady state to another, as quickly as possible or with as little energy expended as possible. Our results show that protocols identified by a flexible neural-network ansatz, which allows the optimization of multiple control parameters and the emergence of sharp features, are more efficient than protocols derived recently by constrained analytical methods. Our learning scheme is straightforward to use in experiment, suggesting a way of designing protocols for the efficient manipulation of active matter in the laboratory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence

Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

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

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

[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↗
NASA NTRS2002

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

Hillslope-Channel Transitions and the Role of Water Tracks in a Changing Permafrost Landscape

The Arctic is experiencing rapid climate change, and the effect on hydrologic processes and resulting geomorphic changes to hillslopes and channels is unclear because we lack quantitative models and theory for rapid changes resulting from thawing permafrost. Here, the presence of permafrost modulates water flow and the stability of soil-mantled slopes, implying that there should be a signature of permafrost processes, including warming-driven disturbance, in channel network extent. To inform understanding of hillslope-channel dynamics under changing climates, we examined soil-mantled hillslopes within a ~300 km 2 area of the Seward Peninsula, western Alaska, where discontinuous permafrost is particularly susceptible to thaw and rapid landscape change. In this study, we pair high-resolution topographic and satellite data to multi-annual observations of InSAR-derived surface displacement over a 5-year period to quantify spatial variations in topographic change across an upland landscape. We find that neither the basin slope nor the presence of knickzones controls the magnitude of recent surface displacements within the study basin, as may be expected under conceptual models of temperate hillslope evolution. Rather, the highest displacement magnitudes tended to occur at the broad hillslope-channel transition zone. In this study area, this zone is occupied by water tracks, which are zero-order ecogeomorphic features that concentrate surface and subsurface flow paths. Our results suggest that water tracks, which appear to occupy hillslope positions between saturation and incision thresholds, are vulnerable to warming-induced subsidence and incision. We hypothesize that gullying within water tracks will outpace infilling by hillslope processes, resulting in the growth of the channel network under future warming.

58 GEOSCIENCES↗
NASA NTRS1991

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

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

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

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

Unsteady turbulent shear flow in shock tube discontinuities

A pressure-ruptured shock tube and an arc driven shock tube, have been used to study the evolution of turbulent fluctuations at contact surfaces with N2O4-2NO2 mixtures and at ionizing shock fronts in argon. The study has focused on point density diagnostics derived from crossed light beam correlations and electric probes. Turbulent bursts are found for which dynamical and spectral analyses suggest a particle-like evolution of fluctuation segments with a unique and characteristic frequency, independent of flow history and overall flow conditions.

Johnson, J. A., III↗
DOE OSTI2021

High-throughput screening of multifunctional nanocoatings based on combinations of polyphenols and catecholamines

Biomimetic surface coatings based on plant polyphenols and catecholamines have been used broadly in a variety of applications. However, the lack of a rational cost-effective platform for screening these coatings and their properties limits the true potential of these functional materials to be unleashed. Here, we investigated the oxidation behavior and coating formation ability of a library consisting of 45 phenolic compounds and catecholamines. UV–vis spectroscopy demonstrated significant acceleration of oxidation and polymerization under UV irradiation. We discovered that several binary mixtures resulted in non-additive behavior (synergistic or antagonistic effect) yielding much thicker or thinner coatings than individual compounds measured by ellipsometry. To investigate the properties of coatings derived from new combinations, we used a miniaturized high-throughput strategy to screen 2,532 spots coated with single, binary, and ternary combinations of coating precursors in one run. We evaluated the use of machine learning models to learn the relation between the chemical structure of the precursors and the thickness of the nanocoatings. Formation and stability of nanocoatings were investigated in a high-throughput manner via discontinuous dewetting. 30 stable combinations (hits) were used to tune the surface wettability and to form water droplet microarray and spot size gradients of water droplets on the coated surface. No toxicity was observed against eukaryotic HeLa cells and Pseudomonas aeruginosa (strain PA30) bacteria after 24 h incubation at 37 C. The strategy introduced here for high-throughput screening of nanocoatings derived from combinations of coating precursors enables the discovery of new functional materials for various applications in science and technology in a cost-effective miniaturized manner.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Factors Controlling a Synthetic Aperture Radar (SAR) Derived Root-Zone Soil Moisture Product over The Seward Peninsula of Alaska

Root-zone soil moisture exerts a fundamental control on vegetation, energy balance, and the carbon cycle in Arctic ecosystems, but it is still not well understood in vast, remote, and understudied regions of discontinuous permafrost. The root-zone soil moisture product (30 m resolution) used in this analysis was retrieved from a time-series P-Band (420–440 MHz) synthetic aperture radar (SAR) backscatter observations (August 2017 & October 2017). While similar approaches have been taken to retrieve surface (0 cm to 5 cm) soil moisture from L-Band (1.2 GHz) SAR backscatter, this is one of the first known attempts at reaching the root-zone in permafrost regions. Here, we analyze secondary factors (excluding primary factors, such as precipitation) controlling summer (August) soil moisture at depths of 6 cm, 12 cm, and 20 cm over a 4500 km2 area on the Seward Peninsula of Alaska. Using a random forest model, we quantify the impact of topography, vegetation, and meteorological factors on soil moisture distributions. In developing the random forest model, we explore a variety of feature scales (30 m, 60 m, 90 m, 120 m, 180 m, and 240 m), tune hyperparameters (the structure of individual decision trees making up the ensemble including the number and depth of trees), and perform the final feature selection using cross-validated recursive feature elimination. Results suggest that root-zone soil moisture on the Seward Peninsula is primarily controlled by vegetation at 6 cm, but deeper in the soil column topography and meteorological factors, such as predominant winter wind direction and summer insolation, play a larger role. The random forest model accounts for 40% to 60% of the variation observed (R2 = 0.44 at 6 cm, R2 = 0.52 at 12 cm, R2 = 0.58 at 20 cm). These results indicate that vegetation is the dominant control on soil moisture shallow in the soil column, but the impact of vegetation does not extend to deeper layers retrieved from P-Band SAR backscatter.

Dann, Julian↗
DOE OSTI2025

A Wind-driven Origin for the Firework Morphology of the Supernova Remnant Pa 30

Pa 30—the likely remnant of the Galactic Type Iax supernova of 1181 AD—displays an unusual, firework-like morphology, consisting of radial filaments extending from a common center, where a white dwarf (WD) currently drives a very fast wind (speed ≳10 4 km s −1 ). We propose the filaments arose from the Rayleigh–Taylor-unstable nature of the interface between the circumstellar medium (CSM) and the shocked wind launched by the natal WD; the filaments then elongated intact due to the Kelvin–Helmholtz-stable nature of the large initial density contrast between the wind and CSM, supplemented by the slowly declining wind density profile (relative to homologously expanding ejecta). To support this interpretation, we present two-dimensional hydrodynamical simulations and derive the filament properties, including their speed, density, and temperature, all of which are consistent with observations. We suggest the filaments elongate until the wind and CSM densities become comparable at the contact discontinuity, which occurs within 1–10 yr, and then truncate because the Rayleigh–Taylor instability (RTI) halts. The subsequent Kelvin–Helmholtz instability growth timescale across the current width of the filaments is longer than the age of Pa 30, so they remain intact. The filament-less central region in Pa 30 is therefore more likely a consequence of the finite timescale over which the RTI operates, rather than a wind termination shock. In general, firework-like filaments may form in other systems, provided there is a sufficiently large density contrast between the ejecta and its surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An Accurate SUPG-stabilized Continuous Galerkin Discretization for Anisotropic Heat Flux in Magnetic Confinement Fusion

We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov–Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗