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

Field Emission Analysis in SRF Cavities for PIP-II Using GEANT4

Field emission (FE) remains a significant hurdle for achieving optimal performance and reliability in superconducting radiofrequency (SRF) cavities used in accelerator cryomodules. A thorough understanding of the generation and propagation of FE-induced radiation is therefore essential to mitigate this problem. The absence of standardized measurement protocols further complicates the comparison of radiation data across different testing phases and facilities. This highlights the need for a precise quantitative method to diagnose and analyze FE-induced radiation. Such efforts could prove beneficial for improving cavity preparation and cleanroom assembly techniques during the prototype and production stages of Fermilab's Proton Improvement Plan-II (PIP-II) project. This study presents the initial steps of detailed Geant4 simulations aimed at analyzing FE-induced radiation in the low-beta 650 MHz 5-cell elliptical (LB650) cavity. Our goal is to combine these results with radiation diagnostics to enhance diagnostic accuracy and optimize detector positioning. This integrated approach ultimately aims to improve the preparation, assembly, and testing procedures for PIP-II SRF cavities, ensuring the delivery of FE-free cryomodules.

Wijethunga, Sajini A.K. [Fermilab] (ORCID:00000002↗
DOE OSTI2025

Field Emission Analysis in SRF Cavities for PIP-II Using GEANT4

Field emission (FE) remains a significant hurdle for achieving optimal performance and reliability in super-conducting radiofrequency (SRF) cavities used in accelerator cryomodules. A thorough understanding of the generation and propagation of FE-induced radiation is therefore essential to mitigate this problem. The absence of standardized measurement protocols further complicates the comparison of radiation data across different testing phases and facilities. This highlights the need for a precise quantitative method to diagnose and analyze FE-induced radiation. Such efforts could prove beneficial for improving cavity preparation and cleanroom assembly techniques during the prototype and production stages of Fermilab's Proton Improvement Plan-II (PIP-II) project. This study presents the initial steps of detailed Geant4 simulations aimed at analyzing FE-induced radiation in the low-beta 650 MHz 5-cell elliptical (LB650) cavity. Our goal is to combine these results with radiation diagnostics to enhance diagnostic accuracy and optimize detector positioning. This integrated approach ultimately aims to improve the preparation, assembly, and testing procedures for PIP-II SRF cavities, ensuring the delivery of FE-free cryomodules.

Wijethunga, S. A.K. [Fermilab] (ORCID:000000028765↗
DOE OSTI2015

HIERARCHICAL PROBABILISTIC INFERENCE OF COSMIC SHEAR

Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing surveys. The first part of our statistical framework requires specification of a likelihood function for the pixel data in an imaging survey given parameterized models for the galaxies in the images. We derive the lensing shear posterior by marginalizing over all intrinsic galaxy properties that contribute to the pixel data (i.e., not limited to galaxy ellipticities) and learn the distributions for the intrinsic galaxy properties via hierarchical inference with a suitably flexible conditional probabilitiy distribution specification. We use importance sampling to separate the modeling of small imaging areas from the global shear inference, thereby rendering our algorithm computationally tractable for large surveys. With simple numerical examples we demonstrate the improvements in accuracy from our importance sampling approach, as well as the significance of the conditional distribution specification for the intrinsic galaxy properties when the data are generated from an unknown number of distinct galaxy populations with different morphological characteristics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Quickest Detection for Changes in Maximal kNN Coherence of Random Matrices

This paper addresses the problem of quickest detection of a change in the maximal coherence between columns of a n p random matrix based on a sequence of matrix observations having a single unknown change point. The random matrix is assumed to have identically distributed rows and the maximal coherence is defined as the largest of the ( 2 p) correlation coefficients associated with any row. Likewise the κ nearest neighbor (kNN) coherence is defined as the κ-th largest of these correlation coefficients. The forms of the preand postchange distributions of the observed matrices are assumed to belong to the family of elliptically contoured densities with sparse dispersion matrices but are otherwise unknown. A nonparametric stopping rule is proposed that is based on the maximal kNN sample coherence between columns of each observed random matrix. This is a summary statistic that is related to a test of existence of a hub vertex in a sample correlation graph having degree at least κ. Performance bounds on the delay and false alarm performance of the proposed stopping rule are obtained in the purely high-dimensional regime where p → ∞ and n is fixed. When the prechange dispersion matrix is diagonal it is shown that, among all functions of the proposed summary statistic, the proposed stopping rule is asymptotically optimal under a minimax quickest change detection model as the stopping threshold approaches infinity. Furthermore, the theory developed also applies to sequential hypothesis testing and fixed sample size tests.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2018

Qualitative Properties of Solutions to Elliptic Equations with Nonpower Nonlinearities in ℝ{sub n}

We consider a class of anisotropic elliptic equations with nonpower nonlinearities. We establish the existence of a solution in the local Sobolev–Orlicz spaces without conditions on the data growth at infinity. We find sufficient uniqueness conditions, obtain estimates characterizing the behavior of the solution at infinity, and prove the continuous dependence of the solution on the right-hand side of the equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Explicit Weierstrass traveling wave solutions and bifurcation analysis for dissipative Zakharov–Kuznetsov modified equal width equation

We consider the dissipative Zakharov–Kuznetsov modified equal width (ZK-MEW) equation and discuss the effect of the dissipation on the existence and nature of traveling wave solutions of the equation. We use Lyapunov function and dynamical system theory in order to show that when viscosity term is added to the ZK-MEW equation, yet there exists still bounded traveling wave solutions in certain regions. Subsequently, we obtain general solution of the ZK-MEW equation in the presence and absence of viscosity in terms of Weirstrass ℘ functions and Jacobi elliptic functions. We also derive a new form of kink-type solution by exploring a factorization method based on functional transformation and the Abel’s first order nonlinear equation. Finally, we use the phase plane analysis and examine the stability of the viscous waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Carbon nanotube biconvex microcavities

Developing highly efficient microcavities with predictive narrow-band resonance frequencies using the least amount of material will allow the applications in nonlinear photonic devices. We have developed a microcavity array that comprised multi-walled carbon nanotubes (MWCNT) organized in a biconvex pattern. The finite element model allowed designing microcavity arrays with predictive transmission properties and assessing the effects of the microarray geometry. The microcavity array demonstrated negative index and produced high Q factors. 2–3 μm tall MWCNTs were patterned as biconvex microcavities, which were separated by 10 μm in an array. The microcavity was iridescent and had optical control over the diffracted elliptical patterns with a far-field pattern, whose properties were predicted by the model. It is anticipated that the MWCNT biconvex microcavities will have implications for the development of highly efficient lenses, metamaterial antennas, and photonic circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A heterogeneous stochastic FEM framework for elliptic PDEs

We introduce a new concept of sparsity for the stochastic elliptic operator −div(a(x,ω)∇(⋅)), which reflects the compactness of its inverse operator in the stochastic direction and allows for spatially heterogeneous stochastic structure. This new concept of sparsity motivates a heterogeneous stochastic finite element method (HSFEM) framework for linear elliptic equations, which discretizes the equations using the heterogeneous coupling of spatial basis with local stochastic basis to exploit the local stochastic structure of the solution space. We also provide a sampling method to construct the local stochastic basis for this framework using the randomized range finding techniques. The resulting HSFEM involves two stages and suits the multi-query setting: in the offline stage, the local stochastic structure of the solution space is identified; in the online stage, the equation can be efficiently solved for multiple forcing functions. An online error estimation and correction procedure through Monte Carlo sampling is given. Numerical results for several problems with high dimensional stochastic input are presented to demonstrate the efficiency of the HSFEM in the online stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

RISE: robust iterative surface extension for sub-nanometer X-ray mirror fabrication

Precision optics have been widely required in many advanced technological applications. X-ray mirrors, as an example, serve as the key optical components at synchrotron radiation and free electron laser facilities. They are rectangular silicon or glass substrates where a rectangular Clear Aperture (CA) needs to be polished to sub-nanometer Root Mean Squared (RMS) to keep the imaging capability of the incoming X-ray wavefront at the diffraction limit. The convolutional polishing model requires a CA to be extended with extra data, from which the dwell time is calculated via deconvolution. However, since deconvolution is very sensitive to boundary errors and noise, the existing surface extension methods can hardly fulfill the sub-nanometer requirement. On one hand, the figure errors in a CA were improperly modeled during the extension, leading to continuity issues along the boundary. On the other hand, uncorrectable high-frequency errors and noise were also extended. In this study, we propose a novel Robust Iterative Surface Extension (RISE) method that resolves these problems with a data fitting strategy. RISE models the figure errors in a CA with orthogonal polynomials and ensures that only correctable errors are fit and extended. Combined with boundary conditions, an iterative refinement of dwell time is then proposed to compensate the errors brought by the extension and deconvolution, which drastically reduces the estimated figure error residuals in a CA while the increase of total dwell time is negligible. To our best knowledge, RISE is the first data fitting-based surface extension method and is the first to optimize dwell time based on iterative extension. An experimental verification of RISE is given by fabricating two elliptic cylinders (10 mm × 80 mm CAs) starting from a sphere with a radius of curvature around 173 m using ion beam figuring. The figure errors in the two CAs greatly improved from 204.96 nm RMS and 190.28 nm RMS to 0.62 nm RMS and 0.71 nm RMS, respectively, which proves that RISE is an effective method for sub-nanometer level X-ray mirror fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Theoretical study of symmetry of flux onto a capsule

An analytic model to describe the flux asymmetry onto a capsule based on the viewfactor approximation is developed and verified with numerical simulations. By using a nested spheres technique to represent the various sources of flux asymmetry, the model can treat spherically and cylindrically symmetric hohlraums, e.g., cylinder, elliptic, and rugby. This approach includes the more realistic case of frequency-dependent flux asymmetry compared with the more standard frequency-integrated or single-frequency approaches [D. W. Phillion and S. M. Pollaine, Phys. Plasmas 1, 2963 (1994)]. Correspondingly, the approach can be used to assess x-ray preheat asymmetry generated from localized laser absorption in the high-Z hohlraum wall. For spherical hohlraums with 4, 6, or 8 laser entrance holes (LEHs), an optimal configuration of LEHs, laser spot placement, and angle-of-incidence of the single-ringed laser beams is defined. An analogy between minimizing the flux asymmetry onto a capsule and the Thomson problem of point charge placement on a sphere for minimized energy is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

On the resolvent of multidimensional operators with frequently changing boundary conditions in the case of the homogenized Dirichlet condition

We consider an elliptic operator in a multidimensional domain with frequently changing boundary conditions in the case when the homogenized operator contains the Dirichlet boundary condition. We prove the uniform resolvent convergence of the perturbed operator to the homogenized operator and obtain estimates for the rate of convergence. A complete asymptotic expansion is constructed for the resolvent when it acts on sufficiently smooth functions. Bibliography: 41 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Upgrade of the Fermilab Spoke Test Cryostat for Testing of PIP-II 650 MHz 5-Cell Elliptical Cavities

Design of the high beta 650 MHz prototype cryomodule for PIP-II is currently undergoing at Fermilab. The cryomodule includes six 5-cell elliptical SRF cavities with accelerating voltage up to 20 MV and low heat dissipation (Q₀ > 3·10^{z}EhNZeHn). Characterization of performance of fully integrated jacketed cavities with high power coupler and tuner is crucial for the project. Such a characterization of jacketed cavity requires a horizontal test cryostat. Existing horizontal testing facilities at Fermilab, Horizontal Test Stand (HTS) and Spoke Test Cryostat (STC), are not large enough to accommodate jacketed 650 MHz 5-cell cavity. An upgrade of the STC is proposed to install extension to the cryostat and modify cryogenic connections and RF infrastructure to provide testing of 650 MHz cavities. In this paper we describe STC upgrade and commissioning of the upgraded facility. We discuss mitigation of issues and problems specific for testing of high Q₀ 650 MHz cavities, which require low residual magnetic field and low acoustic and mechanical vibrations environment.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

A fast Huygens sweeping method for capturing paraxial multi-color optical self-focusing in nematic liquid crystals

We propose a numerically efficient algorithm for simulating the multi-color optical self-focusing phenomena in nematic liquid crystals. The propagation of the nematicon is modeled by a parabolic wave equation coupled with a nonlinear elliptic partial differential equation governing the angle between the crystal and the direction of propagation. Numerically, the paraxial parabolic wave equation is solved by a fast Huygens sweeping method, while the nonlinear elliptic PDE is handled by the alternating direction explicit (ADE) method. The overall algorithm is shown to be numerically efficient for computing high frequency beam propagations. - Highlights: • A simple strategy can already improve the fast Huygens sweeping method for high frequency asymptotic solutions. • The Alternating Direction Explicit method can significant improve the computational efficiency. • Numerical solutions match the Fredericks transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Local parametric instability near elliptic points in vortex flows under shear deformation

The dynamics of two point vortices embedded in an oscillatory external flow consisted of shear and rotational components is addressed. The region associated with steady-state elliptic points of the vortex motion is established to experience local parametric instability. The instability forces the point vortices with initial positions corresponding to the steady-state elliptic points to move in spiral-like divergent trajectories. This divergent motion continues until the nonlinear effects suppress their motion near the region associated with the steady-state separatrices. The local parametric instability is then demonstrated not to contribute considerably to enhancing the size of the chaotic motion regions. Instead, the size of the chaotic motion region mostly depends on overlaps of the nonlinear resonances emerging in the perturbed system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Reconstructing high-dimensional Hilbert-valued functions via compressed sensing

We present and analyze a novel sparse polynomial technique for approximating high-dimensional Hilbert-valued functions, with application to parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our theoretical framework treats the function approximation problem as a joint sparse recovery problem, where the set of jointly sparse vectors is possibly infinite. To achieve the simultaneous reconstruction of Hilbert-valued functions in both parametric domain and Hilbert space, we propose a novel mixed-norm based `1 regularization method that exploits both energy and sparsity. Our approach requires extensions of concepts such as the restricted isometry and null space properties, allowing us to prove recovery guarantees for sparse Hilbert-valued function reconstruction. We complement the enclosed theory with an algorithm for Hilbert-valued recovery, based on standard forward-backward algorithm, meanwhile establishing its strong convergence in the considered infinite-dimensional setting. Finally, we demonstrate the minimal sample complexity requirements of our approach, relative to other popular methods, with numerical experiments approximating the solutions of high-dimensional parameterized elliptic PDEs.

Dexter, Nick↗
DOE OSTI2021

Hierarchical Inference with Bayesian Neural Networks: An Application to Strong Gravitational Lensing

In the past few years, approximate Bayesian Neural Networks (BNNs) have demonstrated the ability to produce statistically consistent posteriors on a wide range of inference problems at unprecedented speed and scale. However, any disconnect between training sets and the distribution of real-world objects can introduce bias when BNNs are applied to data. This is a common challenge in astrophysics and cosmology, where the unknown distribution of objects in our universe is often the science goal. In this work, we incorporate BNNs with flexible posterior parameterizations into a hierarchical inference framework that allows for the reconstruction of population hyperparameters and removes the bias introduced by the training distribution. We focus on the challenge of producing posterior PDFs for strong gravitational lens mass model parameters given Hubble Space Telescope–quality single-filter, lens-subtracted, synthetic imaging data. We show that the posterior PDFs are sufficiently accurate (statistically consistent with the truth) across a wide variety of power-law elliptical lens mass distributions. We then apply our approach to test data sets whose lens parameters are drawn from distributions that are drastically different from the training set. We show that our hierarchical inference framework mitigates the bias introduced by an unrepresentative training set's interim prior. Simultaneously, we can precisely reconstruct the population hyperparameters governing our test distributions. Our full pipeline, from training to hierarchical inference on thousands of lenses, can be run in a day. The framework presented here will allow us to efficiently exploit the full constraining power of future ground- and space-based surveys (https://github.com/swagnercarena/ovejero).

79 ASTRONOMY AND ASTROPHYSICS↗