Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scattering matrix simulations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurements of Zγ +jets differential cross sections in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

Differential cross-section measurements of $Zγ$ production in association with hadronic jets are presented, using the full 139 fb -1 dataset of $\sqrt{s}$ = 13 TeV proton–proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the $Z$ boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLO PS , as well as next-to-leading-order predictions from MADGRAPH5_AMC@NLO and SHERPA.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

DESI mock challenge: constructing DESI galaxy catalogues based on FastPM simulations

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large-scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a halo occupation distribution (HOD) model upon the FastPM simulations, such that the resulting galaxy clustering reproduces high-resolution N-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements – monopole and quadrupole – to a precision required by the DESI Year 1 emission line galaxy sample down to non-linear scales, i.e. $k\lt 0.5\, h\, \mathrm{Mpc}^{-1}$ or $s\gt 10\, \mathrm{Mpc}\, h^{-1}$. Furthermore, we compute covariance matrices based on the resulting FastPM galaxy clustering – monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate (e.g. power spectrum) on the covariance matrix of the Fourier counterpart (e.g. correlation function). We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between 2 and 10 per cent. Consequently, using the current pipeline we can precisely reproduce the clustering of N-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

Probabilistic Simulation of Progressive Fracture in Bolted-Joint Composite Laminates

This report describes computational methods to probabilistically simulate fracture in bolted composite structures. An innovative approach that is independent of stress intensity factors and fracture toughness was used to simulate progressive fracture. The effect of design variable uncertainties on structural damage was also quantified. A fast probability integrator assessed the scatter in the composite structure response before and after damage. Then the sensitivity of the response to design variables was computed. General-purpose methods, which are applicable to bolted joints in all types of structures and in all fracture processes-from damage initiation to unstable propagation and global structure collapse-were used. These methods were demonstrated for a bolted joint of a polymer matrix composite panel under edge loads. The effects of the fabrication process were included in the simulation of damage in the bolted panel. Results showed that the most effective way to reduce end displacement at fracture is to control both the load and the ply thickness. The cumulative probability for longitudinal stress in all plies was most sensitive to the load; in the 0 deg. plies it was very sensitive to ply thickness. The cumulative probability for transverse stress was most sensitive to the matrix coefficient of thermal expansion. In addition, fiber volume ratio and fiber transverse modulus both contributed significantly to the cumulative probability for the transverse stresses in all the plies.

Minnetyan, L.↗
DOE OSTI2026

Efficient many-jet event generation with flow matching

We apply for the first time, to the best of our knowledge, the flow matching method to the problem of phase-space sampling for event generation in high-energy collider physics. By training the model to remap the random numbers used to generate the momenta and helicities of the scattering matrix elements as implemented in the portable partonic event generator pepper, we find substantial efficiency improvements in the studied processes. We focus our study on the highest final-state multiplicities in Drell-Yan and top-antitop pair production used in simulated samples for the Large Hadron Collider, which computationally are the most relevant ones. We find that the unweighting efficiencies improve by factors of 184 and 25, respectively, when compared to the standard approach of using a vegas-based optimization. We also compare continuous normalizing flows trained with flow matching against the previously studied normalizing flows based on coupling layers and find that the former leads to better results, faster training and a better scaling behavior across the studied multiplicity range, while the latter evaluate faster. When combining the advantages of both methods using the regflow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest final-state multiplicities.

Bothmann, E. [CERN; Gottingen U.] (ORCID:000000016↗
NASA NTRS2009

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗
DOE OSTI2024

Spin squeezing with itinerant dipoles: A case for shallow lattices

Entangled spin squeezed states generated via dipolar interactions in lattice models provide unique opportunities for quantum enhanced sensing and are now within reach of current experiments. A critical question in this context is which parameter regimes offer the best prospects under realistic conditions. Light scattering in deep lattices can induce significant decoherence and strong Stark shifts, while shallow lattices face motional decoherence as a fundamental obstacle. Here we analyze the interplay between motion and spin squeezing in itinerant fermionic dipoles in one dimensional chains using exact matrix product state simulations. We demonstrate that shallow lattices can achieve more than 5 dB of squeezing, outperforming deep lattices by up to more than 3 dB, even in the presence of low filling, loss, and decoherence. We relate this finding to SU(2)-symmetric superexchange interactions, which keep spins aligned and protect collective correlations. We show that the optimal regime is achieved for small repulsive off-site interactions, with a trade-off between maximal squeezing and optimal squeezing time. Published by the American Physical Society 2024

Wellnitz, David (ORCID:0000000349788938)↗
NASA NTRS2005

Role of Probabilistic Micromechanics Modeling in Establishing Design Allowables in Composites

One of the major challenges in designing with any new material, and particularly with advanced composite materials, is the fidelity of material design allowables. In the case of composite materials, the concern arises from the inherent nature of these materials, i.e., their heterogeneous make-up and the various factors that affect their properties in a specific design environment. Composites have various scales - micro, macro, laminate and structural, as well as numerous other fabrication related parameters. Many advanced composites in aerospace applications involve complex two- and three-dimensional fiber architectures and requires high-temperature processing. Since there are uncertainties associated with each of these, the observed behavior of composite materials shows scatter. Evaluating the effect of each of these variables on the observed scatter in composite properties solely by teSting is cost and time prohibitive. One alternative is to evaluate these effects by computational simulation. The authors have developed probabilistic composite micromechanics techniques by combining woven composite micromechanics and Fast Probability Integration (FPI) techniques to address these issues. In this paper these techniques will be described and demonstrated through selected examples. Results in the form of cumulative distribution functions (CDF) of the composite properties of a MI (melt-infiltrated) SiC/SiC (silicon carbide fiber in a silicon carbide matrix) Composite will be presented. A CDF is a relationship defined by the value of the property (the response variable) with respect to the cumulative probability of occurrence. Furthermore, input variables causing scatter are identified and ranked based upon their sensitivity magnitude. Sensitivity information is very valuable in quality control. How these results can be utilized to develop design allowables so that these materials may be used by structural analysts/designers will also be discussed.

Mital, Subodh K.↗
DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗
DOE OSTI2023

Resolution of Virtual Depth Sectioning from Four-Dimensional Scanning Transmission Electron Microscopy

Abstract One approach to three-dimensional structure determination using the wealth of scattering data in four-dimensional (4D) scanning transmission electron microscopy (STEM) is the parallax method proposed by Ophus et al. (2019. Advanced phase reconstruction methods enabled by 4D scanning transmission electron microscopy, Microsc Microanal25, 10–11), which determines the scattering matrix and uses it to synthesize a virtual depth-sectioning reconstruction of the sample structure. Drawing on an equivalence with a hypothetical confocal imaging mode, we derive contrast transfer and point spread functions for this parallax method applied to weakly scattering objects, showing them identical to earlier depth-sectioning STEM modes when only bright field signal is used, but that improved depth resolution is possible if dark field signal can be used. Through a simulation-based study of doped Si, we show that this depth resolution is preserved for thicker samples, explore the impact of shot noise on the parallax reconstructions, discuss challenges to making use of dark field signal, and identify cases where the interpretation of the parallax reconstruction breaks down.

3D imaging↗
NASA NTRSDate not supplied

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗
NASA NTRS1992

Radiative transfer theory for polarimetric remote sensing of pine forest

The radiative transfer theory is applied to interpret polarimetric radar backscatter from pine forest with clustered vegetation structures. To take into account the clustered structures with the radiative transfer theory, the scattering function of each cluster is calculated by incorporating the phase interference of scattered fields from each component. Subsequently, the resulting phase matrix is used in the radiative transfer equations to evaluate the polarimetric backscattering coefficients from random medium layers embedded with vegetation clusters. Upon including the multi-scale structures, namely, trunks, primary and secondary branches, as well as needles, we interpret and simulate the polarimetric radar responses from pine forest for different frequencies and looking angles. The preliminary results are shown to be in good agreement with the measured backscattering coefficients at the Landes maritime pine forest during the MAESTRO-1 experiment.

Hsu, C. C.↗
NASA NTRS1991

Applications of the conjugate gradient FFT method in scattering and radiation including simulations with impedance boundary conditions

The theoretical and computational aspects related to the application of the Conjugate Gradient FFT (CGFFT) method in computational electromagnetics are examined. The advantages of applying the CGFFT method to a class of large scale scattering and radiation problems are outlined. The main advantages of the method stem from its iterative nature which eliminates a need to form the system matrix (thus reducing the computer memory allocation requirements) and guarantees convergence to the true solution in a finite number of steps. Results are presented for various radiators and scatterers including thin cylindrical dipole antennas, thin conductive and resistive strips and plates, as well as dielectric cylinders. Solutions of integral equations derived on the basis of generalized impedance boundary conditions (GIBC) are also examined. The boundary conditions can be used to replace the profile of a material coating by an impedance sheet or insert, thus, eliminating the need to introduce unknown polarization currents within the volume of the layer. A general full wave analysis of 2-D and 3-D rectangular grooves and cavities is presented which will also serve as a reference for future work.

Barkeshli, Kasra↗
DOE OSTI2024

3D Heisenberg universality in the van der Waals antiferromagnet NiPS 3

Van der Waals (vdW) magnetic materials are comprised of layers of atomically thin sheets, making them ideal platforms for studying magnetism at the two-dimensional (2D) limit. These materials are at the center of a host of novel types of experiments, however, there are notably few pathways to directly probe their magnetic structure. We confirm the magnetic order within a single crystal of NiPS 3 and show it can be accessed with resonant elastic X-ray diffraction along the edge of the vdW planes in a carefully grown crystal by detecting structurally forbidden resonant magnetic X-ray scattering. We find the magnetic order parameter has a critical exponent of β ~ 0.36, indicating that the magnetism of these vdW crystals is more adequately characterized by the three-dimensional (3D) Heisenberg universality class. We verify these findings with first-principles density functional theory, Monte-Carlo simulations, and density matrix renormalization group calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2012

Efficient Implementation of the Invariant Imbedding T-Matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.

Invariant imbedding method↗
NASA NTRS1990

Electromagnetic scattering from two-dimensional thick material junctions

The problem of the plane wave diffraction is examined by an arbitrary symmetric two dimensional junction, where Generalized Impedance Boundary Conditions (GIBCs) and Generalized Sheet Transition Conditions (GSTCs) are employed to simulate the slabs. GIBCs and GSTCs are constructed for multilayer planar slabs of arbitrary thickness and the resulting GIBC/GSTC reflection coefficients are compared with exact counterparts to evaluate the GIBCs/GSTCs. The plane wave diffraction by a multilayer material slab recessed in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering Matrix Formulation (GDMF) in conjunction with the dual integral equation approach. Various scattering patterns are computed and validated with exact results where possible. The diffraction by a material discontinuity in a thick dielectric/ferrite slab is considered by modelling the constituent slabs with GSTCs. A non-unique solution in terms of unknown constants is obtained, and these constants are evaluated for the recessed slab geometry by comparison with the solution obtained therein. Several other simplified cases are also presented and discussed. An eigenfunction expansion method is introduced to determine the unknown solution constants in the general case. This procedure is applied to the non-unique solution in terms of unknown constants; and scattering patterns are presented for various slab junctions and compared with alternative results where possible.

Ricoy, M. A.↗
NASA NTRS2011

Effect of Thin Cirrus Clouds on Dust Optical Depth Retrievals From MODIS Observations

The effect of thin cirrus clouds in retrieving the dust optical depth from MODIS observations is investigated by using a simplified aerosol retrieval algorithm based on the principles of the Deep Blue aerosol property retrieval method. Specifically, the errors of the retrieved dust optical depth due to thin cirrus contamination are quantified through the comparison of two retrievals by assuming dust-only atmospheres and the counterparts with overlapping mineral dust and thin cirrus clouds. To account for the effect of the polarization state of radiation field on radiance simulation, a vector radiative transfer model is used to generate the lookup tables. In the forward radiative transfer simulations involved in generating the lookup tables, the Rayleigh scattering by atmospheric gaseous molecules and the reflection of the surface assumed to be Lambertian are fully taken into account. Additionally, the spheroid model is utilized to account for the nonsphericity of dust particles In computing their optical properties. For simplicity, the single-scattering albedo, scattering phase matrix, and optical depth are specified a priori for thin cirrus clouds assumed to consist of droxtal ice crystals. The present results indicate that the errors in the retrieved dust optical depths due to the contamination of thin cirrus clouds depend on the scattering angle, underlying surface reflectance, and dust optical depth. Under heavy dusty conditions, the absolute errors are comparable to the predescribed optical depths of thin cirrus clouds.

Feng, Qian↗
DOE OSTI2022

Quantitative assessment of fitting errors associated with streak camera noise in Thomson scattering data analysis

Thomson scattering measurements in high energy density experiments are often recorded using optical streak cameras. In the low-signal regime, noise introduced by the streak camera can become an important and sometimes the dominant source of measurement uncertainty. In this paper, we present a formal method of accounting for the presence of streak camera noise in our measurements. We present a phenomenological description of the noise generation mechanisms and present a statistical model that may be used to construct the covariance matrix associated with a given measurement. This model is benchmarked against simulations of streak camera images. We demonstrate how this covariance may then be used to weight fitting of the data and provide quantitative assessments of the uncertainty in the fitting parameters determined by the best fit to the data and build confidence in the ability to make statistically significant measurements in the low-signal regime, where spatial correlations in the noise become apparent. These methods will have general applicability to other measurements made using optical streak cameras.

47 OTHER INSTRUMENTATION↗