Search NASA⌕ Search

SEARCH · Search NASA

Results for “polarization curves”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

(Invited) Modeling Electrolyzers: Exploring the Applied Voltage Breakdown

To commercialize successfully polymer-electrolyte electrolyzers and optimize their performance, one requires a detailed understanding of the underlying physics and phenomena. Mathematical modeling is ideally suited to explore such intricacies. In this Tutorial, the modeling equations and approaches towards both proton- and hydroxide-conducting polymer electrolyzers will be detailed. This includes the introduction of the applied voltage breakdown that separates the overall polarization curve into its constitutive parts such that the limiting mechanisms can be ascertained. Throughout, different case studies will be explored including impacts of bubble generation, liquid versus vapor feed, and pH changes due to carbonate and hydroxide solution feeds. A key focus of the Tutorial will be in exploring and describing the various transport and competing phenomena within the electrolyzer cell.

42 ENGINEERING↗
DOE OSTI2020

Mathematical Modeling of Hydroxide-Exchange-Membrane Water Electrolyzer

Water electrolyzers can transform intermittent renewable energy like solar energy and wind energy into the chemical energy of hydrogen with zero greenhouse-gas emissions. The hydroxide-exchange membrane electrolyzer (HEME) combines the capability to produce pressurized hydrogen with the advantage of being able to use low or non-platinum group metal (PGM) electrocatalysts in the alkaline environment.1 Hydroxide salts, for example, KOH, are added to the HEME water feed on both anode and cathode to improve its performance. However, the specific mechanism of performance improvement still needs to be further understood. In addition, at high current densities, bubble evolution can result in mass-transport limitations, a less well studied phenomena. Mathematical modeling is ideal to explore these issues as it is cost and time efficient and can deconvolute the physics, processes, and observed phenomena and study the applied-voltage breakdown. In this work, we extend our previously developed 1D two-phase continuum model2 to study the varies processes in the HEME and provide insights on performance optimizations. First, the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) kinetics at different hydroxide concentrations have been studied by rotating disk electrodes (RDE) and implemented in the model. Then, the model is calibrated and validated against experimental HEME polarization curves for different KOH concentrations as a liquid electrolyte. The model clearly shows a performance increase with increasing KOH concentrations, which is consistent with the experimental results. The reduced ohmic resistance and increased electrochemical active surface area (ECSA) are the two main reasons for performance increase. The large amount of hydroxide in the liquid electrolyte not only helps to distribute the reactant hydroxide throughout the catalyst layer (CL), which reduces ohmic loss, but also enables reaction at the interface between the liquid electrolyte and electrocatalyst, which increases the ECSA. Applied-voltage breakdown demonstrates that the electrolyzer performance is dominated by anode kinetics and ohmic loss. A comparison with the DI water feed shows a more uniform current distribution in the anode CL when KOH is added, which indicates a higher utilization of the CL. Second, we present modeling on the effects of bubble coverage. As gas evolves, part of the ECSA is minimized due to bubble coverage. To account for this effect, an empirical relationship between the fractional bubble coverage and the current density is implemented in the model.3 The model shows this bubble coverage effect is more pronounced at large current densities with DI water feed. Acknowledgements This work was funded under the HydroGEN Consortium by the Energy Efficiency and Renewable Energy, Hydrogen and Fuel Cell Technologies Office, of the U. S. Department of Energy under contract number DE-AC02-05CH11231. References R. Abbasi, B. P. Setzler, S. Lin, J. Wang, Y. Zhao, H. Xu, B. Pivovar, B. Tian, X. Chen, G. Wu and Y. Yan, 31, 1805876 (2019). L. N. Stanislaw, M. R. Gerhardt and A. Z. Weber, ECS Transactions, 92, 767 (2019). H. Vogt and R. J. Balzer, Electrochimica Acta, 50, 2073 (2005).

Liu, Jiangjin↗
DOE OSTI2024

Reduction of Methane Leaks through Corrosion Mitigation Pre-treatments for Pipelines with Field Applied Coatings

Corrosion of buried, coated steel pipelines transporting natural gas is a significant source of methane emissions, from pipeline venting required for maintenance and repairs and from pipeline leaks and incidents. Corrosion of steel under field applied coatings is an important safety concern for the pipeline industry. This project investigates the application of a field applied alloy over girth welds to mitigate external corrosion of buried coated steel pipelines. Various metallic coating options were considered, which were required to meet several criteria: (1) it must resist corrosion under open-circuit or mild cathodic protection conditions, (2) it must protect the substrate steel, and (3) it must not negatively affect the adhesion of the field coating. Finite element models and lab testing were performed of alloy coating compositions to identify promising alloy types underneath disbonded coatings. Polarization curves of coating alloys were generated to provide the boundary conditions for the COMSOL model to compute potential and current distributions around coated areas. Sacrificial and corrosion-resistant metal alloy coatings were evaluated and optimized using corrosion modeling and laboratory electrochemical testing, where aluminum alloy 5356 (5% Mg) and steel alloy B9 (9% Cr) were selected. Corrosion test coupons were designed and fabricated using thermal spray aluminum 5356 and welded B9 steel overlays on API 5L grade X42 line pipe steel. The corrosion test coupons, with simulated pipe coating damage, were tested in a laboratory soil box and a field pipeline site in Texas for 3-months. Corrosion test coupons were then tested for 6-months at field pipeline sites in Texas and Tennessee to quantify corrosion rates and performance of the aluminum and steel alloys under polyethylene tape and 2-part epoxy coatings, various coating holidays, and with and without cathodic protection.

03 NATURAL GAS↗
DOE OSTI2026

Full-polarization millimeter wavelength variability of Sagittarius A * during the 2018 EHT campaign

Context. Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims. We investigated the variability and polarization properties of Sgr A* using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results. We find low variability in total intensity (σ/μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions. Our results offer new constraints on the physical mechanisms driving variability in Sgr A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

Galaxy: center↗
DOE OSTI2023

Effective Landau-type model of a Hf x Zr 1 - x O 2 -graphene nanostructure

Here, to describe the charge-polarization coupling in the nanostructure formed by a thin Hf x Zr 1-x O 2 film with a single-layer graphene as a top electrode, we develop the “effective” Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf 1-x Zr x O 2 films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the Hf 1-x Zr x O 2 film polar properties on the film thickness, grain size and surface energy. The proposed Landau free energy with five “effective” expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf 1-x Zr x O 2 film influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 – 25) nm thick Hf 1-x Zr x O 2 films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 ≤ x ≤ 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs, because the working voltages applied to the Hf 1-x Zr x O 2 film (which acts as a gate) can be relatively low (less than 2 V). These low voltages are sufficient to induce the pronounced hysteresis of ferroelectric polarization in the Hf 1-x Zr x O 2 gate, which, due to the strong electric coupling, induces the hysteresis of the graphene charge.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Opportunities and challenges in process modeling and simulation of electrochemical systems

Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.

09 BIOMASS FUELS↗
DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗
DOE OSTI2025

Durable, Ultra-Thin Diaphragms for Liquid Alkaline Water Electrolysis

W. L. Gore & Associates assembled a world class team (De Nora Tech, LLC; New York University; and thyssenkrupp nucera) to develop durable, ultra-low resistance diaphragms that dramatically improve system efficiency, dynamic response, and operational flexibility. To enable the use of ultra-thin diaphragms, an optimized electrode package was chosen as part of an electrode down select study. Ultra-thin diaphragms based on microporous polymer scaffolds imbibed with zirconia nanoparticles were developed and tested. Performance and durability on par with incumbent materials was demonstrated. A transparent electrolyzer to directly observe in-situ diaphragm-bubble interactions was designed and built, and simultaneous optical and electrochemical performance tests were performed. A Pressurized Alkaline Water Electrolysis (PAWE) single cell test station was designed to support advanced R&D in high-pressure electrolysis, particularly for evaluating novel electrode and separator configurations under elevated temperature and pressure operating conditions. Due to early project termination, the original scope could not be fully achieved.

08 HYDROGEN↗
DOE OSTI2026

Hydrogen Infrastructure Modeling and Optimization

HyET, established in 2008, is a small business developing electro-chemical hydrogen compressors. HyET’s goals are for a low-cost, efficient, high-pressure compressor ready for market and installation into hydrogen infrastructure stations servicing mobile fuel cell applications. NLR has been working on hydrogen infrastructure research for over a decade and includes electrolyzer production via electrolysis, storage, compression, dispensing, safety, and performance validation. Accelerating technologies into the marketplace is one goal of NLR’s hydrogen research. Therefore, HyET and NLR have teamed up to provide one year of support for HyET’s modeling and small-scale experiments optimization efforts for HyET’s next-generation electro-chemical compressor.

08 HYDROGEN↗
DOE OSTI2021

Quasi-periodic Oscillations in the TESS Light Curve of TX Col, a Diskless Intermediate Polar on the Precipice of Forming an Accretion Disk

One of the fundamental properties of an intermediate polar is the dynamical nature of the accretion flow as it encounters the white dwarf’s (WD’s) magnetosphere. Many works have presumed a dichotomy between disk-fed accretion, in which the WD accretes from a Keplerian disk, and stream-fed accretion, in which the matter stream from the donor star directly impacts the WD’s magnetosphere without forming a disk. However, there is also a third, poorly understood regime in which the accretion flow consists of a torus of diamagnetic blobs that encircles the WD. This mode of accretion is expected to exist at mass-transfer rates below those observed during disk-fed accretion, but above those observed during pure stream-fed accretion. We invoke the diamagnetic-blob regime to explain the exceptional Transiting Exoplanet Survey Satellite light curve of the intermediate polar TX Col, which transitioned into and out of states of enhanced accretion during Cycles 1 and 3. Power-spectral analysis reveals that the accretion was principally stream fed. However, when the mass-transfer rate spiked, large-amplitude quasi-periodic oscillations (QPOs) abruptly appeared and dominated the light curve for weeks. The QPOs have two striking properties: they appear in a stream-fed geometry at elevated accretion rates, and they occur preferentially within a well-defined range of frequencies (∼10–25 cycle day{sup −1}). We propose that during episodes of enhanced accretion, a torus of diamagnetic blobs forms near the binary’s circularization radius and that the QPOs are beats between the white dwarf’s spin frequency and unstable blob orbits within the WD’s magnetosphere. We discuss how such a torus could be a critical step in producing an accretion disk in a formerly diskless system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Chiral Recognition: A Spin-Driven Process in Chiral Oligothiophene. A Chiral-Induced Spin Selectivity (CISS) Effect Manifestation

In this paper it is experimentally demonstrated that the electron-spin/molecular-handedness interaction plays a fundamental role in the chiral recognition process. This conclusion is inferred comparing current versus potential (I-V) curves recorded using chiral electrode surfaces, which are obtained via chemisorption of an enantiopure thiophene derivative: 3,3'-bibenzothiophene core functionalized with 2,2'-bithiophene wings (BT 2 T 4 ). The chiral recognition capability of these chiral-electrodes is probed via cyclic voltammetry measurements, where, Ag nanoparticles (AgNPs) capped with enantiopure (BT 2 T 4 ) (BT 2 T 4 )@AgNP) are used as the chiral redox probe. Then, the interface handedness is explored by recording spin-polarized (I-V) curves in spin-dependent electrochemistry (SDE) and magnetic-conductive atomic force microscopy (mc-AFM) experiments. The quality of the interfaces is thoroughly cross-checked using X-ray photoemission spectroscopy, Raman, electrodesorption measurements, which further substantiate the metal(electrode)-sulfur(thiophene) central role in the chemisorption process. Spin-polarization values of about 15% and 30% are obtained in the case of SDE and mc-AFM experiments, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Chemical Compensation Challenges in Processing Antiferroelectric PbZrO 3 Thin Films

The electric field-induced antipolar-to-polar phase transition in antiferroelectric materials is accompanied by large changes in the structural and functional response, making them attractive for many applications ranging from pulsed energy capacitors to high-strain-high (blocking-)force actuators, solid-state heating and cooling systems, and optoelectronic devices. PbZrO 3 , as an end member of the PZT, (x)PbZr 1–x Ti x O 3 , solid solution, has been the subject of many studies. However, processing of perovskite PbZrO 3 is extremely challenging due to phase instabilities induced by Pb loss at the temperatures necessary for perovskite crystallization. Here, we discuss the challenges associated with Pb loss, as well as its compensation through bulk Pb overstoichiometry and interfacial PbO additions, in chemical solution processed, highly oriented PbZrO 3 thin films. For both 042 o - and 001 o -oriented (o-orthorhombic distortion) PbZrO 3 thin films, the crystallization interfaces are the most important contributors to Pb loss and off-stoichiometric outcomes, and no single approach was sufficient to address these challenges. Pb-rich and Pb-deficient nonperovskite phases including ZrO x were only observable through microscopic characterization, and X-ray diffraction alone could not rule out the presence of such secondary phases. “Ideal” macroscopic antiferroelectric polarization-electric field double hysteresis curves were observed despite (at times) the large presence of secondary phases. However, reduced saturation polarization and increased phase transition electric fields were observed concomitantly with the presence of secondary phase(s) and tentatively assigned to the voltage drop across the ZrO x nanocrystals. Based on these results, it is imperative that the presence of secondary phases always be addressed (beyond X-ray diffraction spectra) in order to correctly evaluate the properties of antiferroelectric films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

In Situ UV–vis Quantification of Chlorobasicity in Binary Molten Chlorides

Optical basicity measurements provide a basis for determining the relative chlorobasicity of a given molten salt composition. This study uses Bi3+ as a spectroscopic probe to investigate optical basicity across three binary molten chloride systems (LiCl-KCl, MgCl2–KCl, and ZnCl2–KCl) through in situ UV–vis spectroscopy measurements of the 1S0→3P1 transition. LiCl-KCl systems exhibit the highest optical basicity with minimal compositional dependence, while both divalent systems show stronger, yet nonlinear, dependence that is reminiscent of acid–base titration curves. The influence of cation electronegativity and polarizing power (Li+ < Mg2+ < Zn2+) directly reflects the basicity curves, with ZnCl2–KCl displaying the steepest gradients and lowest basicity values. Nonlinear dependencies in divalent systems indicate formation of chloride-bridged network structures with coordination states dictated by the amount of KCl in a salt composition, which was further investigated and validated by molecular dynamics simulations of polymeric structures formed in ZnCl2 and MgCl2-based salts using machine learning force fields. The temperature dependence of chlorobasicity in terms of optical basicity was investigated, and a loose, positive correlation between the two was revealed. These findings establish Bi3+ as a sensitive probe for understanding electronic environments in molten chloride systems and provide insights for predicting the chemical behavior of molten salts for use in high-temperature environments.

Buttice, Kailee [Department of Nuclear Engineering↗
DOE OSTI2023

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Post maximum light and late time optical imaging polarimetry of type I superluminous supernova 2020znr

ABSTRACT Optical imaging polarimetry was conducted on the hydrogen poor superluminous supernova 2020znr during three phases after maximum light (≈ +34 d, +288 d, and +289 d). After instrumental and interstellar polarization correction, all measurements are consistent with null-polarization detection. Modelling the light curve with a magnetar spin-down model shows that SN2020znr has similar magnetar and ejecta parameters to other SLSNe. A comparison of the best-fitting values discussed in the literature on SN 2017egm and SN 2015bn, two hydrogen poor SLSNe showing an increase of polarization after maximum light, suggests that SN 2020znr has higher mass ejecta that may prevent access to the geometry of the inner ejecta with optical polarimetry. The combined information provided by spectroscopy and light-curve analysis of type I SLSNe may be an interesting avenue to categorize the polarization properties of this class of transients. This approach would require to expand the sample of SLSNe polarimetry data currently available with early and late time epochs new measurements.

2015bn↗
DOE OSTI2020

Experimental validation of solid oxide fuel cell polarization modeling: An LSM-YSZ/YSZ/Ni-YSZ case study

Polarization modeling based on the measured I-V curves is a useful tool for analyzing SOFC performance and diagnosing performance limitations. Using the polarization model, the overpotential of a cell is usually separated into different polarization contributions, including area specific ohmic polarization, activation polarization, and anodic and cathodic concentration polarizations. In this study, we aim to experimentally validate the accuracy of the polarization model. As a case study, anode-supported cells with LSM-YSZ/YSZ/Ni-YSZ configuration are used. By curve-fitting the I-V curves of a cell tested under different cathode oxygen partial pressures into the polarization model, the contributions of different polarization losses are quantified. To validate the polarization model, independent experiments/studies including electrochemical impedance spectroscopy, cathode kinetics study, ex-situ diffusivity measurement, and tortuosity measurement using 3D reconstructed anode are conducted. The results from polarization modeling and from validation studies are in good agreement, thereby validating the polarization model for SOFC performance analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Ferroelectric Switching in Hybrid Molecular-Beam-Epitaxy-Grown BaTiO 3 Films

Molecular beam epitaxy (MBE) is a promising synthesis technique for both heterostructure growth and epitaxial integration of ferroelectric BaTiO 3 . However, direct measurement of the remnant polarization (P r ) has not been previously reported in MBE-grown BaTiO 3 films. We report the in situ growth of an all-epitaxial SrRuO 3 /BaTiO 3 /SrRuO 3 heterostructure on Nb-doped SrTiO 3 (001) substrates by hybrid MBE using metal–organic precursors. This capacitor structure consisting of 16 nm SrRuO 3 /40 nm BaTiO 3 /16 nm SrRuO 3 shows hysteretic polarization–electric field (P–E) curves with P r ∼ 15 μC cm –2 at frequencies ranging from 500 Hz to 20 kHz, after isolating the intrinsic ferroelectric response from non-ferroelectric contributions using the Positive-Up-Negative-Down (PUND) method. We hypothesize that the asymmetry in switching behavior and current leakage has origins in structural defects. Furthermore, this work opens the door to defect-engineered ferroelectric BaTiO 3 -based heterostructures grown by hybrid MBE for future electronic, photonic and spintronic applications.

SrRuO3 electrode↗
DOE OSTI2021

A Lattice Boltzmann Scheme for Polarized Radiative Transfer in Planetary Atmospheres

Quantitative calculation of reflected polarized light from planets is of great significance to interpret the polarimetric observations of planetary atmospheres. In this work, we propose a lattice Boltzmann (LB) scheme for polarized radiative transfer in planetary atmospheres. The LB scheme obtains the Stokes vector through the vector LB equation performing simple collision and streaming processes. The vector radiative transfer equation is rigorously derived from the vector LB equation via the Maxwell iteration technique. Polarized radiative transfer of a single point on the planetary surface with Rayleigh and Mie scattering atmospheres are first accurately solved by our LB scheme. Afterward, we systematically investigate the disk-integrated polarization of finite and semi-infinite, conservative and nonconservative Rayleigh scattering planetary atmospheres. The disk-integrated phase curves of the Stokes parameters and degree of polarization and spherical and geometric albedos are produced, which are in good agreement with the benchmark results for different cases. The numerical results indicate that our LB scheme is efficient and accurate for polarized radiative transfer in planetary atmospheres. Our LB scheme is expected to provide a competitive numerical tool to interpret the polarimetric observations of planetary atmospheres.

47 OTHER INSTRUMENTATION↗