Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Asymptotic preserving methods”

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.

28 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗
NASA NTRS1985

An experimental study of plane mixing layer development

Mean flow and turbulence was measured in the near field of two plane mixing layers in air with a maximum velocity of 21 m/sec. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Gases with both tripped and untripped initial boundary layers were studied. It was found that the two stream layer developed to the self preserving state in a distance much shorter than did the single stream layer which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two stream layer and the development of the single stream layer showed agreement with existing data. It is suggested that the two stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single stream mixing layer.

Mehta, R. D.↗
DOE OSTI2023

A finite element method for angular discretization of the radiation transport equation on spherical geodesic grids

Discrete ordinate (S N ) and filtered spherical harmonics (FP N ) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. Here, the resulting method is then compared with both S N and FP N schemes using four test problems and is found to perform well when one of the other methods fail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Private Tabular Survey Data Products through Synthetic Microdata Generation

We propose two synthetic microdata approaches to generate private tabular survey data products for public release. We adapt a pseudo posterior mechanism that downweights by-record likelihood contributions with weights ∈[0,1] based on their identification disclosure risks to producing tabular products for survey data. Our method applied to an observed survey database achieves an asymptotic global probabilistic differential privacy guarantee. Our two approaches synthesize the observed sample distribution of the outcome and survey weights, jointly, such that both quantities together possess a privacy guarantee. The privacy-protected outcome and survey weights are used to construct tabular cell estimates (where the cell inclusion indicators are treated as known and public) and associated standard errors to correct for survey sampling bias. Through a real data application to the Survey of Doctorate Recipients public use file and simulation studies motivated by the application, we demonstrate that our two microdata synthesis approaches to construct tabular products provide superior utility preservation as compared to the additive noise approach of the Laplace Mechanism. Moreover, our approaches allow the release of microdata to the public, enabling additional analyses at no extra privacy cost.

Mathematical Methods In Social Sciences↗
NASA NTRS1984

Near-field turbulence properties of single- and two-stream plane mixing layers

Detailed mean-flow and turbulence measurements have been made in the near-field of two plane mixing layers in air with a maximum velocity of 21 m/s. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Cases with both tripped and untripped initial boundary layers were studied. In all cases, it was found that the two-stream layer developed to the self-preserving state in a distance much shorter than the single-stream layer which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two-stream layer, and the development of the single-stream layer showed agreement with existing data. The results suggest that the two-stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single-stream mixing layer.

Mehta, R. D.↗
NASA NTRS2024

Tracing Presolar Grain Populations Through Matrix Alteration Histories

Primitive meteorites preserve remnants of evolved stars as presolar grains, from asymptotic and red giant branch stars, novae and supernovae. These grains retain records of the formation mechanisms of the stellar environment from which they condensed, as well as physical and chemical processes in the interstellar medium (ISM), solar nebula and the parent bodies in which they are found. Consequently, their study can refine our understanding of the solar system and external stellar and planetary systems. However, presolar grains and their records can be altered, re-equilibrated with the surrounding matrix, and eventually lost through aqueous and thermal processing on the meteorite parent body. NanoSIMS is the most rapid and effective method for in situ identification of presolar grains via their highly anomalous isotopic ratios. Previous studies have measured variations in presolar grain populations between different matrix regions within the same meteorite thin section. However, most of these studies do not identify if there are any ties between these variations in presolar grain populations and the chemical and mineralogical differences in measured matrix regions. Therefore, the goal of this study is to determine if there is measurable variations in the degree of alteration of different matrix regions within the same thin section, and if mineralogical and chemical changes in matrix regions can be used as indicators for presolar grain preservation. Such a correlation could improve the efficiency of presolar grain identification during searches across large surface areas. To do this, we carried out coordinated scanning electron microscope (SEM) and NanoSIMS analyses, using a new SEM method which measures the ppm concentration of 100 nm to 5 μm sized Mg-, S-, Al-, C- and Ca-rich matrix grains.

presolar grains↗
DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗