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.

240 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Oblique ion-acoustic cnoidal waves in two temperature superthermal electrons magnetized plasma

A study is presented for the oblique propagation of ion acoustic cnoidal waves in a magnetized plasma consisting of cold ions and two temperature superthermal electrons modelled by kappa-type distributions. Using the reductive perturbation method, the nonlinear Korteweg de-Vries equation is derived, which further gives the solutions with a special type of cnoidal elliptical functions. Both compressive and rarefactive structures are found for these cnoidal waves. Nonlinear periodic cnoidal waves are explained in terms of plasma parameters depicting the Sagdeev potential and the phase curves. It is found that the density ratio of hot electrons to ions μ significantly modifies compressive/refractive wave structures. Furthermore, the combined effects of superthermality of cold and hot electrons κ{sub c},κ{sub h}, cold to hot electron temperature ratio σ, angle of propagation and ion cyclotron frequency ω{sub ci} have been studied in detail to analyze the height and width of compressive/refractive cnoidal waves. The findings in the present study could have important implications in understanding the physics of electrostatic wave structures in the Saturn's magnetosphere where two temperature superthermal electrons are present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Partial differential equations preconditioner resilient to soft and hard faults

We present a domain-decomposition-based preconditioner for the solution of partial differential equations (PDEs) that is resilient to both soft and hard faults. The algorithm reformulates the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to both soft and hard faults. This reformulation allows us to recast the problem as a set of independent tasks, and exploit data locality to reduce global communication. We discuss two different parallel implementations: (a) a single program multiple data (SPMD) version based on a one-to-one mapping between subdomain and MPI processes responsible for both state and computation; and (b) an asynchronous server–client implementation where all state information is held by the servers and clients are designed solely as computational units. We present a scalability comparison of both implementations under nominal conditions, showing efficiency within ~80% for up to 12,000 cores. We present a resilience analysis under different fault scenarios based on the server–client implementation. This framework provides resiliency to hard faults such that if a client crashes, it stops asking for work, and the servers simply distribute the work among all of the other clients alive. Erroneous subdomain solves (e.g. due to soft faults) appear as corrupted data, which is either rejected if that causes a task to fail, or is seamlessly filtered out during the regression stage through a suitable noise model. Three different types of faults are modeled: hard faults modeling nodes (or clients) crashing; soft faults occurring during the communication of the tasks between server and clients; and soft faults occurring during task execution. We demonstrate the resiliency of the approach for a 2D elliptic PDE, and explore the effect of the faults at various failure rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

On the cluster physics of Sunyaev–Zel’dovich and X-ray surveys. IV. Characterizing density and pressure clumping due to infalling substructures

Understanding the outskirts of galaxy clusters at the virial radius (R{sub 200}) and beyond is critical for an accurate determination of cluster masses, structure growth, and to ensure unbiased cosmological parameter estimates from cluster surveys. This problem has drawn renewed interest due to recent determinations of gas mass fractions beyond R{sub 200}, which appear to be considerably larger than the cosmic mean. Here, we use a large suite of cosmological hydrodynamical simulations to study the inhomogeneity of the intra-cluster medium and employ different variants of simulated physics, including radiative gas physics and thermal feedback by active galactic nuclei. We find that density and pressure clumping closely trace each other as a function of radius, but the bias on density remains on average <20% within R{sub 200}. At larger radii, clumping increases steeply due to the continuous infall of coherent structures that have not yet passed the accretion shock. Density and pressure clumping increase with cluster mass and redshift, which probes on average dynamically younger objects that are still in the process of assembling. The angular power spectra of gas density and pressure show that the clumping signal is dominated by large-scale cosmic filaments that reach from the cosmic web into the clusters, signaling the presence of gravitationally driven “super clumping.” While the prolateness of the gravitational halo potential implies an approximate radial correlation of these gaseous large-scale structures, gas density and pressure lose coherence on small scales across different radii due to dissipative gas physics. In contrast, the angular power spectrum of dark matter shows an almost uniform size distribution due to unimpeded subhalos. We provide a synopsis of the radial dependence of the clusters’ non-equilibrium measures (kinetic pressure support, ellipticity, and clumping) that all increase sharply beyond R{sub 200}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

CLASSIFICATION OF STELLAR ORBITS IN AXISYMMETRIC GALAXIES

It is known that two supermassive black holes (SMBHs) cannot merge in a spherical galaxy within a Hubble time; an emerging picture is that galaxy geometry, rotation, and large potential perturbations may usher the SMBH binary through the critical three-body scattering phase and ultimately drive the SMBH to coalesce. We explore the orbital content within an N-body model of a mildly flattened, non-rotating, SMBH-embedded elliptical galaxy. When used as the foundation for a study on the SMBH binary coalescence, the black holes bypassed the binary stalling often seen within spherical galaxies and merged on gigayear timescales. Using both frequency-mapping and angular momentum criteria, we identify a wealth of resonant orbits in the axisymmetric model, including saucers, that are absent from an otherwise identical spherical system and that can potentially interact with the binary. We quantified the set of orbits that could be scattered by the SMBH binary, and found that the axisymmetric model contained nearly six times the number of these potential loss cone orbits compared to our equivalent spherical model. In this flattened model, the mass of these orbits is more than three times that of the SMBH, which is consistent with what the SMBH binary needs to scatter to transition into the gravitational wave regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Sum-Frequency Generation from a Thin Spherical Layer: II. Analysis of Solution

We have analyzed the influence of the angle between the wave vectors of two incident waves (the opening angle) and of the type of the anisotropy of a nonlinear layer on the shape of the directivity pattern of a sum-frequency harmonic generated by the two plane elliptically polarized electromagnetic waves from a thin spherical optically nonlinear layer deposited on the surface of a dielectric spherical particle placed in a dielectric. Our analysis has shown that, if the radius of the thin spherical nonlinear layer is small, the shape of the directivity pattern will change significantly with increasing opening angle only for some types of anisotropy: the main lobes shift toward the direction that is opposite to the direction of the sum of the wave vectors of the incident waves. For three types of anisotropy, the directivity patterns have similar shapes. We have also found that, for one of the types of the anisotropy, the shape of the directivity pattern remains unchanged upon a change in the opening angle. The mathematical properties of functions describing the spatial distribution of the generated harmonic have been determined. In particular, it has been found that, upon incidence of linearly polarized waves on a thin nonlinear spherical layer that possesses either solely chiral or solely nonchiral nonlinear properties, linearly polarized radiation of the sum-frequency harmonic is generated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Second-harmonic generation in shear wave beams with different polarizations

A coupled pair of nonlinear parabolic equations was derived by Zabolotskaya [1] that model the transverse components of the particle motion in a collimated shear wave beam propagating in an isotropic elastic solid. Like the KZK equation, the parabolic equation for shear wave beams accounts consistently for the leading order effects of diffraction, viscosity and nonlinearity. The nonlinearity includes a cubic nonlinear term that is equivalent to that present in plane shear waves, as well as a quadratic nonlinear term that is unique to diffracting beams. The work by Wochner et al. [2] considered shear wave beams with translational polarizations (linear, circular and elliptical), wherein second-order nonlinear effects vanish and the leading order nonlinear effect is third-harmonic generation by the cubic nonlinearity. The purpose of the current work is to investigate the quadratic nonlinear term present in the parabolic equation for shear wave beams by considering second-harmonic generation in Gaussian beams as a second-order nonlinear effect using standard perturbation theory. In order for second-order nonlinear effects to be present, a broader class of source polarizations must be considered that includes not only the familiar translational polarizations, but also polarizations accounting for stretching, shearing and rotation of the source plane. It is found that the polarization of the second harmonic generated by the quadratic nonlinearity is not necessarily the same as the polarization of the source-frequency beam, and we are able to derive a general analytic solution for second-harmonic generation from a Gaussian source condition that gives explicitly the relationship between the polarization of the source-frequency beam and the polarization of the second harmonic.

42 ENGINEERING↗