Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuous fields”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗
NASA NTRS1974

Electric field measurements across the Harang discontinuity

The Harang discontinuity, the area separating the positive and negative bay regions in the midnight sector of the auroral zone, is a focal point for changes in behavior of many phenomena. Through this region the electric field, in a frame corotating with the earth, rotates through the west from a basically northward field in the positive bay region to a basically southward field in the negative bay region, appearing as a reversal in a single-axis measurement of the north-south component. Thirty-two of these reversals have been identified in the Ogo 6 data from November and December 1969. The discontinuity is dynamic in nature, moving southward and steepening its latitudinal profile as magnetic activity is increased. As activity decreases, it relaxes poleward and spreads out in latitudinal width. It occurs over several hours of magnetic local time.

Maynard, N. C.↗
DOE OSTI2024

Physics-informed neural networks for heterogeneous poroelastic media

This study presents a novel physics-informed neural network (PINN) framework for modeling poroelasticity in heterogeneous media with material interfaces. The approach introduces a composite neural network (CoNN) where separate neural networks predict displacement and pressure variables for each material. While sharing identical activation functions, these networks are independently trained for all other parameters. To address challenges posed by heterogeneous material interfaces, the CoNN is integrated with the Interface-PINNs (I-PINNs) framework (Sarma et al., Comput. Methods Appl. Mech. Eng. 429: 117135, 2024), allowing different activation functions across material interfaces. Further, this ensures accurate approximation of discontinuous solution fields and gradients. Performance and accuracy of this combined architecture were evaluated against the conventional PINNs approach, a single neural network (SNN) architecture, and the eXtended PINNs (XPINNs) framework through two one-dimensional benchmark examples with discontinuous material properties. The results show that the proposed CoNN with I-PINNs architecture achieves an RMSE that is two orders of magnitude better than the conventional PINNs approach and is at least 40 times faster than the SNN framework. Compared to XPINNs, the proposed method achieves an RMSE at least one order of magnitude better and is 40% faster.

42 ENGINEERING↗
NASA NTRS1991

Generation of Pc 1 waves by the ion temperature anisotropy associated with fast shocks caused by sudden impulses

The high correlation of Pc 1 events with magnetospheric compressions is known. A mechanism is proposed which leads to the generation of Pc 1 waves. The interaction of a dynamic pressure pulse with the earth's bow shock leads to the formation of a weak fast-mode shock propagating into the magnetoshealth. The shock wave can pass right through a tangential discontinuity (magnetopause) and into the magnetosphere, without disturbing either of the structures. In a quasiperpendicular geometry, the shock wave exhibits anisotropic heating. This anisotropy drives unstable ion-cyclotron waves which can contribute to the generation of the Pc 1 waves which are detected. The viability of the mechanism is demonstrated with simulations. This mechanism could explain the peak in the occurrence of observed Pc 1 waves in the postnoon sector where a field-aligned discontinuity in the solar wind would most often be parallel to the magnetopause surface due to the average Parker-spiral magnetic-field configuration.

Mandt, M. E.↗
DOE OSTI2018

An updated Lagrangian discontinuous Galerkin hydrodynamic method for gas dynamics

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for gas dynamics. The new method evolves conserved unknowns in the current configuration, which obviates the Jacobi matrix that maps the element in a reference coordinate system or the initial coordinate system to the current configuration. The density, momentum, and total energy (ρ, ρu, E ) are approximated with conservative higher-order Taylor expansions over the element and are limited toward a piecewise constant field near discontinuities using a limiter. Two new limiting methods are presented for enforcing the bounds on the primitive variables of density, velocity, and specific internal energy (ρ, u, e). The nodal velocity, and the corresponding forces, are calculated by solving an approximate Riemann problem at the element nodes. An explicit second-order method is used to temporally advance the solution. This new Lagrangian DG hydrodynamic method conserves mass, momentum, and total energy. 1D Cartesian coordinates test problem results are presented to demonstrate the accuracy and convergence order of the new DG method with the new limiters.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Discontinuities in Jovian sulphur plasma

The radial distribution of the Jovian sulfur plasma is discussed. Spectra of the ionized sulfur in the Jovian magnetosphere indicate a sharp discontinuity in visible sulfur emission at a radial distance of six Jupiter radii (Io orbit), outside of which emission is greatly reduced and the plasma is observed not to be in corotation with the planet. Possible explanations of this phenomenon include an effect of Io on the electron temperature outside its orbit, leading to the suppression of S II emission by second ionization, and a current system in Jupiter's ionosphere which prevents plasma from diffusing outward and decouples it from the corotation field. This discontinuity also requires reconciliation with Pioneer observations of the radial and pitch-angle diffusion of energetic charged particles.

Mekler, Y.↗
DOE OSTI2016

RELATIONSHIP BETWEEN DISTRIBUTION OF MAGNETIC DECAY INDEX AND FILAMENT ERUPTIONS

The decay index n of a horizontal magnetic field is considered to be an important parameter in judging the stability of a flux rope. However, the spatial distribution of this parameter has not been extensively explored so far. In this paper, we present a delineative study of the three-dimensional maps of n for two eruptive events, in which filaments underwent asymmetrical eruptions. The corresponding n -distributions are both found to show that the filaments tend to erupt at abnormal regions (dubbed ABN regions) of n . These ABN regions appear to be divided into two subregions, with larger and smaller n . Moreover, an analysis of the magnetic topological configuration of the ABN regions has been also performed. The results indicate that these ABN regions are associated with a kind of special quasi-separatrix layer across which the connectivity of magnetic field is discontinuous. The presented observations and analyses strongly suggest that the torus instability in ABN regions may play a crucial role for the triggering of an asymmetrical eruption. Additionally, our investigation can provide a way of forecasting how a filament might erupt, and predicting the location for an asymmetrically eruptive filament to be split through analyzing the spatial structure of n .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Small-scale energy storage and release as the cause of the stellar X-ray corona

It has yet to be established why ordinary stars possess X-ray coronas. Detailed observations of the Sun make it doubtful that the X-ray corona is produced by the dissipation of Alfven waves. The uniform X-ray brightness over all scales from 10,000 to 100,000 km simply does not look like a resonance phenomenon. It is demonstrated that the arbitrary winding patterns introduced into the bipolar magnetic fields of X-ray corona produce discontinuities within the field. It is suggested that the random motions of the footprints of the field cause the accumulation of internal strains in the field, which dissipate through neutral reconnection across the associated discontinuities to provide the primary heat source for the X-ray corona. It is emphasized, however, that the essential high resolution observations of the Sun have yet to be carried out, and until the theory is firmly established for the Sun, it cannot be known how to interpret the X-ray emission of the Sun or of the other stars.

Parker, E. N.↗
DOE OSTI2017

A multispacecraft event study of Pc5 ultralow-frequency waves in the magnetosphere and their external drivers

We investigate a quiet time event of magnetospheric Pc5 ultralow-frequency (ULF) waves and their likely external drivers using multiple spacecraft observations. Enhancements of electric and magnetic field perturbations in two narrow frequency bands, 1.5–2 mHz and 3.5–4 mHz, were observed over a large radial distance range from r ~ 5 to 11 RE. During the first half of this event, perturbations were mainly observed in the transverse components and only in the 3.5–4 mHz band. In comparison, enhancements were stronger during the second half in both transverse and compressional components and in both frequency bands. No indication of field line resonances was found for these magnetic field perturbations. Perturbations in these two bands were also observed in the magnetosheath, but not in the solar wind dynamic pressure perturbations. For the first interval, good correlations between the flow perturbations in the magnetosphere and magnetosheath and an indirect signature for Kelvin-Helmholtz (K-H) vortices suggest K-H surface waves as the driver. For the second interval, good correlations are found between the magnetosheath dynamic pressure perturbations, magnetopause deformation, and magnetospheric waves, all in good correspondence to interplanetary magnetic field (IMF) discontinuities. The characteristics of these perturbations can be explained by being driven by foreshock perturbations resulting from these IMF discontinuities. This event shows that even during quiet periods, K-H-unstable magnetopause and ion foreshock perturbations can combine to create a highly dynamic magnetospheric ULF wave environment

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Graphic statics and their potential for digital design and fabrication with concrete

The paper describes the influence of the lower bound theorem of theory of plasticity on the design of structures during the last 70 years. Starting with the fundamental theorems developed in the middle of the last century, the influence of rigid plastic (discontinuous) stress fields on the design of reinforced beam, wall, slab and shell structures is worked out. The potential of graphic statics for digital design and fabrication is illustrated by the results of several research projects on form finding and design of reinforced concrete structures carried out during the last decade. It is shown that the design methods allow for simultaneous considerations about structure and architecture. This leads to leaner structure and reduced material consumption. However, while a current limitation of manufacturing objects in these ''nonstandard'' shapes is presently cost limiting, digital fabrication may change this and open a much greater scope for the use and application of graphic statics in a near future.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interface PINNs (I-PINNs): A physics-informed neural networks framework for interface problems

Here, we present a novel physics-informed neural networks (PINNs) framework for modeling interface problems, termed Interface PINNs (I-PINNs). I-PINNs uses different neural networks for any two subdomains separated by a sharp interface such that the neural networks differ only through their activation functions while the other parameters remain identical. The performance of I-PINNs, conventional PINNs, and other existing domain-decomposition PINNs methods such as extended PINNs (XPINNs) and multi-domain PINN (M-PINN) is compared through several one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems. The results demonstrate that I-PINNs provides a root-mean-square-error accuracy, at least two orders of magnitude better than conventional PINNs and XPINNs at approximately one-tenth of the computational cost of conventional PINNs and half the cost of XPINNs. Additionally, while I-PINNs and M-PINN provide comparable accuracies, M-PINN is found to be approximately 50% more expensive.

42 ENGINEERING↗
DOE OSTI2016

Modeling Crack Propagation in Polycrystalline Microstructure Using Variational Multiscale Method

Crack propagation in a polycrystalline microstructure is analyzed using a novel multiscale model. The model includes an explicit microstructural representation at critical regions (stress concentrators such as notches and cracks) and a reduced order model that statistically captures the microstructure at regions far away from stress concentrations. Crack propagation is modeled in these critical regions using the variational multiscale method. In this approach, a discontinuous displacement field is added to elements that exceed the critical values of normal or tangential tractions during loading. Compared to traditional cohesive zone modeling approaches, the method does not require the use of any special interface elements in the microstructure and thus can model arbitrary crack paths. The capability of the method in predicting both intergranular and transgranular failure modes in an elastoplastic polycrystal is demonstrated under tensile and three-point bending loads.

36 MATERIALS SCIENCE↗
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 NTRS1993

A multilevel adaptive projection method for unsteady incompressible flow

There are two main requirements for practical simulation of unsteady flow at high Reynolds number: the algorithm must accurately propagate discontinuous flow fields without excessive artificial viscosity, and it must have some adaptive capability to concentrate computational effort where it is most needed. We satisfy the first of these requirements with a second-order Godunov method similar to those used for high-speed flows with shocks, and the second with a grid-based refinement scheme which avoids some of the drawbacks associated with unstructured meshes. These two features of our algorithm place certain constraints on the projection method used to enforce incompressibility. Velocities are cell-based, leading to a Laplacian stencil for the projection which decouples adjacent grid points. We discuss features of the multigrid and multilevel iteration schemes required for solution of the resulting decoupled problem. Variable-density flows require use of a modified projection operator--we have found a multigrid method for this modified projection that successfully handles density jumps of thousands to one. Numerical results are shown for the 2D adaptive and 3D variable-density algorithms.

Howell, Louis H.↗