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 289 records · Page 16

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↗

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↗

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↗

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↗

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↗

Non-Grassmann mechanical model of the Dirac equation

We construct a new example of the spinning-particle model without Grassmann variables. The spin degrees of freedom are described on the base of an inner anti-de Sitter space. This produces both {Gamma}{sup {mu}} and {Gamma}{sup {mu}{nu}}-matrices in the course of quantization. Canonical quantization of the model implies the Dirac equation. We present the detailed analysis of both the Lagrangian and the Hamiltonian formulations of the model and obtain the general solution to the classical equations of motion. Comparing Zitterbewegung of the spatial coordinate with the evolution of spin, we ask on the possibility of space-time interpretation for the inner space of spin. We enumerate similarities between our analogous model of the Dirac equation and the two-body system subject to confining potential which admits only the elliptic orbits of the order of de Broglie wavelength. The Dirac equation dictates the perpendicularity of the elliptic orbits to the direction of center-of-mass motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Exact solitary solution of Schamel equation in plasmas with negative ions

A theoretical investigation is carried out for understanding the properties of the solitary solution in plasmas with negative ions. Schamel equation for a plasma consisting of electron, positive ions and negative ions has been derived by using the reductive perturbation method. The effects of negative ions and the density on the properties of the solitary solution is discussed. We make use of the extended mapping method and auxiliary equation to obtain the solution of Schamel equation. This solution includes the Jacobi elliptic function solutions, hyperbolic function solutions, rational solutions, and periodic wave solutions. Furthermore, we show that the incorporate negative ion effects in the nonlinear propagation of ion acoustic waves that are controlled by trapped electrons and the results of the solitary solution in plasmas with negative ions model are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NONLINEAR COLOR-METALLICITY RELATIONS OF GLOBULAR CLUSTERS. III. ON THE DISCREPANCY IN METALLICITY BETWEEN GLOBULAR CLUSTER SYSTEMS AND THEIR PARENT ELLIPTICAL GALAXIES

One of the conundrums in extragalactic astronomy is the discrepancy in observed metallicity distribution functions (MDFs) between the two prime stellar components of early-type galaxies-globular clusters (GCs) and halo field stars. This is generally taken as evidence of highly decoupled evolutionary histories between GC systems and their parent galaxies. Here we show, however, that new developments in linking the observed GC colors to their intrinsic metallicities suggest nonlinear color-to-metallicity conversions, which translate observed color distributions into strongly peaked, unimodal MDFs with broad metal-poor tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved field stars in nearby elliptical galaxies and those produced by chemical evolution models of galaxies. The GC MDF shape, characterized by a sharp peak with a metal-poor tail, indicates a virtually continuous chemical enrichment with a relatively short timescale. The characteristic shape emerges across three orders of magnitude in the host galaxy mass, suggesting a universal process of chemical enrichment among various GC systems. Given that GCs are bluer than field stars within the same galaxy, it is plausible that the chemical enrichment processes of GCs ceased somewhat earlier than that of the field stellar population, and if so, GCs preferentially trace the major, vigorous mode of star formation events in galactic formation. We further suggest a possible systematic age difference among GC systems, in that the GC systems in more luminous galaxies are older. This is consistent with the downsizing paradigm whereby stars of brighter galaxies, on average, formed earlier than those of dimmer galaxies; this additionally supports the similar nature shared by GCs and field stars. Although the sample used in this study (the Hubble Space Telescope Advanced Camera for Surveys/Wide Field Channel, WFPC2, and WFC3 photometry for the GC systems in the Virgo galaxy cluster) confines our discussion to R {approx}< R{sub e} for giant ellipticals and {approx}<10 R{sub e} for normal ellipticals, our findings suggest that GC systems and their parent galaxies have shared a more common origin than previously thought, and hence greatly simplify theories of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Interoperable Technologies for Advanced Petascale Simulations

Our final report on the accomplishments of ITAPS at Stony Brook during period covered by the research award includes component service, interface service and applications. On the component service, we have designed and implemented a robust functionality for the Lagrangian tracking of dynamic interface. We have migrated the hyperbolic, parabolic and elliptic solver from stage-wise second order toward global second order schemes. We have implemented high order coupling between interface propagation and interior PDE solvers. On the interface service, we have constructed the FronTier application programer's interface (API) and its manual page using doxygen. We installed the FronTier functional interface to conform with the ITAPS specifications, especially the iMesh and iMeshP interfaces. On applications, we have implemented deposition and dissolution models with flow and implemented the two-reactant model for a more realistic precipitation at the pore level and its coupling with Darcy level model. We have continued our support to the study of fluid mixing problem for problems in inertial confinement fusion. We have continued our support to the MHD model and its application to plasma liner implosion in fusion confinement. We have simulated a step in the reprocessing and separation of spent fuels from nuclear power plant fuel rods. We have implemented the fluid-structure interaction for 3D windmill and parachute simulations. We have continued our collaboration with PNNL, BNL, LANL, ORNL, and other SciDAC institutions.

97 MATHEMATICS AND COMPUTING↗

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↗

ORIGIN AND GROWTH OF NUCLEAR STAR CLUSTERS AROUND MASSIVE BLACK HOLES

The centers of stellar spheroids less luminous than {approx}10{sup 10} L {sub Sun} are often marked by the presence of nucleated central regions, called 'nuclear star clusters' (NSCs). The origin of NSCs is still unclear. Here we investigate the possibility that NSCs originate from the migration and merger of stellar clusters at the center of galaxies where a massive black hole (MBH) may sit. We show that the observed scaling relation between NSC masses and the velocity dispersion of their host spheroids cannot be reconciled with a purely 'in situ' dissipative formation scenario. On the other hand, the observed relation appears to be in agreement with the predictions of the cluster merger model. A dissipationless formation model also reproduces the observed relation between the size of NSCs and their total luminosity, R{proportional_to}{radical}(L{sub NSC}). When an MBH is included at the center of the galaxy, such dependence becomes substantially weaker than the observed correlation, since the size of the NSC is mainly determined by the fixed tidal field of the MBH. We evolve through dynamical friction a population of stellar clusters in a model of a galactic bulge taking into account dynamical dissolution due to two-body relaxation, starting from a power-law cluster initial mass function and adopting an initial total mass in stellar clusters consistent with the present-day cluster formation efficiency of the Milky Way (MW). The most massive clusters reach the center of the galaxy and merge to form a compact nucleus; after 10{sup 10} years, the resulting NSC has properties that are consistent with the observed distribution of stars in the MW NSC. When an MBH is included at the center of a galaxy, globular clusters are tidally disrupted during inspiral, resulting in NSCs with lower densities than those of NSCs forming in galaxies with no MBHs. We suggest this as a possible explanation for the lack of NSCs in galaxies containing MBHs more massive than {approx}10{sup 8} M {sub Sun }. Finally, we investigate the orbital evolution of globular clusters in giant elliptical galaxies which are believed to always host an MBH at their center rather than an NSC. In these systems an additional mechanism can prevent an NSC from forming: the time for globular clusters to reach the center of the galaxy is much longer than the Hubble time.

79 ASTRONOMY AND ASTROPHYSICS↗

Breatherlike electromagnetic wave propagation in an antiferromagnetic medium with Dzyaloshinsky-Moriya interaction

We investigate the nature of propagation of electromagnetic waves (EMWs) in an antiferromagnetic medium with Dzyaloshinsky-Moriya (DM) interaction environment. The interplay of bilinear and DM exchange spin coupling with the magnetic field component of the EMW has been studied by solving Maxwell's equations coupled with a nonlinear spin equation for the magnetization of the medium. We made a nonuniform expansion of the magnetization and magnetic field along the direction of propagation of EMW, in the framework of reductive perturbation method, and the dynamics of the system is found to be governed by a generalized derivative nonlinear Schroedinger (DNLS) equation. We employ the Jacobi-elliptic function method to solve the DNLS equation, and the electromagnetic wave propagation in an antiferromagnetic medium is governed by the breatherlike spatially and temporally coherent localized modes under the influence of DM interaction parameter.

36 MATERIALS SCIENCE↗

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↗

Nonlinear periodic waves solutions of the nonlinear self-dual network equations

The new classes of periodic solutions of nonlinear self-dual network equations describing the breather and soliton lattices, expressed in terms of the Jacobi elliptic functions have been obtained. The dependences of the frequencies on energy have been found. Numerical simulations of soliton lattice demonstrate their stability in the ideal lattice and the breather lattice instability in the dissipative lattice. However, the lifetime of such structures in the dissipative lattice can be extended through the application of ac driving terms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗