Search NASA⌕ Search

SEARCH · Search NASA

Results for “activity coefficient model”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Aeroservoelastic modeling and applications using minimum-state approximations of the unsteady aerodynamics

The theoretical basis of flexible-aircraft modeling techniques encompassing aerodynamic, control, and elastic-structure effects is investigated analytically, with a focus on methods which employ minimum-state approximations for the unsteady aerodynamics. Rational-function approximations to generalized aerodynamic forces are reviewed; constraints and lag-coefficient optimization are explained; the problem of physical weighting in the minimum-state equations of motion is examined; and results of typical analyses from the NASA Active Flexible Wing project (Perry et al., 1988) are presented in extensive tables and graphs and discussed in detail. The minimum-state approach is shown to produce accurate models at significantly reduced computation costs.

Tiffany, Sherwood H.↗
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermodynamic properties of the Nd-Bi system via emf measurements, $\mathrm{DFT}$ calculations, machine learning, and $\mathrm{CALPHAD}$ modeling

Thermodynamic properties of the Nd-Bi system were investigated using a combination of experimental measurements, first-principles calculations based on density functional theory (DFT), data mining and machine learning (DM + ML) predictions, and calculation of phase diagrams (CALPHAD) modeling. The electromotive force (emf) of Nd-Bi alloys in molten LiCl-KCl-NdCl 3 at 773–973 K was measured via coulometric titration of Nd into Bi for the determination of thermochemical properties such as activity coefficients and solubilities of Nd in Bi. A new peritectic reaction of [liquid + NdBi 2 = Nd 3 Bi 7 ] at 774 K was confirmed using differential scanning calorimetry, structural (X-ray diffraction), and microstructural (scanning electron microscopy) analyses. The unknown crystal structure of NdBi2 was suggested to be a mixture of the anti-La 2 Sb configuration and the La 2 Te-type configuration based on ML predictions for over 26,000 data-mined AB 2 -type configurations together with DFT-based verifications. Using the newly acquired experimental data and DFT-based calculations, the thermodynamic description of the Nd-Bi system was remodeled, and a more complete Nd-Bi phase diagram was calculated, including the Nd 3 Bi 7 compound, invariant transition reactions, and liquidus temperatures.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Solar Sail Torque Model Characterization for the Near Earth Asteroid Scout Mission

Near Earth Asteroid Scout (NEA Scout) was a mission to test solar sail propulsion for orbital transfer from cislunar space to flyby and image an asteroid. Had it succeeded, one of the mission goals was to characterize the solar torque on the sail to ensure successful attitude control for the orbit transfer and imaging the asteroid. The simulation used to develop the flight attitude control software uses the generalized model for solar sails, a tensor equation of the forces and torques on sails of arbitrary shape. Rios-Reyes and Scheeres developed a general process to update the torque tensor coefficients using estimates of sail torque over a range of directions to the sun. Their process was adapted and implemented for the specific case of NEA Scout using spacecraft telemetry collected during sail characterization maneuvers in combination with simulation models and parameters. The NEA Scout maneuvers were limited to the operating range of the mission and constraints of the control hardware and allowed safe testing of each attitude before proceeding to the next. The NEA Scout reaction wheel speeds are used to measure accumulated momentum, while the Active Mass Translator (AMT) position is used to subtract out the torque from the center of mass crossed with the sail force and isolate the torque from only the sail shape. The process was tested by running attitude control simulations of the characterization maneuvers, generating simulated telemetry, estimating the solar torques, then using a least squares estimating the solar torque coefficients using least-squares and then performing a least-squares fit to the solar torque tensor coefficients. These estimated coefficients were tested by evaluating the solar torques under the same conditions as the simulated telemetry and comparing to the true simulated torques. Solar force model updates can be performed separately by observing the effect of the sail on the trajectory, and the torque model can be refined using those solar force updates. This process met the needs of the NEA Scout mission and can be adapted to characterize the solar torque for other missions with different sails.

solar sail↗
NASA NTRSDate not supplied

Solar Sail Torque Model Characterization for the Near Earth Asteroid Scout Mission

Near Earth Asteroid Scout (NEA Scout) was a mission to test solar sail propulsion for orbital transfer from cislunar space to flyby and image an asteroid. One of the goals of the mission was to characterize the solar torque on the sail to ensure successful attitude control for the orbit transfer and imaging the asteroid. The simulation used to develop the flight attitude control software uses the generalized model for solar sails, a tensor equation of the forces and torques on sails of arbitrary shape. Rios-Reyes and Scheeres developed a general process to update the torque tensor coefficients using estimates of sail torque over a range of directions to the sun. Their process was adapted and implemented for the specific case of NEA Scout using spacecraft telemetry collected during sail characterization maneuvers in combination with simulation models and parameters. The NEA Scout maneuvers were limited to the operating range of the mission and constraints of the control hardware and allowed safe testing of each attitude before proceeding to the next. The NEA Scout reaction wheel speeds are used to measure accumulated momentum, while the Active Mass Translator (AMT) position is used to subtract out the torque from the center of mass crossed with the sail force and isolate the torque from only the sail shape. The process was tested by running attitude control simulations of the characterization maneuvers, generating simulated telemetry, estimating the solar torques, then using a least squares estimating the solar torque coefficients using least-squares and then performing a least-squares fit to the solar torque tensor coefficients. These estimated coefficients were tested by evaluating the solar torques under the same conditions as the simulated telemetry and comparing to the true simulated torques. Solar force model updates can be performed separately by observing the effect of the sail on the trajectory, and the torque model can be refined using those solar force updates. This process met the needs of the NEA Scout mission and can be adapted to characterize the solar torque for other missions with different sails.

solar sail↗
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients: Impact of Particle Deposition in ICRP 66 Human Respiratory Tract Model

Inhaled radioactive materials can pose a long-term health concern, as the material can be incorporated into the body’s metabolic pathways and remain in organs and tissues for extended durations. During the retention period, the radioactive material may localize in a source organ and irradiate adjacent target organs and tissues. Distribution of these materials changes over time, requiring biokinetic modeling to evaluate their movement through various tissues and organs. The evolving distribution depends on multiple inputs characterizing the inhaled material, such as particle size and size distribution, particle density, aspect ratio, specific radionuclide, the chemical form, and solubility. In addition, biological parameters such as breathing rate, breathing type (nasal or nasal/oral), respiratory system morphometry, tidal volume, functional residual capacity, and anatomical dead space all influence material transport. These aerosol properties and physiological characteristics of the respiratory tract jointly define a range of initial conditions that influence the time-dependent distribution of radioactive material. To evaluate both uncertainty in the initial conditions of inhalation exposure and the final output (committed effective dose) from biokinetic models, a Python-based software tool, Radiological Exposure Dose Calculator (REDCAL), was developed to propagate uncertainty within the human respiratory tract model. Focusing on deposition fraction uncertainty, the primary objective was to characterize the initial activity distribution across respiratory regions as a function of anticipated particle sizes and distributions. The impact of the deposition fraction uncertainty was propagated to committed effective dose coefficients for selected radionuclides in a companion publication. For each particle size, a lognormal distribution, characterized by its geometric mean as defined within ICRP Publication 66, serves as the basis for introducing uncertainty into the physical processes governing deposition in various lung regions. Finally, this study addresses the deposition process and examines how uncertainty in deposition mechanisms affects activity distribution in the airways, ultimately presenting the expected range and standard deviation of deposited activity as a function of particle size.

International Commission on Radiological Protectio↗
DOE OSTI2025

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High Efficiency 2-Bladed Utility Wind Turbine Enhanced by CoFlow Jet Active Flow Control

The aerodynamic performance and power output of a high efficiency Co-Flow Jet (CFJ) 2-bladed Wind Turbine is studied in this paper. CFJ is a Zero-Net-Mass-Flux (ZNMF) active flow control method that dramatically increases airfoil lift coefficient and suppresses flow separation at a low energy expenditure. The 3D Reynolds Averaged Navier-Stokes (RANS) equations with one-equation Spalart-Allmaras (SA) turbulence model are solved to simulate the 3D flows of the wind turbines. The CFJ-Wind Turbine in this paper is modified from the NREL 3-bladed 3.4 MW reference wind turbine with the same RPM, blade length, distribution of chord, twist, and bending. The CFJ 2-bladed wind turbine utilizes the CFJ-NACA 64 series airfoils at most of the blade span from 31.88\% to the tip. The CFJ injection and suction slots are implemented along the full blade span to enhance the effectiveness and power output. The study shows that the 2-bladed CFJ turbine significantly increases power output compared with the 3-bladed baseline turbine at all wind speeds by 15.8\% at the wind speed of 9.8 m/s and by 25.9\% at 4 m/s respectively. The efficiency enhancement is substantially more at the low wind speed than at high wind speed. This would be very beneficial to increase the annual power output or capacity factor as wind turbines work at the speeds lower than the rated speed for most of the time. The results show that the flow field around the blade surface presents a significant enhanced suction peak near the blade leading edge, which results in a very high tangential force and thus increased wind turbine net power coefficient. Consistent with the previous study for the small NREL Phase VI research wind turbine, this paper demonstrates that CFJ 2-bladed wind turbine is an effective technology to extract more power from the wind at all speeds. This study is just an initial effort to demonstrate the advantage of the CFJ blade with no design optimization. It is believed that there is more potential to further enhance the efficiency improvement.

active flow control↗
DOE OSTI2025

Effects of threshold resummation for large-x PDF in large momentum effective theory

Parton distribution functions (PDFs) at large x are challenging to extract from experimental data, yet they are essential for understanding hadron structure and searching for new physics beyond the Standard Model. Within the framework of the large momentum P z expansion of lattice quasi-PDFs, we investigate large x PDFs, where the matching coefficient is factorized into the hard kernel, related to the active quark momentum xP z , and the threshold soft function, associated with the spectator momentum (1 − x)P z . The renormalization group equation of the soft function enables the resummation of the threshold double logarithms α k ln 2k (1 − x), which is crucial for a reliable and controllable calculation of large x PDFs. Our analysis with pion valence PDFs indicates that perturbative matching breaks down when the spectator momentum (1 − x)P z approaches Λ QCD , but remains valid when both xP z and (1 − x)P z are much larger than Λ QCD . Additionally, we incorporate leading renormalon resummation within the threshold framework, demonstrating good perturbative convergence in the region where both spectator and active quark momenta are perturbative scales.

factorization↗
DOE OSTI2021

Raman Spectroscopy Coupled with Chemometric Analysis for Speciation and Quantitative Analysis of Aqueous Phosphoric Acid Systems

Complex chemical systems that exhibit varied and matrix-dependent speciation are notoriously difficult to monitor and characterize on-line and in real-time. Optical spectroscopy is an ideal tool for in situ characterization of chemical species that can enable quantification as well as species identification. Chemometric modeling, a multivariate method, has been successfully paired with optical spectroscopy to enable measurement of analyte concentrations even in complex solutions where univariate methods such as Beer’s law analysis fail. Here, Raman spectroscopy is used to quantify the concentration of phosphoric acid and its three deprotonated forms during a titration. In this system, univariate approaches would be difficult to apply due to multiple species being present simultaneously within the solution as pH is varied. Locally Weighted Regression (LWR) modeling was used to determine phosphate concentration from spectral signature. LWR results, in tandem with Multivariate Curve Resolution modeling, provide direct measurement of the concentration of each phosphate species using only the Raman signal. Furthermore, results are presented within the context of fundamental solution chemistry, including Pitzer equations to compensate for activity coefficients and non-idealities associated with high ionic strength systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1992

Modeling and measurement of microwave emission and backscattering from bare soil surfaces

A multifrequency ground-based radiometer-scatterometer system working at frequencies between 3.0 GHz and 11.0 GHz has been used to study the effect of soil moisture and roughness on microwave emission and backscattering. The freezing and thawing effect of the soil surface and the changes of the surface roughness due to rain and erosion are reported. To analyze the combined active and passive data, a scattering model based on physical optics approximation for the low frequency and geometrical optics approximation for high frequency has been developed. The model is used to calculate the bistatic scattering coefficients from the surface. By considering the conservation of energy, the result has been integrated over a hemisphere above the surface to calculate the emissivity. The backscattering and emission model has been coupled with the observed data in order to extract soil moisture and surface roughness.

Saatchi, S.↗
NASA NTRS1989

Atmospheric drag model calibrations for spacecraft lifetime prediction

Although solar activity prediction uncertainty normally dominates decay prediction error budget for near-Earth spacecraft, the effect of drag force modeling errors for given levels of solar activity needs to be considered. Two atmospheric density models, the modified Harris-Priester model and the Jacchia-Roberts model, to reproduce the decay histories of the Solar Mesosphere Explorer (SME) and Solar Maximum Mission (SMM) spacecraft in the 490- to 540-kilometer altitude range were analyzed. Historical solar activity data were used in the input to the density computations. For each spacecraft and atmospheric model, a drag scaling adjustment factor was determined for a high-solar-activity year, such that the observed annual decay in the mean semimajor axis was reproduced by an averaged variation-of-parameters (VOP) orbit propagation. The SME (SMM) calibration was performed using calendar year 1983 (1982). The resulting calibration factors differ by 20 to 40 percent from the predictions of the prelaunch ballistic coefficients. The orbit propagations for each spacecraft were extended to the middle of 1988 using the calibrated drag models. For the Jaccia-Roberts density model, the observed decay in the mean semimajor axis of SME (SMM) over the 4.5-year (5.5-year) predictive period was reproduced to within 1.5 (4.4) percent. The corresponding figure for the Harris-Priester model was 8.6 (20.6) percent. Detailed results and conclusions regarding the importance of accurate drag force modeling for lifetime predictions are presented.

Binebrink, A. L.↗
DOE OSTI2024

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1989

Selection of optimum median-filter-based ambiguity removal algorithm parameters for NSCAT

The NASA Scatterometer, NSCAT, is an active spaceborne radar designed to measure the normalized radar backscatter coefficient (sigma0) of the ocean surface. These measurements can, in turn, be used to infer the surface vector wind over the ocean using a geophysical model function. Several ambiguous wind vectors result because of the nature of the model function. A median-filter-based ambiguity removal algorithm will be used by the NSCAT ground data processor to select the best wind vector from the set of ambiguous wind vectors. This process is commonly known as dealiasing or ambiguity removal. The baseline NSCAT ambiguity removal algorithm and the method used to select the set of optimum parameter values are described. An extensive simulation of the NSCAT instrument and ground data processor provides a means of testing the resulting tuned algorithm. This simulation generates the ambiguous wind-field vectors expected from the instrument as it orbits over a set of realistic meoscale wind fields. The ambiguous wind field is then dealiased using the median-based ambiguity removal algorithm. Performance is measured by comparison of the unambiguous wind fields with the true wind fields. Results have shown that the median-filter-based ambiguity removal algorithm satisfies NSCAT mission requirements.

Shaffer, Scott↗
NASA NTRS2007

Effects on Spacecraft Radiometric Data at Superior Solar Conjunction

During 2006, three ESA interplanetary spacecraft, Rosetta, Mars Express (MEX) and Venus Express (VEX), passed through superior solar conjunction. For all three spacecraft, the noise in the post-fit range-rate residuals from the orbit determination was analysed. At small Sun-Earth-Probe (SEP) angles the level was almost two orders of magnitude higher than normal. The main objective was to characterize the Doppler (rangerate) noise as a function of SEP angle. At least then the range-rate data can be appropriately weighted within the orbit determination so that the solution uncertainties are realistic. For VEX, some intervals of particularly noisy Doppler data could be correlated with unusual solar activity. For Rosetta, the biases in the range data residuals were analysed with the aim of improving the model used for calibrating the signal delay due to the solar plasma. The model, which originally had fixed coefficients, was adjusted to achieve better fits to the data. Even the relatively small Doppler biases were well represented. Using the improved model, the electron density at 20 solar radii was compared with earlier results obtained by radio science studies using Voyager 2 and Ulysses radiometric data. There is some evidence for a dependency of the density on the phase within the 11 years solar cycle.

Morley, Trevor↗
DOE OSTI2021

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

A parameterization of eddy transfer coefficients for two-level seasonal statistical dynamical zonally averaged models

This paper examines a parameterization of a quasi-geostrophic eddy transport that takes into account the time variation of eddy transfer coefficients according to Green's (1970) theory. It was found that, in the original eddy transfer relationship of Green, connecting the integral of the northward eddy entropy flux through midlatitudes with the second power of the difference in 500-mb entropy across the region of baroclinic activity, a value of 4 for the exponent is obtained when the temperature gradients at 500 mb are used. When the gradients at 1000 mb are used, an exponent of 1.5 is obtained. The differences in the powers in the eddy transfer relation were explored in a two-level zonally averaged model. It was found that an appropriate choice of power may be of special importance if the model is devised to simulate the seasonal climate cycle or to test astronomical changes inducing different seasonalities.

Neeman, Binyamin U.↗