Search NASA⌕ Search

SEARCH · Search NASA

Results for “ASYMPTOTIC SOLUTIONS”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗
DOE OSTI2020

Horizon constraints on holographic Green’s functions

We explore a new class of general properties of thermal holographic Green’s functions that can be deduced from the near-horizon behaviour of classical perturbations in asymptotically anti-de Sitter spacetimes. We show that at negative imaginary Matsubara frequencies and appropriate complex values of the wavenumber the retarded Green’s functions of generic operators are not uniquely defined, due to the lack of a unique ingoing solution for the bulk perturbations. From a boundary perspective these ‘pole-skipping’ points correspond to locations in the complex frequency and momentum planes at which a line of poles of the retarded Green’s function intersects with a line of zeroes. As a consequence the dispersion relations of collective modes in the boundary theory at energy scales ω ~ T are directly constrained by the bulk dynamics near the black-brane horizon. For the case of conserved U(1) current and energy-momentum tensor operators we give examples where the dispersion relations of hydrodynamic modes pass through a succession of poleskipping points as real wavenumber is increased. We discuss implications of our results for transport, hydrodynamics and quantum chaos in holographic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Small instantons in weakly-gauged holographic models

Small instantons can play an important role in Yang-Mills theories whose gauge couplings are sizeable at small distances. An interesting class of theories where this could occur is in weakly-gauged holographic models (dual to Yang-Mills theories interacting with strongly-coupled CFTs), since gauge couplings are indeed enhanced towards the UV boundary of the 5D AdS space. However, contrary to expectations, we show that small instantons in these non-asymptotically-free models are highly suppressed and ineffective. This is due to the conservation of topological charge that forbids instantons to be localized near the UV boundary. Despite this fact we find non-trivial UV localized instanton-anti-instanton solutions of the Yang-Mills equations where the topological charges annihilate in the AdS bulk. These analytic solutions arise from a 5D conformal transformation of the uplifted 4D instanton. Our analysis therefore reveals unexpected nonperturbative configurations of Yang-Mills theories when they interact with strongly-coupled CFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗
DOE OSTI2020

Asymptotic growth of the 4d $ \mathcal{N} $ = 4 index and partially deconfined phases

We study the Cardy-like asymptotics of the 4d N = 4 index and demonstrate the existence of partially deconfined phases where the asymptotic growth of the index is not as rapid as in the fully deconfined case. We then take the large-N limit after the Cardy-like limit and make a conjecture for the leading asymptotics of the index. While the Cardy-like behavior is derived using the integral representation of the index, we demonstrate how the same results can be obtained using the Bethe ansatz type approach as well. In doing so, we discover new non-standard solutions to the elliptic Bethe ansatz equations including continuous families of solutions for SU(N) theory with N ≥ 3. Here, we argue that the existence of both standard and continuous non-standard solutions has a natural interpretation in terms of vacua of N = 1* theory on R 3 × S 1 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A mathematical model of asynchronous data flow in parallel computers *

Abstract We present a simplified model of data flow on processors in a high-performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton–Jacobi equation associated with this conservation law for which the existence and uniqueness of solutions can be proven. We then present the results of numerical experiments for both discrete and continuum models; these show a qualitative agreement between the two and the effect of variations in the computing environment’s processing capabilities on the progress of the modelled computation.

Barnard, Richard C.↗
DOE OSTI2023

On triality defects in 2d CFT

We consider the triality fusion category discovered in the c = 1 Kosterlitz-Thouless theory. We analyze this fusion category using the tools from the group theoretical fusion category and compute the simple lines, fusion rules and F-symbols. We then studied the physical implication of this fusion category including deriving the spin selection rule, computing the asymptotic density of states of irreducible representations of the fusion category symmetries, and analyzing its anomaly and constraints under the renormalization group flow. There is another set of F-symbols for the fusion categories with the same fusion rule known in the literature. We find these two solutions are different as they lead to different spin selection rules. This gives a complete list of the fusion categories with the same fusion rule by the classification result in.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Stability Analysis of Parameter Varying Genetic Toggle Switches Using Koopman Operators

The genetic toggle switch is a well known model in synthetic biology that represents the dynamic interactions between two genes that repress each other. The mathematical models for the genetic toggle switch that currently exist have been useful in describing circuit dynamics in rapidly dividing cells, assuming fixed or time-invariant kinetic rates. There is a growing interest in being able to model and extend synthetic biological function for growth conditions such as stationary phase or during nutrient starvation. As cells transition from one growth phase to another, kinetic rates become time-varying parameters. In this paper, we propose a novel class of parameter varying nonlinear models that can be used to describe the dynamics of genetic circuits, including the toggle switch, as they transition from different phases of growth. We show that there exists unique solutions for this class of systems, as well as for a class of systems that incorporates the microbial phenomena of quorum sensing. Further, we show that the domain of these systems, which is the positive orthant, is positively invariant. We also showcase a theoretical control strategy for these systems that would grant asymptotic monostability of a desired fixed point. We then take the general form of these systems and analyze their stability properties through the framework of time-varying Koopman operator theory. A necessary condition for asymptotic stability is also provided as well as a sufficient condition for instability. A Koopman control strategy for the system is also proposed, as well as an analogous discrete time-varying Koopman framework for applications with regularly sampled measurements.

Harrison, Jamiree (ORCID:0000000265351579)↗
DOE OSTI2022

Formation of Point Shocks for 3D Compressible Euler

We consider the 3D isentropic compressible Euler equations with the ideal gas law. We provide a constructive proof of the formation of the first point shock from smooth initial datum of finite energy, with no vacuum regions, with nontrivial vorticity present at the shock, and under no symmetry assumptions . We prove that for an open set of Sobolev‐class initial data that are a small L ∞ perturbation of a constant state, there exist smooth solutions to the Euler equations which form a generic stable shock in finite time. The blowup time and location can be explicitly computed, and solutions at the blowup time are smooth except for a single point , where they are of cusp‐type with Hölder C 1/3 regularity. Our proof is based on the use of modulated self‐similar variables that are used to enforce a number of constraints on the blowup profile, necessary to establish global existence and asymptotic stability in self‐similar variables. © 2022 Wiley Periodicals LLC.

Mathematics↗
DOE OSTI2023

Accelerating cosmology from Λ < 0 gravitational effective field theory

A large class of Λ < 0 cosmologies have big-bang/big crunch spacetimes with time-symmetric backgrounds and asymptotically AdS Euclidean continuations suggesting a possible holographic realization. We argue that these models generically have time-dependent scalar fields, and these can lead to realistic cosmologies at the level of the homogeneous background geometry, with an accelerating phase prior to the turnaround and crunch. We first demonstrate via explicit effective field theory examples that models with an asymptotically AdS Euclidean continuation can also exhibit a period of accelerated expansion without fine tuning. We then show that certain significantly more tuned examples can give predictions arbitrarily close to a ΛCDM model. Finally, we demonstrate via an explicit construction that the potentials of interest can arise from a superpotential, thus suggesting that these solutions may be compatible with an underlying supersymmetric theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Finite-difference multiple fluid solution for source-driven rotation in highly magnetized linear plasma device

The rotation profile of a magnetized plasma cylinder composed of multiple fluids is investigated analytically, expanding on previous results. The analytic steady-state solution is used as a benchmark for a time-dependent multiple-fluid finite-difference code, MITNS: Multiple-Ion Transport Numerical Solver. Magnetic field evolution is taken into account, both analytically and numerically. Its details are shown to be of importance when particles are allowed out of the domain. MITNS reproduces the asymptotic expansion results for a small parameter $δ<<<1$. For $\sqrt{m_elm_i} ~ δ << 1$ , a slightly different regime, dominated by viscosity-induced transport of ions, is found numerically and analytically. Finally, this verification supports the use of this code for more complex time-dependent calculations in the future. Additionally, we derive the angular velocity profile of each species due to radial particle and charge fluxes of various strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Optical tweezer measurements of asymptotic nonlinearities in complex fluids

This article presents micro-medium-amplitude oscillatory shear (μMAOS), a method to measure the frequency-dependent micromechanical properties of soft materials in the asymptotically nonlinear regime using optical tweezers. For this work, we have developed a theoretical framework to extract these nonlinear mechanical properties of the material from experimental measurements and also proposed a physical interpretation of the third-order nonlinearities measured in single-tone oscillatory tests. We validate the method using a well-characterized surfactant solution of wormlike micelles, and subsequently employ this technique to demonstrate that the cytoplasm of a living cell undergoes strain softening and shear thinning when locally subjected to weakly nonlinear oscillatory deformations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Interrogating organo-mineral interactions by Py-GC/MS and ToF-SIMS

Soil organic matter (SOM) is crucial to soil fertility and agricultural sustainability in the face of climate change, and organo-mineral interactions are a key mechanism for stabilizing SOM. The stability of mineral-associated organic matter (MAOM) can be approximated by the strength of bonds between biomolecules and mineral surfaces. We hypothesized that factors such as mineral reactivity and solution pH would impact bonding strength. To test these factors, sorption isotherms were performed on minerals differing in particle size, morphology, and reactivity (quartz, goethite, gibbsite, kaolinite, montmorillonite), under two pH conditions (4 and 6), and across a range of tannic acid concentration (0-1400 ppm). The quantity and quality of organo-mineral bonds was assessed by quantifying post-sorption solution compounds using an absorbance plate reader and assessing the thermal energy needed to break bonds using evolved gas analysis (EGA) on a pyrolysis gas chromatograph-mass spectrometer (Py-GC/MS). Organo-mineral surfaces were also interrogated using Time-of-flight secondary ion mass spectrometry (ToF-SIMS). Sorption experiments revealed an asymptotic relationship between sorbed tannic acid and increasing initial concentrations which suggests our concentration range saturated available mineral surfaces. Results indicated that mineral identity, chiefly surface area and native pH, impacted the amount of compound sorbed to mineral surfaces. Evolved gas analyses of organo-mineral samples uncovered increased thermal stability of sorbed compounds compared to our tannic acid standard particularly on the iron oxide, goethite. ToF-SIMS revealed that goethite surface ions were covered by increasing concentrations of adsorbed tannic acid with the replacement of chlorine and sulfate ions by tannic acid related ions. Our results indicate that Py-GC/MS and ToF-SIMS are useful techniques to determine the nature and strength of organo-mineral bonds and are steppingstones to understanding controls on carbon sequestration in terrestrial systems.

58 GEOSCIENCES↗
DOE OSTI2022

Conformal wavefunctions for graviton amplitudes

The extended-BMS algebra of asymptotically flat spacetime contains an SO(3, 1) subgroup that acts by conformal transformations on the celestial sphere. It is of interest to study the representations of this subgroup associated with gravitons. To reduce the equation of motion to a Schrodinger-like equation it is necessary to impose a non-covariant gauge condition. Using these solutions, leading-order gauge invariant Weyl scalars are then computed and decomposed into families of unitary principal series representations. An invertible holographic mapping is constructed between these unitary principal series operators and massless spin-2 perturbations of flat spacetime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗