Search NASA⌕ Search

SEARCH · Search NASA

Results for “Elliptic problems”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An Improved Discrete Least-Squares/Reduced-Basis Method for Parameterized Elliptic PDEs

It is shown that the computational efficiency of the discrete least-squares (DLS) approximation of solutions of stochastic elliptic PDEs is improved by incorporating a reduced-basis method into the DLS framework. In particular, we consider stochastic elliptic PDEs with an affine dependence on the random variable. The goal is to recover the entire solution map from the parameter space to the finite element space. To this end, first, a reduced-basis solution using a weak greedy algorithm is constructed, then a DLS approximation is determined by evaluating the reduced-basis approximation instead of the full finite element approximation. The main advantage of the new approach is that one only need apply the DLS operator to the coefficients of the reduced-basis expansion, resulting in huge savings in both the storage of the DLS coefficients and the online cost of evaluating the DLS approximation. In addition, the recently developed quasi-optimal polynomial space is also adopted in the new approach, resulting in superior convergence rates for a wider class of problems than previous analyzed. Numerical experiments are provided that illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Validation of general purpose mathematical efficiency modeling with ISOCS - 15579

Gamma spectroscopy is a commonly used method for radioactive waste characterization and activity determination in in-situ objects. The gamma rays emitted by radionuclides in the waste container matrix are highly penetrating and can be measured using one or more gamma ray detectors, thus allowing for a passive nondestructive assay. One problem is that a source based calibration of full energy peak efficiency can be difficult or impractical because of the large size and/or an unusual shape of the sample geometry. One attractive solution is to use mathematical efficiency modeling to calculate the full energy peak efficiency. The In-Situ Object Calibration Software (ISOCS) has successfully been used for more than a decade for mathematical efficiency modeling of template based geometries. Templates include the most common objects in waste characterization and D and D, including but not limited to: boxes, drums, and pipes with hotspots, beams, tanks and walls. Scenarios exist where these templates cannot model the desired measurement geometry. Canberra has developed a general purpose version of ISOCS where primitive objects can be combined into practically any geometry. The primitive objects are right circular cylinder, right elliptical cylinder, right circular cone, right elliptical cone, sphere, spheroid, ellipsoid, box (rectangular parallelepiped), torus, and beaker (any cylindrical symmetric object). Multiple primitive objects can be combined using the priority scheme and cut by planes into almost any desired shape. In addition, multiple objects can be defined in the same scene, for example cylinders next to each other or a pipe and a box. Collimators can be constructed as easily as in template based ISOCS. The ray tracing algorithm in ISOCS has been extended to support the new primitive objects for attenuation calculation between the radioactive source and the detector end cap. The new version uses the same vacuum efficiency response, or characterization, in the energy range 10 - 7000 keV of individual detectors to be able to calculate the efficiency as accurately as possible. MCNP has been shown to be able to calculate detection efficiencies for any geometry and the new version of ISOCS have been validated to MCNP. The vacuum efficiency is based on an MCNP model and this model was the base for constructing an identical model in MCNP as in ISOCS. The calculated efficiencies from the two codes have been compared and the agreement is within 3% except for the lowest energies for one geometry. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Towards Using Eshelby Calculations to Enhance Kinetic Model for Zirconium Hydride Precipitation

A C++ library (called Eshelby ) was implemented in fiscal year 2015 based upon the formulas documented in this report. The library implements a generalized version of Eshelby's inclusion problem. The library was written as a set of functions which can be called from another program; the principle intended use cases are kinetic models of precipitate formation in zirconium claddings where use of the Eshelby library provides needed elastic energy density calculations, as well as calculations of stress and strain in and around precipitates; it is intended that the library will be made open source. For isotropic inclusions in the form of oblate and prolate ellipsoids, the Eshelby library can be used for nearly any relevant/appropriate shape parameters to calculate strains, stresses and energy density at interior and exterior points. The Eshelby library uses a combination of analytical formulas and numerical routines making it very extensible. For example, the library can can easily be extended to include inclusions such as spheres since analytical expressions exist for the required elliptic integrals; similarly, general ellipsoids do not have analytical results for the required elliptic integrals but those integrals can be numerically evaluated and thus fit naturally into the Eshelby library. This report documents all formulas implemented in the Eshelby library and presents some demonstration calculations relevant to the intended application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Large-scale assembly bias of dark matter halos

We present precise measurements of the assembly bias of dark matter halos, i.e. the dependence of halo bias on other properties than the mass, using curved 'separate universe' N-body simulations which effectively incorporate an infinite-wavelength matter overdensity into the background density. This method measures the LIMD (local-in-matter-density) bias parameters b {sub n} in the large-scale limit. We focus on the dependence of the first two Eulerian biases b {sup E} {sup {sub 1}} and b {sup E} {sup {sub 2}} on four halo properties: the concentration, spin, mass accretion rate, and ellipticity. We quantitatively compare our results with previous works in which assembly bias was measured on fairly small scales. Despite this difference, our findings are in good agreement with previous results. We also look at the joint dependence of bias on two halo properties in addition to the mass. Finally, using the excursion set peaks model, we attempt to shed new insights on how assembly bias arises in this analytical model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Comparison of multigrid algorithms for high‐order continuous finite element discretizations

Summary We present a comparison of different multigrid approaches for the solution of systems arising from high‐order continuous finite element discretizations of elliptic partial differential equations on complex geometries. We consider the pointwise Jacobi, the Chebyshev‐accelerated Jacobi, and the symmetric successive over‐relaxation smoothers, as well as elementwise block Jacobi smoothing. Three approaches for the multigrid hierarchy are compared: (1) high‐order h ‐multigrid, which uses high‐order interpolation and restriction between geometrically coarsened meshes; (2) p ‐multigrid, in which the polynomial order is reduced while the mesh remains unchanged, and the interpolation and restriction incorporate the different‐order basis functions; and (3) a first‐order approximation multigrid preconditioner constructed using the nodes of the high‐order discretization. This latter approach is often combined with algebraic multigrid for the low‐order operator and is attractive for high‐order discretizations on unstructured meshes, where geometric coarsening is difficult. Based on a simple performance model, we compare the computational cost of the different approaches. Using scalar test problems in two and three dimensions with constant and varying coefficients, we compare the performance of the different multigrid approaches for polynomial orders up to 16. Overall, both h ‐multigrid and p ‐multigrid work well; the first‐order approximation is less efficient. For constant coefficients, all smoothers work well. For variable coefficients, Chebyshev and symmetric successive over‐relaxation smoothing outperform Jacobi smoothing. While all of the tested methods converge in a mesh‐independent number of iterations, none of them behaves completely independent of the polynomial order. When multigrid is used as a preconditioner in a Krylov method, the iteration number decreases significantly compared with using multigrid as a solver. Copyright © 2015 John Wiley & Sons, Ltd.

Sundar, Hari↗
DOE OSTI2017

Optical properties of magnetic photonic crystals with an arbitrary magnetization orientation

We have studied the peculiarities of diffraction of light in magnetic photonic crystals at large values of magnetooptical activity parameter and modulation depth. We have considered the case of an arbitrary angle between the directions of the external static magnetic field and the normal to the layer. The problem has been solved by the modified Ambartsumyan layer summation method. It has been shown that the given system is nonreciprocal with respect to not only circular, but linear polarizations also. In this case, a new type of nonreciprocity is observed (namely, the relation R(α) ≠ R(–α) holds, where R is the reflection coefficient and α is the angle of incidence). It has been demonstrated that in the case of oblique incidence, there appears a new photonic forbidden band that is not selective relative to the polarization of incident light. We have detected strong dependences of reflectance, absorbance, transmittance nonreciprocity, and other characteristics on the angle between the direction of the external static magnetic field and the normal to the layer boundary. Such a system can be used as a controllable polarization filter and a mirror, as well as a source of circular (elliptic) polarization, a controllable optical diode, and so on.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Dynamical Friction and the Evolution of Supermassive Black Hole Binaries: The Final Hundred-parsec Problem

The supermassive black holes originally in the nuclei of two merging galaxies will form a binary in the remnant core. The early evolution of the massive binary is driven by dynamical friction before the binary becomes “hard” and eventually reaches coalescence through gravitational-wave emission. We consider the dynamical friction evolution of massive binaries consisting of a secondary hole orbiting inside a stellar cusp dominated by a more massive central black hole. In our treatment, we include the frictional force from stars moving faster than the inspiralling object, which is neglected in the standard Chandrasekhar treatment. We show that the binary eccentricity increases if the stellar cusp density profile rises less steeply than ρ∝r{sup −2}. In cusps shallower than ρ∝r{sup −1}, the frictional timescale can become very long due to the deficit of stars moving slower than the massive body. Although including fast stars increases the decay rate, low mass-ratio binaries (q≲10{sup −3}) in sufficiently massive galaxies have decay timescales longer than one Hubble time. During such minor mergers, the secondary hole stalls on an eccentric orbit at a distance of order one-tenth the influence radius of the primary hole (i.e., ≈10--100 pc for massive ellipticals). We calculate the expected number of stalled satellites as a function of the host galaxy mass and show that the brightest cluster galaxies should have ≳1 of such satellites orbiting within their cores. Our results could provide an explanation for a number of observations, which include multiple nuclei in core ellipticals, off-center AGNs, and eccentric nuclear disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Solving Critical Problems of the Muon Collider Higgs Factory: Optics, Magnets and their Protection, Detector Backgrounds

A low-energy medium-luminosity Muon Collider (MC) is being studied as a possible Higgs Factory (HF). Electrons from muon decays will deposit more than 300 kW in superconducting magnets of the HF collider ring. This imposes significant challenges to superconducting (SC) magnets used in the MC storage ring (SR) and interaction regions (IR). Magnet designs are proposed which provide high operating gradient and magnetic field in a large aperture to accommodate the large size of muon beams (due to low beta*), as well as a cooling system to intercept the large heat deposition from the showers induced by decay electrons. The distribution of heat deposition in the MC SR lattice elements requires large-aperture magnets in order to accommodate thick high-Z absorbers to protect the SC coils. Based on the developed MARS15 model and intensive simulations, a sophisticated radiation protection system was designed for the collider SR and IR to bring the peak power density in the superconducting coils below the quench limit and reduce the dynamic heat deposition in the cold mass by a factor of 100. The system consists of tight tungsten masks in the magnet interconnect regions and elliptical tungsten liners in the magnet aperture optimized individually for each magnet. These also reduce the background particle fluxes in the collider detector.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Generation of circularly polarized XUV and soft-x-ray high-order harmonics by homonuclear and heteronuclear diatomic molecules subject to bichromatic counter-rotating circularly polarized intense laser fields

Recently, studies of bright circularly polarized high-harmonic beams from atoms in the soft-x-ray region as a source for x-ray magnetic circular dichroism measurement in a tabletop-scale setup have received considerable attention. In this paper, we address the problem with molecular targets and perform a detailed quantum study of H 2 +, CO, and N 2 molecules in bichromatic counter-rotating circularly polarized laser fields where we adopt wavelengths (1300 and 790 nm) and intensities (2 × 10 14 W/cm 2 ) reported in a recent experiment [Proc. Natl. Acad. Sci. USA 112, 14206 (2015)]. Our treatment of multiphoton processes in homonuclear and heteronuclear diatomic molecules is nonperturbative and based on the time-dependent density-functional theory for multielectron systems. The calculated radiation spectrum contains doublets of left and right circularly polarized harmonics with high-energy photons in the XUV and soft-x-ray ranges. Overall, our results reveal intriguing and substantially different nonlinear optical responses for homonuclear and heteronuclear diatomic molecules subject to circularly polarized intense laser fields. We study in detail the below- and above-threshold harmonic regions and analyze the ellipticity and phase of the generated harmonic peaks.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

A scalable geometric multigrid solver for nonsymmetric elliptic systems with application to variable-density flows

A geometric multigrid algorithm is introduced for solving nonsymmetric linear systems resulting from the discretization of the variable density Navier–Stokes equations on nonuniform structured rectilinear grids and high-Reynolds number flows. The restriction operation is defined such that the resulting system on the coarser grids is symmetric, thereby allowing for the use of efficient smoother algorithms. To achieve an optimal rate of convergence, the sequence of interpolation and restriction operations are determined through a dynamic procedure. A parallel partitioning strategy is introduced to minimize communication while maintaining the load balance between all processors. To test the proposed algorithm, we consider two cases: 1) homogeneous isotropic turbulence discretized on uniform grids and 2) turbulent duct flow discretized on stretched grids. Testing the algorithm on systems with up to a billion unknowns shows that the cost varies linearly with the number of unknowns. This O (N) behavior confirms the robustness of the proposed multigrid method regarding ill-conditioning of large systems characteristic of multiscale high-Reynolds number turbulent flows. The robustness of our method to density variations is established by considering cases where density varies sharply in space by a factor of up to 104, showing its applicability to two-phase flow problems. Strong and weak scalability studies are carried out, employing up to 30,000 processors, to examine the parallel performance of our implementation. Excellent scalability of our solver is shown for a granularity as low as 104 to 105 unknowns per processor. At its tested peak throughput, it solves approximately 4 billion unknowns per second employing over 16,000 processors with a parallel efficiency higher than 50%.

Computer Science↗
DOE OSTI2017

CONTRIBUTIONS OF THE ELASTICITY TO THE PRECESSION OF A TWO-LAYER EARTH MODEL

We focus on the updating of a specific contribution to the precession of the equator in longitude, usually named as “second order.” It stems from the crossing of certain terms of the lunisolar gravitational potential. The IAU2006 precession theory assigns it the value of −46.8 mas/cy that was derived for a rigid Earth model. Instead of that model, we consider a two-layer Earth composed of an elastic mantle and a liquid core, working out the problem within the Hamiltonian framework developed by Getino and Ferrándiz. The targeted effect is obtained without further simplifying assumptions through Hori’s canonical perturbation method applied up to the second order of perturbation. On account of using a more realistic Earth model, the revised value of the second-order contribution is significantly changed and reaches −55.29 mas/cy. That variation of the second-order contribution is larger than other contributions included in IAU2006. It must be compensated with an increase of −8.51 mas/cy in the value of the lunisolar first-order component p{sub A}{sup ′} of the precession of the equator rate, which is derived from the total rate by subtracting the remaining contributions accounted for in IAU2006 precession. The updating of the second-order contribution implies that the p{sub A}{sup ′} parameter has to be changed, from 5040684.593 to 5040693.104 mas/cy in absence of potential revisions of other contributions. It entails a proportional variation of Earth’s dynamical ellipticity H{sub d}, for which the estimation associated with IAU2006, 0.00327379448, should be updated to 0.00327380001, about 1.7 ppm larger.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Electromagnetically induced absorption and electromagnetically induced transparency for optical transitions F{sub g} → F{sub e} in the field of elliptically polarized waves

Nonlinear laser spectroscopy is considered in the scheme of two collinear waves with arbitrary elliptical polarizations. Emphasis is placed on investigating the nonlinear corrections in the absorption spectrum of one of the waves. The spontaneous transfer of low-frequency Zeeman coherence is shown to affect the sign of the subnatural-width resonance. For a closed transition, the direction of the resonance profile has been found to depend only on the angular momenta F{sub e} and F{sub g}. On this basis, a classification has been developed for various transitions by the direction of the subnatural-width resonance profile.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Efficient Bayesian inference with latent Hamiltonian neural networks in No-U-Turn Sampling

When sampling for Bayesian inference, one popular approach in the computational field is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS), which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice when relying on computationally expensive forward models. We propose Latent Hamiltonian neural networks (L-HNNs) with HMC and NUTS for solving Bayesian inference problems. Once trained, L-HNNs do not require numerical gradients of the target density during sampling, and hence numerous evaluations of the forward computational model. Moreover, L-HNNs satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well-suited for use within HMC and NUTS because stationarity can be shown. We also propose the integration of L-HNNs in an online error monitoring scheme, in which numerical gradients of the target density are used for a few samples whenever the L-HNNs prediction errors are large. This online error monitor scheme prevents sample degeneracy in regions of low probability density and ensures robust uncertainty quantification. We demonstrate L-HNNs in NUTS with online error monitoring on several analytical examples involving complex, heavy-tailed, and high-local-curvature probability densities. We then demonstrate the applicability of L-HNNs in NUTS to two computational case studies, namely the Allen-Cahn stochastic partial differential equation and an elliptic partial differential equation with 25 and 50 inference parameters, respectively. Overall, the L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. In conclusion, compared to traditional NUTS, L-HNNs in NUTS with online error monitoring required 1–2 orders of magnitude fewer numerical gradients of the target density and improved the effective sample size (ESS) per gradient (which is a measure of both the sampling quality and the computational expense) by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Vortex equations governing the fractional quantum Hall effect

An existence theory is established for a coupled non-linear elliptic system, known as “vortex equations,” describing the fractional quantum Hall effect in 2-dimensional double-layered electron systems. Via variational methods, we prove the existence and uniqueness of multiple vortices over a doubly periodic domain and the full plane. In the doubly periodic situation, explicit sufficient and necessary conditions are obtained that relate the size of the domain and the vortex numbers. For the full plane case, existence is established for all finite-energy solutions and exponential decay estimates are proved. Quantization phenomena of the magnetic flux are found in both cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

On the Post-Keplerian Corrections to the Orbital Periods of a Two-body System and Their Application to the Galactic Center

We perform detailed numerical analyses of the orbital motion of a test particle around a spinning primary, with the aim of investigating the possibility of using the post-Keplerian (pK) corrections to the orbiter’s periods (draconitic, anomalistic, and sidereal) as a further opportunity to perform new tests of post-Newtonian gravity. As a specific scenario, the S-stars orbiting the massive black hole (MBH) supposedly lurking in Sgr A* at the center of the Galaxy are adopted. We first study the effects of the pK Schwarzchild, Lense–Thirring, and quadrupole moment accelerations experienced by a target star for various possible initial orbital configurations. It turns out that the results of the numerical simulations are consistent with the analytical ones in the small eccentricity approximation for which almost all the latter ones were derived. For highly elliptical orbits, the sizes of the three pK corrections considered turn out to increase remarkably. The periods of the observed S2 and S0-102 stars as functions of the MBH’s spin axis orientation are considered as well. The pK accelerations lead to corrections of the orbital periods of the order of 1–100 days (Schwarzschild), 0.1–10 hr (Lense–Thirring), and 1–10{sup 3} s (quadrupole) for a target star with a = 300–800 au and e ≈ 0.8, which could be measurable with future facilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Confocal Annular Josephson Tunnel Junctions with Large Eccentricity

Confocal annular Josephson tunnel junctions (CAJTJs), which are the natural generalization of the circular annular Josephson tunnel junctions, have a rich nonlinear phenomenology due to the intrinsic non-uniformity of their planar tunnel barrier delimited by two closely spaced confocal ellipses. In the presence of a uniform magnetic field in the barrier plane, the periodically changing width of the elliptical annulus generates a asymmetric double-well for a Josephson vortex trapped in a long and narrow CAJTJ. The preparation and readout of the vortex pinned in one of the two potential minima, which are important for the possible realization of a vortex qubit, have been numerically and experimentally investigated for CAJTJs with the moderate aspect ratio 2:1. In this work, we focus on the impact of the annulus eccentricity on the properties of the vortex potential profile and study the depinning mechanism of a fluxon in more eccentric samples with aspect ratio 4:1. We also discuss the effects of the temperature-dependent losses as well as the influence of the current and magnetic noise.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

On quasiconformal maps and semilinear equations in the plane

Assume that Ω is a domain in the complex plane ℂ and A(z) is a symmetric 2×2 matrix function with measurable entries, detA = 1; and such that 1/K|ξ|{sup 2} ≤ 〈A(z)ξ, ξ〉 ≤ K|ξ|{sup 2}, ξ ∈ ℝ{sup 2}, 1 ≤ K < ∞ . In particular, for semilinear elliptic equations of the form div (A(z)∇u(z)) = f(u(z)) in Ω; we prove a factorization theorem that asserts that every weak solution u to the above equation can be expressed as the composition u = To𝜔; where 𝜔 : Ω → G stands for a K−quasiconformal homeomorphism generated by the matrix function A(z); and T(w) is a weak solution of the semilinear equation ∇T(w) = J(w)f(T(w)) in G: Here, the weight J(w) is the Jacobian of the inverse mapping 𝜔{sup −1}: Similar results hold for the corresponding nonlinear parabolic and hyperbolic equations. Some applications of these results to anisotropic media are given.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Cosmic-Ray Feedback Heating of the Intracluster Medium

Active galactic nuclei (AGNs) play a central role in solving the decades-old cooling-flow problem. Although there is consensus that AGNs provide the energy to prevent catastrophically large star formation, one major problem remains: How is the AGN energy thermalized in the intracluster medium (ICM)? We perform a suite of three-dimensional magnetohydrodynamical adaptive mesh refinement simulations of AGN feedback in a cool core cluster including cosmic rays (CRs). CRs are supplied to the ICM via collimated AGN jets and subsequently disperse in the magnetized ICM via streaming, and interact with the ICM via hadronic, Coulomb, and streaming instability heating. We find that CR transport is an essential model ingredient at least within the context of the physical model considered here. When streaming is included, (i) CRs come into contact with the ambient ICM and efficiently heat it, (ii) streaming instability heating dominates over Coulomb and hadronic heating, (iii) the AGN is variable and the atmosphere goes through low-/high-velocity dispersion cycles, and, importantly, (iv) CR pressure support in the cool core is very low and does not demonstrably violate observational constraints. However, when streaming is ignored, CR energy is not efficiently spent on the ICM heating and CR pressure builds up to a significant level, creating tension with the observations. Overall, we demonstrate that CR heating is a viable channel for the AGN energy thermalization in clusters and likely also in ellipticals, and that CRs play an important role in determining AGN intermittency and the dynamical state of cool cores.

79 ASTRONOMY AND ASTROPHYSICS↗