Search NASA⌕ Search

SEARCH · Search NASA

Results for “radiation transfer”

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 163 records · Page 9

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. II. Connection with γ-Rays

Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LAT γ -ray bands as well as closely correlated optical and γ -ray flares. The swings in optical polarization angle are also accompanied by γ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast γ -ray flares or orphan γ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

79 ASTRONOMY AND ASTROPHYSICS↗

Separating radiative forcing by aerosol–cloud interactions and rapid cloud adjustments in the ECHAM–HAMMOZ aerosol–climate model using the method of partial radiative perturbations

Using the method of offline radiative transfer modeling within the partial radiative perturbation (PRP) approach, the effective radiative forcing by aerosol–cloud interactions (ERF aci ) in the ECHAM–HAMMOZ aerosol climate model is decomposed into a radiative forcing by anthropogenic cloud droplet number change and adjustments of the liquid water path and cloud fraction.The simulated radiative forcing by anthropogenic cloud droplet number change and liquid water path adjustment are of approximately equal magnitude at –0.52 and –0.53 W m –2 , respectively, while the cloud-fraction adjustment is somewhat weaker at –0.31 W m –2 (constituting 38 %, 39 %, and 23 %of the total ERFaci, respectively); geographically, all three ERF aci components in the simulation peak over China, the subtropical eastern ocean boundaries,the northern Atlantic and Pacific oceans, Europe, and eastern North America (in order of prominence). Spatial correlations indicate that the temporal-mean liquid water path adjustment is proportional to the temporal-mean radiative forcing, while the relationship between cloud-fraction adjustment and radiative forcing is less direct. While the estimate of warm-cloud ERF aci is relatively insensitive to the treatment of ice and mixed-phase cloud overlying warm cloud, there are indications that more restrictive treatments of ice in the column result in a low bias in the estimated magnitude of the liquid water path adjustment and a high bias in the estimated magnitude of the droplet number forcing. Since the present work is the first PRP decomposition of the aerosol effective radiative forcing into radiative forcing and rapid cloud adjustments, idealized experiments are conducted to provide evidence that the PRP results are accurate. The experiments show that using low-frequency(daily or monthly) time-averaged model output of the cloud property fields underestimates the ERF, but 3-hourly mean output is sufficiently frequent.

54 ENVIRONMENTAL SCIENCES↗

Beyond the visible: Accounting for ultraviolet and far‐red radiation in vegetation productivity and surface energy budgets

Photosynthetically active radiation (PAR) is typically defined as light with a wavelength within 400–700 nm. However, ultra-violet (UV) radiation within 280–400 nm and far-red (FR) radiation within 700–750 nm can also excite photosystems, though not as efficiently as PAR. Vegetation and land surface models (LSMs) typically do not explicitly account for UV's contribution to energy budgets or photosynthesis, nor FR's contribution to photosynthesis. However, whether neglecting UV and FR has significant impacts remains unknown. Here, we explored how canopy radiative transfer (RT) and photosynthesis are impacted when explicitly implementing UV in the canopy RT model and accounting for UV and FR in the photosynthesis models within a next-generation LSM that can simulate hyperspectral canopy RT. We validated our improvements using photosynthesis measurements from plants under different light sources and intensities and surface reflection from an eddy-covariance tower. Our model simulations suggested that at the whole plant level, after accounting for UV and FR explicitly, chlorophyll content, leaf area index (LAI), clumping index, and solar radiation all impact the modeling of gross primary productivity (GPP). At the global scale, mean annual GPP within a grid would increase by up to 7.3% and the increase is proportional to LAI; globally integrated GPP increases by 4.6 PgC year −1 (3.8% of the GPP without accounting for UV + FR). Further, using PAR to proxy UV could overestimate surface albedo by more than 0.1, particularly in the boreal forests. Our results highlight the importance of improving UV and FR in canopy RT and photosynthesis modeling and the necessity to implement hyperspectral or multispectral canopy RT schemes in future vegetation and LSMs.

energy budget↗

Radiative heating rate profiles over the southeast Atlantic Ocean during the 2016 and 2017 biomass burning seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96° S, 14.35° W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d -1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the PV Resource Dataset in Central and North America: Preprint

A recent effort from the National Renewable Energy Laboratory (NREL) has led to a new radiative transfer model, Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT), to efficiently compute plane-of-array (POA) irradiance. This model has been implemented by the National Solar Radiation Data Base (NSRDB) to provide photovoltaic (PV) resource over Central and North America in both narrow- and broad wavelength bands. This study conducts a comprehensive evaluation of the PV resource dataset using surface-based observations in various climate zones. The results demonstrate the PV resource has an excellent agreement with the cloudy-sky observations from thermopiles and reference cells.

41 EE - Solar Energy Technologies Office (EE-4S)↗

CF + excitation in the interstellar medium

The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Vegetation clumping modulates global photosynthesis through adjusting canopy light environment

Abstract The spatial dispersion of photoelements within a vegetation canopy, quantified by the clumping index (CI), directly regulates the within‐canopy light environment and photosynthesis rate, but is not commonly implemented in terrestrial biosphere models to estimate the ecosystem carbon cycle. A few global CI products have been developed recently with remote sensing measurements, making it possible to examine the global impacts of CI. This study deployed CI in the radiative transfer scheme of the Community Land Model version 5 (CLM5) and used the revised CLM5 to quantitatively evaluate the extent to which CI can affect canopy absorbed radiation and gross primary production (GPP), and for the first time, considering the uncertainty and seasonal variation of CI with multiple remote sensing products. Compared to the results without considering the CI impact, the revised CLM5 estimated that sunlit canopy absorbed up to 9%–15% and 23%–34% less direct and diffuse radiation, respectively, while shaded canopy absorbed 3%–18% more diffuse radiation across different biome types. The CI impacts on canopy light conditions included changes in canopy light absorption, and sunlit–shaded leaf area fraction related to nitrogen distribution and thus the maximum rate of Rubisco carboxylase activity ( V cmax ), which together decreased photosynthesis in sunlit canopy by 5.9–7.2 PgC year −1 while enhanced photosynthesis by 6.9–8.2 PgC year −1 in shaded canopy. With higher light use efficiency of shaded leaves, shaded canopy increased photosynthesis compensated and exceeded the lost photosynthesis in sunlit canopy, resulting in 1.0 ± 0.12 PgC year −1 net increase in GPP. The uncertainty of GPP due to the different input CI datasets was much larger than that caused by CI seasonal variations, and was up to 50% of the magnitude of GPP interannual variations in the tropical regions. This study highlights the necessity of considering the impacts of CI and its uncertainty in terrestrial biosphere models.

59 BASIC BIOLOGICAL SCIENCES↗

Tree-based solvers for adaptive mesh refinement code $\scriptsize{FLASH}$ – IV. An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

X-ray radiation, in particular radiation between 0.1 and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the $\scriptsize{TREERAY}$ reverse ray trace algorithm implemented in the flash magnetohydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV < E γ < 100 keV, dubbed $\scriptsize{XRAYTHESPOT}$. $\scriptsize{XRAYTHESPOT}$ allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross-sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multibin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disc and simulations of molecular clouds with active chemistry, radiation pressure, and protostellar radiation feedback from infrared to X-ray radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of Aerosol Embedded in Shallow Cumulus Cloud Fields on the Surface Solar Irradiance

Ubiquitous shallow cumulus clouds are associated with complex variability in surface solar irradiance (SSI). Aerosol embedded in the cloud field typically has a much smaller overall radiative effect, but can significantly perturb the shape of the SSI probability density function (PDF). These perturbations have important implications for several applications that utilize SSI, but are poorly quantified and are the subject of this study. Multiple cases of shallow cumulus cloud fields with embedded aerosol are simulated at the Southern Great Plains Atmospheric Observatory using large eddy simulation (LES). The LES-derived cloud and aerosol fields are then ingested into Monte Carlo three-dimensional (3D) radiative transfer to simulate SSI. We find a variety of perturbations to the SSI PDF that depend on aerosol presence and optical properties. The processes leading to these perturbations include extinction of the direct beam that often increases from the clear-sky region toward cloud edge due to aerosol hygroscopic growth, and scattering of radiation by aerosol into cloud shadows. The ability to predict the SSI PDF in the presence of aerosol is assessed by adding three representative aerosol optical properties into an existing machine learning framework. We show that machine learning accurately predicts the SSI PDF across a wide range of conditions with negligible computational expense. Importance metrics reveal the relatively high influence of aerosol optical properties in making the predictions. These new findings highlight the important role that aerosol plays in SSI variability for highly 3D cloud-aerosol environments and provides a computationally efficient route forward for its simulation.

3D radiative transfer↗

Radiative impact of record-breaking wildfires from integrated ground-based data

The radiative effects of wildfires have been traditionally estimated by models using radiative transfer calculations. Assessment of model-predicted radiative effects commonly involves information on observation-based aerosol optical properties. However, lack or incompleteness of this information for dense plumes generated by intense wildfires reduces substantially the applicability of this assessment. Here we introduce a novel method that provides additional observational constraints for such assessments using widely available ground-based measurements of shortwave and spectrally resolved irradiances and aerosol optical depth (AOD) in the visible and near-infrared spectral ranges. We apply our method to quantify the radiative impact of the record-breaking wildfires that occurred in the Western US in September 2020. For our quantification we use integrated ground-based data collected at the Atmospheric Measurements Laboratory in Richland, Washington, USA with a location frequently downwind of wildfires in the Western US. We demonstrate that remarkably dense plumes generated by these wildfires strongly reduced the solar surface irradiance (up to 70% or 450 Wm -2 for total shortwave flux) and almost completely masked the sun from view due to extremely large AOD (above 10 at 500 nm wavelength). We also demonstrate that the plume-induced radiative impact is comparable in magnitude with those produced by a violent volcano eruption occurred in the Western US in 1980 and continental cumuli.

54 ENVIRONMENTAL SCIENCES↗

An Evaluation of the Spectral Irradiance Data from the NSRDB

A recent effort by the National Renewable Energy Laboratory (NREL) has led to a new radiative transfer model, the Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT), to efficiently compute spectral irradiances in the plane-of-array (POA). This model has been implemented in the National Solar Radiation Data Base (NSRDB) to provide photovoltaic (PV) resource in both narrowband and broadband wavelengths. This study conducts an evaluation of the spectral irradiances from the PV Resource dataset using surface-based observations at NREL’s Solar Radiation Research Laboratory (SRRL) and University of Oregon (UO). The results demonstrate that the PV resource has a generally good agreement with the long-term observations in both clear-sky and cloudy-sky conditions. Further research is needed to reduce the overestimation of visible irradiances in clear sky conditions and underestimation of near-infrared irradiances in cloudy-sky conditions.

14 SOLAR ENERGY↗

Profiles of Radiative Fluxes at ENA

Profiles of radiative fluxes observed at the Atmospheric Radiation Measurement (ARM)’s Eastern North Atlantic (ENA) observatory along with the ancillary measurements are reported. The below-cloud drizzle properties were derived by combining the data from the ceilometer and Ka-band ARM Zenith Radar (KAZR) following the technique explained by Ghate et al. (2021 JAMC). The cloud and drizzle water path values were derived from the brightness temperatures reported by the microwave radiometer following the technique of Cadeddu et al. (2020 AMT). The cloud water path was then scaled to the KAZR-reported radar reflectivity to calculate profiles of liquid water content (LWC). Following the analysis from Ghate et al. (2023 JGR), cloud droplet effective radius was calculated using the number concentration value of 100 cm-3. The cloud properties, along with the thermodynamic properties, served as an input to the Rapid Radiative Transfer Model (RRTM) to yield profiles of radiative fluxes at a 1-minute temporal and 50-m vertical resolution. The fluxes were then averaged to hourly temporal resolution for analysis. In Mitra et al. (2025 JClim), the calculated profiles were compared against those derived from the satellite measurements (SYN1deg). Flux profiles from the SYN1deg and the thermodynamic and cloud properties used for deriving them are also reported here. Both all-sky and clear-sky radiative flux profiles were calculated. Due to the large data volume, the surface and top-of-atmosphere (TOA) radiative fluxes for the six-year period, and the hourly profiles of the radiative fluxes for January 2018, are submitted here. Full profiles of radiative fluxes calculated from the thermodynamic and cloud properties measured at the ENA site at 1-minute temporal and 50-m vertical resolution for a six-year period are available from the authors. Six files here correspond to the following data: 1_ENARAD_CERES_with_cld_amount_timeseries.nc: Time-series of hourly values of RRTM-simulated values of upwelling and downwelling fluxes at the surface and TOA, observed boundary-layer cloud fractions, and upwelling and downwelling fluxes from the SYN1deg from July 2015 to January 2022. 2_CERES_2018_at_CERES_levels.nc: SYN1deg radiative fluxes at six levels for the year 2018. 3_ENARad_2018_at_CERES_levels.nc: RRTM calculated fluxes at the SYN1deg vertical levels for the year 2018. 4_ENARad_rrtminputs_hourly_201801.nc: Thermodynamic and cloud properties used as an input to the RRTM for January 2018. 5_CERES_inputs_hourly_201801.nc: Thermodynamic and cloud properties utilized by SYN1deg algorithm for January 2018. 6_ENARAD_hourly_201801.nc: Full profiles of hourly averaged radiative fluxes from the RRTM simulations for January 2018.

Atmosphere↗

Discretization Writeup for Grey Flux-Limited Radiation Diffusion

This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗

Mapping foliar photosynthetic capacity in sub-tropical and tropical forests with UAS-based imaging spectroscopy: Scaling from leaf to canopy

Accurate understanding of the variability in foliar physiological traits across landscapes is critical to improve parameterization and evaluation of terrestrial biosphere models (TBMs) that seek to represent the response of terrestrial ecosystems to a changing climate. Numerous studies suggest imaging spectroscopy can characterize foliar biochemical and morphological traits at the canopy scale, but there is only limited evidence for retrieving canopy photosynthetic capacity (e.g., maximum carboxylation rate, V c,max and maximum electron transport rate, J max ). Moreover, the effect of canopy structure within forest communities on scaling up spectra-trait relationships from leaf to canopy level is not well known. To advance the spectra-trait approach and enable the estimation of key traits using remote sensing, we collected imaging spectroscopy data from an Unoccupied Aerial System (UAS) platform over two forest sites in China (a subtropical forest in Mt. Dinghu and a tropical rainforest in Xishuangbanna). At these sites, we also collected ground measurements of leaf spectra and traits, including biochemical (leaf nitrogen, phosphorus, chlorophyll, and water content), morphological (leaf mass per area, LMA) and physiological (V c,max25 and J max25 ) traits (n=135 tree-crowns from 42 species across two sites). Using a partial least-squares regression (PLSR) approach, we built and tested spectra-trait models with repeated cross-validation. The spectral models developed with leaf spectra were directly transferred to canopy spectra to evaluate the effect of canopy structure. Here we further applied canopy spectral models to map these traits at individual tree-crown scale. The results demonstrate that (1) UAS-based canopy spectra can be used to estimate V c,max (R 2 =0.55, nRMSE=11.79%), Jmax (R 2 =0.54, nRMSE=12.34%), and five additional foliar traits (R 2 =0.38-0.60, nRMSE=10.11-13.56%) at the tree-crown scale with demonstrated generalizability across two sites; (2) canopy structure strongly affects the spectratrait relationships from leaf to canopy level, but the effects vary considerably across foliar traits and cannot be well captured by the 4SAIL canopy radiative transfer model. UAS-based imaging spectroscopy maps large variability in all foliar traits (including physiological traits) with spatially explicit information, reproducing the field-observed inter- and intra-specific variations. These results demonstrate the capability of using UAS-based imaging spectroscopy for characterizing the variability of foliar physiological traits at individual tree-crown scale over forest landscapes and highlight the similar generalizability but different biophysical mechanisms underlying spectra-trait relationships at leaf and canopy levels.

54 ENVIRONMENTAL SCIENCES↗

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of the finite slab geometry in the NORA code

In this work we develop the complete linearization formalism for the radiation transfer equation in slab geometry. This development was motivated by an observed convergence instability arising from the NORA code’s inner boundary while modeling a suite of simple test problems in slab geometry presented at the Line Transfer Workshop held in LLNL in May of 2019. The convergence issue was seen to emerge early on in Newton-Raphson iterations. These test problems required parameters in the former NORA inner boundary to achieve the asymptotic limit of zero incoming intensity. In this work we derive the limiting case directly to gain insight into the former algorithm’s pathology.

97 MATHEMATICS AND COMPUTING↗