Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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 1,117 records · Page 62

A Single-column Model Ensemble Approach Applied to the TWP-ICE Experiment

Single-column models (SCM) are useful test beds for investigating the parameterization schemes of numerical weather prediction and climate models. The usefulness of SCM simulations are limited, however, by the accuracy of the best estimate large-scale observations prescribed. Errors estimating the observations will result in uncertainty in modeled simulations. One method to address the modeled uncertainty is to simulate an ensemble where the ensemble members span observational uncertainty. This study first derives an ensemble of large-scale data for the Tropical Warm Pool International Cloud Experiment (TWP-ICE) based on an estimate of a possible source of error in the best estimate product. These data are then used to carry out simulations with 11 SCM and two cloud-resolving models (CRM). Best estimate simulations are also performed. All models show that moisture-related variables are close to observations and there are limited differences between the best estimate and ensemble mean values. The models, however, show different sensitivities to changes in the forcing particularly when weakly forced. The ensemble simulations highlight important differences in the surface evaporation term of the moisture budget between the SCM and CRM. Differences are also apparent between the models in the ensemble mean vertical structure of cloud variables, while for each model, cloud properties are relatively insensitive to forcing. The ensemble is further used to investigate cloud variables and precipitation and identifies differences between CRM and SCM particularly for relationships involving ice. This study highlights the additional analysis that can be performed using ensemble simulations and hence enables a more complete model investigation compared to using the more traditional single best estimate simulation only.

numerical weather forecasting↗

Investigation of Primary Mirror Segment's Residual Errors for the Thirty Meter Telescope

The primary mirror segment aberrations after shape corrections with warping harness have been identified as the single largest error term in the Thirty Meter Telescope (TMT) image quality error budget. In order to better understand the likely errors and how they will impact the telescope performance we have performed detailed simulations. We first generated unwarped primary mirror segment surface shapes that met TMT specifications. Then we used the predicted warping harness influence functions and a Shack-Hartmann wavefront sensor model to determine estimates for the 492 corrected segment surfaces that make up the TMT primary mirror. Surface and control parameters, as well as the number of subapertures were varied to explore the parameter space. The corrected segment shapes were then passed to an optical TMT model built using the Jet Propulsion Laboratory (JPL) developed Modeling and Analysis for Controlled Optical Systems (MACOS) ray-trace simulator. The generated exit pupil wavefront error maps provided RMS wavefront error and image-plane characteristics like the Normalized Point Source Sensitivity (PSSN). The results have been used to optimize the segment shape correction and wavefront sensor designs as well as provide input to the TMT systems engineering error budgets.

warping harness↗

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

Modeling and Prediction of the Noise from Non-Axisymmetric Jets

The new source model was combined with the original sound propagation model developed for rectangular jets to produce a new version of the rectangular jet noise prediction code. This code was validated using a set of rectangular nozzles whose geometries were specified by NASA. Nozzles of aspect ratios two, four and eight were studied at jet exit Mach numbers of 0.5, 0.7 and 0.9, for a total of nine cases. Reynolds-averaged Navier-Stokes solutions for these jets were provided to the contactor for use as input to the code. Quantitative comparisons of the predicted azimuthal and polar directivity of the acoustic spectrum were made with experimental data provided by NASA. The results of these comparisons, along with a documentation of the propagation and source models, were reported in a journal article publication (Ref. 4). The complete set of computer codes and computational modules that make up the prediction scheme, along with a user's guide describing their use and example test cases, was provided to NASA as a deliverable of this task. The use of conformal mapping, along with simplified modeling of the mean flow field, for noise propagation modeling was explored for other nozzle geometries, to support the task milestone of developing methods which are applicable to other geometries and flow conditions of interest to NASA. A model to represent twin round jets using this approach was formulated and implemented. A general approach to solving the equations governing sound propagation in a locally parallel nonaxisymmetric jet was developed and implemented, in aid of the tasks and milestones charged with selecting more exact numerical methods for modeling sound propagation, and developing methods that have application to other nozzle geometries. The method is based on expansion of both the mean-flowdependent coefficients in the governing equation and the Green's function in series of orthogonal functions. The method was coded and tested on two analytically prescribed mean flows which were meant to represent noise reduction concepts being considered by NASA. Testing (Ref. 5) showed that the method was feasible for the types of mean flows of interest in jet noise applications. Subsequently, this method was further developed to allow use of mean flow profiles obtained from a Reynolds-averaged Navier-Stokes (RANS) solution of the flow. Preliminary testing of the generalized code was among the last tasks completed under this contract. The stringent noise-reduction goals of NASA's Fundamental Aeronautics Program suggest that, in addition to potentially complex exhaust nozzle geometries, next generation aircraft will also involve tighter integration of the engine with the airframe. Therefore, noise generated and propagated by jet flows in the vicinity of solid surfaces is expected to be quite significant, and reduced-order noise prediction tools will be needed that can deal with such geometries. One important source of noise is that generated by the interaction of a turbulent jet with the edge of a solid surface (edge noise). Such noise is generated, for example, by the passing of the engine exhaust over a shielding surface, such as a wing. Work under this task supported an effort to develop a RANS-based prediction code for edge noise based on an extension of the classical Rapid Distortion Theory (RDT) to transversely sheared base flows (Refs. 6 and 7). The RDT-based theoretical analysis was applied to the generic problem of a turbulent jet interacting with the trailing edge of a flat plate. A code was written to evaluate the formula derived for the spectrum of the noise produced by this interaction and results were compared with data taken at NASA Glenn for a variety of jet/plate configurations and flow conditions (Ref. 8). A longer-term goal of this task was to work toward the development of a high-fidelity model of sound propagation in spatially developing non-axisymmetric jets using direct numerical methods for solving the relevant equations. Working with NASA Glenn Acoustics Branch personnel, numerical methods and boundary conditions appropriate for use in a high-resolution calculation of the full equations governing sound propagation in a steady base flow were identified. Computer codes were then written (by NASA) and tested (by OAI) for an increasingly complex set of flow conditions to validate the methods. The NASA-supplied codes were ported to the High-End Computing resources of the NASA Advanced Supercomputing facility for testing and validation against analytical (where possible) and independent numerical solutions. The cases which were completed during the course of this contract were solutions of the two-dimensional linearized Euler equations with no mean flow, a uniform mean flow and a nonuniform mean flow representative of a parallel flow jet.

Leib, Stewart J.↗

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗

Structural dynamics of the renewable energy economy: A longitudinal input-output insights for a resilient transition

As countries accelerate their energy transitions, understanding how renewable energy (RE) systems structurally integrate into national economies is essential. This study presents a longitudinal economic input-output (EIO) analysis of the renewable energy sector in South Korea from 2016 to 2022. We develop a novel EIO-based framework that disaggregates the RE sector both by energy source (thermal, hydro, nuclear and renewable) and by industrial function (manufacturing, generation, and services), allowing for a detailed assessment of production dynamics, value-added creation, and import dependency. By quantifying backward and forward linkages and induced economic effects, the analysis reveals persistent structural vulnerabilities in renewable manufacturing and increasing sectoral interdependencies. Results reveal that while the renewable energy sector's production and value-added shares have increased, critical segments remain highly import-dependent, particularly in equipment manufacturing. The analysis highlights systemic gaps in domestic supply chain resilience and offers sector-specific insights for reducing vulnerability and enhancing energy security. Although applied to South Korea as a case study, the proposed framework is designed to be transferable to other national contexts where renewable energy planning requires economic structural insights. The findings offer policy-relevant guidance for enhancing domestic energy resilience and aligning industrial strategy with long-term decarbonization goals.

Economic linkage↗

Constraining Upper Troposphere/Lower Stratosphere Aerosol Physical Processes with High-Altitude Aircraft Measurements

Interest in a more complete understanding of the sources, composition and microphysics of stratospheric aerosol particles has intensified during recent years for several reasons: (1) small volcanic eruptions have been recognized as a driver of short-term changes in climate forcing; (2) emissions of sulfur dioxide (SO2) and other aerosol precursors have shifted to south Asia and other low latitude regions with intense vertical transport; (3) organic material has been recognized as a key contributor to lower stratospheric aerosol mass; and (4) interest in possible solar radiation management (geoengineering) through significant enhancements in stratospheric aerosols has intensified. To address stratospheric aerosol science issues, we are proposing a NASA Earth Ventures mission to NASA to provide extensive high-altitude aircraft measurements of critical gas-phase and aerosol properties at multiple locations across the planet. In this presentation, we will discuss the objectives of the proposed campaign, the measurements provided, the sampling strategy, and the modeling and analysis approaches that would be used to address specific science questions.

Jensen, Eric↗

Steps Toward Improved Integration, Search, and Analysis of Heterogeneous Data in the Astrobiology Habitable Environments Database

The Astrobiology Habitable Environments Database (AHED) is a new data system being developed as a long-term, open-access repository for astrobiology data. AHED is intended to store user-contributed results from NASA or externally-funded research in astrobiology, and to encourage sharing and synergy within the astrobiology community. However, the interdisciplinary nature of astrobiology presents some specific challenges to data management, integration, and analysis within AHED. In some disciplines (e.g., genomics), open databases thrive because the contributed products are fairly uniform and standardized (e.g., sequence data). In astrobiology, each investigation produces a unique set of data products; this makes it difficult to search across different datasets to find similar data, or to combine results from separate investigations. With AHED, we are taking steps to ensure there is adequate metadata - both at the dataset and record levels - to facilitate search, integration, and analysis. At the dataset level, we are developing a new metadata standard for describing astrobiology datasets, with detailed information about content, funding source, and scientific relevance, along with a set of topical keywords for characterizing datasets. At the record level, we are encouraging users to provide more structured content and finer-grained metadata. In many user-contributed science data repositories, few restrictions are placed on the uploaded data format, and minimal or no record-level metadata is required; thus users are unburdened when it comes to data preparation. The tradeoff is that deep integration and search across datasets is almost impossible without standardized structures and metadata. Although AHED users are free to upload minimally-described datasets, they will be encouraged to use database authoring tools (supplied by the underlying platform - Open Data Repository's Data Publisher) plus a set of customizable astrobiology-specific templates to help structure their data and provide standardized metadata. In reward for their extra effort, AHED will be able to deliver enhanced search, discovery, and analysis capabilities.

astrobiology↗

REDOTHERM: a thermodynamic modeling framework for redox-based thermochemical processes

Two-step thermochemical redox cycles are being developed as a potential pathway for the production of hydrogen and syngas. While there are many possible reactor and system configurations, moving oxide systems are considered promising in terms of the redox thermodynamics, due to the potential implementation of a countercurrent system that can achieve higher performance compared to other configurations. There is a lack of a robust thermodynamic modeling framework in the field, with multiple models incorporating incorrect thermodynamic assumptions that violate the second law of thermodynamics. We present in this work REDOTHERM, an open-source system model for moving oxides that incorporates the correct thermodynamic limits, as well as various options for the system auxiliary units including product separation, heat recovery, and oxygen removal. The model is agnostic to the energy source, and could be used for solar thermal or other configurations. We highlight the uses of this model, presenting some of the tradeoffs and challenges in redox-active material selection and how they affect the entire thermochemical hydrogen production process. This model could be easily adapted and used for material exploration, system/reactor design, and technoeconomic analysis.

08 HYDROGEN↗

Modeling Holocene Peatland Carbon Accumulation in North America

Peatlands are a large carbon reservoir. Yet the quantification of their carbon stock still has a large uncertainty due to lacking observational data and well-tested peatland biogeochemistry models. Here, a process-based peatland model was calibrated using long-term peat carbon accumulation data at multiple sites in North America. The model was then applied to quantify the peat carbon accumulation rates and stocks within North America over the last 12,000 years. We estimated that 85–174 Pg carbon was accumulated in North American peatlands over the study period including 0.37–0.76 Pg carbon in subtropical peatlands. During the period from 10,000 to 8,000 years ago, the warmer and wetter conditions might have played an important role in stimulating peat carbon accumulation by enhancing plant photosynthesis. Enhanced peat decomposition due to warming slowed the carbon accumulation through the rest of the Holocene. While recent modeling studies indicate that the northern peatlands will continue to act as a carbon sink in this century, our studies suggest that future enhanced peat decomposition accompanied by peatland areal changes induced by permafrost degradation and other disturbances shall confound the sink and source analysis.

54 ENVIRONMENTAL SCIENCES↗

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained (i.e., all complicated source code for the derivative calculations is constructed quickly with accuracy). The basic equations for second-order sensitivity derivatives are presented; four methods are compared. Each scheme requires that large systems are solved first for the first-order derivatives and, in all but one method, for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent flow example; both geometric-shape and flow-condition design variables are considered. Several methods are tested; results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method; for six independent variables, it is at least two to four times faster than central finite differences and requires only 60 percent more memory than the original code; the performance is expected to improve further in the future.

Taylor, Arthur C., III↗

The Development of Models for Carbon Dioxide Reduction Technologies for Spacecraft Air Revitalization

Through the respiration process, humans consume oxygen (O2) while producing carbon dioxide (CO2) and water (H2O) as byproducts. For long term space exploration, CO2 concentration in the atmosphere must be managed to prevent hypercapnia. Moreover, CO2 can be used as a source of oxygen through chemical reduction serving to minimize the amount of oxygen required at launch. Reduction can be achieved through a number of techniques. The National Aeronautics and Space Administration (NASA) is currently exploring the Sabatier reaction, the Bosch reaction, and co-electrolysis of CO2 and H2O for this process. Proof-of-concept experiments and prototype units for all three processes have proven capable of returning useful commodities for space exploration. While all three techniques have demonstrated the capacity to reduce CO2 in the laboratory, there is interest in understanding how all three techniques would perform at a system-level within a spacecraft. Consequently, there is an impetus to develop predictive models for these processes that can be readily re-scaled and integrated into larger system models. Such analysis tools provide the ability to evaluate each technique on a comparable basis with respect to processing rates. This manuscript describes the current models for the carbon dioxide reduction processes under parallel developmental e orts. Comparison to experimental data is provided were available for veri cation purposes.

Swickrath, Michael J.↗

Geospatial Cluster Analysis of the State, Duration and Severity of Drought Over Paraíba State, Northeastern Brazil

Droughts threaten water resources, agriculture, socio-economic activities and the population at the global and regional level, so identifying areas with homogeneous drought behaviors is an important consideration in improving the management of water resources. The objective of this study is to identify homogenous zones over Paraíba State in relation to the state, duration and severity of droughts that have occurred over the last 20 years (1998–2017) using hierarchical cluster analysis based on both gauge-measured and Tropical Rainfall Measuring Mission (TRMM) estimated rainfall data (TMPA 3B42). The drought series were calculated using the Standardized Precipitation Index (SPI) based on eight time scales and were grouped according to drought state, duration and severity time series. The integrated results of state, duration and severity of droughts indicate that there is a basis for dividing Paraíba State into two major regions (a) Zone I, formed by Mata Paraibana and Agreste Paraibano, and (b) Zone II, composed by Borborema and Sertão Paraibano. This division is evident when assessing short-term droughts, but in the case of long-term droughts, Paraíba State has a high similarity in terms of drought state, duration, and severity. Factors such as proximity to the ocean, active climatic systems, and the local relief configuration were identified as influencing the drought regime. Finally, it is concluded that TMPA rainfall estimates represent a valuable source of data to regionalize and identify drought patterns over this part of Brazil and that other studies of this type should be carried out to monitor these phenomena based on other satellite-based rainfall data, including the Global Precipitation Mission (GPM).

SPI↗

Comparison of Single-Event Transients in an Epitaxial Silicon Diode Resulting from Heavy Ion-, Focused X-Ray-, and Pulsed Laser-Induced Charge Generation

Heavy ion, focused x-ray, and pulsed laser single event transient experiments are performed on a silicon epitaxial diode. Collected charge, transient rise times, and transient fall times are calculated and compared between the different sources. It is observed that these transient shape characteristics depend on the source (ion, x-ray, or laser), even when similar amounts of charge are generated. Finally, the observed differences are examined and explained in terms of basic charge collection mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Real-space texture and pole-figure analysis using the 3D pair distribution function on a platinum thin film

An approach is described for studying texture in nanostructured materials. The approach implements the real-space texture pair distribution function (PDF), txPDF, laid out by Gong & Billinge {(2018). arXiv:1805.10342 [cond-mat]}. It is demonstrated on a fiber-textured polycrystalline Pt thin film. The approach uses 3D PDF methods to reconstruct the orientation distribution function of the powder crystallites from a set of diffraction patterns, taken at different tilt angles of the substrate with respect to the incident beam, directly from the 3D PDF of the sample. A real-space equivalent of the reciprocal-space pole figure is defined in terms of interatomic vectors in the PDF and computed for various interatomic vectors in the Pt film. Furthermore, it is shown how a valid isotropic PDF may be obtained from a weighted average over the tilt series, including the measurement conditions for the best approximant to the isotropic PDF from a single exposure, which for the case of the fiber-textured film was in a nearly grazing incidence orientation of ∼10°. Finally, an open-source Python software package, FouriGUI , is described that may be used to help in studies of texture from 3D reciprocal-space data, and indeed for Fourier transforming and visualizing 3D PDF data in general.

3D pair distribution functions↗

Exploration of the polarization angle variability of the Crab Nebula with POLARBEAR and its application to the search for axionlike particles

The Crab Nebula, also known as Tau A, is a polarized astronomical source at millimeter wavelengths. It has been used as a stable light source for polarization angle calibration in millimeter-wave astronomy. However, it is known that its intensity and polarization vary as a function of time at a variety of wavelengths. Thus, it is of interest to verify the stability of the millimeter-wave polarization. If detected, polarization variability may be used to better understand the dynamics of Tau A, and for understanding the validity of Tau A as a calibrator. One intriguing application of such observation is to use it for the search of axionlike particles (ALPs). Ultralight ALPs couple to photons through a Chern-Simons term, and induce a temporal oscillation in the polarization angle of linearly polarized sources. After assessing a number of systematic errors and testing for internal consistency, we evaluate the variability of the polarization angle of the Crab Nebula using 2015 and 2016 observations with the 150 GHz P instrument. We place a median 95% upper bound of polarization oscillation amplitude A < 0.06 5 ° over the oscillation frequencies from 0.75 year − 1 to 0.66 hour − 1 . Assuming that no sources other than ALP are causing Tau A’s polarization angle variation, that the ALP constitutes all the dark matter, and that the ALP field is a stochastic Gaussian field, this bound translates into a median 95% upper bound of ALP-photon coupling g a γ γ < 2.16 × 10 − 12 GeV − 1 × ( m a / 10 − 21 eV ) in the mass range from 9.9 × 10 − 23 eV to 7.7 × 10 − 19 eV . This demonstrates that this type of analysis using bright polarized sources is as competitive as those using the polarization of cosmic microwave background in constraining ALPs. Published by the American Physical Society 2024

Adachi, Shunsuke (ORCID:0000000204007555)↗