Search NASA⌕ Search

SEARCH · Search NASA

Results for “RIEMANN FUNCTION”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

A new criterion for bar-forming instability in rapidly rotating gaseous and stellar systems. 2: Nonaxisymmetric form

We have previously introduced the parameter alpha as an indicator of stability to m = 2 nonaxisymmetric modes in rotating, self-gravitating, axisymmetric, gaseous (alpha less than or approximately equal to 0.34) and stellar (alpha less than or approximately equal to 0.25) systems. This parameter can be written as alpha = (ft/2)(exp 1/2), where t is defined as T/(absolute value of W), T is the total rotational kinetic energy, W is the total gravitational potential energy, and f is a function characteristic of the topology/connectedness and the geometric shape of a system. In this paper, we extend the stability criterion to nonaxisymmetric equilibrium systems by determining empirically the appropriate form of the function f for ellipsoids and elliptical disks and cylinders. We test the validity of this extension of the stability indicator alpha by considering its predictions for previously published, gaseous and stellar, nonaxisymmetric models. The formulation and critical values account accurately for the stability properties of m =2 modes in gaseous Riemann S-type ellipsoids (including the Jacobi and Dedekind ellipsoids) and elliptical Riemann disks as well as in stellar elliptical Freeman disks and cylinders: all these systems are dynamically stable except the stellar elliptical Freeman disks that exhibit a relatively small region of m = 2 dynamical instability. A partial disagreement in the case of stellar Freeman ellipsoids in maximum rotation may be due to the fact that the region of instability has not been previously determined with sufficient accuracy.

Christodoulou, Dimitris M.↗
DOE OSTI2016

Warped $AdS_6 × S^2$ in Type IIB supergravity I: local solutions

We investigate the existence of solutions with 16 residual supersymmetries to Type IIB supergravity on a space-time of the form $AdS_6×S^2$ warped over a two-dimensional Riemann surface $Σ$. The $SO(2,5)×SO(3)$ isometry extends to invariance under the exceptional Lie superalgebra $F$(4). In the present paper, we construct the general Ansatz compatible with these symmetries, derive the corresponding reduced BPS equations, and obtain their complete local solution in terms of two locally holomorphic functions $A_±$ on $Σ$, subject to certain positivity and regularity conditions. Globally, ($A_+,A_-$) are allowed to be multiple-valued on $Σ$ and be holomorphic sections of a holomorphic bundle over $Σ$ with structure group contained in $SU$(1,1)×C. Globally regular solutions are expected to provide the near-horizon geometry of (p,q) 5-brane and 7-brane webs which are holographic duals to five-dimensional conformal field theories. A preliminary analysis of the positivity and regularity conditions will be presented here, leaving the construction of globally regular solutions to a subsequent paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1990

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗
NASA NTRS1987

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗
NASA NTRS2017

Feature Detection and Curve Fitting Using Fast Walsh Transforms for Shock Tracking: Applications

Walsh functions form an orthonormal basis set consisting of square waves. Square waves make the system well suited for detecting and representing functions with discontinuities. Given a uniform distribution of 2p cells on a one-dimensional element, it has been proven that the inner product of the Walsh Root function for group p with every polynomial of degree < or = (p - 1) across the element is identically zero. It has also been proven that the magnitude and location of a discontinuous jump, as represented by a Heaviside function, are explicitly identified by its Fast Walsh Transform (FWT) coefficients. These two proofs enable an algorithm that quickly provides a Weighted Least Squares fit to distributions across the element that include a discontinuity. The detection of a discontinuity enables analytic relations to locally describe its evolution and provide increased accuracy. Time accurate examples are provided for advection, Burgers equation, and Riemann problems (diaphragm burst) in closed tubes and de Laval nozzles. New algorithms to detect up to two C0 and/or C1 discontinuities within a single element are developed for application to the Riemann problem, in which a contact discontinuity and shock wave form after the diaphragm bursts.

Gnoffo, Peter A.↗
DOE OSTI2020

BPS spectra and 3-manifold invariants

We provide a physical definition of new homological invariants H a (M 3 ) of 3-manifolds (possibly, with knots) labeled by abelian flat connections. The physical system in question involves a 6d fivebrane theory on M 3 times a 2-disk, D 2 , whose Hilbert space of BPS states plays the role of a basic building block in categorification of various partition functions of 3d N = 2 theory T[M 3 ]: D 2 × S 1 half-index, S 2 × S 1 superconformal index, and S 2 × S 1 topologically twisted index. The first partition function is labeled by a choice of boundary condition and provides a refinement of Chern–Simons (WRT) invariant. A linear combination of them in the unrefined limit gives the analytically continued WRT invariant of M 3 . The last two can be factorized into the product of half-indices. Here, we show how this works explicitly for many examples, including Lens spaces, circle fibrations over Riemann surfaces, and plumbed 3-manifolds.

3-manifold↗
NASA NTRS1991

Development of a grid-independent approximate Riemannsolver

A grid-independent approximate Riemann solver for use with the Euler and Navier-Stokes equations was introduced and explored. The two-dimensional Euler and Navier-Stokes equations are described in Cartesian and generalized coordinates, as well as the traveling wave form of the Euler equations. The spatial and temporal discretization are described for both explicit and implicit time-marching schemes. The grid-aligned flux function of Roe is outlined, while the 5-wave grid-independent flux function is derived. The stability and monotonicity analysis of the 5-wave model are presented. Two-dimensional results are provided and extended to three dimensions. The corresponding results are presented.

Rumsey, Christopher Lockwood↗
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗
DOE OSTI2020

Riemann-Hilbert correspondence and blown up surface defects

The relationship of two dimensional quantum field theory and isomonodromic deformations of Fuchsian systems has a long history. Recently four-dimensional N = 2 gauge theories joined the party in a multitude of roles. In this paper we study the vacuum expectation values of intersecting half-BPS surface defects in SU(2) theory with N f = 4 fundamental hypermultiplets. We show they form a horizontal section of a Fuchsian system on a sphere with 5 regular singularities, calculate the monodromy, and define the associated isomonodromic tau-function. Using the blowup formula in the presence of half-BPS surface defects, initiated in the companion paper, we obtain the GIL formula, establishing an unexpected relation of the topological string/free fermion regime of supersymmetric gauge theory to classical integrability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗
DOE OSTI2018

Discrete Hahn polynomials for numerical solution of two-dimensional variable-order fractional Rayleigh–Stokes problem

The main aim of this paper is to find the numerical solutions of 2D Rayleigh–Stokes problem with the variable-order fractional derivatives in the Riemann–Liouville sense. The presented method is based on collocation procedure in combination with the new operational matrix of the variable-order fractional derivatives, in the Caputo sense, for the discrete Hahn polynomials. The main advantage of the proposed method is obtaining a global approximation for spatial and temporal discretizations, and it reduced the problem to an algebraic system, which is easier to solve. Also, the profit of approximating a continuous function by Hahn polynomials is that for computing the coefficients of the expansion, we only have to compute a summation and the calculation of coefficients is exact. The error bound for the approximate solution is estimated. Finally, we evaluate results of the presented method with other numerical methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Coaction and double-copy properties of configuration-space integrals at genus zero

We investigate configuration-space integrals over punctured Riemann spheres from the viewpoint of the motivic Galois coaction and double-copy structures generalizing the Kawai-Lewellen-Tye (KLT) relations in string theory. For this purpose, explicit bases of twisted cycles and cocycles are worked out whose orthonormality simplifies the coaction. We present methods to efficiently perform and organize the expansions of configuration-space integrals in the inverse string tension α' or the dimensional-regularization parameter ϵ of Feynman integrals. Generating-function techniques open up a new perspective on the coaction of multiple polylogarithms in any number of variables and analytic continuations in the unintegrated punctures. We present a compact recursion for a generalized KLT kernel and discuss its origin from intersection numbers of Stasheff polytopes and its implications for correlation functions of two-dimensional conformal field theories. We find a non-trivial example of correlation functions in (p, 2) minimal models, which can be normalized to become uniformly transcendental in the p → ∞ limit.

Conformal Field Theory↗
DOE OSTI2020

Convex Programs for Minimal-Area Problems

The minimal-area problem that defines string diagrams in closed string field theory asks for the metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length at least 2π. This is an extremal length problem in conformal geometry as well as a problem in systolic geometry. Here, we consider the analogous minimal-area problem for homology classes of curves and, with the aid of calibrations and the max flow-min cut theorem, formulate it as a local convex program. We derive an equivalent dual program involving maximization of a concave functional. These two programs give new insights into the form of the minimal-area metric and are amenable to numerical solution. We explain how the homology problem can be modified to provide the solution to the original homotopy problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

On a criterion of conformal parabolicity of a Riemannian manifold

The paper relates to the circle of problems concerning the connection between the conformal type of a Riemannian manifold and the canonical form of its isoperimetric function. Two special examples of 2-manifolds are constructed, which explain the meaning, role and importance of the conditions involved in the criterion, previously obtained by the author, which decides whether a noncompact Riemannian n-manifold is conformally parabolic. Bibliography: 8 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

A thermodynamically compatible splitting procedure in hyperelasticity

A material is hyperelastic if the stress tensor is obtained by variation of the stored energy function. The corresponding 3D mathematical model of hyperelasticity written in the Eulerian coordinates represents a system of 14 conservative partial differential equations submitted to stationary differential constraints. A classical approach for numerical solving of such a 3D system is a geometrical splitting: the 3D system is split into three 1D systems along each spatial direction and solved then by using a Godunov type scheme. Each 1D system has 7 independent eigenfields corresponding to contact discontinuity, longitudinal waves and shear waves. The construction of the corresponding Riemann solvers is not an easy task even in the case of isotropic solids. Indeed, for a given specific energy it is extremely difficult, if not impossible, to check its rank-one convexity which is a necessary and sufficient condition for hyperbolicity of the governing equations. In this paper, we consider a particular case where the specific energy is a sum of two terms. The first term is the hydrodynamic energy depending only on the density and the entropy, and the second term is the shear energy which is unaffected by the volume change. In this case a very simple criterion of hyperbolicity can be formulated. We propose then a new splitting procedure which allows us to find a numerical solution of each 1D system by solving successively three 1D sub-systems. Each sub-system is hyperbolic, if the full system is hyperbolic. Moreover, each sub-system has only three waves instead of seven, and the velocities of these waves are given in explicit form. The last property allows us to construct reliable Riemann solvers. Numerical 1D tests confirm the advantage of the new approach. A multi-dimensional extension of the splitting procedure is also proposed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

A nominally second-order cell-centered Lagrangian scheme for simulating elastic–plastic flows on two-dimensional unstructured grids

In this paper, we describe a cell-centered Lagrangian scheme devoted to the numerical simulation of solid dynamics on two-dimensional unstructured grids in planar geometry. This numerical method, utilizes the classical elastic-perfectly plastic material model initially proposed by Wilkins [M.L. Wilkins, Calculation of elastic–plastic flow, Meth. Comput. Phys. (1964)]. In this model, the Cauchy stress tensor is decomposed into the sum of its deviatoric part and the thermodynamic pressure which is defined by means of an equation of state. Regarding the deviatoric stress, its time evolution is governed by a classical constitutive law for isotropic material. The plasticity model employs the von Mises yield criterion and is implemented by means of the radial return algorithm. The numerical scheme relies on a finite volume cell-centered method wherein numerical fluxes are expressed in terms of sub-cell force. The generic form of the sub-cell force is obtained by requiring the scheme to satisfy a semi-discrete dissipation inequality. Sub-cell force and nodal velocity to move the grid are computed consistently with cell volume variation by means of a node-centered solver, which results from total energy conservation. The nominally second-order extension is achieved by developing a two-dimensional extension in the Lagrangian framework of the Generalized Riemann Problem methodology, introduced by Ben-Artzi and Falcovitz [M. Ben-Artzi, J. Falcovitz, Generalized Riemann Problems in Computational Fluid Dynamics, Cambridge Monogr. Appl. Comput. Math. (2003)]. Finally, the robustness and the accuracy of the numerical scheme are assessed through the computation of several test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The Effective Dynamics of the Volume Preserving Mean Curvature Flow

We consider the dynamics of small closed submanifolds (‘bubbles’) under the volume preserving mean curvature flow. We construct a map from (n+1)-dimensional Euclidean space into a given (n+1)-dimensional Riemannian manifold which characterizes the existence, stability and dynamics of constant mean curvature submanifolds. This is done in terms of a reduced area function on the Euclidean space, which is given constructively and can be computed perturbatively. This allows us to derive adiabatic and effective dynamics of the bubbles. The results can be mapped by rescaling to the dynamics of fixed size bubbles in almost Euclidean Riemannian manifolds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

An a 0 resonance in strongly coupled π η , K K ¯ scattering from lattice QCD

Here, we present the first calculation of coupled-channel meson-meson scattering in the isospin $=1$, $G$-parity negative sector, with channels $\pi \eta$, $K\overline{K}$ and $\pi \eta'$, in a first-principles approach to QCD. From the discrete spectrum of eigenstates in three volumes extracted from lattice QCD correlation functions we determine the energy dependence of the $S$-matrix, and find that the $S$-wave features a prominent cusp-like structure in $\pi \eta \to \pi \eta$ close to $K\overline{K}$ threshold coupled with a rapid turn on of amplitudes leading to the $K\overline{K}$ final-state. This behavior is traced to an $a_0(980)$-like resonance, strongly coupled to both $\pi \eta$ and $K\overline{K}$, which is identified with a pole in the complex energy plane, appearing on only a single unphysical Riemann sheet. Consideration of $D$-wave scattering suggests a narrow tensor resonance at higher energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗