Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Appendices”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Modeling fuels and fire effects in 3D: Model description and applications

Scientists and managers critically need ways to assess how fuel treatments alter fire behavior, yet few tools currently exist for this purpose. In this work, we present a spatially-explicit-fuel-modeling system, FuelManager, which models fuels, vegetation growth, fire behavior (using a physics-based model, FIRETEC), and fire effects. FuelManager's flexible approach facilitates modeling fuels across a wide range of detail. Large trees or shrubs with specific coordinates are modeled as individual “Plants”, while understory plants are modeled as collections of plants called “LayerSets”. Both Plants and LayerSets contain various fuel particles (leaves, needles, twigs) with various properties including shape, size and surface area to volume ratio. A wide range of vegetation and treatments can be modeled, analyzed quantitatively and visualized in a 3D viewer. We describe the modeling approach and demonstrate fuel modeling at different levels of detail, fuel treatment and fire effects capabilities. Lastly, detailed model equations are provided in the Appendices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A modified Villain formulation of fractons and other exotic theories

We reformulate known exotic theories (including theories of fractons) on a Euclidean spacetime lattice. We write them using the Villain approach and then we modify them to a convenient range of parameters. The new lattice models are closer to the continuum limit than the original lattice versions. In particular, they exhibit many of the recently found properties of the continuum theories including emergent global symmetries and surprising dualities. Also, these new models provide a clear and rigorous formulation to the continuum models and their singularities. In appendices, we use this approach to review well-studied lattice models and their continuum limits. These include the XY-model, the $\mathbb{Z}_N$ clock-model, and various gauge theories in diverse dimensions. Here, this presentation clarifies the relation between the condensed-matter and the high-energy views of these systems. It emphasizes the role of symmetries associated with the topology of field space, duality, and various anomalies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A transient near to far field transformation method and verification benchmarking procedure

The numerical calculation of electromagnetic far fields in the time-domain requires a near to far field transformation (NTFF) method. While time-domain NTFF methods for popular finite-difference time-domain (FDTD) approaches are well established, there is little discourse on NTFF methods for finite-element time-domain (FETD) codes. Here, this work is concerned with the development of an NTFF method for the Empire FETD code, which utilizes curl and divergence conforming elements. This discretization presents a difficulty in obtaining the equivalent electric current for the NTFF. Straightforward finite element interpolation of the fields is shown to give poor accuracy. Alternative interpolation methods are recommended. An expanding magnetic quadrupole pulse benchmark problem, which is fully developed in the appendices, provides the basis for quantitative comparison.

FETD↗
DOE OSTI2016

Adaptive Sequential Monte Carlo for Multiple Changepoint Analysis

Process monitoring and control requires detection of structural changes in a data stream in real time. This paper introduces an efficient sequential Monte Carlo algorithm designed for learning unknown changepoints in continuous time. The method is intuitively simple: new changepoints for the latest window of data are proposed by conditioning only on data observed since the most recent estimated changepoint, as these observations carry most of the information about the current state of the process. The proposed method shows improved performance over the current state of the art. Another advantage of the proposed algorithm is that it can be made adaptive, varying the number of particles according to the apparent local complexity of the target changepoint probability distribution. This saves valuable computing time when changes in the changepoint distribution are negligible, and enables re-balancing of the importance weights of existing particles when a significant change in the target distribution is encountered. The plain and adaptive versions of the method are illustrated using the canonical continuous time changepoint problem of inferring the intensity of an inhomogeneous Poisson process, although the method is generally applicable to any changepoint problem. Performance is demonstrated using both conjugate and non-conjugate Bayesian models for the intensity. Lastly, appendices to the article are available online, illustrating the method on other models and applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

RADIATIVELY EFFICIENT MAGNETIZED BONDI ACCRETION

We have carried out a numerical study of the effect of large-scale magnetic fields on the rate of accretion from a uniform, isothermal gas onto a resistive, stationary point mass. Only mass, not magnetic flux, accretes onto the point mass. The simulations for this study avoid complications arising from boundary conditions by keeping the boundaries far from the accreting object. Our simulations leverage adaptive refinement methodology to attain high spatial fidelity close to the accreting object. Our results are particularly relevant to the problem of star formation from a magnetized molecular cloud in which thermal energy is radiated away on timescales much shorter than the dynamical timescale. Contrary to the adiabatic case, our simulations show convergence toward a finite accretion rate in the limit in which the radius of the accreting object vanishes, regardless of magnetic field strength. For very weak magnetic fields, the accretion rate first approaches the Bondi value and then drops by a factor of {approx}2 as magnetic flux builds up near the point mass. For strong magnetic fields, the steady-state accretion rate is reduced by a factor of {approx}0.2 {beta}{sup 1/2} compared to the Bondi value, where {beta} is the ratio of the gas pressure to the magnetic pressure. We give a simple expression for the accretion rate as a function of the magnetic field strength. Approximate analytic results are given in the Appendices for both time-dependent accretion in the limit of weak magnetic fields and steady-state accretion for the case of strong magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE UNIVERSAL INITIAL MASS FUNCTION IN THE EXTENDED ULTRAVIOLET DISK OF M83

We report deep Subaru H{alpha} observations of the extended ultraviolet (XUV) disk of M83. These new observations enable the first complete census of very young stellar clusters over the entire XUV disk. Combining Subaru and Galaxy Evolution Explorer data with a stellar population synthesis model, we find that (1) the standard, but stochastically sampled, initial mass function (IMF) is preferred over the truncated IMF because there are low-mass stellar clusters (10{sup 2-3} M{sub Sun }) that host massive O-type stars; (2) the standard Salpeter IMF and a simple aging effect explain the counts of far-UV (FUV)-bright and H{alpha}-bright clusters with masses >10{sup 3} M{sub Sun }; and (3) the H{alpha}-to-FUV flux ratio over the XUV disk supports the standard IMF. To reach conclusion (2), we assumed instantaneous cluster formation and a constant cluster formation rate over the XUV disk. The Subaru Prime Focus Camera covers a large area even outside the XUV disk-far beyond the detection limit of the H I gas. This enables us to statistically separate the stellar clusters in the disk from background contamination. The new data, model, and previous spectroscopic studies provide overall consistent results with respect to the internal dust extinction (A{sub V} {approx} 0.1 mag) and low metallicity ({approx}0.2 Z{sub Sun }) using the dust extinction curve of the Small Magellanic Cloud. The minimum cluster mass for avoiding the upper IMF incompleteness due to stochastic sampling and the spectral energy distributions of O, B, and A stars are discussed in the Appendices.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Vorticity, shocks, and magnetic fields in subsonic, ICM-like turbulence

We analyze data from high resolution simulations of the generation of compressible, MHD turbulence with properties chosen to resemble conditions in galaxy clusters. In particular, the flow is driven to have turbulence Mach number M{sub t}∼1/2 in an isothermal medium with an initially very weak, uniform seed magnetic field (β=P{sub g}/P{sub B}=10{sup 6}). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal (∇⋅δu=0), while in the other it is entirely compressive (∇×δu=0). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates, and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy–Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn–Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree–Fock and Bethe ansatz local density approximation) to show their successes and limitations. Here, we also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

CMB lensing and primordial squeezed non-gaussianity

Squeezed primordial non-Gaussianity can strongly constrain early-universe physics, but it can only be observed on the CMB after it has been gravitationally lensed. We give a new simple non-perturbative prescription for accurately calculating the effect of lensing on any squeezed primordial bispectrum shape, and test it with simulations. We give the generalization to polarization bispectra, and discuss the effect of lensing on the trispectrum. We explain why neglecting the lensing smoothing effect does not significantly bias estimators of local primordial non-Gaussianity, even though the change in shape can be ∼>10%. We also show how τ{sub NL} trispectrum estimators can be well approximated by much simpler CMB temperature modulation estimators, and hence that there is potentially a ∼ 10–30% bias due to very large-scale lensing modes, depending on the range of modulation scales included. Including dipole sky modulations can halve the τ{sub NL} error bar if kinematic effects can be subtracted using known properties of the CMB temperature dipole. Lensing effects on the g{sub NL} trispectrum are small compared to the error bar. In appendices we give the general result for lensing of any primordial bispectrum, and show how any full-sky squeezed bispectrum can be decomposed into orthogonal modes of distinct angular dependence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

On the importance of direct detection combined limits for spin independent and spin dependent dark matter interactions

In this work we show how the inclusion of dark matter (DM) direct detection upper bounds in a theoretically consistent manner can affect the allowed parameter space of a DM model. Traditionally, the limits from DM direct detection experiments on the elastic scattering cross section of DM particles as a function of their mass are extracted under simplifying assumptions. Relaxing the assumptions related to the DM particle nature, such as the neutron to proton ratio of the interactions, or the possibility of having similar contributions from the spin independent (SI) and spin dependent (SD) interactions can vary significantly the upper limits. Furthermore, it is known that astrophysical and nuclear uncertainties can also affect the upper bounds. To exemplify the impact of properly including all these factors, we have analysed two well motivated and popular DM scenarios: neutralinos in the NMSSM and a Z' portal with Dirac DM. We have found that the allowed parameter space of these models is subject to important variations when one includes both the SI and SD interactions at the same time, realistic neutron to proton ratios, as well as using different self-consistent speed distributions corresponding to popular DM halo density profiles, and distinct SD structure functions. Finally, we provide all the necessary information to include the upper bounds of SuperCDMS and LUX taking into account all these subtleties in the investigation of any particle physics model. The data for each experiment and example codes are available at this site http://goo.gl/1CDFYi, and their use is detailed in the appendices of this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Dark matter in the Sun: scattering off electrons vs nucleons

The annihilation of dark matter (DM) particles accumulated in the Sun could produce a flux of neutrinos, which is potentially detectable with neutrino detectors/telescopes and the DM elastic scattering cross section can be constrained. Although the process of DM capture in astrophysical objects like the Sun is commonly assumed to be due to interactions only with nucleons, there are scenarios in which tree-level DM couplings to quarks are absent, and even if loop-induced interactions with nucleons are allowed, scatterings off electrons could be the dominant capture mechanism. We consider this possibility and study in detail all the ingredients necessary to compute the neutrino production rates from DM annihilations in the Sun (capture, annihilation and evaporation rates) for velocity-independent and isotropic, velocity-dependent and isotropic and momentum-dependent scattering cross sections for DM interactions with electrons and compare them with the results obtained for the case of interactions with nucleons. Moreover, we improve the usual calculations in a number of ways and provide analytical expressions in three appendices. Interestingly, we find that the evaporation mass in the case of interactions with electrons could be below the GeV range, depending on the high-velocity tail of the DM distribution in the Sun, which would open a new mass window for searching for this type of scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Interacting dark sector and precision cosmology

We consider a recently proposed model in which dark matter interacts with a thermal background of dark radiation. Dark radiation consists of relativistic degrees of freedom which allow larger values of the expansion rate of the universe today to be consistent with CMB data (H 0 -problem). Scattering between dark matter and radiation suppresses the matter power spectrum at small scales and can explain the apparent discrepancies between ΛCDM predictions of the matter power spectrum and direct measurements of Large Scale Structure LSS (σ 8 -problem). We go beyond previous work in two ways: 1. we enlarge the parameter space of our previous model and allow for an arbitrary fraction of the dark matter to be interacting and 2. we update the data sets used in our fits, most importantly we include LSS data with full k-dependence to explore the sensitivity of current data to the shape of the matter power spectrum. We find that LSS data prefer models with overall suppressed matter clustering due to dark matter - dark radiation interactions over ΛCDM at 3-4 σ. However recent weak lensing measurements of the power spectrum are not yet precise enough to clearly distinguish two limits of the model with different predicted shapes for the linear matter power spectrum. In two Appendices we give a derivation of the coupled dark matter and dark radiation perturbation equations from the Boltzmann equation in order to clarify a confusion in the recent literature, and we derive analytic approximations to the solutions of the perturbation equations in the two physically interesting limits of all dark matter weakly interacting or a small fraction of dark matter strongly interacting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The reconstructed power spectrum in the Zeldovich approximation

Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. Here, we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the ${\tt DarkSky}$ simulation, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical findings for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Transport-driven toroidal rotation with general viscosity profile

Abstract Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up↗
DOE OSTI2021

Mapping the magnetic field in the Taurus/B211 filamentary cloud with SOFIA HAWC + and comparing with simulation

ABSTRACT Optical and infrared polarization mapping and recent Planck observations of the filametary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to the long axis of the cloud. We use the HAWC + polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of 36°, about five times that measured on a larger scale by Planck. Applying the Davis–Chandrasekhar–Fermi (DCF) method with velocity information obtained from Institut de Radioastronomie Millimétrique 30 m C18O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfv$\acute{\rm e}$nic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfv$\acute{\rm e}$nic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Millimeter observational signatures of flares in magnetically arrested black hole accretion models

In general relativistic magnetohydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes in very high resolution GRMHD simulations. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating γ ≳ 100 (T e ≳ 10 11 K) electrons. For models in which the electrons are heated to T e ∼ 10 11 K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multiwavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulating image coaddition with the Nancy Grace Roman Space Telescope – I. Simulation methodology and general results

The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near-infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired ‘target’ point spread function (PSF). We present here an initial application of imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the ‘leakage’), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper (‘Paper II’) explores the implications for weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗