SEARCH · Engineering Papers

Results for “asymptotic analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A comparison of locally adaptive multigrid methods: LDC, FAC and FIC

This study is devoted to a comparative analysis of three 'Adaptive ZOOM' (ZOom Overlapping Multi-level) methods based on similar concepts of hierarchical multigrid local refinement: LDC (Local Defect Correction), FAC (Fast Adaptive Composite), and FIC (Flux Interface Correction)--which we proposed recently. These methods are tested on two examples of a bidimensional elliptic problem. We compare, for V-cycle procedures, the asymptotic evolution of the global error evaluated by discrete norms, the corresponding local errors, and the convergence rates of these algorithms.

Khadra, Khodor↗

DOE OSTI2024

Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P∝ k m , where k is the wavenumber. The hydrodynamic condition is the same as in the $\theta$-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The $\theta$-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ~1.86 shock and gamma-law = 5/3 fluids with Atwood number A = 0.5. The bubble amplitude h B is found to grow in two stages. Initially, h B undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, h B grows as a power law $~t^{\theta}$ when $\langle{k}\rangle{|h_B|}$ > O(1), where $\langle{k}\rangle$ is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Quantum Foundations of Classical Reversible Computing

The reversible computation paradigm aims to provide a new foundation for general classical digital computing that is capable of circumventing the thermodynamic limits to the energy efficiency of the conventional, non-reversible digital paradigm. However, to date, the essential rationale for, and analysis of, classical reversible computing (RC) has not yet been expressed in terms that leverage the modern formal methods of non-equilibrium quantum thermodynamics (NEQT). In this paper, we begin developing an NEQT-based foundation for the physics of reversible computing. We use the framework of Gorini-Kossakowski-Sudarshan-Lindblad dynamics (a.k.a. Lindbladians) with multiple asymptotic states, incorporating recent results from resource theory, full counting statistics and stochastic thermodynamics. Important conclusions include that, as expected: (1) Landauer’s Principle indeed sets a strict lower bound on entropy generation in traditional non-reversible architectures for deterministic computing machines when we account for the loss of correlations; and (2) implementations of the alternative reversible computation paradigm can potentially avoid such losses, and thereby circumvent the Landauer limit, potentially allowing the efficiency of future digital computing technologies to continue improving indefinitely. We also outline a research plan for identifying the fundamental minimum energy dissipation of reversible computing machines as a function of speed.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Optimization of Orifice Geometry for Cross-Flow Mixing in a Cylindrical Duct

Mixing of gaseous jets in a cross-flow has significant applications in engineering, one example of which is the dilution zone of a gas turbine combustor. Despite years of study, the design of jet injection in combustors is largely based on practical experience. A series of experiments was undertaken to delineate the optimal mixer orifice geometry. A cross-flow to core-flow momentum-flux ratio of 40 and a mass flow ratio of 2.5 were selected as representative of an advanced design. An experimental test matrix was designed around three variables: the number of orifices, the orifice aspect ratio (long-to-short dimension), and the orifice angle. A regression analysis was performed on the data to arrive at an interpolating equation that predicted the mixing performance of orifice geometry combinations within the range of the test matrix parameters. Results indicate that mixture uniformity is a non-linear function of the number of orifices, the orifice aspect ratio, and the orifice angle. Optimum mixing occurs when the asymptotic mean jet trajectories are in the range of 0.35 less than r/R less than 0.5 (where r = 0 is at the mixer wall) at z/R = 1.0. At the optimum number of orifices, the difference between shallow-angled slots with large aspect ratios and round holes is minimal and either approach will lead to good mixing performance. At the optimum number of orifices, it appears possible to have two local optimums where one corresponds to an aspect ratio of 1.0 and the other to a high aspect ratio.

Sowa, W. A.↗

DOE OSTI2024

Asymptotic behavior of the Fourier coefficients and the analytic structure of the QCD equation of state

In this paper we study the universal properties of the baryon chemical potential Fourier coefficients in quantum chromodynamics. We show that by following a well-defined strategy, the Fourier coefficients can be used to locate Yang-Lee edge singularities associated with chiral phase transition (and by extension with the Roberge-Weiss) in the complex chemical potential plane. We comment on the viability of performing this analysis using lattice QCD data. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

NASA NTRS1995

Semianalytic continuum spectra of Type 2 supernovae

We extend the approximate radiative transfer analysis of Hershkowitz, Linder, & Wagoner (1986) to a more general class of supernova model atmospheres, using a simple fit to the effective continuum opacity produced by lines (Wagoner, Perez, & Vasu 1991). At the low densities considered, the populations of the excited states of hydrogen are governed mainly by photoionization and recombination, and scattering dominates absorptive opacity. We match the asymptotic expressions for the spectral energy density J(sub nu) at the photosphere, whose location at each frequency is determined by a first-order calculation of the deviation of J(sub nu) from the Planck function B(sub nu). The emergent spectral luminosity then assumes the form L(sub nu) = 4 pi(squared)r(squared)(sub *) zeta(squared)B(sub nu)(T(sub p)), where T(sub p)(nu) is the photospheric temperature zeta is the dilution factor, and r(sub *) is a fiducial radius (ultimately taken to be the photospheric radius r(sub p)(nu)). The atmosphere is characterized by an effective temperature T(sub e) (varies as L(sup 1/4)r(sup -1/2)(sub *)) and hydrogen density n(sub H) = dependence of zeta on frequency nu and the parameters T(sub p), r(sub p), and alpha. The resulting understanding of the dependence of the spectral luminosity on observable parameters which characterize the relevant physical conditions will be of particular use in assessing the reliability of the expanding photosphere method of distance determination. This is particularly important at cosmological distances, where no information about the progenitor star will be available. This technique can also be applied to other low-density photosphere.

Montes, Marcos J.↗

DOE OSTI2020

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Resummation of QED radiative corrections in a strong constant crossed field

By considering radiative corrections of up to 3rd-loop order, Ritus and Narozhny conjectured that the proper expansion parameter for QED in a strong constant crossed field is $g=\alpha\chi^{2/3}$, where the dynamical quantum parameter $\chi=e\sqrt{-(Fp)^2}/m^3$ combines the particle momentum $p$ with the external field strength tensor $F$. Here we present and discuss the first non-perturbative result in this context, the resummed bubble-type polarization corrections to the electron self-energy in a constant crossed field. Our analysis confirms the relevance of the scaling parameter $g$ to the enhancement of bubble-type radiative corrections. This parameter actually represents the characteristic value of the ratio of the 1-loop polarization bubble to the photon virtuality. After an all-order resummation we identify and discuss two contributions to the self-energy with different formation regions and asymptotic behavior for $g\gg1$. Whereas the breakdown of perturbation theory occurs already for $g\gtrsim1$, the leading-order result remains dominant until the asymptotic regime $g\gg 1$ is reached. However, the latter is specific to processes like elastic scattering or photon emission and does not have to remain true for general higher-order QED processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1973

A new approach to the design of time varying control systems with application to the space shuttle booster

A general approach toward the analysis and design of closed loop control for slowly time-varying linear systems was investigated. The approach is by the method of generalized multiple scales, in which slow and fast dynamics are systematically separated by employing different clocks which measure time at varying rates. The clocks, which are necessarily nonlinear functions of time, so that the system dynamics are asymptotically invariant with respect to the new time scale. A transfer function relating the output to the input for general linear slowly time-varying systems is developed and represents the actual system under certain conditions. The clock functions are shown to satisfy an algebraic characteristic equation and can be determined in general in terms of the coefficients.

Ramnath, R. V.↗

NASA NTRS1982

The optimization of convergence for Chebyshev polynomial methods in an unbounded domain

Grosch and Orszag (1977) have performed a numerical analysis of the problem of solving differential equations in a semiinfinite or infinite domain using Chebyshev polynomials. The principal limitation of the conducted study was that it was entirely empirical. Various differential equations were solved in different ways and the numbers were compared. The present investigation has the objective to extend the studies conducted by Grosch and Orszag by deriving asymptotic approximations to the Chebyshev coefficients of simple model functions. This approach makes it possible to conduct more systematic comparisons of different methods, extend the range of comparisons, and, perhaps most important, give simple analytic formulas for choosing the optimum domain size or mapping parameter L for various situations.

Boyd, J. P.↗

NASA NTRS1992

The prediction of noise radiation from supersonic elliptic jets

This paper describes the prediction of noise radiation from supersonic elliptic jets. The noise is associated with the large scale structures in the jet mixing layer. These structures are described as instability waves. The local characteristics of the instability waves are determined from a compressible, linear, analysis. The jet mean velocity and density are described in elliptic cylindrical coordinates. The local eigensolution for the instability waves is determined from a finite difference solution of the non-separable boundary value problem. This inner solution which is formulated in terms of the method of multiple scales is matched with the radiated field using the method of matched asymptotic expansions. The form of the far-field directivity is derived. Predictions are presented for the noise radiation by the several modes of instability in the elliptic jet. The radiated field is not axisymmetric and certain modes radiate strongly in the directions of the major and minor axes of the jet. The extension of the present work to other geometries and flow fields is discussed.

Morris, Philip J.↗

NASA NTRS1992

Rapid time variations and radiative instabilities of astrophysical masers

The time-varying intensities are obtained for astrophysical masers that are radiatively unstable. Numerical integrations of the time-dependent, nonlinear equations of radiative transfer are performed with the usual approximation of a linear maser. At long times after changes in the physical conditions, the intensity of maser radiation reaches an asymptotic behavior and oscillates permanently in these idealized calculations with a period that is related to the length of the maser divided by the speed of light. The intensity varies by more than a factor of 10. These intensities depend upon the same four parameters as we originally found to determine the regime for radiative instabilities based on a stability analysis of the steady state. A detailed comparison is made between the predictions of the stability analysis and the time variations of the intensities. Calculations are performed for interacting pairs as well as isolated, individual masers.

Scappaticci, Gerardo A.↗

NASA NTRS1999

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗

NASA NTRS2006

Computational Analysis of the Effect of Porosity on Shock Cell Strength at Cruise

A computational flow field analysis is presented of the effect of core cowl porosity on shock cell strength for a modern separate flow nozzle at cruise conditions. The goal of this study was to identify the primary physical mechanisms by which the application of porosity can reduce shock cell strength and hence the broadband shock associated noise. The flow is simulated by solving the asymptotically steady, compressible, Reynoldsaveraged Navier-Stokes equations on a structured grid using an implicit, up-wind, flux-difference splitting finite volume scheme. The standard two-equation k - epsilon turbulence model with a linear stress representation is used with the addition of a eddy viscosity dependence on total temperature gradient normalized by local turbulence length scale. Specific issues addressed in this study were the optimal area required to weaken a shock impinging on the core cowl surface and the optimal level of porosity and placement of porous areas for reduction of the overall shock cell strength downstream. Two configurations of porosity were found to reduce downstream shock strength by approximately 50%.

Massey, Steven J.↗

NASA NTRS2008

Transient Effects in Planar Solidification of Dilute Binary Alloys

The initial transient during planar solidification of dilute binary alloys is studied in the framework of the boundary integral method that leads to the non-linear Volterra integral governing equation. An analytical solution of this equation is obtained for the case of a constant growth rate which constitutes the well-known Tiller's formula for the solute transient. The more physically relevant, constant ramping down temperature case has been studied both numerically and analytically. In particular, an asymptotic analytical solution is obtained for the initial transient behavior. A numerical technique to solve the non-linear Volterra equation is developed and the solution is obtained for a family of the governing parameters. For the rapid solidification condition, growth rate spikes have been observed even for the infinite kinetics model. When recirculating fluid flow is included into the analysis, the spike feature is dramatically diminished. Finally, we have investigated planar solidification with a fluctuating temperature field as a possible mechanism for frequently observed solute trapping bands.

Mazuruk, K.↗

NASA NTRS2008

Transient Effects in Planar Solidification of Dilute Binary Alloys

Mazuruk, Konstantin↗

NASA NTRS1986

NLTE-analysis of three extremely helium-rich O-type subdwarfs

Three extremely helium-rich sdO stars (LSE 153, LSE 259 and LSE 263) were analyzed spectroscopically by means of detailed NLTE model atmospheres. These stars are very hot, with effective temperatures ranging from 70,000 to 75,000 K and gravities between log g = 4.4 and 4.9. Upper limits for the hydrogen abundance were also derived. The evolutionary status of the sdO stars is discussed and it is concluded that they evolve from the asymptotic giant branch towards the white dwarf stage. A possible evolutionary link between these hot stars and the extreme helium stars of spectral type B is discussed.

Husfeld, D.↗

DOE OSTI2024

First simultaneous global QCD analysis of dihadron fragmentation functions and transversity parton distribution functions

We perform a comprehensive study within quantum chromodynamics (QCD) of dihadron observables in electron-positron annihilation, semi-inclusive deep-inelastic scattering, and proton-proton collisions, including recent cross section data from Belle and azimuthal asymmetries from STAR. We extract simultaneously for the first time π + π − dihadron fragmentation functions (DiFFs) and the nucleon transversity distributions for up and down quarks as well as antiquarks. For the transversity distributions we impose their small- x asymptotic behavior and the Soffer bound. In addition, we utilize a new definition of DiFFs that has a number density interpretation to then calculate expectation values for the dihadron invariant mass and momentum fraction. Furthermore, we investigate the compatibility of our transversity results with those from single-hadron fragmentation (from a transverse momentum dependent/collinear twist-3 framework) and the nucleon tensor charges computed in lattice QCD. We find a universal nature to all of this available information. Future measurements of dihadron production can significantly further this research, especially, as we show, those that are sensitive to the region of large parton momentum fractions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗