Search NASA⌕ Search

SEARCH · Search NASA

Results for “SCATTERING MATRIX”

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 451 records · Page 25

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↗

A New Galerkin Quadrature Method Not Requiring a Matrix Inverse

We derive a new Galerkin quadrature (GQ) method for S 𝑛 calculations that differs from the two methods preceding it in that a matrix inverse for an 𝑁 𝑑 × 𝑁 𝑑 matrix, where 𝑁𝑑 is the number of directions in the quadrature set, is no longer required. Galerkin quadrature methods are designed for calculations with highly anisotropic scattering. Such methods are not simply special angular quadratures but also are methods for representing the S 𝑛 scattering source that offers several advantages relative to the standard scattering source representation when highly truncated Legendre cross-section expansions must be used. Galerkin quadrature methods are also useful when the scattering is moderately anisotropic, but the quadrature being used is not sufficiently accurate for the order of the scattering source expansion that is required. Furthermore, we derive the new method and present computational results showing that its performance for two challenging problems is comparable to those of the two GQ methods that preceded it.

Galerkin quadrature↗

Creep fatigue life prediction for engine hot section materials (isotropic)

The activities performed during the first year of the NASA HOST Program, Creep Fatigue Life Prediction for Engine Hot Section Materials (Isotropic), being conducted by Pratt & Whitney Aircraft are summarized. The program is a 5 year, two part effort aimed at improving the high temperature crack initiation prediction technology for gas turbine hot section components. Significant results of the program produced thus far are discussed. Cast B1900 + Hf and wrought IN 718 were selected as the base and alternate materials, respectively. A single heat of B1900 + Hf was obtained and test specimens fabricated. The material was characterized with respect to grain size, gamma prime size, carbide distribution, and dislocation density. Monotonic tensile and creep testing has shown engineering properties within anticipated scatter for this material. Examination of the tensile tests has shown a transition from inhomogeneous planar slip within the grains at lower temperatures to more homogeneous matrix deformation. Examination of the creep tests has shown a transgranular failure mode at 1400 F and an intergranular failure mode at 1600 F and 1800 F.

Moreno, V.↗

Polarimetric signatures of a coniferous forest canopy based on vector radiative transfer theory

Complete polarization signatures of a coniferous forest canopy are studied by the iterative solution of the vector radiative transfer equations up to the second order. The forest canopy constituents (leaves, branches, stems, and trunk) are embedded in a multi-layered medium over a rough interface. The branches, stems and trunk scatterers are modeled as finite randomly oriented cylinders. The leaves are modeled as randomly oriented needles. For a plane wave exciting the canopy, the average Mueller matrix is formulated in terms of the iterative solution of the radiative transfer solution and used to determine the linearly polarized backscattering coefficients, the co-polarized and cross-polarized power returns, and the phase difference statistics. Numerical results are presented to investigate the effect of transmitting and receiving antenna configurations on the polarimetric signature of a pine forest. Comparison is made with measurements.

Karam, M. A.↗

MOM3D method of moments code theory manual

MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.

Shaeffer, John F.↗

On the Retrieval of Ice Cloud Particle Shapes from POLDER Measurements

Shapes of ice crystals can significantly affect the radiative transfer in ice clouds. The angular distribution of the polarized reflectance over ice clouds strongly depends on ice crystal shapes. Although the angular-distribution features of the total or polarized reflectance over ice clouds implies a possibility of retrieving ice cloud particle shapes by use of remote sensing data, the accuracy of the retrieval must be evaluated. In this study, a technique, which applies single ice crystal habit and multidirectional polarized radiance to retrieve ice cloud particle shapes, is assessed. Our sensitivity studies show that the retrieved particle shapes &om this algorithm can be considered good approximations to those in actual clouds in calculation of the phase matrix elements. Although a fractal poly-crystal shape or an inhomogeneous hexagonal column may also produce this type of phase functions, more representative single-scattering properties from combinations of natural particle shapes and size distributions may still be necessary in accurate retrieval of other cloud properties such as optical thickness and particle size.

Sun, W.↗

The prediction of pseudo-resonance positions in the Schwinger variational principle

The Schwinger Variational Principle is appled to s-wave electron-hydrogen atom scattering. It is shown computationally, that, consistent with a paper by B. Apagyi, P. Levay, and K. Ladanyi, there are pseudo-resonances at the static exchange level of approximation, but not at the static level. The T-matrix was employed, as well as, the K-matrix version of the Schwinger Principle, with a real Slater basis, and obtained the same results in both. The origin of the pseudo-resonances as resulting from singularities in the separable potential that is effectively employed in the Lippmann-Schwinger equation from which the Schwinger Variational Principle can be derived. The determination of the pseudo-resonance parameters from the separable potential is computationally inexpensive and may be used to predict the pseudo-resonance parameters for the scattering calculations so that they may be avoided.

Weatherford, Charles A.↗

Properties of states near 𝐸 𝑥 = 6 MeV in 18 Ne through 17 F +𝑝 scattering

The rate of energy production in the hot-CNO cycle and breakout to the rapid-proton capture process in Type I X-ray bursts is strongly related to the 14 O⁡(𝛼,𝑝)⁢ 17 F reaction rate. The properties of states in 18 Ne near 𝐸 𝑥 = 6.1–6.3 MeV are important for understanding the 14 O⁡(𝛼,𝑝)⁢ 17 F reaction rate. In order to study 18 Ne resonances around this energy region, the RESOLUT radioactive-ion beam facility at Florida State University was used to perform 17 F⁡(𝑝,𝑝)⁢ 17 F elastic scattering on a polypropylene target under inverse kinematics. Scattered protons were detected in a silicon-strip detector array while recoiling 17 F ions were detected in coincidence in a gas ionization detector. An 𝑅-matrix analysis of measured cross sections was conducted along with a reanalysis of data from previous 17 F +𝑝 measurements. All the data analyzed are well described by a consistent set of parameters with a 1 − assignment for a state at 6.14(1) MeV. A second comparable solution is also found with a 3 − assignment for the 6.14(1) MeV state. There is no conclusive evidence supporting one of the two solutions over the other. The rates of the 14 O⁡(𝛼,p)⁢ 17 F reaction that are determined from the two solutions differ by up to an order of magnitude.

6 ≤ A ≤ 19↗

Scattering neutrinos, spin models, and permutations

We consider a class of Heisenberg all-to-all coupled spin models inspired by neutrino interactions in a supernova with N degrees of freedom. These models are characterized by a coupling matrix that is relatively simple in the sense that there are only a few, relative to N , nontrivial eigenvalues, in distinction to the classic Heisenberg spin-glass models, leading to distinct behavior in both the high-temperature and low-temperature regimes. When the momenta of the neutrinos are uniform and random in directions, we can calculate the large- N partition function for the quantum Heisenberg model. In particular, the high-temperature partition function predicts a non-Gaussian density of states, providing interesting counterexamples showing the limits of general theorems on the density of states for quantum spin models. We can repeat the same argument for classical Heisenberg models, also known as rotor models, and we find the high-temperature expansion is completely controlled by the eigenvalues of the coupling matrix, and again predicts non-Gaussian behavior for the density of states as long as the number of eigenvalues does not scale linearly with N . Indeed, we derive the amusing fact that these partition functions are essentially the generating function for counting permutations in the high-temperature regime. Finally, for the case relevant to neutrinos in a supernova, we identify the low-temperature phase as a unique state with the direction of the momenta of the neutrino dictating its coherent state in flavor-space, a state we dub the “flavor-momentum-locked” state. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Weak bosons as partons below 10 TeV partonic center-of-momentum

We investigate the modeling of weak boson number densities for leptons and hadrons in practical calculations in the Standard Model. Working in the framework of the Effective $W$ Approximation (EWA) and in $R_\xi$ and axial gauges, we derive the unrenormalized, tree-level parton number density functions for weak bosons from massless leptons at next-to-leading power in the collinear expansion. Corrections exhibit a number of properties, including those conjectured but never universally derived. Parallels with heavy quark factorization are also found. We avoid pathologies through a novel set of kinematical consistency conditions. When satisfied, good agreement between the full and approximated matrix elements is achieved. Findings suggest that the EWA may be testable at the Large Hadron Collider with $450$ fb$^{-1}$ luminosity of same-sign $WW$ scattering data at $\sqrt{s}=13.6$ TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Probing quantum entanglement with generalized parton distributions at the Electron-Ion Collider

Within the collinear factorization framework based on generalized parton distributions, we calculate the spin density matrix of exclusively produced quark and antiquark pairs $𝑢\bar{𝑢}$, $𝑑⁢\bar{d}$, $𝑠⁢\bar{𝑠}$, $𝑐\bar{⁢𝑐}$, $𝑏\bar{𝑏}$ in electron-proton scattering. The presence of both real and imaginary parts in the scattering amplitudes leads to a rich pattern of entanglement between the quark and the antiquark. We map out kinematical regions where the pairs exhibit entanglement, Bell nonlocality, and nonstabilizerness. We also predict that massive quarks and antiquarks are transversely polarized, similar to the well-known transverse hyperon polarization in unpolarized collisions. In strangeness, charm, and bottom productions, the polarization can reach 50%–80% in certain kinematic regions in the low-energy runs of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculations of polarization and radiance in the atmosphere

Two different methods for the solution of radiative transfer problems, matrix operator and Monte Carlo, are discussed. As an example of the use of the method, results are given for the radiance and polarization of the radiation scattered from haze layers as well as from models of the real atmosphere. The variation of the radiance, polarization, and ellipticity with the aerosol amount in the real atmosphere model is presented. The Monte Carlo method is applied to the problem of calculation of the radiance and polarization of the photons in the atmosphere when there is an ocean as the lower surface. It is found that the turbidity of the ocean can be deduced from the upwelling radiance measured at some height in the atmosphere. The downwelling radiation just beneath the ocean surface is elliptically polarized at those angles where it is derived from the total internal reflection of the upwelling radiation at the ocean surface.

Plass, G. N.↗

Mechanical Properties of Aerogels

Aerogels are extremely low density solids that are characterized by a high porosity and pore sizes on the order of nanometers. Their low thermal conductivity and sometimes transparent appearance make them desirable for applications such as insulation in cryogenic vessels and between double paned glass in solar architecture. An understanding of the mechanical properties of aerogels is necessary before aerogels can be used in load bearing applications. In the present study, the mechanical behavior of various types of fiber-reinforced silica aerogels was investigated with hardness, compression, tension and shear tests. Particular attention was paid to the effects of processing parameters, testing conditions, storage environment, and age on the aerogels' mechanical response. The results indicate that the addition of fibers to the aerogel matrix generally resulted in softer, weaker materials with smaller elastic moduli. Furthermore, the testing environment significantly affected compression results. Tests in ethanol show an appreciable amount of scatter, and are not consistent with results for tests in air. In fact, the compression specimens appeared to crack and begin to dissolve upon exposure to the ethanol solution. This is consistent with the inherent hydrophobic nature of these aerogels. In addition, the aging process affected the aerogels' mechanical behavior by increasing their compressive strength and elastic moduli while decreasing their strain at fracture. However, desiccation of the specimens did not appreciably affect the mechanical properties, even though it reduced the aerogel density by removing trapped moisture. Finally, tension and shear test results indicate that the shear strength of the aerogels exceeds the tensile strength. This is consistent with the response of brittle materials. Future work should concentrate on mechanical testing at cryogenic temperatures, and should involve more extensive tensile tests. Moreover, before the mechanical response of reinforced aerogels can be fully understood, more tests of unreinforced aerogels are necessary. Unreinforced aerogels are of particular use because their birefringent nature allows for visual determination of stress fields during mechanical testing. The success of any future tests depends on the availability of a large supply of quality specimens with well-documented preparation and storage histories.

Parmenter, Kelly E.↗

Scattering Properties of Heterogeneous Mineral Particles with Absorbing Inclusions

We analyze the results of numerically exact computer modeling of scattering and absorption properties of randomly oriented poly-disperse heterogeneous particles obtained by placing microscopic absorbing grains randomly on the surfaces of much larger spherical mineral hosts or by imbedding them randomly inside the hosts. These computations are paralleled by those for heterogeneous particles obtained by fully encapsulating fractal-like absorbing clusters in the mineral hosts. All computations are performed using the superposition T-matrix method. In the case of randomly distributed inclusions, the results are compared with the outcome of Lorenz-Mie computations for an external mixture of the mineral hosts and absorbing grains. We conclude that internal aggregation can affect strongly both the integral radiometric and differential scattering characteristics of the heterogeneous particle mixtures.

minerals↗

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗

Dynamic Interfacial Design in Adaptive Hybrid Materials Enables Reversible and Tunable Mechano-Optic Smart Responses

Next-generation polymeric materials are shifting toward adaptive and interactive behaviors of living systems; however, designing materials that can reversibly modulate optical properties under mechanical deformation while maintaining mechanical robustness remains a key challenge. Here, we report a mechanically robust vitrimer-based adaptive hybrid material (AHM) that exhibits a stretch-induced reversible transparency-to-opacity transition, enabled by the integration of dynamic interactions at the polymer–silica nanoparticle interface and controlled nanoparticle self-assembly. The AHM combines boronic ester–functionalized polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (S-Bpin) with diol-functionalized silica nanoparticles (diol-SiNPs) to form a hybrid network hosting both dynamic boronic ester and hydrogen-bonding interactions. These reversible linkages facilitate controlled nanoparticle self-assembly and enable strain-induced nanoparticle alignment/aggregation. Upon stretching, SiNP-rich domains align and aggregate within the polymer matrix, while local modulus mismatch between stiff aggregated SiNP/borylated-styrene-rich regions and the softer elastomeric midblock induces surface microwrinkle formation. These internal aggregates and surface wrinkles cooperatively enhance light scattering, producing the opaque state under strain. Furthermore, the tailored AHM exhibits high toughness, thermomechanical stability, reprocessability, and programmable shape-memory behavior. This work presents a dynamic interfacial design strategy for mechanically robust, optically reconfigurable, and reusable soft materials for adaptive optics, smart windows, sensing, soft robotics, and circular smart-material platforms.

adaptive hybrid materials↗

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Collisional quenching of O(1D) by O(3P)

Metastable O(1D) atoms may be quenched in collisions with ground state O(3P) atoms by transitions in the avoided crossing regions of the three lowest 3Pi(g) states of O2 of which the lowest separates to O(3P) + O(3P) and the two upper to O(3P) + O(1D). Quantal calculations of the adiabatic potential energy curves of the 3Pi(g) states are carried out with particular attention to an avoided crossing region in the lowest two states around a nuclear separation of 3.2a(0). Diabatic potential matrix elements are constructed from the adiabatic curves by imposing the requirement that they be smooth everywhere. A multi-state diabatic formulation is used to describe the scattering and the cross-sections for the collision-induced quenching of O(1D) atoms are calculated.

Yee, J.-H.↗