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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Exploring the Interplay of Resilience and Energy Consumption for a Task-Based Partial Differential Equations Preconditioner

We discuss algorithm-based resilience to silent data corruption (SDC) in a task- based domain-decomposition preconditioner for partial differential equations (PDEs). The algorithm exploits a reformulation of the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to SDC. The implementation is based on a server-client model where all state information is held by the servers, while clients are designed solely as computational units. Scalability tests run up to ~ 51 K cores show a parallel efficiency greater than 90%. We use a 2D elliptic PDE and a fault model based on random single bit-flip to demonstrate the resilience of the application to synthetically injected SDC. We discuss two fault scenarios: one based on the corruption of all data of a target task, and the other involving the corruption of a single data point. We show that for our application, given the test problem considered, a four-fold increase in the number of faults only yields a 2% change in the overhead to overcome their presence, from 7% to 9%. We then discuss potential savings in energy consumption via dynamics voltage/frequency scaling, and its interplay with fault-rates, and application overhead.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Exploring the interplay of resilience and energy consumption for a task-based partial differential equations preconditioner

Here, we discuss algorithm-based resilience to silent data corruptions (SDCs) in a task-based domain-decomposition preconditioner for partial differential equations (PDEs). The algorithm exploits a reformulation of the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to SDCs. The implementation is based on a server-client model where all state information is held by the servers, while clients are designed solely as computational units. Scalability tests run up to ~51K cores show a parallel efficiency greater than 90%. We use a 2D elliptic PDE and a fault model based on random single and double bit-flip to demonstrate the resilience of the application to synthetically injected SDC. We discuss two fault scenarios: one based on the corruption of all data of a target task, and the other involving the corruption of a single data point. We show that for our application, given the test problem considered, a four-fold increase in the number of faults only yields a 2% change in the overhead to overcome their presence, from 7% to 9%. We then discuss potential savings in energy consumption via dynamic voltage/frequency scaling, and its interplay with fault-rates, and application overhead.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Trigonometric and hyperbolic functions method for constructing analytic solutions to nonlinear plane magnetohydrodynamics equilibrium equations

In this paper, we present a solution method for constructing exact analytic solutions to magnetohydrodynamics (MHD) equations. The method is constructed via all the trigonometric and hyperbolic functions. The method is applied to MHD equilibria with mass flow. Applications to a solar system concerned with the properties of coronal mass ejections that affect the heliosphere are presented. Some examples of the constructed solutions which describe magnetic structures of solar eruptions are investigated. Moreover, the constructed method can be applied to a variety classes of elliptic partial differential equations which arise in plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A DLM/FD/IB method for simulating compound vesicle motion under creeping flow condition

In this article we present first a new distributed Lagrange multiplier/fictitious domain (DLM/FD) method for simulating fluid–particle interaction in Stokes flow. A conjugate gradient method driven by both pressure and distributed Lagrange multiplier, called one-shot method, has been developed to solve the discrete Stokes problem while enforcing the rigid body motion within the region occupied by the particle. The methodology is validated by comparing the numerical results of a neutrally buoyant particle of either a circular or elliptic shape with the associated Jeffery's solutions. We have successively combined the above methodology with an immersed boundary (IB) method and an elastic membrane modeled by a spring network to simulate the dynamics of a compound vesicle. In simple shear flow under creeping flow condition, the results are consistent with those obtained in literature. In Poiseuille flow, the compound vesicle motion is dominated by the motion of the vesicle membrane as expected and stays in the central region of the channel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mixed gauge-global symmetries, elliptic modes, and black hole thermodynamics in Hořava-Lifshitz gravity

In Hořava-Lifshitz gravity, a putative consistent theory of quantum gravity for which there is evidence for both black hole thermodynamics and a holographic construction, spacetime is endowed with a preferred dynamical spacelike foliation. The theory has a leaf reparameterization symmetry that is neither global nor local gauge, hyperbolic and elliptic equations of motion, a lack of splittability, and universal horizon black hole solutions. The reparameterization symmetry is “mixed”: it is a local symmetry in one coordinate yet global on each leaf. More broadly it is an example of both unfree and projectable gauge symmetries. The mixed symmetry and associated charge has not yet been accounted for in calculations of universal horizon thermodynamics in Hořava-Lifshitz gravity. This has led to problems, in particular the failure of the first law in a class of asymptotically AdS solutions where the normal to the leaves of the foliation is not aligned with the time translation Killing vector at infinity. We show how the dynamics of the charge corresponding to this symmetry coupled with the other features above resolves this issue. We then briefly comment how this mixed symmetry, the corresponding charge, and the elliptic equations of motion also conspire to evade recent holographic arguments for only local gauge fields in consistent theories of quantum gravity due to the lack of splittability of the elliptic equation and associated mode.

Global Symmetries↗
DOE OSTI2016

Blade Load Unsteadiness and Turbulence Statistics in an Actuator-Line Computed Turbine–Turbine Interaction Problem

The objective of this study is to investigate how different volumetric projection techniques used in actuator-line modeling affect the unsteady blade loads and wake turbulence statistics. The two techniques for the body-force projection radius are based on either (i) the grid spacing or (ii) the combination of grid spacing and an equivalent elliptic blade planform. An array of two National Renewable Energy Laboratory 5-MW turbines separated by seven rotor diameters is simulated for 2000s (about rotor 300 revolutions) within a large-eddy simulation (LES) solver of the neutral and moderately convective atmospheric boundary layer (ABL). The statistics of sectional angle of attack (AOA), blade loads, and turbine power histories are quantified. Moreover, the degree of unsteadiness of sectional blade loads in response to atmospheric and wake turbulence is computed via a reduced frequency based on the rate-of-change in sectional AOA. Finally, the goal of this work is to make the wind energy community aware of the uncertainties associated with actuator-line modeling approaches.

17 WIND ENERGY↗
DOE OSTI2020

Statistical Treatment of Inverse Problems Constrained by Differential Equations-based Models with Stochastic Terms

This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations. In both cases we show that the variogram and the hybrid scores produce better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Summary Report: Multigrid for Systems of Elliptic PDEs

We are interested in determining if multigrid can be effectively applied to the system. The conclusion that I seem to be drawn to is that it is impossible to develop a blackbox multigrid solver for these general systems. Analysis of the system of PDEs must be conducted first to determine pre-processing procedures on the continuous problem before applying a multigrid method. Determining this pre-processing is currently not incorporated in black-box multigrid strategies. Nevertheless, we characterize some system features that will make the system more amenable to multigrid approaches, techniques that may lead to more amenable systems, and multigrid procedures that are generally more appropriate for these systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Sum-Frequency Generation from a Thin Spherical Layer: I. Analytical Solution

In the tensor form, using the generalized Rayleigh–Gans–Debye approximation, the problem of the sum-frequency generation from a thin nonlinear layer deposited on a dielectric spherical particle placed in a dielectric medium is solved. The second-order nonlinear dielectric susceptibility tensor is chosen in a general form containing chiral components. In the vector and tensor forms, expressions are obtained that describe the spatial distribution of the sum-frequency radiation field generated by two plane electromagnetic elliptically polarized waves. Limiting expressions that describe the spatial distribution of the sum-frequency harmonic at small and large radii of the spherical layer are obtained. It is revealed that, at small radii of the spherical layer, the radiation due to the chiral anisotropy coefficients makes a dominant contribution to the generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Energy levels of a quantum particle on a cylindrical surface with non-circular cross-section in electric and magnetic fields

We calculate the energy levels of a quantum particle on a cylindrical surface with non-circular cross-section in uniform electric and magnetic fields. Using separation of variables method and a change of independent variable, we show that the problem can be reduced to a one-dimensional Schrödinger equation for a periodic potential. The effects of varying the shape of the cross-section while keeping the same perimeter and the strengths of the electric and magnetic fields are investigated for elliptical, corrugated, and nearly-rectangular tubes with radial dimensions of the order of a nanometer. The geometric potential has minima at the angular positions where there is a significant amount of curvature. For the elliptical and corrugated tubes, it is shown that as the tube departs from the circular shape of cross-section the double-degeneracy between the energy levels is lifted. For the nearly-rectangular tube, it is shown that energy level crossings occur as the horizontal dimension of the tube is varied while keeping the same perimeter and radius of circular corners. The interplay between the curvature and the strength of the electric and magnetic fields determines the overall behavior of the energy levels. As the strength of the electric field increases, the overall potential gets skewed creating a potential well on the side corresponding to the more negative electric potential. The energy levels of the first few excited states approach more positive values while the ground state energy level approaches a more negative value. For large electric fields, all bound state energy levels tend to more negative values. The contribution of weak magnetic fields to the overall potential behaves in the same way as the electric field contribution but with its sign depending on the direction of the component of the momentum parallel to the cylindrical axis. Large magnetic fields lead to pairing of energy levels reminiscent of 2D Landau levels for the elliptical and nearly-rectangular tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Image feature extraction and galaxy classification: a novel and efficient approach with automated machine learning

ABSTRACT In this work, we explore the possibility of applying machine learning methods designed for 1D problems to the task of galaxy image classification. The algorithms used for image classification typically rely on multiple costly steps, such as the point spread function deconvolution and the training and application of complex Convolutional Neural Networks of thousands or even millions of parameters. In our approach, we extract features from the galaxy images by analysing the elliptical isophotes in their light distribution and collect the information in a sequence. The sequences obtained with this method present definite features allowing a direct distinction between galaxy types. Then, we train and classify the sequences with machine learning algorithms, designed through the platform Modulos AutoML. As a demonstration of this method, we use the second public release of the Dark Energy Survey (DES DR2). We show that we are able to successfully distinguish between early-type and late-type galaxies, for images with signal-to-noise ratio greater than 300. This yields an accuracy of $86{{\ \rm per\ cent}}$ for the early-type galaxies and $93{{\ \rm per\ cent}}$ for the late-type galaxies, which is on par with most contemporary automated image classification approaches. The data dimensionality reduction of our novel method implies a significant lowering in computational cost of classification. In the perspective of future data sets obtained with e.g. Euclid and the Vera Rubin Observatory, this work represents a path towards using a well-tested and widely used platform from industry in efficiently tackling galaxy classification problems at the peta-byte scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Rigid-Mode Limit of the Yokoya Matrix Formalism and the Burov-Lebedev Dispersion Equation

Transverse single-bunch instabilities of space-charge-dominated coasting beams with round and flat transverse geometries are studied using a unified dispersion-relation framework. The analysis combines the Burov-Lebedev formalism, which captures space-charge tune spread, Landau damping, and instability threshold behavior, with Yokoya’s projection method for representing coherent transverse mode structure and its dependence on beam aspect ratio. In the rigid-beam limit, the formulation reduces to a scalar dispersion relation of Burov-Lebedev paper. For non-rigid transverse oscillations, truncation of Yokoya’s Hermite-based expansion yields a finite-dimensional matrix eigenvalue problem in which space-charge and coupling impedance effects enter through Burov-Lebedev–type denominators. This approach provides a consistent basis for comparing rigid and non-rigid instability behavior in round and flat beams and for assessing the role of beam ellipticity in modifying coherent mode structure and stability thresholds.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Novel techniques and devices for in-situ film coatings of long, small diameter tubes or elliptical and other surface contours

In this study, devices and techniques that can, via physical vapor deposition,coat various surface contours or very long small aperture pipes, are described. Recently, a magnetron mole was developed in order to in-situ coat accelerator tube sections of the Brookhaven National Lab relativistic heavy ion collider that have 7.1 cm diameter with access points that are 500 m apart, for copper coat the accelerator vacuum tube in order to alleviate the problems of unacceptable ohmic heating and of electron clouds. A magnetron with a 50 cm long cathode was designed fabricated and successfully operated to copper coat a whole assembly containing a full-size, stainless steel, cold bore, of the accelerator magnet tubing connected to two types bellows, to which two additional pipes made of accelerator tubing were connected. The magnetron is mounted on a carriage with spring loaded wheels that successfully crossed bellows and adjusted for variations in vacuum tube diameter, while keeping the magnetron centered. Electrical power and cooling water were fed through a cable bundle. The umbilical cabling system, which is enclosed in a flexible braided metal sleeve, is driven by a motorized spool. To increase cathode lifetime, movable magnet package was developed, and thickest possible cathode was made, with a rather challenging target to substrate distance of less than 1.5 cm. Optimized process to ensure excellent adhesion was developed. Coating thickness of 10 μm Cu passed all industrial tests and even exceeded maximum capability of a 12 kg pull test fixture. Room temperature radio frequency (RF) resistivity measurement indicated that 10 μm Cu coated stainless steel accelerator tube has conductivity close to copper tubing. Work is in progress to repeat the RF resistivity measurement at cryogenic temperatures. Over 20 years ago, a device using multi axis robotic manipulators controlling separate robotic assemblies resulted in nine-axes of motion combined with conformal shape of the cathodes that can adapt to various curved surface contours was developed and successfully used for depositing optical coating on aircraft canopies. The techniques can be utilized for in situ coating of elliptical and other surface contour RF cavities and long beam pipes with thick superconducting films. Plans are to incorporate ion assisted deposition in those techniques for attaining dense, adherent and defect free coatings.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗
DOE OSTI2015

Interaction of elliptically polarised cross-degenerate cnoidal waves in an isotropic gyrotropic medium with spatial dispersion of cubic nonlinearity

Three unusual classes of particular analytical solutions to a system of four nonlinear equations are found for slowly varying complex amplitudes of circularly polarised components of the electric field. The system describes the self-action and interaction of two elliptically polarised plane waves collinearly propagating in an isotropic medium with second-order frequency dispersion and spatial dispersion of cubic nonlinearity. The solutions correspond to self-consistent combinations of two elliptically polarised cnoidal waves whose mutually orthogonal polarisation components vary in accordance with pairwise identical laws during propagation. At the same time, the amplitudes of the component with the same circular polarisation are proportional to two different elliptic Jacobi functions with the same periods. (nonlinear optical phenomena)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Image segmentation for analyzing galaxy-galaxy strong lensing systems

The goal of this Letter is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17–25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as 10 8.5 M ⊙ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.

79 ASTRONOMY AND ASTROPHYSICS↗