Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 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 OSTI2014

Selection of polynomial chaos bases via Bayesian model uncertainty methods with applications to sparse approximation of PDEs with stochastic inputs

Generalized polynomial chaos (gPC) expansions allow us to represent the solution of a stochastic system using a series of polynomial chaos basis functions. The number of gPC terms increases dramatically as the dimension of the random input variables increases. When the number of the gPC terms is larger than that of the available samples, a scenario that often occurs when the corresponding deterministic solver is computationally expensive, evaluation of the gPC expansion can be inaccurate due to over-fitting. We propose a fully Bayesian approach that allows for global recovery of the stochastic solutions, in both spatial and random domains, by coupling Bayesian model uncertainty and regularization regression methods. It allows the evaluation of the PC coefficients on a grid of spatial points, via (1) the Bayesian model average (BMA) or (2) the median probability model, and their construction as spatial functions on the spatial domain via spline interpolation. The former accounts for the model uncertainty and provides Bayes-optimal predictions; while the latter provides a sparse representation of the stochastic solutions by evaluating the expansion on a subset of dominating gPC bases. Moreover, the proposed methods quantify the importance of the gPC bases in the probabilistic sense through inclusion probabilities. We design a Markov chain Monte Carlo (MCMC) sampler that evaluates all the unknown quantities without the need of ad-hoc techniques. The proposed methods are suitable for, but not restricted to, problems whose stochastic solutions are sparse in the stochastic space with respect to the gPC bases while the deterministic solver involved is expensive. We demonstrate the accuracy and performance of the proposed methods and make comparisons with other approaches on solving elliptic SPDEs with 1-, 14- and 40-random dimensions.

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 OSTI2014

The Effect of Various Actuator-Line Modeling Approaches on Turbine-Turbine Interactions and Wake-Turbulence Statistics in Atmospheric Boundary-Layer Flow

When using an actuator-line representation of a wind turbine for computational fluid dynamics, it is common practice to volumetrically project the line force onto the flow field to create a body force in the fluid momentum equation. The objective of this study is to investigate how different volumetric projection techniques of the body force created by an actuator-line wind turbine rotor model affect the generated wake characteristics and blade loads in a turbine-turbine interaction problem. Two techniques for the body-force projection width are used , and they 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 NREL 5-MW turbines separated by seven rotor diameters is simulated within a large-eddy simulation solver subject to offshore neutral and moderately-convective atmospheric boundary-layer inflow. Power, thrust,and bending moment histories of both turbines, the statistics of angle of attack and blade loads over 2000 sec, variations in the mean and fluctuating velocity components, and turbulent kinetic energy and selected Reynolds stresses along vertical and spanwise sampling locations in the wake are analyzed. Comparisons for the different techniques of determining the body-force projection width of the actuator-line method are made and their effect on different physical quantities are assessed.

Atmospheric boundary layer, wind turbine wake↗
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 OSTI2015

PLANET FORMATION IN BINARIES: DYNAMICS OF PLANETESIMALS PERTURBED BY THE ECCENTRIC PROTOPLANETARY DISK AND THE SECONDARY

Detections of planets in eccentric, close (separations of ∼20 AU) binary systems such as α Cen or γ Cep provide an important test of planet formation theories. Gravitational perturbations from the companion are expected to excite high planetesimal eccentricities, resulting in destruction rather than growth of objects with sizes of up to several hundred kilometers in collisions of similar-sized bodies. It was recently suggested that the gravity of a massive axisymmetric gaseous disk in which planetesimals are embedded drives rapid precession of their orbits, suppressing eccentricity excitation. However, disks in binaries are themselves expected to be eccentric, leading to additional planetesimal excitation. Here we develop a secular theory of eccentricity evolution for planetesimals perturbed by the gravity of an elliptical protoplanetary disk (neglecting gas drag) and the companion. For the first time, we derive an expression for the disturbing function due to an eccentric disk, which can be used for a variety of other astrophysical problems. We obtain explicit analytical solutions for planetesimal eccentricity evolution neglecting gas drag and delineate four different regimes of dynamical excitation. We show that in systems with massive (≳ 10{sup –2} M {sub ☉}) disks, planetesimal eccentricity is usually determined by the gravity of the eccentric disk alone, and is comparable to the disk eccentricity. As a result, the latter imposes a lower limit on collisional velocities of solids, making their growth problematic. In the absence of gas drag, this fragmentation barrier can be alleviated if the gaseous disk rapidly precesses or if its own self-gravity is efficient at lowering disk eccentricity.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2013

Multibump solutions for quasilinear elliptic equations with critical growth

The current paper is concerned with constructing multibump solutions for a class of quasilinear Schrödinger equations with critical growth. This extends the classical results of Coti Zelati and Rabinowitz [Commun. Pure Appl. Math. 45, 1217–1269 (1992)] for semilinear equations as well as recent work of Liu, Wang, and Guo [J. Funct. Anal. 262, 4040–4102 (2012)] for quasilinear problems with subcritical growth. The periodicity of the potentials is used to glue ground state solutions to construct multibump bound state solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Distance correlation methods for discovering associations in large astrophysical databases

High-dimensional, large-sample astrophysical databases of galaxy clusters, such as the Chandra Deep Field South COMBO-17 database, provide measurements on many variables for thousands of galaxies and a range of redshifts. Current understanding of galaxy formation and evolution rests sensitively on relationships between different astrophysical variables; hence an ability to detect and verify associations or correlations between variables is important in astrophysical research. In this paper, we apply a recently defined statistical measure called the distance correlation coefficient, which can be used to identify new associations and correlations between astrophysical variables. The distance correlation coefficient applies to variables of any dimension, can be used to determine smaller sets of variables that provide equivalent astrophysical information, is zero only when variables are independent, and is capable of detecting nonlinear associations that are undetectable by the classical Pearson correlation coefficient. Hence, the distance correlation coefficient provides more information than the Pearson coefficient. We analyze numerous pairs of variables in the COMBO-17 database with the distance correlation method and with the maximal information coefficient. We show that the Pearson coefficient can be estimated with higher accuracy from the corresponding distance correlation coefficient than from the maximal information coefficient. For given values of the Pearson coefficient, the distance correlation method has a greater ability than the maximal information coefficient to resolve astrophysical data into highly concentrated horseshoe- or V-shapes, which enhances classification and pattern identification. These results are observed over a range of redshifts beyond the local universe and for galaxies from elliptical to spiral.

79 ASTRONOMY AND ASTROPHYSICS↗