Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian simulation”

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 127 records · Page 7

Electron cyclotron emission quasi-optical transmission system on the HL-3 tokamak

A new quasi-optical (QO) Electron Cyclotron Emission (ECE) transmission system has been established on the HL-3 tokamak, which includes a focusing QO mirror combination and a long-distance transmission line. This system was developed to meet the requirements for poloidal spatial resolution and the high signal-to-noise ratio needed for magnetohydrodynamic (MHD) instability studies using ECE on the HL-3. The QO mirror combination was installed inside the vacuum chamber for focusing. Laboratory test results, theoretical calculations, and synthetic ECE simulation results indicate that the Gaussian beam can meet the spatial resolution requirements for the accurate measurement of the MHD instability on the q = 1/2/3 surfaces, corresponding to the poloidal mode numbers m = 3/6/9. This includes good diagnostic poloidal spatial resolution for the important 2/1 and 3/2 modes. At the front end of the transmission line, a high-efficiency mode converter was designed to transition the TE 10 mode to the HE 11 mode for input into the transmission line, with an insertion loss of less than 1.5 dB. A 30 m long-distance corrugated oversized waveguide was constructed, with transmission losses ranging from 6 to 10 dB in the 60–120 GHz range. Polarization adjustment results show that the polarization offset and geometric spatial polarization angle change consistently, which can provide a reference for polarization adjustment in other complex structured transmission lines. As a result, the newly established ECE QO transmission system will provide strong support for future physics research involving ECE on the HL-3.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

Impact of hydrogenation on the stability and mechanical properties of amorphous boron nitride

Abstract Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g cm −3 ), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of aBN upon hydrogen doping. With the introduction of a Gaussian approximation potential for atomistic simulations, we investigate the changing morphology of aBN with varying H doping concentrations. We found that for 8 at% of H doping, the concentration ofsp 3 -hybridized atoms reaches to a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of aBN films for numerous applications.

Materials Science↗

Cosmological constraints from the Subaru Hyper Suprime-Cam year 1 shear catalogue lensing convergence probability distribution function

Here we utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) year 1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude S 8 = σ 8 $\sqrt {Ω_{m}/0.3}$. The large-scale structure’s statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF—a complementary higher-order statistic. By defining the PDF over the standard-deviation-normalized convergence map, we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find S 8 = $0.860^{+0.066}_{–0.109}$ from the PDF alone and S 8 = 0.798$^{+0.029}_{–0.042}$ from the combination of the PDF and power spectrum (68% confidential level (CL)). The PDF improves the power-spectrum-only constraint by about 10%.

79 ASTRONOMY AND ASTROPHYSICS↗

Multifrequency Models of Black Hole Photon Rings from Low-luminosity Accretion Disks

Images of black holes encode both astrophysical and gravitational properties. Detecting highly lensed features in images can differentiate between these two effects. We present an accretion disk emission model coupled to the Adaptive Analytical Ray Tracing (AART) code that allows a fast parameter space exploration of black hole photon ring images produced from synchrotron emission from 10 to 670 GHz. As an application, we systematically study several disk models and compute their total flux density, average radii, and optical depth. The model parameters are chosen around fiducial values calibrated to general relativistic magnetohydrodynamic (GRMHD) simulations and observations of M87*. For the parameter space studied, we characterize the transition between optically thin and thick regimes and the frequency at which the first photon ring is observable. Our results highlight the need for careful definitions of photon ring radius in the image domain, as in certain models the highly lensed photon ring is dimmer than the direct emission at certain angles. We find that at low frequencies the ring radii are set by the electron temperature, while at higher frequencies the magnetic field strength plays a more significant role, demonstrating how multifrequency analysis can also be used to infer plasma parameters. Lastly, we show how our implementation can qualitatively reproduce multifrequency black hole images from GRMHD simulations when adding time variability to our disk model through Gaussian random fields. This approach provides a new method for simulating observations from the Event Horizon Telescope and the proposed Black Hole Explorer space mission.

79 ASTRONOMY AND ASTROPHYSICS↗

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey. II. Anisotropic Large-scale Coherence in Hot Gas, Galaxies, and Dark Matter

Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev–Zel’dovich (tSZ) maps to measure the anisotropic distribution of hot gas 2.5–40 Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of ∼15 Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a ∼3× more significant measurement of the extended gas compared to the technique’s proof-of-concept. Decomposing stacks into cosine multipoles of order m , we detect a dipole ( m = 1) and quadrupole ( m = 2) at 8σ–10σ, as well as evidence for m = 4 signal at up to 6σ, indicating sensitivity to late-time non-Gaussianity. We compare to Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of photon correlation techniques for processing of laser velocimeter signals

The objective was to provide the theory and a system design for a new type of photon counting processor for low level dual scatter laser velocimeter (LV) signals which would be capable of both the first order measurements of mean flow and turbulence intensity and also the second order time statistics: cross correlation auto correlation, and related spectra. A general Poisson process model for low level LV signals and noise which is valid from the photon-resolved regime all the way to the limiting case of nonstationary Gaussian noise was used. Computer simulation algorithms and higher order statistical moment analysis of Poisson processes were derived and applied to the analysis of photon correlation techniques. A system design using a unique dual correlate and subtract frequency discriminator technique is postulated and analyzed. Expectation analysis indicates that the objective measurements are feasible.

Mayo, W. T., Jr.↗

A simulation of air pollution model parameter estimation using data from a ground-based LIDAR remote sensor

One way to obtain estimates of the unknown parameters in a pollution dispersion model is to compare the model predictions with remotely sensed air quality data. A ground-based LIDAR sensor provides relative pollution concentration measurements as a function of space and time. The measured sensor data are compared with the dispersion model output through a numerical estimation procedure to yield parameter estimates which best fit the data. This overall process is tested in a computer simulation to study the effects of various measurement strategies. Such a simulation is useful prior to a field measurement exercise to maximize the information content in the collected data. Parametric studies of simulated data matched to a Gaussian plume dispersion model indicate the trade offs available between estimation accuracy and data acquisition strategy.

Kibler, J. F.↗

Modeling laser velocimeter signals as triply stochastic Poisson processes

Previous models of laser Doppler velocimeter (LDV) systems have not adequately described dual-scatter signals in a manner useful for analysis and simulation of low-level photon-limited signals. At low photon rates, an LDV signal at the output of a photomultiplier tube is a compound nonhomogeneous filtered Poisson process, whose intensity function is another (slower) Poisson process with the nonstationary rate and frequency parameters controlled by a random flow (slowest) process. In the present paper, generalized Poisson shot noise models are developed for low-level LDV signals. Theoretical results useful in detection error analysis and simulation are presented, along with measurements of burst amplitude statistics. Computer generated simulations illustrate the difference between Gaussian and Poisson models of low-level signals.

Mayo, W. T., Jr.↗

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗

Attitude determination and parameter estimation using vector observations - Application

This paper presents tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and earth-pointing attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. In the majority of tests performed, the two methods converge to the same estimates in the same number of iterations, but the new algorithm requires about 60 percent more computational effort. Some cases were found in which the new method converges in fewer iterations than the differential correction, and some for which the differential correction requires fewer iterations.

Markley, F. Landis↗

Measured and computed performance of a microstrip filter composed of semi-insulating GaAs on a fused quartz substrate

The performance of a microstrip hammerhead filter that has been fabricated on an electrically thin layer of semiinsulating GaAs backed by a fused quartz substrate was measured and compared to results of a three-dimensional finite-difference time-domain (FD-TD) program used to calculate the response of the filter both with and without the GaAs layer. The program, presented by Sheen et al. (1990), discretizes the entire structure and then simulates the propagation of a Gaussian pulse through the filter. The microstrip filter is intended for applications involving ultrathin lifted-off or etched-back GaAs containing both active devices and passive microstrip circuitry backed by a much thicker mechanically rigid low-loss, low-dielectric-constant substrate. The low-pass characteristics of the hammerhead filter with the intermediate GaAs layer are compared with those of the same filter on quartz alone. Both the measured and computed data show a significant shift in cutoff frequency (about 10 percent at the 3 dB points) for a GaAs layer that is 0.007 wavelengths thick at 4 GHz.

Siegel, Peter H.↗

Bit-wise arithmetic coding for data compression

This article examines the problem of compressing a uniformly quantized independent and identically distributed (IID) source. We present a new compression technique, bit-wise arithmetic coding, that assigns fixed-length codewords to the quantizer output and uses arithmetic coding to compress the codewords, treating the codeword bits as independent. We examine the performance of this method and evaluate the overhead required when used block-adaptively. Simulation results are presented for Gaussian and Laplacian sources. This new technique could be used as the entropy coder in a transform or subband coding system.

Kiely, A. B.↗

Turbulence Model Validation Through Joint Experimental /Computational Studies of Separated Flow Over A Three-Dimensional Tapered Bump: Part I - Experimental Investigation

The results of a three-year experimental campaign aimed at comprehensively documenting the separated flow over a three-dimensional bump are presented with the purpose of generating a benchmark experimental database useful in validating computational fluid dynamics (CFD) flow simulations. The bump model geometry was designed to provide well-defined and repeatable smooth-body flow separation conditions that were suitable for both experiments and simulations. The bump had a Gaussian streamwise profile with a constant height equal to 8.5% of its width over the central 60% of the test section width. The remaining 40% were outboard spanwise portions that gradually taper to zero using an error function profile to minimize side-wall boundary layer effects. The model was immersed in a turbulent boundary layer that was developed on a suspended flat plate in the Notre Dame Mach 0.6 Wind Tunnel. In order to document the effect of the incoming boundary layer thickness on the flow separation, the bump model could be located at two streamwise positions. The mean velocity and turbulence intensity of the wind tunnel freestream flow field and approaching turbulent boundary layer were fully documented. The measurements of the flow separation region included surface visualization, wall shear stress using oil-film interferometry, mean and dynamic surface pressure, and planar and stereoscopic particle image velocimetry. The experiments were conducted over a range of Mach numbers from 0.05 to 0.2 corresponding to a range of Reynolds numbers based on the test section spanwise dimension (0.914m) of 1.0 × 10 6 ≤ Re L = U ∞ L/ν ≤ 4.0 × 10 6 . The bulk of the results are presented for the higher Mach number conditions of 0.1 and 0.2 with Re L = 2.0 × 10 6 and 4.0 × 10 6 , respectively. Extensive uncertainty analysis of the data was performed. The data is archived in the NASA Langley Turbulence Modeling Resource website at https://turbmodels.larc.nasa.gov/Other_exp_Data/speedbump_sep_exp.html. In addition to the experiments, a computational effort was made in parallel by the CFD group at Boeing Research & Technology highlighting the usefulness of the data set, which is outlined in the accompanying CFD report.

Flow separation↗

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. II. cosmological results

Here, we present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years of the Dark Energy Survey. We implement the following: (1) second and third moments; (2) wavelet phase harmonics; (3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven 𝑤 Cold Dark Matter (𝑤⁢ CDM) cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the figure of merit (𝑆 8 , Ω m ) is improved by 70% (Λ ⁢CDM) and 90% (𝑤 ⁢CDM). When all the summary statistics are combined, we achieve a 2% constraint on the amplitude of fluctuations parameter 𝑆 8 ≡ 𝜎 8 ⁢(Ω m /0.3) 0.5 , obtaining 𝑆 8 = 0.794 ±0.017 (Λ⁢ CDM) and 𝑆 8 = 0.817 ±0.021 (𝑤 ⁢CDM), and a ∼10% constraint on Ω m , obtaining Ω m =0.259 ±0.025 (Λ ⁢CDM) and Ω m = 0.273 ±0.029 (𝑤⁢ CDM). In the context of the 𝑤⁢ CDM scenario, these statistics also strengthen the constraints on the parameter 𝑤, obtaining 𝑤 <−0.72. The constraints from different statistics are shown to be internally consistent (with a 𝑝-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the first three years of the Dark Energy Survey data, finding good consistency; we also compare with results from external datasets, such as planck constraints from the cosmic microwave background, finding statistical agreement, with discrepancies no greater than <2.2⁢𝜎.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear estimation theory applied to the interplanetary orbit determination problem.

Martingale theory and appropriate smoothing properties of Loeve (1953) have been used to develop a modified Gaussian second-order filter. The performance of the filter is evaluated through numerical simulation of a Jupiter flyby mission. The observations used in the simulation are on-board measurements of the angle between Jupiter and a fixed star taken at discrete time intervals. In the numerical study, the influence of each of the second-order terms is evaluated. Five filter algorithms are used in the simulations. Four of the filters are the modified Gaussian second-order filter and three approximations derived by neglecting one or more of the second-order terms in the equations. The fifth filter is the extended Kalman-Bucy filter which is obtained by neglecting all of the second-order terms.

Tapley, B. D.↗

Microstructure-Based Modeling of Laser Beam Shaping During Additive Manufacturing

Recent experimental studies suggest the use of spatially extended laser beam profiles as a strategy to control the melt pool during laser powder bed fusion (LPBF) additive manufacturing. However, linkages connecting laser beam profiles to thermal fields and resultant microstructures have not been established. Herein, we employ a coupled thermal transport-Monte Carlo model to predict the evolution of temperature fields and grain microstructures during LPBF using Gaussian, ring, and Bessel beam profiles. In addition, simulation results reveal that the ring-shaped beam yields lower temperatures compared with the Gaussian beam. Owing to the small melt pool size when using the Bessel beam, the grains are smaller in size and more equiaxed compared to those using the Gaussian and ring beams. Our approach provides future avenues to predict the impact of laser beam shaping on microstructure development during LPBF.

36 MATERIALS SCIENCE↗

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗