Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete dipole approximation”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Single-Scattering Properties of Melting Precipitation for GPM Passive Microwave and Radar Remote Sensing Applications

Over the past two decades, detailed computational simulations of the intricate three dimensional structures of ice-phase crystals and aggregates of those crystals have been developed. The microwave single-scattering properties of these simulated ice particles have been computed and used to improve quantitative estimates of snow rates and to better deRne the vertical structure of snow water contents in deep convective systems, as derived from satellite-borne passive microwave and/or radar remote sensing measurements. The same icephase particles have more recently been used as the starting point for simulations of melting precipitation using computational melting methods. In the current study, a heuristic melting method, as well as a physically-based melting procedure based on smoothed-particle hydrodynamics, are applied to ice particle models to describe the full evolution of the particles from dry snow to liquid drops. The discrete dipole approximation is utilized to calculate the single-scattering properties of the mixed-phase particles throughout the melting process. Then, the properties of the particles are “mapped” into simpliRed microphysical simulations of particle size spectra in the melting layers of stratiform, precipitating clouds. The bulk single-scattering properties of the melting layers and the sensitivity of their properties to modeling assumptions are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
NASA NTRS2021

Towards a Mass-Consistent Methodology for Realistic Melting Hydrometeor Retrieval

To address the acute challenge posed by the melting layer to accurate surface precipitation retrievals from space, we en-sure the compositional consistency in ice, liquid, and total masses of synthetic melting hydrometeors with a method of stochastic compensation. The method is applied to simulated melting hydrometeors prior to calculating their scattering properties using the discrete dipole approximation (DDA). We investigate the impact of this stochastic compensation to calculated scattering properties by contrasting it with a naïve approach and report our findings.

precipitation remote sensing↗
NASA NTRS2021

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗
NASA NTRSDate not supplied

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
NASA NTRS2023

The Polarimetric Radar Scattering Properties of Oriented Aggregates

The scattering properties of aggregates are studied herein. Early aggregates (<7 monomers) of branched planar crystals and mature aggregates (up to 100 monomers) of columns are randomly generated with varying assumptions about the monomer attachment processes and the orientation behavior during collection. The resulting physical properties of the aggregates correspond well with prior in situ and retrieved sizes and shapes. Assumed azimuthally uniform orientations during collection and monomer pivoting upon attachment resulted in flatter and denser aggregates. The column aggregates had lower density and more spherical shapes than the branched planar crystal aggregates. The scattering properties were calculated using the discrete dipole approximation for a set of orientation angles and transformed to spectral coefficients representing modes of orientation angle variability. The zeroth- and second-order coefficients dominate this variability, with the zeroth-order coefficients representing the scattering properties for randomly oriented particles. The second-order coefficients for backscatter showed differences between horizontal and vertical polarization increasing with density, and these coefficients for specific differential phase increase with both mass and density. Similarly, coefficients for the copolar covariance decreased with density. Rapid changes in the contributions to the radar moments from the second-order coefficients from low to moderate density were observed, likely due to the increasing presence of horizontally aligned monomers in the aggregate structure. Differences in how differential reflectivity and correlation coefficient evolve with the orientation distribution parameters suggest that these measurements, along with specific differential phase and reflectivity, provide complementary information about aggregate sizes, shapes, and orientation distributions.

Cloud microphysics↗
NASA NTRS2023

Developing a Radar Signal Simulator for the Community Radiative Transfer Model

Active radar instruments provide vertically resolved clouds and precipitation measurements that cannot be provided by the passive instruments. These active measurements are not conventionally assimilated into the data assimilation systems because of the lack of fast forward radiative transfer models and also difficulties in the error modelling of the measurements. This paper describes the development, evaluation, and sensitivity analysis for a forward radar model implemented in the Community Radiative Transfer Model (CRTM). The scattering properties required by the forward model are provided by the hydrometeor lookup tables that were generated using the discrete dipole approximation. The model is able to calculate both the reflectivity and the attenuated reflectivity for any given radar instrument at any given zenith angles as long as CRTM instrument specific coefficients are available. The evaluation using CloudSat measurements shows a very good agreement between the simulations and measurements as long as the input profiles of hydrometeors are consistent with the measured reflectivity profiles. Major sources contributing to the differences between the measured and simulated reflectivities are input hydrometeor profiles, scattering lookup tables, lack of melting layer in the forward model, CRTM scattering solvers, and attenuation calculations. In addition to the forward model, both Tangent Linear and Adjoint of the model are also implemented and tested within CRTM. These components may be required by some data assimilation systems for the assimilation of radar measurements.

radar↗
NASA NTRSDate not supplied

A Comprehensive Forward Model for Spaceborne Radar Instruments

We present the development and validation of a comprehensive forward model designed to enhance remote sensing capabilities of spaceborne radar instruments. To overcome limitations in existing models, we integrated a Discrete Dipole Approximation (DDA) cloud scattering database into our Radiative Transfer Model (RTM), focusing on microwave frequencies. By simulating the optical properties of non-spherical frozen hydrometeors, the DDA technique effectively reduced discrepancies between simulated and observed values, surpassing traditional Mie tables. The evaluation of DDA lookup tables involved comparisons with a collocated dataset comprising short-term forecasts and satellite microwave data, providing evidence of their superiority. Additionally, we address the challenges of assimilating active radar measurements, which offer vertically resolved insights into clouds and precipitation. We explored the assimilation of spaceborne radar measurements in Numerical Weather Prediction (NWP) models by integrating a forward radar model, along with its adjoint and tangent linear, into the data assimilation system. Evaluation using CloudSat measurements demonstrated promising agreement between simulations and observations, particularly when the input hydrometeor profiles aligned with the measured reflectivity profiles, showcasing the potential of the developed forward radar model. Moreover, we discuss other challenges in radar measurement assimilation within NWP models, including potential observation errors and biases.

Isaac Moradi↗
NASA NTRSDate not supplied

Assimilation of Active MW and Radar Observations in the NWP Models

Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. This talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗
NASA NTRSDate not supplied

Simulating Microwave and Radar Signals in Severe Weather Conditions

The development of radiative transfer simulators for radar and microwave signals, spanning a range of frequencies from 10 to 800 GHz, is paramount for enhancing weather forecasting accuracy, particularly for severe weather events. These tools facilitate the assimilation of microwave and radar observations into numerical weather prediction models, thereby improving accuracy of weather forecasts. Additionally, they directly simulate radar signals crucial for autonomous vehicle operation, with frequencies commonly used in radar systems such as 24, 74, 77, and 79 GHz. However, these frequencies are susceptible to weather phenomena like severe rain and snow, which can significantly impact vehicle safety and performance. These simulation tools, including radiative transfer simulators, serve a dual purpose. Firstly, they enhance forecasts for severe weather events, contributing to autonomous vehicle safety by providing early warnings and risk assessments. Secondly, they enable the direct simulation of radar signals in autonomous vehicle driving systems, allowing researchers and engineers to evaluate radar system performance under various weather conditions. In addition to the radar signal simulator, this abstract discusses the incorporation of advanced scattering properties developed using the discrete dipole approximation (DDA). The DDA technique enhances scattering calculations for frozen hydrometeors at microwave frequencies, thereby improving the accuracy of radar signal simulations and enabling more accurate assessments of radar system performance in adverse weather conditions. In summary, this abstract explores the development and utilization of comprehensive simulation tools, emphasizing their significance in simulating microwave and radar signals and improving weather forecasts. Special attention is given to the simulation of radar signals at critical frequencies for autonomous vehicle sensing and navigation, addressing challenges posed by severe weather phenomena and their effects on signal propagation and detection

Isaac Moradi↗
DOE OSTI2024

Finite-Temperature Quantum Matter with Rydberg or Molecule Synthetic Dimensions

Synthetic-dimension platforms offer unique pathways for engineering quantum matter. We compute the phase diagram of a many-body system of ultracold atoms (or polar molecules) with a set of Rydberg states (or rotational states) as a synthetic dimension, where the particles are arranged in real space in optical microtrap arrays and interact via dipole-dipole exchange interaction. Using mean-field theory, we find three ordered phases—two are localized in the synthetic dimension, predicted as zero-temperature ground states, and one is a delocalized phase. We characterize them by identifying the spontaneously broken discrete symmetries of the Hamiltonian. We also compute the phase diagram as a function of temperature and interaction strength for both signs of the interaction. For system sizes with more than six synthetic sites and attractive interactions, we find that the thermal phase transitions can be first or second order, which leads to a tricritical point on the phase boundary. Furthermore, by examining the dependence of the tricritical point and other special points of the phase boundary on the synthetic dimension size, we shed light on the physics for thermodynamically large synthetic dimension.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Visualization of Plasmonic Couplings Using Ultrafast Electron Microscopy

Hybrids of graphene and metal plasmonic nanostructures are promising building blocks for applications in optoelectronics, surface-enhanced scattering, biosensing, and quantum information. An understanding of the coupling mechanism in these hybrid systems is of vital importance to its applications. Previous efforts in this field mainly focused on spectroscopic studies of strong coupling within the hybrids with no spatial resolution. In this study, we report direct imaging of the local plasmonic coupling between single Au nanocapsules and graphene step edges at the nanometer scale by photon-induced near-field electron microscopy in an ultrafast electron microscope for the first time. The proximity of a step in the graphene to the nanocapsule causes asymmetric surface charge density at the ends of the nanocapsules. Computational electromagnetic simulations confirm the experimental observations. The results reported here indicate that this hybrid system could be used to manipulate the localized electromagnetic field on the nanoscale, enabling promising future plasmonic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1993

Scattering Properties of Aggregate Particles

We are engaged in modeling the scattering and emitting properties of irregular particles with the goal of understanding the composition and structure of cometary dust particles.

discrete dipole approximation non-spherical monome↗
NASA NTRS1982

Theoretical studies of photoexcitation and ionization in H2O

Theoretical studies using Franck-Condon and static-exchange approximations are reported for the complete dipole excitation and ionization spectrum in H2O, where (1) large Cartesian Gaussian basis sets are used to represent the required discrete and continuum electronic eigenfunctions at the ground state equilibrium geometry, and (2) previously devised moment-theory techniques are employed in constructing the continuum oscillator-strength densities from the calculated spectra. Comparisons are made of the calculated excitation and ionization profiles with recent experimental photoabsorption studies and corresponding spectral assignments, electron impact-excitation cross sections, and dipole and synchrotron-radiation studies of partial-channel photoionization cross sections. The calculated partial-channel cross sections are found to be atomic-like, and dominated by 2p-kd components. It is suggested that the latter transition couples with the underlying 1b(1)-kb(1) channel, accounting for a prominent feature in recent synchrotron-radiation measurements.

Diercksen, G. H. F.↗
NASA NTRS1981

Photoexcitation and ionization in carbon dioxide - Theoretical studies in the separated-channel static-exchange approximation

Vertical-electronic static-exchange photoexcitation and ionization cross sections are reported which provide a first approximation to the complete dipole spectrum of CO2. Separated-channel static-exchange calculations of vertical-electronic transition energies and oscillator strengths, and Stieltjes-Chebyshev moment methods were used in the development. Detailed comparisons were made of the static-exchange excitation and ionization spectra with photoabsorption, electron-impact excitation, and quantum-defect estimates of discrete transition energies and intensities, and with partial-channel photoionization cross sections obtained from fluorescence measurements and from tunable-source and (e, 2e) photoelectron spectroscopy. Results show that the separate-channel static-exchange approximation is generally satisfactory in CO2.

Padial, N.↗
NASA NTRS2010

Charge-Spot Model for Electrostatic Forces in Simulation of Fine Particulates

The charge-spot technique for modeling the static electric forces acting between charged fine particles entails treating electric charges on individual particles as small sets of discrete point charges, located near their surfaces. This is in contrast to existing models, which assume a single charge per particle. The charge-spot technique more accurately describes the forces, torques, and moments that act on triboelectrically charged particles, especially image-charge forces acting near conducting surfaces. The discrete element method (DEM) simulation uses a truncation range to limit the number of near-neighbor charge spots via a shifted and truncated potential Coulomb interaction. The model can be readily adapted to account for induced dipoles in uncharged particles (and thus dielectrophoretic forces) by allowing two charge spots of opposite signs to be created in response to an external electric field. To account for virtual overlap during contacts, the model can be set to automatically scale down the effective charge in proportion to the amount of virtual overlap of the charge spots. This can be accomplished by mimicking the behavior of two real overlapping spherical charge clouds, or with other approximate forms. The charge-spot method much more closely resembles real non-uniform surface charge distributions that result from tribocharging than simpler approaches, which just assign a single total charge to a particle. With the charge-spot model, a single particle may have a zero net charge, but still have both positive and negative charge spots, which could produce substantial forces on the particle when it is close to other charges, when it is in an external electric field, or when near a conducting surface. Since the charge-spot model can contain any number of charges per particle, can be used with only one or two charge spots per particle for simulating charging from solar wind bombardment, or with several charge spots for simulating triboelectric charging. Adhesive image-charge forces acting on charged particles touching conducting surfaces can be up to 50 times stronger if the charge is located in discrete spots on the particle surface instead of being distributed uniformly over the surface of the particle, as is assumed by most other models. Besides being useful in modeling particulates in space and distant objects, this modeling technique is useful for electrophotography (used in copiers) and in simulating the effects of static charge in the pulmonary delivery of fine dry powders.

Walton, Otis R.↗
DOE OSTI2026

Non-Markovian relaxation spectroscopy of fluxonium qubits

Recent studies have shown that parasitic two-level systems (TLS) in superconducting qubits, which are a leading source of decoherence, can have relaxation times longer than the qubits themselves. However, the standard techniques used to characterize qubit relaxation is only valid for measuring T 1 under the Born-Markov approximation and could mask environmental memory effects in practice. Here, we introduce two-timescale relaxometry, a technique to probe the qubit and environment relaxation simultaneously and efficiently. We apply it to high-coherence fluxonium qubits over a frequency range of 0.1-0.4 GHz, and reveal a discrete spectrum of TLS with millisecond lifetimes. Our analysis of the spectrum is consistent with a random distribution of TLS in the aluminum oxide tunnel barrier of the Josephson junction chain of the fluxonium, with a spectral and volumetric density and average electric dipole similar to previous TLS studies at much higher frequencies. Our study suggests that investigating and mitigating TLS in the junction chain is crucial to the development of various types of noise-protected qubits in circuit QED.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Primary $\gamma$-ray intensities and $\gamma$-strength functions from discrete two-step $\gamma$-ray cascades in radiative proton-capture experiments

Background: Reaction rates of radiative capture reactions can play a crucial role in the nucleosynthesis of heavy nuclei in explosive stellar environments. These reaction rates depend strongly on $\gamma$-ray decay widths in the reaction products, which are, for nonresonant capture reactions at high excitation energies, derived from the $\gamma$-ray strength function and the nuclear level density. Recently, the ratio method was applied to primary $\gamma$ rays observed from ($\textit{d, p}$) reactions and nuclear resonance fluorescence measurements to extract the dipole strength in atomic nuclei and to test the generalized Brink-Axel hypothesis. Purpose: The purpose of this work is to apply the ratio method to primary $\gamma$-ray intensities of the 63,65 Cu($\textit{p},\gamma$) reactions to extract $\gamma$-ray strength information on the nuclei 64,66 Zn. Here, the impact of spin distribution, total $\gamma$-ray decay widths, level densities, and width fluctuations on the application of the ratio method will be discussed. Additionally, by comparing the relative $\gamma$-ray strength at different excitation energies, conclusions on the validity of the generalized Brink-Axel hypothesis can be made. Method: The radiative proton capture reaction measurements have been performed at the HORUS $\gamma$-ray spectrometer of the University of Cologne at one excitation energy for each reaction. Primary $\gamma$-ray intensities have been determined by normalizing secondary $\gamma$-ray transitions in two-step cascades using their absolute branching ratio. The ratio method was applied to the measured primary $\gamma$-ray intensities as well as to previous measurements by Erlandsson et al. at different excitation energies. Results: The relative strength function curve for 64 Zn from our measurement shows no significant deviation from the previous measurement at a different excitation energy. The same is true for 66 Zn where both measurements were at almost the same excitation energy. Absolute $\gamma$-strength function values have been obtained by normalizing the relative curves to quasiparticle random phase approximation calculations because of the absence of experimental data in the respective energy region. Conclusion: The generalized Brink-Axel hypothesis, i.e., the independence of the strength function on the excitation energy, seems to hold in the studied energy region and nuclei. The method to obtain primary $\gamma$-ray intensities from two-step cascade spectra was shown to be a valuable and sensitive tool although its uncertainties are connected to the knowledge of the low-energy level scheme of the investigated nucleus. The scaling in the ratio method should be taken with care, because the relative strength is not a simple sum of $f_{E1}$ and $f_{M1}$ but a somewhat complex linear combination dependent on the excitation energy of the nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗