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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗
DOE OSTI2019

In Quest of Calibration Density and Consistency in Hydrologic Modeling: Distributed Parameter Calibration against Streamflow Characteristics

Conventional basin-by-basin approaches to calibrate hydrologic models are limited to gauged basins and typically result in spatially discontinuous parameter fields. Moreover, the consequent low calibration density in space falls seriously behind the need from present-day applications like high resolution river hydrodynamic modeling. In this study we calibrated three key parameters of the Variable Infiltration Capacity (VIC) model at every 1/8° grid-cell using machine learning-based maps of four streamflow characteristics for the conterminous United States (CONUS), with a total of 52,663 grid-cells. This new calibration approach, as an alternative to parameter regionalization, applied to ungauged regions too. A key difference made here is that we tried to regionalize physical variables (streamflow characteristics) instead of model parameters whose behavior may often be less well understood. The resulting parameter fields no longer presented any spatial discontinuities and the patterns corresponded well with climate characteristics, such as aridity and runoff ratio. The calibrated parameters were evaluated against observed streamflow from 704/648 (calibration/validation period) small-to-medium-sized catchments used to derive the streamflow characteristics, 3941/3809 (calibration/validation period) small-to-medium-sized catchments not used to derive the streamflow characteristics as well as five large basins. Comparisons indicated marked improvements in bias and Nash-Sutcliffe efficiency. Model performance was still poor in arid and semiarid regions, which is mostly due to both model structural and forcing deficiencies. Although the performance gain was limited by the relative small number of parameters to calibrate, the study and results here served as a proof-of-concept for a new promising approach for fine-scale hydrologic model calibrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Spherical compression of an applied magnetic field in inertial confinement fusion

Applying an external magnetic field to laser-driven inertial confinement fusion implosions is a promising approach for enhancing fusion yield. The field is compressed with the plasma, producing a magnetized hotspot that anisotropically suppresses thermal losses and traps alpha particles, making performance sensitive to the compressed field orientation. Here, we derive a simple, readily applicable analytic model that enables rapid evaluation of the compressed field topology and show that ablation into the hotspot amplifies the central field, while the ablated ice near the hotspot edge develops a decaying, radially bent field, with a discontinuity in the field direction. The radially bent field renders thermal insulation at the hotspot edge negligible and largely independent of the applied field strength, whereas insulation in the hotspot core still depends strongly on the applied field. Applying the model to non-axial initial field configurations, we find that an initially applied mirror field provides the greatest suppression, followed by the standard axial field.

Physics - Plasma physics↗
DOE OSTI2017

Validation of radiative transfer computation with Monte Carlo method for ultra-relativistic background flow

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation. Although we examined the effect of the shock structure on the emitted spectrum with a large number of cells, it is hard to employ so many computational cells per dimension in multi-dimensional simulations. Therefore, a further investigation with a smaller number of cells is required for obtaining realistic high-energy photons with multi-dimensional computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An analytical approach to evaluating magnetic-field closure and topological changes in FRC devices

Here, we describe mathematical methods based on optimizing a modified non-linear flux function (MFF) to evaluate whether odd-parity perturbations affect the local closure of magnetic field lines in field-reversed configurations. Using the MFF methodology, quantitative formulas are derived that provide the shift of the field minimum and the threshold for field-line opening, a discontinuous change in field topology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
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↗
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↗
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↗
DOE OSTI2018

Analysis of dynamics, stability, and flow fields' structure of an accelerated hydrodynamic discontinuity with interfacial mass flux by a general matrix method

We develop a general matrix method to analyze from a far field the dynamics of an accelerated interface between incompressible ideal fluids of different densities with interfacial mass flux and with negligible density variations and stratification. We rigorously solve the linearized boundary value problem for the dynamics conserving mass, momentum, and energy in the bulk and at the interface. We find a new hydrodynamic instability that develops only when the acceleration magnitude exceeds a threshold. This critical threshold value depends on the magnitudes of the steady velocities of the fluids, the ratio of their densities, and the wavelength of the initial perturbation. The flow has potential velocity fields in the fluid bulk and is shear-free at the interface. The interface stability is set by the interplay of inertia and gravity. For weak acceleration, inertial effects dominate, and the flow fields experience stable oscillations. For strong acceleration, gravity effects dominate, and the dynamics is unstable. For strong accelerations, this new hydrodynamic instability grows faster than accelerated Landau-Darrieus and Rayleigh-Taylor instabilities. For given values of the fluids' densities and their steady bulk velocities, and for a given magnitude of acceleration, we find the critical and maximum values of the initial perturbation wavelength at which this new instability can be stabilized and at which its growth is the fastest. The quantitative, qualitative, and formal properties of the accelerated conservative dynamics depart from those of accelerated Landau-Darrieus and Rayleigh-Taylor dynamics. New diagnostic benchmarks are identified for experiments and simulations of unstable interfaces.

Flow dynamics↗
DOE OSTI2021

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Investigation of the discontinuous precipitation of U-Nb alloys via thermodynamic analysis and phase-field modeling

U-Nb’s discontinuous precipitation, $\gamma^{bcc}_{matrix}$, $\alpha^{orth}_{cellular}$, $\gamma\prime^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $\gamma\prime$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb’s discontinuous precipitation hypothesized that the energy of $\alpha$ forms an intermediate common tangent with the first potential of the double-well energy of $\gamma$ at the $\gamma\prime$ composition. In this work, we examine different possible mechanisms by which the discontinuous precipitation product in the U-Nb system can be stabilized. We put forward a mechanism by which the bulk free energy of the $\gamma$ can develop a non-monotonic dependency with composition. Additionally we posit that local contributions due to lattice mismatch between the $\alpha$ and $\gamma$ phases may be responsible for the generation of metastable states that may stabilized by thermodynamics as well as by kinetics. Our work suggests that local misfit strain tends to play a crucial role in the growth of the discontinuous precipitation product. Depending on the magnitude of strain developed at the $\alpha$/$\gamma\prime$ interfaces, either an increasing $\gamma\prime$ composition or a random distribution of $\gamma\prime$ compositions around the equiatomic value with respect to increasing temperature could be expected. Moreover, we show how it is possible to stabilize the discontinuous precipitation front through highly anisotropic and fast interface diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Modeling strong discontinuities in the material point method using a single velocity field

We propose a new way to include strong velocity discontinuities in the material point method (MPM). We use the recent idea of damage gradient partitioning to dynamically determine whether or not deformed locations of material points are on the same side of a discontinuity as each node of a finite element discretization of the Eulerian velocity field. Velocity basis functions associated with the finite element nodes are then altered to ensure that their contributions to discontinuous components of the velocity are supported only at points assigned to the same side of the discontinuity. The partitioning of material points relative to each finite element node eliminates the need for indexing partitions, and accommodates discontinuities that do not partition material into disjoint bodies (e.g., propagating cracks). By building velocity discontinuities directly into the approximation space, sliding contact, separation, and multi-body interaction are possible within a single velocity field. The methodology avoids the need to track contact surfaces, construct and evolve explicit cracks, or use multiple velocity fields. Several numerical examples are computed to support these claims. The examples are also compared to the conventional MPM, demonstrating the method’s superiority in dealing with sliding contact and separation.

42 ENGINEERING↗