Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLOSURES”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Mapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing

A novel mapping closure approximation (MCA) technique is developed to construct a model for the conditional dissipation rate (CDR) of a scalar in homogeneous turbulence. It is shown that the CDR model from amplitude mapping closure is incorrect in asymptotic behavior for unsymmetric binary mixing. The correct asymptotic behavior can be described by the CDR model formulated by the MCA technique. The MCA approach is outlined for constructing successive approximation to probability density function (PDF) and conditional moment.

He, Gouwei↗
NASA NTRS2000

Clear-Sky Closure Studies of Lower Tropospheric Aerosol and Water Vapor During ACE-2 Using Airborne Sunphotometer, Airborne In-Situ, Space-Borne, and Ground-Based Measurements

We report on clear-sky column closure experiments (CLEARCOLUMN) performed in the Canary Islands during the second Aerosol Characterization Experiment (ACE-2) in June/July 1997. We present CLEARCOLUMN results obtained by combining airborne sunphotometer and in-situ (optical particle counter, nephelometer, and absorption photometer) measurements taken aboard the Pelican aircraft, space-borne NOAA/AVHRR data and ground-based lidar and sunphotometer measurements. During both days discussed here, vertical profiles flown in cloud-free air masses revealed 3 distinctly different layers: a marine boundary layer (MBL) with varying pollution levels, an elevated dust layer, and a very clean layer between the MBL and the dust layer. A key result of this study is the achievement of closure between extinction or layer aerosol optical depth (AOD) computed from continuous in-situ aerosol size-distributions and composition and those measured with the airborne sunphotometer. In the dust, the agreement in layer AOD (lambda = 380-1060 nm) is 3-8%. In the MBL there is a tendency for the in-situ results to be slightly lower than the sunphotometer measurements (10-17% at lambda = 525 nm), but these differences are within the combined error bars of the measurements and computations.

Schmid, Beat↗
NASA NTRS1999

In Situ Aerosol-Size Distributions and Clear-Column Radiative Closure During ACE-2

As part of the second Aerosol Characterization Experiment (ACE-2) during June and July of 1997, aerosol-size distributions were measured on board the CIRPAS Pelican aircraft through the use of a Differential Mobility Analyzer (DMA) and 2 Optical Particle Counters (OPCs). During the campaign, the boundary-layer aerosol typically possessed characteristics representative of a background marine aerosol or a continentally influenced aerosol, while the free-tropospheric aerosol was characterized by the presence or absence of a Saharan dust layer. A range of radiative closure comparisons were made using the data obtained during vertical profiles flown on 4 missions. Of particular interest here are the comparisons made between the optical properties as determined through the use of measured aerosol-size distributions and those measured directly by an airborne 14-wavelength sunphotometer and 3 nephelometers. Variations in the relative humidity associated with each of the direct measurements required consideration of the hygroscopic properties of the aerosol for size-distribution-based calculations. Simultaneous comparison with such a wide range of directly-measured optical parameters not only offers evidence of the validity of the physicochemical description of the aerosol when closure is achieved, but also provides insight into potential sources of error when some or all of the comparisons result in disagreement. Agreement between the derived and directly-measured optical properties varied for different measurements and for different cases. Averaged over the 4 case studies, the derived extinction coefficient at 525 nm exceeded that measured by the sunphotometer by 2.5% in the clean boundary layer, but underestimated measurements by 13% during pollution events. For measurements within the free troposphere, the mean derived extinction coefficient was 3.3% and 17% less than that measured by the sunphotometer during dusty and non-dusty conditions, respectively. Likewise, averaged discrepancies between the derived and measured scattering coefficient were -9.6%, +4.7%, +17%, and -41% for measurements within the clean boundary layer, polluted boundary layer, free troposphere with a dust layer, and free troposphere without a dust layer, respectively. Each of these quantities, as well as the majority of the more than 100 individual comparisons from which they were averaged, were within estimated uncertainties.

Collins, R.↗
NASA NTRS2006

Second Moment Closure Near the Two-component Limit

The purpose of this paper is to explore some wider implications of the two-component limit for both single point turbulence models and spectral closure theories. Although the two-component limit arises most naturally in inhomogeneous problems like wall-bounded turbulence, the analysis will be restricted to homogeneous turbulence. But since homogeneous turbulence is the crucial case for realizability, the conclusions will nevertheless be applicable to modeling. Th essential point of our argument is that whereas the evolution of the stochastic velocity field is Markovian because it is governed by the Navier-Stokes equations, the exact stress evolution equation is not Markovian because it is unclosed. This property of moment evolution has been stressed by Kraichnan (1959). We will show that modeling stress evolution at the two-component limit with a closure that is Markovian in the stresses alone leads to basic inconsistencies in single-point modeling and, perhaps surprisingly, in spectral modes as well.

Rubinstein, Robert↗
NASA NTRS2008

ST5 Observations of the Imbalance of Region 1 and 2 Field-Aligned Currents and its Implication to Ionospheric Closure Currents

A major unsolved question in the physics of ionosphere-magnetosphere coupling is how field-aligned currents (FACs) close. In order to maintain the divergence free condition, overall downward FACs (carried mainly by upward electrons) must eventually balance the overall upward FACs associated with the precipitating electrons through ionospheric Pedersen currents. Although much of the current closure may take place via local Pedersen currents flowing between Region 1 (R1) and Region 2 (R2) FACs, there is a generally an imbalance, i.e., more currents in R1 than in R2, in total currents between them. The net currents may be closed within R1 via cross-polar cap Pedersen currents. In this study, we use the magnetic field observations from Space Technology 5 mission to quantify the imbalance of R1 and R2 currents. We will determine the net R1-R2 currents under various solar wind conditions and discuss the implication of such imbalance to the ionospheric closure currents.

Le, G.↗
NASA NTRS2003

Toward Aerosol/Cloud Condensation Nuclei (CCN) Closure during CRYSTAL-FACE

During July 2002, measurements of cloud condensation nuclei were made in the vicinity of southwest Florida as part of the Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) field campaign. These observations, at supersaturations of 0.2 and 0.85%, are presented here. The performance of each of the two CCN counters was validated through laboratory calibration and an in situ intercomparison. The measurements indicate that the aerosol sampled during the campaign was predominantly marine in character: the median concentrations were 233 cm-3 (at S = 0.2%) and 371 cm(sup -3) (at S = 0.85%). Three flights during the experiment differed from this general trend; the aerosol sampled during the two flights on 18 July was more continental in character, and the observations on 28 July indicate high spatial variability and periods of very high aerosol concentrations. This study also includes a simplified aerosol/CCN closure analysis. Aerosol size distributions were measured simultaneously with the CCN observations, and these data are used to predict a CCN concentration using Kohler theory. For the purpose of this analysis, an idealized composition of pure ammonium sulfate was assumed. The analysis indicates that in this case, there was good general agreement between the predicted and observed CCN concentrations: at S = 0.2%, N(sub predicted)/N(sub observed)= 1.047 (R(sup 2)= 0.911)); at S = 0.85%, N(sub predicted)/N(sub observed)=1.201 (R(sup 2)= 0.835)). The impacts of the compositional assumption and of including in-cloud data in the analysis are addressed. The effect of removing the data from the 28 July flight is also examined; doing so improves the result of the closure analysis at S = 0.85%. When omitting that atypical flight, N(sub predicted)/N(sub observed) = 1.085 (R(sup 2) = 0.770) at S = 0.85%.

VanReken, Timothy M.↗
NASA NTRS2010

Toward Realistic Simulation of low-Level Clouds Using a Multiscale Modeling Framework With a Third-Order Turbulence Closure in its Cloud-Resolving Model Component

This study presents preliminary results from a multiscale modeling framework (MMF) with an advanced third-order turbulence closure in its cloud-resolving model (CRM) component. In the original MMF, the Community Atmosphere Model (CAM3.5) is used as the host general circulation model (GCM), and the System for Atmospheric Modeling with a first-order turbulence closure is used as the CRM for representing cloud processes in each grid box of the GCM. The results of annual and seasonal means and diurnal variability are compared between the modified and original MMFs and the CAM3.5. The global distributions of low-level cloud amounts and precipitation and the amounts of low-level clouds in the subtropics and middle-level clouds in mid-latitude storm track regions in the modified MMF show substantial improvement relative to the original MMF when both are compared to observations. Some improvements can also be seen in the diurnal variability of precipitation.

Xu, Kuan-Man↗
NASA NTRS2011

A Rocket-Base Study of Auroral Electrodynamics Within the Current Closure Ionosphere

The Auroral Current and Electrodynamics Structure (ACES) mission consisted of two sounding rockets launched nearly simultaneously from Poker Flat Research Range, AK on January 29, 2009 into a dynamic multiple-arc aurora. The ACES rocket mission, in conjunction with the PFISR Radar, was designed to observe the three-dimensional current system of a stable auroral arc system. ACES utilized two well instrumented payloads flown along very similar magnetic field footprints, at various altitudes with small temporal separation between both payloads. ACES High, the higher altitude payload (apogee 360 km), took in-situ measurements of the plasma parameters above the current closure region to provide the input signature into the lower ionosphere. ACES Low, the low-altitude payload (apogee 130 km), took similar observations within the current closure region, where cross-field currents can flow. We present results comparing observations of the electric fields, magnetic fields, electron flux, and the electron temperature at similar magnetic footpoints between both payloads. We further present data from all-sky imagers and PFISR detailing the evolution of the auroral event as the payloads traversed regions connected by similar magnetic footpoints. Current measurements derived from the magnetometers on both payloads are further compared. We examine data from both PFISR and observations on the high-altitude payload which we interpreted as a signature of electron acceleration by means of Alfv n waves. We further examine all measurements to understand ionospheric conductivity and how energy is being deposited into the ionosphere through Joule heating. Data from ACES is compared against models of Joule heating to make inferences regarding the effect of collisions at various altitudes.

Kaeppler, Stephen R.↗
NASA NTRSDate not supplied

Concept for a Far-infrared Outgoing Radiation Closure Experiment – Antarctica (FORCE-A)

The next decade promises to be an incredibly exciting time in climate science. There are two new space flight missions, PREFIRE and FORUM, that will open the far-infrared spectrum to direct, accurate observations for the first time. PREFIRE is planned to operate between 2022 and 2024 and FORUM will launch in late 2025 or early 2026. The TICFIRE instrument is also a candidate for the NASA A-CCP mission to be launched in the 2028 timeframe. A key focus of these missions and instruments is improved understanding of polar climates. In support of these missions we present a concept for a radiative closure experiment to be conducted in Antarctica during the PREFIRE mission lifetime and then again during the operational FORUM and TICFIRE/A-CCP missions. The main component of the campaigns would be a long-duration balloon flight launched from McMurdo Station with the potential of 1-2 months aloft. Candidate balloon flight instrumentation includes a far-IR Fourier transform spectrometer and far-IR radiometers. Ground based instrumentation includes zenith viewing infrared and far-infrared spectrometers, lidars, and microwave radiometers. The objective of the FORCE-A campaign is to demonstrate radiative closure in the infrared with the multiple campaign instruments combined with the numerous relevant satellite instruments that pass overhead every 30 minutes (AIRS, CrIS, IASI, MODIS, VIIRS, CERES, BBR, Libera). The campaign will serve to advance radiation sciences as well as to provide the means for validation of the new far-infrared observations.

Martin G Mlynczak↗
NASA NTRS2019

Closure Relations of Gamma-Ray Bursts in High Energy Emission

The synchrotron external shock model predicts the evolution of the spectral (β) and temporal (α) indices during the gamma-ray burst (GRB) afterglow for different environmental density profiles, electron spectral indices, electron cooling regimes, and regions of the spectrum. We study the relationship between α and β, the so-called “closure relations” with GRBs detected by Fermi Large Area Telescope (Fermi-LAT) from 2008 August to 2018 August. The spectral and temporal indices for the >100 MeV emission from the Fermi-LAT as determined in the Second Fermi-LAT Gamma-Ray Burst Catalog (2FLGC) are used in this work. We select GRBs whose spectral and temporal indices are well constrained (58 long-duration GRBs and 1 short-duration GRBs) and classify each GRB into the best-matched relation. As a result, we found that a number of GRBs require a very small fraction of the total energy density contained in the magnetic field (òB <~10−7). The estimated mean and standard deviation of electron spectral index p are 2.40 and 0.44, respectively. The GRBs satisfying a closure relation of the slow cooling tend to have a softer p value compared to those of the fast cooling. Moreover, the Kolmogorov–Smirnov test of the two p distributions from the fast and slow coolings rejects a hypothesis that the two distributions are drawn from the single reference distribution with a significance of 3.2σ. Lastly, the uniform density medium is preferred over the medium that decreases like the inverse of distance squared for long-duration GRBs.

Gamma-ray bursts↗
NASA NTRS2017

Global Energetics of Solar Flares. V. Energy Closure in Flares and Coronal Mass Ejections

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles (E(sub nt)), the energy of direct heating (E(sub dir)), and the energy in coronal mass ejections (E(sub CME)), which are the primary energy dissipation processes in a flare, is found to have a ratio of (E(sub nt) + E(sub dir) + E(sub CME))/E(sub mag) = 0.87±0.18, compared with the dissipated magnetic free energy E(sub mag), which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: 0.51±0.17 in nonthermal energy of ≥ 6 keV electrons, 0.17± 0.17 in nonthermal ≥ 1 MeV ions, 0.07 ± 0.14 in CMEs, and 0.07 ± 0.17 in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction ≈ 0.03 of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at e(sub c) ≈ 6 keV, based on the mean differential emission measure (DEM) peak temperature of T(sub e) = 8.6 MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar activity↗
NASA NTRS2017

Global Energetics of Solar Flares. V. Energy Closure in Flares and Coronal Mass Ejections

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles (E(sub nt)), the energy of direct heating (E(sub dir)), and the energy in coronal mass ejections (E(sub CME)), which are the primary energy dissipation processes in a flare, is found to have a ratio of (E(sub nt)+E(sub dir)+E(sub CME))/E(sub mag) = 0.87±0.18, compared with the dissipated magnetic free energy E(sub mag), which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: 0.51±0.17 in nonthermal energy of ≥ 6 keV electrons, 0.17± 0.17 in nonthermal ≥ 1 MeV ions, 0.07 ± 0.14 in CMEs, and 0.07 ± 0.17 in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction ≈ 0.03 of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at ec ≈ 6 keV, based on the mean differential emission measure (DEM) peak temperature of T(sub e) = 8.6 MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar activity↗
NASA NTRS2017

Global Energetics of Solar Flares and CMEs: V. Energy Closure

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles (E(sub nt)), the energy of direct heating (E(sub dir)), and the energy in coronal mass ejections (E(sub CME)), which are the primary energy dissipation processes in a flare, is found to have a ratio of (E(sub nt)+E(sub dir)+E(sub CME))/E(sub mag) = 0.87±0.18, compared with the dissipated magnetic free energy E(sub mag), which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: 0.51±0.17 in nonthermal energy of ≥ 6 keV electrons, 0.17± 0.17 in nonthermal ≥ 1 MeV ions, 0.07 ± 0.14 in CMEs, and 0.07 ± 0.17 in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction ≈ 0.03 of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at e(sub c) ≈ 6 keV, based on the mean differential emission measure (DEM) peak temperature of T(sub e) = 8.6 MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar activity↗
NASA NTRS2024

Assessment of Dust Size Retrievals Based on AERONET: A Case Study of Radiative Closure From Visible-Near-Infrared to Thermal Infrared

Super-coarse dust particles (diameters >10 μm) are evidenced to be more abundant in the atmosphere than model estimates and contribute significantly to the dust climate impacts. Since super-coarse dust accounts for less dust extinction in the visible-to-near-infrared (VIS-NIR) than in the thermal infrared (TIR) spectral regime, they are suspected to be underestimated by remote sensing instruments operates only in VIS-NIR, including Aerosol Robotic Networks (AERONET), a widely used data set for dust model validation. In this study, we perform a radiative closure assessment using the AERONET-retrieved size distribution in comparison with the collocated Atmospheric Infrared Sounder (AIRS) TIR observations with comprehensive uncertainty analysis. The consistently warm bias in the comparisons suggests a potential underestimation of super-coarse dust in the AERONET retrievals due to the limited VIS-NIR sensitivity. An extra super-coarse mode included in the AERONET-retrieved size distribution helps improve the TIR closure without deteriorating the retrieval accuracy in the VIS-NIR.

thermal infrared (TIR)↗
DOE OSTI2012

Closure End States for Facilities, Waste Sites, and Subsurface Contamination

The United States (U.S.) Department of Energy (DOE) manages the largest groundwater and soil cleanup effort in the world. DOE’s Office of Environmental Management (EM) has made significant progress in its restoration efforts at sites such as Fernald and Rocky Flats. However, remaining sites, such as Savannah River Site, Oak Ridge Site, Hanford Site, Los Alamos, Paducah Gaseous Diffusion Plant, Portsmouth Gaseous Diffusion Plant, and West Valley Demonstration Project possess the most complex challenges ever encountered by the technical community and represent a challenge that will face DOE for the next decade. Closure of the remaining 18 sites in the DOE EM Program requires remediation of 75 million cubic yards of contaminated soil and 1.7 trillion gallons of contaminated groundwater, deactivation & decommissioning (D&D) of over 3000 contaminated facilities and thousands of miles of contaminated piping, removal and disposition of millions of cubic yards of legacy materials, treatment of millions of gallons of high level tank waste and disposition of hundreds of contaminated tanks. The financial obligation required to remediate this volume of contaminated environment is estimated to cost more than 7% of the to-go life-cycle cost. Critical in meeting this goal within the current life-cycle cost projections is defining technically achievable end states that formally acknowledge that remedial goals will not be achieved for a long time and that residual contamination will be managed in the interim in ways that are protective of human health and environment. Formally acknowledging the long timeframe needed for remediation can be a basis for establishing common expectations for remedy performance, thereby minimizing the risk of re-evaluating the selected remedy at a later time. Once the expectations for long-term management are in place, remedial efforts can be directed towards near-term objectives (e.g., reducing the risk of exposure to residual contamination) instead of focusing on long-term cleanup requirements. An acknowledgement of the long timeframe for complete restoration and the need for long-term management can also help a site transition from the process of pilot testing different remedial strategies to selecting a final remedy and establishing a long-term management and monitoring approach. This approach has led to cost savings and the more efficient use of resources across the Department of Defense complex and at numerous industrial sites across the U.S. Defensible end states provide numerous benefits for the DOE environmental remediation programs including cost-effective, sustainable long-term monitoring strategies, remediation and site transition decision support, and long-term management of closure sites.

environmental remediation↗
DOE OSTI2015

The Cementitious Barriers Partnership Experimental Programs and Software Advancing DOE’s Waste Disposal/Tank Closure Efforts – 15436

The U.S. Department of Energy Environmental Management (DOE-EM) Office of Tank Waste Management-sponsored Cementitious Barriers Partnership (CBP) is chartered with providing the technical basis for implementing cement-based waste forms and radioactive waste containment structures for long-term disposal. DOE needs in this area include the following to support progress in final treatment and disposal of legacy waste and closure of High-Level Waste (HLW) tanks in the DOE complex: long-term performance predictions, flow sheet development and flow sheet enhancements, and conceptual designs for new disposal facilities. The DOE-EM Cementitious Barriers Partnership is producing software and experimental programs resulting in new methods and data needed for end-users involved with environmental cleanup and waste disposal. Both the modeling tools and the experimental data have already benefited the DOE sites in the areas of performance assessments by increasing confidence backed up with modeling support, leaching methods, and transport properties developed for actual DOE materials. In 2014, the CBP Partnership released the CBP Software Toolbox –“Version 2.0” which provides concrete degradation models for 1) sulfate attack, 2) carbonation, and 3) chloride initiated rebar corrosion, and includes constituent leaching. These models are applicable and can be used by both DOE and the Nuclear Regulatory Commission (NRC) for service life and long-term performance evaluations and predictions of nuclear and radioactive waste containment structures across the DOE complex, including future SRS Saltstone and HLW tank performance assessments and special analyses, Hanford site HLW tank closure projects and other projects in which cementitious barriers are required, the Advanced Simulation Capability for Environmental Management (ASCEM) project which requires source terms from cementitious containment structures as input to their flow simulations, regulatory reviews of DOE performance assessments, and Nuclear Regulatory Commission reviews of commercial nuclear power plant (NPP) structures which are part of the overall US Energy Security program to extend the service life of NPPs. In addition, the CBP experimental programs have had a significant impact on the DOE complex by providing specific data unique to DOE sodium salt wastes at Hanford and SRS which are not readily available in the literature. Two recent experimental programs on cementitious phase characterization and on technetium (Tc) mobility have provided significant conclusions as follows: recent mineralogy characterization discussed in this paper illustrates that sodium salt waste form matrices are somewhat similar to but not the same as those found in blended cement matrices which to date have been used in long-term thermodynamic modeling and contaminant sequestration as a first approximation. Utilizing the CBP generated data in long-term performance predictions provides for a more defensible technical basis in performance evaluations. In addition, recent experimental studies related to technetium mobility indicate that conventional leaching protocols may not be conservative for direct disposal of Tc-containing waste forms in vadose zone environments. These results have the potential to influence the current Hanford supplemental waste treatment flow sheet and disposal conceptual design.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Moment Closures on Two-Dimensional Cartesian Grids

Some moment methods for kinetic equations are complicated and take time to develop. Over the course of a couple years, this software was developed to test different closures on standard test problems in the literature. With this software, researchers in the field of moment closures will be able to rapidly test new methods.

Garrett, Charles K.↗
DOE OSTI2017

Kinetic fluid moments closure for a magnetized plasma with collisions

A novel method aimed at a kinetic moments closure for a magnetized plasma with collisions is proposed. The velocity distribution function for each species is expanded in 8 Gaussian Radial Basis Functions (GRBFs) which are essentially shifted Maxwellians at eight representative 3D-velocity points of drift. The vector of 8 fluid moments (for particle density, 3 particle fluxes, total energy density, and 3 energy fluxes) has an 8x8 analytic linear matrix relation to the vector of 8 GRBF density weighs in 3D-real space. The 8 fluid moments with sources for each species are advanced in time while the 8 GRBF weighs are determined from the 8x8 inverse matrix. The two closure moments (for the stress tensor and the energy weighted stress tensor) are linearly determined from the GRBF weighs. Most importantly the velocity moments of the nonlinear Coulomb Fokker-Planck collision operator [Rosenbluth et al, Phys. Rev 107, 1957] are evaluated from the GRBF weights. Generalization from 8 to 12, 16, 20 .., in an energy weighted moment hierarchy is straightforward. The electric field follows from a generalized vorticity (quasi-neutral charge conservation) equation. A strong drift ordering approximation can be applied to eliminate any spuriously unstable high frequency cyclotron motions. A novel weak drift ordering two-time step scheme avoids the vorticity equation by following ion cyclotron motion in time to get the electric field with ion gyroaveraging. Inclusion of low-beta magnetic perturbations is straightforward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗