Search NASA⌕ Search

SEARCH · Search NASA

Results for “Vector calculus”

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

Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network.

Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs [6, 7]. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN’s optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques [8]. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.

97 MATHEMATICS AND COMPUTING↗

Generalized Grad–Shafranov equation for non-axisymmetric MHD equilibria

The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad–Shafranov equation, which is a single elliptic partial differential equation in two space dimensions. By revealing a hidden symmetry, here we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad–Shafranov equation. In contrast to solutions of the classical Grad–Shafranov equation, solutions of the generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad–Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computation of the Biot–Savart line integral with higher-order convergence using straight segments

One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot–Savart field from each straight segment is analytically known. However, if the endpoints of the straight segments are chosen to lie on the coil, then the accuracy of the Biot–Savart computation is generally only the second order in the number of endpoints. In this work, we propose a simple modification: shift each end point of the coil in the outward normal direction by an amount proportional to the local curvature. With this modification, the Biot–Savart accuracy increases to the fourth order and the numerical error is dramatically reduced for a given number of discretization points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A radial basis function Galerkin method for inhomogeneous nonlocal diffusion

We introduce a discretization for a nonlocal diffusion problem using a localized basis of radial basis functions. The stiffness matrix entries are assembled by a special quadrature routine unique to the localized basis. Combining the quadrature method with the localized basis produces a well-conditioned, sparse, symmetric positive definite stiffness matrix. We demonstrate that both the continuum and discrete problems are well-posed and present numerical results for the convergence behavior of the radial basis function method. As a result, we explore approximating the solution to anisotropic differential equations by solving anisotropic nonlocal integral equations using the radial basis function method.

97 MATHEMATICS AND COMPUTING↗

Convective mixing in vertically-layered porous media: The linear regime and the onset of convection

In this paper, we study the effect of permeability heterogeneity on the stability of gravitationally unstable, transient, diffusive boundary layers in porous media. Permeability is taken to vary periodically in the horizontal plane normal to the direction of gravity. In contrast to the situation for vertical permeability variation, the horizontal perturbation structures are multimodal. We therefore use a two-dimensional quasi-steady eigenvalue analysis as well as a complementary initial value problem to investigate the stability behavior in the linear regime, until the onset of convection. We find that thick permeability layers enhance instability compared with thin layers when heterogeneity is increased. On the contrary, for thin layers the instability is weakened progressively with increasing heterogeneity to the extent that the corresponding homogeneous case is more unstable. For high levels of heterogeneity, we find that a small change in the permeability field results in large variations in the onset time of convection, similar to the instability event in the linear regime. However, this trend does not persist unconditionally because of the reorientation of vorticity pairs due to the interaction of evolving perturbation structures with heterogeneity. Consequently, an earlier onset of instability does not necessarily imply an earlier onset of convection. A resonant amplification of instability is observed within the linear regime when the dominant perturbation mode is equal to half the wavenumber of permeability variation. On the other hand, a substantial damping occurs when the perturbation mode is equal to the harmonic and sub-harmonic components of the permeability wavenumber. Finally, the phenomenon of such harmonic interactions influences both the onset of instability as well as the onset of convection.

36 MATERIALS SCIENCE↗

An apparatus for generation and quantitative measurement of homogeneous isotropic turbulence in He ii

The superfluid phase of helium-4, known as He ii, exhibits extremely small kinematic viscosity and may be a useful tool for economically producing and studying high Reynolds number turbulent flow. Such applications are not currently possible because a comprehensive understanding of the complex two-fluid behavior of He ii is lacking. This situation could be remedied by a systematic investigation of simple, well controlled turbulence that can be directly compared with theoretical models. To this end, we have developed a new apparatus that combines flow visualization with second sound attenuation to study turbulence in the wake of a mesh grid towed through a He ii filled channel. One of three mesh grids (mesh number M = 3, 3.75, or 5 mm) can be pulled at speeds between 0.1 and 60 cm/s through a cast acrylic flow channel which has a 16 mm × 16 mm cross section and measures 330 mm long. The motion of solidified deuterium tracer particles, with diameter of the order 1 μm, in the resulting flow is captured by a high speed camera, and a particle tracking velocimetry algorithm resolves the Lagrangian particle trajectories through the turbulent flow field. A pair of oscillating superleak second sound transducers installed in the channel allows complementary measurement of vortex line density in the superfluid throughout the turbulent decay process. Here, success in early experiments demonstrates the effectiveness of both probes, and preliminary analysis of the data shows that both measurements strongly correlate with each other. Further investigations will provide comprehensive information that can be used to address open questions about turbulence in He ii and move toward the application of this fluid to high Reynolds number fluid research

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlocal Convection-Diffusion Problems on Bounded Domains and Finite-Range Jump Processes

In this paper, a nonlocal convection-diffusion model is introduced for the master equation of Markov jump processes in bounded domains. With minimal assumptions on the model parameters, the nonlocal steady and unsteady state master equations are shown to be well-posed in a weak sense. Finally, then the nonlocal operator is shown to be the generator of finite-range nonsymmetric jump processes and, when certain conditions on the model parameters hold, the generators of finite and infinite activity Lévy and Lévy-type jump processes are shown to be special instances of the nonlocal operator.

97 MATHEMATICS AND COMPUTING↗

Formulation analysis and computation of an optimization-based local-to-nonlocal coupling method

We present an optimization-based coupling method for local and nonlocal continuum models. Our approach couches the coupling of the models into a control problem where the states are the solutions of the nonlocal and local equations, the objective is to minimize their mismatch on the overlap of the local and nonlocal problem domains, and the virtual controls are the nonlocal volume constraint and the local boundary condition. We present the method in the context of Local-to-Nonlocal diffusion coupling. Numerical examples illustrate the theoretical properties of the approach.

97 MATHEMATICS AND COMPUTING↗

Differential geometric approaches to momentum-based formulations for fluids [Slides]

This SAND report documents CIS Late Start LDRD Project 22-0311, "Differential geometric approaches to momentum-based formulations for fluids". The project primarily developed geometric mechanics formulations for momentum-based descriptions of nonrelativistic fluids, utilizing a differential geometry/exterior calculus treatment of momentum and a space+time splitting. Specifically, the full suite of geometric mechanics formulations (variational/Lagrangian, Lie-Poisson Hamiltonian and Curl-Form Hamiltonian) were developed in terms of exterior calculus using vector-bundle valued differential forms. This was done for a fairly general version of semi-direct product theory sufficient to cover a wide range of both neutral and charged fluid models, including compressible Euler, magnetohydrodynamics and Euler-Maxwell. As a secondary goal, this project also explored the connection between geometric mechanics formulations and the more traditional Godunov form (a hyperbolic system of conservation laws). Unfortunately, this stage did not produce anything particularly interesting, due to unforeseen technical difficulties. There are two publications related to this work currently in preparation, and this work will be presented at SIAM CSE 23, at which the PI is organizing a mini-symposium on geometric mechanics formulations and structure-preserving discretizations for fluids. The logical next step is to utilize the exterior calculus based understanding of momentum coupled with geometric mechanics formulations to develop (novel) structure-preserving discretizations of momentum. This is the main subject of a successful FY23 CIS LDRD "Structure-preserving discretizations for momentum-based formulations of fluids".

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computer-aided detection of renal calculi from noncontrast CT images using TV-flow and MSER features

Purpose: Renal calculi are common extracolonic incidental findings on computed tomographic colonography (CTC). This work aims to develop a fully automated computer-aided diagnosis system to accurately detect renal calculi on CTC images. Methods: The authors developed a total variation (TV) flow method to reduce image noise within the kidneys while maintaining the characteristic appearance of renal calculi. Maximally stable extremal region (MSER) features were then calculated to robustly identify calculi candidates. Finally, the authors computed texture and shape features that were imported to support vector machines for calculus classification. The method was validated on a dataset of 192 patients and compared to a baseline approach that detects calculi by thresholding. The authors also compared their method with the detection approaches using anisotropic diffusion and nonsmoothing. Results: At a false positive rate of 8 per patient, the sensitivities of the new method and the baseline thresholding approach were 69% and 35% (p < 1e − 3) on all calculi from 1 to 433 mm{sup 3} in the testing dataset. The sensitivities of the detection methods using anisotropic diffusion and nonsmoothing were 36% and 0%, respectively. The sensitivity of the new method increased to 90% if only larger and more clinically relevant calculi were considered. Conclusions: Experimental results demonstrated that TV-flow and MSER features are efficient means to robustly and accurately detect renal calculi on low-dose, high noise CTC images. Thus, the proposed method can potentially improve diagnosis.

60 APPLIED LIFE SCIENCES↗

Fortran mimetic abstraction language (Formal) v0.1.

The Fortran mimetic abstraction language ("Formal") is a domain-specific language (DSL) embedded in Fortran 202Y [1]. Formal provides novel software abstractions for simulating phenomena governed by the partial differential equations (PDEs) of vector and tensor calculus. Such equations model an extremely broad set of physical phenomena, ranging from atmospheric winds to light propagation. Formal's data structures and algorithms mimic in form and behavior continuous functions and operators. Formal supports these mathematical constructs using mimetic discretizations that define a discrete calculus satisfying various tensor calculus theorems, thereby ensuring high-fidelity representations of the physics being modeled. [2] Formal 0.1.0 also lays a foundation for the future use of Fortran 202Y type-safe templates to facilitate the formal verification of tensor contractions in computational physics and artificial intelligence [3]. [1] "Fortran 202Y" is Fortran standard committee's informal designation for the next Fortran revision, which will likely be "Fortran 2028". [2] Corbino, J. and Castillo, J. (2020) Journal of Computational and Applied Mathematics, https://doi.org/10.1016/j.cam.2019.06.042. [3] Haveraaen, M., Järvi, J., & Rouson, D. (2019). Reflecting on Generics for Fortran. https://j3-fortran.org/doc/year/19/19-188.pdf.

Rouson, Damian [Lawrence Berkeley National Laborat↗

MiniTensor v. 1.0

Algebra, calculus, manipulation and optimization for small vectors and tensors.

Mota, Alejandro↗

Reggeometry of lepton- and hadron-induced reactions

We append a simple Pomeron pole amplitude by t and Q{sup 2}, MV dependencies inspired by geometrical ideas. The experimentally transition from soft to hard dynamics is realized by the introduction of a two-component Pomeron with different Q{sup 2} and M{sub V}-dependent residues. A unified description of deeply virtual Compton scattering (DVCS), of exclusive electroproduction of all vector mesons as well as of elastic pp scattering is suggested.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A model for soft high-energy scattering: Tensor pomeron and vector odderon

A model for soft high-energy scattering is developed. The model is formulated in terms of effective propagators and vertices for the exchange objects: the pomeron, the odderon, and the reggeons. The vertices are required to respect standard rules of QFT. The propagators are constructed taking into account the crossing properties of amplitudes in QFT and the power-law ansätze from the Regge model. We propose to describe the pomeron as an effective spin 2 exchange. This tensor pomeron gives, at high energies, the same results for the pp and pp{sup -bar} elastic amplitudes as the standard Donnachie–Landshoff pomeron. But with our tensor pomeron it is much more natural to write down effective vertices of all kinds which respect the rules of QFT. This is particularly clear for the coupling of the pomeron to particles carrying spin, for instance vector mesons. We describe the odderon as an effective vector exchange. We emphasise that with a tensor pomeron and a vector odderon the corresponding charge-conjugation relations are automatically fulfilled. We compare the model to some experimental data, in particular to data for the total cross sections, in order to determine the model parameters. The model should provide a starting point for a general framework for describing soft high-energy reactions. It should give to experimentalists an easily manageable tool for calculating amplitudes for such reactions and for obtaining predictions which can be compared in detail with data. -- Highlights: •A general model for soft high-energy hadron scattering is developed. •The pomeron is described as effective tensor exchange. •Explicit expressions for effective reggeon–particle vertices are given. •Reggeon–particle and particle–particle vertices are related. •All vertices respect the standard C parity and crossing rules of QFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗