Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling coefficient”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

WE-B-207-01: CT Lung Cancer Screening and the Medical Physicist: Background, Findings and Participant Dosimetry Summary of the National Lung Screening Trial (NLST)

The US National Lung Screening Trial (NLST) was a multi-center randomized, controlled trial comparing a low-dose CT (LDCT) to posterior-anterior (PA) chest x-ray (CXR) in screening older, current and former heavy smokers for early detection of lung cancer. Recruitment was launched in September 2002 and ended in April 2004 when 53,454 participants had been randomized at 33 screening sites in equal proportions. Funded by the National Cancer Institute this trial demonstrated that LDCT screening reduced lung cancer mortality. The US Preventive Services Task Force (USPSTF) cited NLST findings and conclusions in its deliberations and analysis of lung cancer screening. Under the 2010 Patient Protection and Affordable Care Act, the USPSTF favorable recommendation regarding lung cancer CT screening assisted in obtaining third-party payers coverage for screening. The objective of this session is to provide an introduction to the NLST and the trial findings, in addition to a comprehensive review of the dosimetry investigations and assessments completed using individual NLST participant CT and CXR examinations. Session presentations will review and discuss the findings of two independent assessments, a CXR assessment and the findings of a CT investigation calculating individual organ dosimetry values. The CXR assessment reviewed a total of 73,733 chest x-ray exams that were performed on 92 chest imaging systems of which 66,157 participant examinations were used. The CT organ dosimetry investigation collected scan parameters from 23,773 CT examinations; a subset of the 75,133 CT examinations performed using 97 multi-detector CT scanners. Organ dose conversion coefficients were calculated using a Monte Carlo code. An experimentally-validated CT scanner simulation was coupled with 193 adult hybrid computational phantoms representing the height and weight of the current U.S. population. The dose to selected organs was calculated using the organ dose library and the abstracted scan parameters. This session will review the results and summarize the individualized doses to major organs and the mean effective dose and CTDIvol estimate for 66,157 PA chest and 23,773 CT examinations respectively, using size-dependent computational phantoms coupled with Monte Carlo calculations. Learning Objectives: Review and summarize relevant NLST findings and conclusions. Understand the scope and scale of the NLST specific to participant dosimetry. Provide a comprehensive review of NLST participant dosimetry assessments. Summarize the results of an investigation providing individualized organ dose estimates for NLST participant cohorts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-B-207-02: CT Lung Cancer Screening and the Medical Physicist: A Dosimetry Summary of CT Participants in the National Lung Cancer Screening Trial (NLST)

The US National Lung Screening Trial (NLST) was a multi-center randomized, controlled trial comparing a low-dose CT (LDCT) to posterior-anterior (PA) chest x-ray (CXR) in screening older, current and former heavy smokers for early detection of lung cancer. Recruitment was launched in September 2002 and ended in April 2004 when 53,454 participants had been randomized at 33 screening sites in equal proportions. Funded by the National Cancer Institute this trial demonstrated that LDCT screening reduced lung cancer mortality. The US Preventive Services Task Force (USPSTF) cited NLST findings and conclusions in its deliberations and analysis of lung cancer screening. Under the 2010 Patient Protection and Affordable Care Act, the USPSTF favorable recommendation regarding lung cancer CT screening assisted in obtaining third-party payers coverage for screening. The objective of this session is to provide an introduction to the NLST and the trial findings, in addition to a comprehensive review of the dosimetry investigations and assessments completed using individual NLST participant CT and CXR examinations. Session presentations will review and discuss the findings of two independent assessments, a CXR assessment and the findings of a CT investigation calculating individual organ dosimetry values. The CXR assessment reviewed a total of 73,733 chest x-ray exams that were performed on 92 chest imaging systems of which 66,157 participant examinations were used. The CT organ dosimetry investigation collected scan parameters from 23,773 CT examinations; a subset of the 75,133 CT examinations performed using 97 multi-detector CT scanners. Organ dose conversion coefficients were calculated using a Monte Carlo code. An experimentally-validated CT scanner simulation was coupled with 193 adult hybrid computational phantoms representing the height and weight of the current U.S. population. The dose to selected organs was calculated using the organ dose library and the abstracted scan parameters. This session will review the results and summarize the individualized doses to major organs and the mean effective dose and CTDIvol estimate for 66,157 PA chest and 23,773 CT examinations respectively, using size-dependent computational phantoms coupled with Monte Carlo calculations. Learning Objectives: Review and summarize relevant NLST findings and conclusions. Understand the scope and scale of the NLST specific to participant dosimetry. Provide a comprehensive review of NLST participant dosimetry assessments. Summarize the results of an investigation providing individualized organ dose estimates for NLST participant cohorts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-CD-206-09: Learning-Based MRI-CT Prostate Registration Using Spare Patch-Deformation Dictionary

Purpose: To enable MRI-guided prostate radiotherapy, MRI-CT deformable registration is required to map the MRI-defined tumor and key organ contours onto the CT images. Due to the intrinsic differences in grey-level intensity characteristics between MRI and CT images, the integration of MRI into CT-based radiotherapy is very challenging. We are developing a learning-based registration approach to address this technical challenge. Methods: We propose to estimate the deformation between MRI and CT images in a patch-wise fashion by using the sparse representation technique. Specifically, we assume that two image patches should follow the same deformation if their patch-wise appearance patterns are similar. We first extract a set of key points in the new CT image. Then, for each key point, we adaptively construct a coupled dictionary from the training MRI-CT images, where each coupled element includes both appearance and deformation of the same image patch. After calculating the sparse coefficients in representing the patch appearance of each key point based on the constructed dictionary, we can predict the deformation for this point by applying the same sparse coefficients to the respective deformations in the dictionary. Results: This registration technique was validated with 10 prostate-cancer patients’ data and its performance was compared with the commonly used free-form-deformation-based registration. Several landmarks in both images were identified to evaluate the accuracy of our approach. Overall, the averaged target registration error of the intensity-based registration and the proposed method was 3.8±0.4 mm and 1.9±0.3 mm, respectively. Conclusion: We have developed a novel prostate MR-CT registration approach based on patch-deformation dictionary, demonstrated its clinical feasibility, and validated its accuracy. This technique will either reduce or compensate for the effect of patient-specific treatment variation measured during the course of radiotherapy, is therefore well-suited for a number of MRI-guided adaptive radiotherapy, and potentially enhance prostate radiotherapy treatment outcome.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Electrochemical Reactivity of Hydrogen Storage Materials: Exploring Borohydride and Hydrazinium Salt Mixtures

Electrochemical characterization of hydrogen storage materials was conducted in a non-aqueous environment to investigate the direct electrochemical release and consumption of hydrogen and the potential for regeneration. We first address the challenge of minimal solubility of the synthetic precursors, sodium borohydride (NaBH 4 ) and hydrazinium bromide (N 2 H 5 Br), in both organic and inorganic solvents. We next determine and calibrate a reference electrode formulation compatible with our non-aqueous media and analytes that demonstrates a stable reference potential. We employ cyclic voltammetry (CV) to characterize the precursors and mixtures thereof. Each CV peak is assigned to a corresponding electrochemical reaction. Using the rate-dependent CV method and Randles–Ševčík equation, we calculate the diffusion coefficient of each chemical (NaBH 4 and N 2 H 5 Br). Analysis of the CVs, coupled with 11B NMR analysis, reveals a room temperature chemical transformation of NaBH 4 and N 2 H 5 Br mixtures into hydrazine borane (N 2 H 4 BH 3 ). These results are particularly significant, considering the limited information available on the electrochemical characterization of metal borohydride and hydrazinium salt in non-aqueous media. This work establishes a foundation for adapting a non-aqueous electrochemical system to further study the borohydride family of chemistries and to design and develop electrochemical devices for direct electrical and chemical energy interconversion with hydrogen storage materials.

08 HYDROGEN↗
DOE OSTI2021

The Cosmological Bootstrap: Spinning correlators from symmetries and factorization

We extend the cosmological bootstrap to correlators involving massless spinning particles, focusing on spin-1 and spin-2. In de Sitter space, these correlators are constrained both by symmetries and by locality. In particular, the de Sitter isometries become conformal symmetries on the future boundary of the spacetime, which are reflected in a set of Ward identities that the boundary correlators must satisfy. We solve these Ward identities by acting with weight-shifting operators on scalar seed solutions. Using this weight-shifting approach, we derive three- and four-point correlators of massless spin-1 and spin-2 fields with conformally coupled scalars. Four-point functions arising from tree-level exchange are singular in particular kinematic configurations, and the coefficients of these singularities satisfy certain factorization properties. We show that in many cases these factorization limits fix the structure of the correlators uniquely, without having to solve the conformal Ward identities. The additional constraint of locality for massless spinning particles manifests itself as current conservation on the boundary. We find that the four-point functions only satisfy current conservation if the s, t, and u-channels are related to each other, leading to nontrivial constraints on the couplings between the conserved currents and other operators in the theory. For spin-1 currents this implies charge conservation, while for spin-2 currents we recover the equivalence principle from a purely boundary perspective. For multiple spin-1 fields, we recover the structure of Yang--Mills theory. Finally, we apply our methods to slow-roll inflation and derive a few phenomenologically relevant scalar-tensor three-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2024

Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence

Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Quantum calculation of inelastic CO collisions with H. II. Pure rotational quenching of high rotational levels

Carbon monoxide is a simple molecule present in many astrophysical environments, and collisional excitation rate coefficients due to the dominant collision partners are necessary to accurately predict spectral line intensities and extract astrophysical parameters. We report new quantum scattering calculations for rotational deexcitation transitions of CO induced by H using the three-dimensional potential energy surface (PES) of Song et al. State-to-state cross sections for collision energies from 10{sup −5} to 15,000 cm{sup −1} and rate coefficients for temperatures ranging from 1 to 3000 K are obtained for CO (v = 0, j) deexcitation from j=1−45 to all lower j′ levels, where j is the rotational quantum number. Close-coupling and coupled-states calculations were performed in full-dimension for j=1−5, 10, 15, 20, 25, 30, 35, 40, and 45 while scaling approaches were used to estimate rate coefficients for all other intermediate rotational states. The current rate coefficients are compared with previous scattering results using earlier PESs. Astrophysical applications of the current results are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

HYDRODYNAMIC TESTING OF ZIRCALOY TO STEEL HUB COUPLINGS FOR USE IN ADVANCED TEST REACTOR I-LOOP TUBE

The ATR I-positions' isolation from the flux trap region provides an opportunity for a test facility that could create conditions representative of commercial light water reactors. Such a facility requires piping couplings and joints that withstand extreme thermal and pressure cycling. Also, transitioning from steel to zircaloy piping inside the reactor vessel allows the neutronic properties of zircaloy to be leveraged while using cost-effective stainless steel outside of the core. Grayloc-style hub couplings will act as transition points between stainless steel piping and zircaloy piping inside the reactor vessel. These couplings have not been sufficiently tested for leakage where there is a transition between materials. This transition is a concern because the thermal expansion coefficients of the materials differ. Seal integrity between the two components may be greatly affected by thermal expansion at the connection. Thus, the effect of thermal cycling must be assessed. Thermal and pressure cycling and hydrodynamic testing will be conducted to verify seal integrity of the steel to zircaloy hub coupling when subjected to ATR conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, we theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O 3 . However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. Herein, to examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients ( trans -aconitic acid (C 6 H 6 O 6 ), maleic acid (C 4 H 4 O 4 ), and sodium nitrite (NaNO 2 )). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O 3 partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O 3 and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O 3 , compared to gas and bulk phases is rapidly maintained and remains constant even as O 3 is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O 3 loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O 3 arrival in the bulk by diffusion. As a result, bulk-phase O 3 becomes depleted from its Henry's law solubility. This bulk-phase O 3 depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k rxn ≈ 10 3 -10 4 M -1 s -1 ), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k rxn ≈ 10 5 M -1 s -1 ) where surface reactions strongly impact the observed decay kinetics.

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

Steady-State Thermal-Hydraulic Analysis and Bowing Reactivity Evaluation Methods Based on Neutron and Gamma Transport Calculations

With the advancement of computer technology, the variational nodal transport code VARIANT is now routinely used for sodium-cooled fast reactor (SFR) design analysis. In particular, VARIANT transport calculations are usually performed for whole core neutronics analyses and for fuel cycle analyses with REBUS-3. However, the steady-state thermal hydraulics analysis with the SE2-ANL coed and the assembly bowing reactivity calculation with the structural analysis with the NUBOW-3D code are still performed with the DIF3D finite difference diffusion theory code. Furthermore, the DIF3D diffusion calculation for these applications is performed with only six triangular meshes per hexagonal assembly, and the SE2-ANL code determines the pin power distributions approximately by assuming that the power distribution within an assembly is separable in the radial and axial directions. In order to overcome the limitations and to improve the accuracy of the existing methods, it was proposed to develop the relevant computational methods based on the transport calculations with the VARIANT and PROTEUS-SN codes. The objectives of this study are to develop a new heating calculation method based on the coupled neutron and gamma transport calculations with VARIANT, to implement the new heating calculation method into the SE2-ANL steady-state thermal-hydraulics analysis code with improved numerical algorithms and additional capabilities for automatic flow allocation calculations, and to develop a method to calculate bowing reactivity coefficients based on the VARIANT core solutions and the PROTEUS-SN assembly solutions, which can be used in NUBOW-3D calculations. A new, coupled neutron and gamma heating calculation procedure has been developed based on VARIANT transport calculations. In the new heating calculation procedure, the neutron flux distribution is first determined by solving an eigenvalue transport problem using the VARIANT transport code. Using the calculated neutron flux, the intra-nodal gamma source distributions are calculated to be consistent with the neutron flux distribution with the GAMSOR code. For this, the GAMSOR code has been updated to generate the intra-nodal gamma source distributions in the form of the VARSRC dataset of VARIANT. With this gamma source distribution, the gamma flux distribution is determined by solving a fixed source gamma transport problem using VARIANT again. With the neutron and gamma flux distributions, the heat generation rates are evaluated at individual fuel pins and assembly duct walls, without assuming separable axial and radial profiles for the intra-assembly power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Tunable two-dimensional interfacial coupling in molecular heterostructures

Two-dimensional van der Waals heterostructures are of considerable interest for the next generation nanoelectronics because of their unique interlayer coupling and optoelectronic properties. Here, we report a modified Langmuir–Blodgett method to organize twodimensional molecular charge transfer crystals into arbitrarily and vertically stacked heterostructures, consisting of bis(ethylenedithio)tetrathiafulvalene (BEDT–TTF)/C 60 and poly (3-dodecylthiophene-2,5-diyl) (P3DDT)/C 60 nanosheets. A strong and anisotropic interfacial coupling between the charge transfer pairs is demonstrated. The van der Waals heterostructures exhibit pressure dependent sensitivity with a high piezoresistance coefficient of -4.4 × 10 -6 Pa -1 , and conductance and capacitance tunable by external stimuli (ferroelectric field and magnetic field). Density functional theory calculations confirm charge transfer between the n-orbitals of the S atoms in BEDT–TTF of the BEDT–TTF/C 60 layer and the π* orbitals of C atoms in C 60 of the P3DDT/C 60 layer contribute to the inter-complex CT. Thus, the two-dimensional molecular van der Waals heterostructures with tunable optical–electronic–magnetic coupling properties are promising for flexible electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Robustness of cluster synchronous patterns in small-world networks with inter-cluster co-competition balance

All edges in the classical Watts and Strogatz's small-world network model are unweighted and cooperative (positive). By introducing competitive (negative) inter-cluster edges and assigning edge weights to mimic more realistic networks, this paper develops a modified model which possesses co-competitive weighted couplings and cluster structures while maintaining the common small-world network properties of small average shortest path lengths and large clustering coefficients. Based on theoretical analysis, it is proved that the new model with inter-cluster co-competition balance has an important dynamical property of robust cluster synchronous pattern formation. More precisely, clusters will neither merge nor split regardless of adding or deleting nodes and edges, under the condition of inter-cluster co-competition balance. Numerical simulations demonstrate the robustness of the model against the increase of the coupling strength and several topological variations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters

We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A DFT analysis of the ground and charge-transfer excited states of Sc 3 N@I h –C 80 fullerene coupled with metal-free and zinc-phthalocyanine

Endohedral metallofullerenes and phthalocyanine derivatives are recognized as excellent active materials in organic photovoltaics (OPVs). The tri-metallic nitride endohedral C 80 fullerenes have greater absorption coefficients in the visible region and electron-accepting abilities similar to C 60 , which can allow for higher efficiencies in OPV devices. Here, we have investigated the ground and charge transfer excited states of two co-facial donor–acceptor (D–A) molecular conjugates formed by the non-covalent coupling of trimetallic nitride endohedral fullerene (Sc 3 N@I h –C 80 ) with metal-free (H 2 Pc) and zinc-phthalocyanine (ZnPc) chromophores using DFT calculations. The charge transfer (CT) excitation energies are calculated using the perturbative delta-SCF method that enforces orthogonality between the ground and excited states. The binding energies calculated using the PBE and DFT-D3 methods indicate that the dispersion effects play an important role in the stabilization of these complexes. The ground state dipole moment of the Sc 3 N@C 80 –H 2 Pc dyad is much larger than that of Sc 3 N@C 80 –ZnPc, but this is reversed in the excited state where the dipole moment of Sc 3 N@C 80 –ZnPc increases significantly. The lowest few excitation energies in the gas phase for the two complexes are very close, in the range of 1.51–2.66 eV for Sc 3 N@C 80 –ZnPc and 1.51–2.71 eV for the Sc 3 N@C 80 –H 2 Pc complex. However, the lower ionization potential and lower exciton binding energy make the Sc 3 N@C 80 –ZnPc dyad a better candidate for OPVs as compared to the Sc 3 N@C 80 –H 2 Pc dyad.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗