Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIZE DETERMINATION”

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 1,009 records · Page 56

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗

Small-Scale Spray Releases: Additional Aerosol Test Results

One of the events postulated in the hazard analysis at the Waste Treatment and Immobilization Plant (WTP) and other U.S. Department of Energy (DOE) nuclear facilities is a breach in process piping that produces aerosols with droplet sizes in the respirable range. The current approach for predicting the size and concentration of aerosols produced in a spray leak involves extrapolating from correlations reported in the literature. These correlations are based on results obtained from small engineered spray nozzles using pure liquids with Newtonian fluid behavior. The narrow ranges of physical properties on which the correlations are based do not cover the wide range of slurries and viscous materials that will be processed in the WTP and across processing facilities in the DOE complex. To expand the data set upon which the WTP accident and safety analyses were based, an aerosol spray leak testing program was conducted by Pacific Northwest National Laboratory (PNNL). PNNL’s test program addressed two key technical areas to improve the WTP methodology (Larson and Allen 2010). The first technical area was to quantify the role of slurry particles in small breaches where slurry particles may plug the hole and prevent high-pressure sprays. The results from an effort to address this first technical area can be found in Mahoney et al. (2012a). The second technical area was to determine aerosol droplet size distribution and total droplet volume from prototypic breaches and fluids, including sprays from larger breaches and sprays of slurries for which literature data are largely absent. To address the second technical area, the testing program collected aerosol generation data at two scales, commonly referred to as small-scale and large-scale. The small-scale testing and resultant data are described in Mahoney et al. (2012b) and the large-scale testing and resultant data are presented in Schonewill et al. (2012). In tests at both scales, simulants were used to mimic the relevant physical properties projected for actual WTP process streams.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Large-Scale Spray Releases: Additional Aerosol Test Results

One of the events postulated in the hazard analysis for the Waste Treatment and Immobilization Plant (WTP) and other U.S. Department of Energy (DOE) nuclear facilities is a breach in process piping that produces aerosols with droplet sizes in the respirable range. The current approach for predicting the size and concentration of aerosols produced in a spray leak event involves extrapolating from correlations reported in the literature. These correlations are based on results obtained from small engineered spray nozzles using pure liquids that behave as a Newtonian fluid. The narrow ranges of physical properties on which the correlations are based do not cover the wide range of slurries and viscous materials that will be processed in the WTP and in processing facilities across the DOE complex. To expand the data set upon which the WTP accident and safety analyses were based, an aerosol spray leak testing program was conducted by Pacific Northwest National Laboratory (PNNL). PNNL’s test program addressed two key technical areas to improve the WTP methodology (Larson and Allen 2010). The first technical area was to quantify the role of slurry particles in small breaches where slurry particles may plug the hole and prevent high-pressure sprays. The results from an effort to address this first technical area can be found in Mahoney et al. (2012a). The second technical area was to determine aerosol droplet size distribution and total droplet volume from prototypic breaches and fluids, including sprays from larger breaches and sprays of slurries for which literature data are mostly absent. To address the second technical area, the testing program collected aerosol generation data at two scales, commonly referred to as small-scale and large-scale testing. The small-scale testing and resultant data are described in Mahoney et al. (2012b), and the large-scale testing and resultant data are presented in Schonewill et al. (2012). In tests at both scales, simulants were used to mimic the relevant physical properties projected for actual WTP process streams.

42 ENGINEERING↗

Structural, magnetic and catalytic properties of cobalt chromite obtained through precursor method

Highlights: • CoCr{sub 2}O{sub 4} was synthesized through the tartarate and gluconate precursor routes. • Both routes led to the formation of the single-phase CoCr{sub 2}O{sub 4}. • The crystallite size was in the range of 14–21 nm. • CoCr{sub 2}O{sub 4} samples presented ferrimagnetic ordering below Currie temperature T{sub c} = 97 K. • CoCr{sub 2}O{sub 4} samples presented catalytic performance in the total oxidation of CH{sub 4}. - Abstract: Cobalt chromite (CoCr{sub 2}O{sub 4}) was synthesized through the precursor method. The precursors: (NH{sub 4}){sub 3}[CoCr{sub 2}(C{sub 4}O{sub 6}H{sub 4}){sub 4}(OH){sub 3}]·4H{sub 2}O, (NH{sub 4}){sub 3}[CoCr{sub 2}(C{sub 6}O{sub 7}H{sub 10}){sub 4}(C{sub 6}O{sub 7}H{sub 9})]·5H{sub 2}O were characterized by elemental chemical analysis, infrared (IR) and ultraviolet–visible (UV–vis) spectroscopy, and thermal analysis. The final oxides were characterized by X-ray diffraction (XRD), scanning and transmission electron microscopy (SEM/TEM), UV–vis, IR, Raman spectroscopy (RS), magnetic measurements, N{sub 2} adsorption–desorption analyses and X-ray photoelectron spectroscopy (XPS). XRD confirmed the cubic CoCr{sub 2}O{sub 4} phase only and determined average crystallite sizes between 14 and 21 nm. Electron microscopy revealed morphology corresponding to the complete crystallization into cubic CoCr{sub 2}O{sub 4}. All the samples presented ferrimagnetic ordering below the Currie temperature (T{sub c}), and a phase transition at T{sub s} ∼26 K attributed to the onset of long-range spiral magnetic order. The CoCr{sub 2}O{sub 4} nanoparticles generated through the gluconate route following calcination at 700 °C for 1 h were found to have the best catalytic activity in the total oxidation of methane.

36 MATERIALS SCIENCE↗

Evaluation of Hydrogen Production Feasibility for a Light Water Reactor in the Midwest

Increased electricity production from renewable energy resources coupled with low natural gas prices has caused existing light-water reactor (LWRs) to experience ever diminishing returns from the electricity market. Via partnership among Idaho National Laboratory (INL), The National Renewable Energy Laboratory (NREL), Argonne National Laboratory (ANL), Exelon, and Fuel Cell Energy a techno-economic analysis on the viability of retrofitting existing pressurized water reactors (PWRs) to produce hydrogen (H2) via high temperature steam electrolysis (HTSE) has been conducted. Such integration would allow nuclear facilities to expand into additional markets that may be more profitable in the long-term. To accommodate such an integration, a detailed analysis of HTSE process operation, requirements, and flexibility was conducted. Technical analysis includes proposed nuclear system control scheme modifications to allow dynamic operation of the HTSE via both thermal and electrical connection to the nuclear plant. High fidelity Modelica simulations showcase the viability of such control schemes. However, due to limited knowledge of solid oxide fuel cell (SOFC) stack degradation due to thermal gradients; thermal cycling of the HTSE was not included. Therefore, the control schemes proposed are only utilized to re-distribute steam at startup and only the portion of electricity utilized in the electrolyzers is cycled. From the detailed analysis of the nuclear integration and the HTSE process design, a comprehensive cost estimation was conducted in the APEA and H2A models to elucidate capital and operational costs associated with the production, compression, and distribution of hydrogen from a nuclear facility. Alongside this costing analysis, market analyses were conducted by NREL and ANL on the electric and hydrogen markets respectively in the PJM interconnect. Utilizing the electricity data market projections in the PJM interconnect from NREL and hydrogen demand/pricing projections from ANL a five variable sweep over component capacities, discount rates, and hydrogen pricing was completed using the stochastic framework RAVEN (Risk Analysis Virtual ENVironment) through its resource dispatch plugin HERON (Heuristic Energy Resource Optimization Network). Each combination of variables was evaluated stochastically four times over a seventeen year timespan from 2026-2042 (inclusive) to determine the most economically advantageous solution. Following the five variable sweep an optimization was conducted about the best sweep point to determine optimal component sizing and setpoints. Results suggest positive gain is achievable at all projected hydrogen market pricing levels and at all discount rates. However, exact component sizing and net returns vary based on these values and if incorrect sizing is selected; major net losses can occur. The optimal result occurred with set points as follows: high hydrogen prices, the largest possible HTSE unit in the sweep set at 7.47 kg/sec [645.4 tpd], a contractual hydrogen market agreement 7.29 kg/sec [629.8 tpd], and a hydrogen storage size 115,188 kg. The analysis suggested that with a discount rate of 8%, a ΔNPV = 1.2 billion over the seventeen year span can be achieved. The results illuminate that by operating in multiple markets the nuclear facility is capable of avoiding sale of electricity during times of low electricity market pricing, while maintaining the ability to capitalize on the high electricity market pricing. It should be noted that the analysis conducted in this report is a differential cash flow analysis and as such does not present profit levels. Instead it highlights the net benefit between building and competing in the hydrogen market utilizing nuclear facilities and conducting business as usual in the electricity markets. Additionally, results presented in this report exhibit conservatism due to four key assumptions: Given the limited knowledge on SOFC stack degradation due to thermal gradients the high temperature steam electrolysis plant is not allowed to thermally cycle. This limitation decreases electricity generation capacity in the nuclear plant, reducing capacity payments and the maximum output to the electrical grid. This analysis considers building a separate hydrogen pipeline for use in the nuclear facility. Research is currently being conducted to determine if existing natural gas pipelines can accommodate direct hydrogen injection. If hydrogen can be integrated into existing natural gas pipelines then a capital savings of ~$19,000,000 per kg/sec [1 kg/sec = 86.4 tpd] of installed HTSE capacity could be realized. Fast transients in the electric grid are typically served by flywheels and electric batteries. With electrolysis only operation the nuclear plant can flex its electricity just as quickly. However, for this analysis, the ancillary services market has been neglected. Should the nuclear plant have the ability to operate in the ancillary services market during periods when it is producing hydrogen, an additional substantial income stream that would increase overall system profitability could be considered. Subsidies for non-emitting generation technologies are not introduced in the cash flow analysis, but instead are completed assuming zero subsidies. Profitability is based solely on electricity pricing. Electric pricing is assumed to never go negative. Instead, it is assumed overall system curtailment by the Independent System Operator (ISO) will occur in such a manner that pricing remains at or above zero. This assumption inflates the “business as usual” scenario, thus decreasing the ΔNPV calculated. If negative pricing scenarios continue to occur as they do today, one would expect a further increase in ΔNPV. Overall, results advocate that, through market diversification, nuclear power plants have the potential to substantially increase current profit margins, increase market penetration, and ultimately solidify their place as a mainstay in energy production.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Additive Construction using Basalt Regolith Fines

Planetary surfaces are often covered in regolith (crushed rock), whose geologic origin is largely basalt. The lunar surface is made of small-particulate regolith and areas of boulders located in the vicinity of craters. Regolith composition also varies with location, reflecting the local bedrock geology and the nature and efficiency of the micrometeorite-impact processes. In the lowland mare areas (suitable for habitation), the regolith is composed of small granules (20 - 100 microns average size) of mare basalt and volcanic glass. Impacting micrometeorites may cause local melting, and the formation of larger glassy particles, and this regolith may contain 10-80% glass. Studies of lunar regolith are traditionally conducted with lunar regolith simulant (reconstructed soil with compositions patterned after the lunar samples returned by Apollo). The NASA Kennedy Space Center (KSC) Granular Mechanics & Regolith Operations (GMRO) lab has identified a low fidelity but economical geo-technical simulant designated as Black Point-1 (BP-1). It was found at the site of the Arizona Desert Research and Technology Studies (RATS) analog field test site at the Black Point lava flow in adjacent basalt quarry spoil mounds. This paper summarizes activities at KSC regarding the utilization of BP-1 basalt regolith and comparative work with lunar basalt simulant JSC-1A as a building material for robotic additive construction of large structures. In an effort to reduce the import or in-situ fabrication of binder additives, we focused this work on in-situ processing of regolith for construction in a single-step process after its excavation. High-temperature melting of regolith involves techniques used in glassmaking and casting (with melts of lower density and higher viscosity than those of metals), producing basaltic glass with high durability and low abrasive wear. Most Lunar simulants melt at temperatures above 1100 C, although melt processing of terrestrial regolith at 1500 C is not uncommon. These temperatures are achievable by laser heating or by using solar concentrators. Similar to volcanic magma, the cooling rate determines the crystallite size - slower cooling develops larger crystals, and rapid quenching can result in fully amorphous glass.

Additive↗

Preparation and analysis of particulate metal deposits

Small particles and clusters of palladium were grown by deposition from the vapor phase under ultrahigh vacuum conditions. Amorphous and crystalline support films of Al2O3 and ultrathin amorphous carbon films were used as substrate materials. The growth of the metal deposit was monitored in situ by scanning transmission diffraction of energy-filtered 100 kV electrons and high resolution transmission electron microscopy (TEM) analysis was performed in a separate instrument. It was established by in situ TEM, however, that the transfer of specimens in this case did not unduly affect the size and distribution of deposit particles. It was found that the cleanness, stoichiometry, crystallinity and structural perfection of the support surface play an essential role in determining the crystalline perfection and structure of the particles. The smallest palladium clusters reproducibly prepared contained not more than six atoms but size determinations below 1 nm average particle diameter are very problematic with conventional TEM. Palladium particles grown on carbon supports feature an impurity-stabilized mosaic structure.

Poppa, H.↗

Black carbon mixing state impacts on cloud microphysical properties: Effects of aerosol plume and environmental conditions

Black carbon (BC) is usually mixed with other aerosol species within individual aerosol particles. This mixture, along with the particles' size and morphology, determines the particles' optical and cloud condensation nuclei properties, and hence black carbon's climate impacts. In this study the particle-resolved aerosol model PartMC-MOSAIC was used to quantify the importance of black carbon mixing state for predicting cloud microphysical quantities. Based on a set of about 100 cloud parcel simulations a process level analysis framework was developed to attribute the response in cloud microphysical properties to changes in the underlying aerosol population ("plume effect") and the cloud parcel cooling rate ("parcel effect"). It shows that the response of cloud droplet number concentration to changes in BC emissions depends on the BC mixing state. When the aerosol population contains mainly aged BC particles an increase in BC emission results in increasing cloud droplet number concentrations ("additive effect"). In contrast, when the aerosol population contains mainly fresh BC particles they act as sinks for condensable gaseous species, resulting in a decrease in cloud droplet number concentration as BC emissions are increased ("competition effect"). Additionally, we quantified the error in cloud microphysical quantities when neglecting the information on BC mixing state, which is often done in aerosol models. Here, the errors ranged from -12% to +45% for the cloud droplet number fraction, from 0% to +1022% for the nucleation-scavenged black carbon (BC) mass fraction, from -12% to +4% for the effective radius, and from -30% to +60% for the relative dispersion.

54 ENVIRONMENTAL SCIENCES↗

Cover crops and poultry litter impact on soil structural stability in dryland soybean production in southeastern United States

Abstract This study explored the efficacy of soil aggregate indices in quantifying soil structural development, utilizing 5‐year field experiment data from the Southeastern United States. The experiment utilized a split‐plot design with cover crops (native vegetation as control, cereal rye (Secale cerealeL.), winter wheat (Triticum aestivum), hairy vetch (Vicia villosa), and mustard (Brassica rapa) plus cereal rye as the main factor and fertilizer source (no fertilizer as control, inorganic fertilizer with phosphorus, potassium, and elemental sulfur, and poultry litter) as the secondary factor. Aggregate size fractions were determined using the wet‐sieving method, and aggregate stability index (ASI), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (FD) were calculated to assess soil structural stability. Main effects results indicated that cereal rye (55.11%) and poultry litter (50.97%) exhibited the highest ASI values. The highest MWD, GMD, and FD were observed under mustard plus cereal rye (1.187 mm), cereal rye (0.462 mm), and hairy vetch (2.573), respectively. Principal component analysis revealed that cover crops significantly improved soil aggregate structure and stability, overcoming limitations of sole fertilization practices. Regression analysis suggested that ASI, MWD, and GWD positively correlated with soil organic carbon, whereas FD negatively correlated with MWD, GMD, and ASI. Principal component analysis exhibited that FD decreased with increasing soil organic carbon, ASI, MWD, and GMD, demonstrating that lower FD values indicate enhanced soil aggregation and structure. Assessed indices, FD included, effectively gauged soil structural stability. These metrics should be prioritized in managerial decisions to support soil productivity and health in agricultural systems.

Agriculture↗

Rotating Disk Slurry Au Electrodeposition at Unsupported Carbon Vulcan XC-72 and Ce 3+ Impregnation for Ethanol Oxidation in Alkaline Media

A robust electrodeposition method consisting of the rotating disk slurry electrode (RoDSE) technique to obtain Au nanoparticles highly dispersed on a conductive carbonaceous support, i.e., Vulcan XC-72R, for ethanol electrooxidation reaction in alkaline media was developed. Ceria was used as a cocatalyst using a Ce(III)-EDTA impregnation method in order to enhance the catalytic activity and improve the catalyst’s overall stability. Furthermore, the RoDSE method used to obtain highly dispersed Au nanoparticles does not require the use of a reducing agent or stabilizing agent, and the noble-metal loading was controlled by the addition and tuning of the metal precursor concentration. Inductively coupled plasma and thermogravimetric analysis indicated that the Au loading in the catalyst was 9 %. We determined the particle size and characteristic Au fcc crystal facets by X-ray diffraction. The morphology of the catalyst was also investigated using electron microscopy techniques. In addition, X-ray absorption spectroscopy was used to corroborate the presence and identify the oxidation state of Ce in the system and to observe if there are any electronic interactions within the 8 % Au/CeO x /C system. Cyclic voltammetry of electrodeposited 9 % Au/C and Ce-promoted 8 % Au/C showed a higher catalytic current density for ethanol oxidation when compared with commercially available catalysts (20 % Au/C) of a higher precious metal loading. Additionally, we report a higher stability toward the ethanol electrooxidation process, which was corroborated by 1 mV/s linear sweep voltammetry and chronoamperometric studies.

36 MATERIALS SCIENCE↗

Design and development of PEGylated liposomal formulation of HER2 blocker Lapatinib for enhanced anticancer activity and diminished cardiotoxicity

Highlights: • Pegylated Liposomal preparation enhances the cytotoxic potential of Lapatinib. • PEGylated liposomal formulation enhances apoptotic potential of Lapatinib on tumor cells. • PEGylated liposomal formulation diminishes cardio and hepatotoxic effect of Lapatinib. Breast cancer is most frequently diagnosed cancer and fifth leading cause of death in women. About 20–30% of all breast cancers overexpress HER2/neu receptors. Lapatinib is a dual tyrosin kinase inhibitor of EGFR and HER2. It exhibits its anticancer effect via blocking intracellular domain of HER2 receptor in breast cancer. Lapatinib belongs to class II of BSC classification due to its poor solubility restricting its clinical application. Due to presence of HER2 receptor on cardiomyocytes, it is associated with generation of cardiotoxicity. The present study was aimed to design a PEGylated liposomal formulation of Lapatinib and evaluate its anticancer potential. Lapatinib liposomes were prepared using lipid layer hydration method and its characterization was done by determining its particle size, zeta potential, entrapment efficiency and in vitro release profiling. The anti-tumor activity of PEGylated liposomal formulation was evaluated in xenografted tumor induced by MDA-MB-453 breast cancer cells in chick embryos. The anti-tumor effect of lapatinib was enhanced by its PEGylated liposomal preparation as it led to the reduction in tumor size to a greater extent compared to the embryos treated with free lapatinib. Flowcytometric analysis and immunoflurescence study using cleaved PARP antibody demonstrated the enhanced apoptotic potential of PEGylated liposomes of lapatinib. SGOT levels, marker for cardiotoxicity and hepatotoxicity, significantly decreased in serum of embryos treated with PEGylated liposomes of lapatinib compared to free drug treated embryos. Hence, the PEGylated liposomal formulation of lapatinib can be used as a therapeutic strategy against HER2 positive breast cancer either alone or in combination with conventional anticancer agents and hormonal therapies.

60 APPLIED LIFE SCIENCES↗

Surface evolution of synthetic bilgewater emulsion

Bilgewater is oil-in-water emulsion emitted from shipboard and it is a source of pollution to the ocean. Wastewater treatment needs fundamental knowledge of the bilgewater surface changes; however, limited information is available. We have reported the first surface characterization of synthetic bilgewater emulsion using time-of-flight secondary ion mass spectrometry (ToF-SIMS) coupled with optical microscopy. A NAVY synthetic bilgewater emulsion consisting of three oils and a detergent mixture is used as the synthetic bilgewater model. Both fresh and aged emulsion samples are analyzed to determine their droplet size distributions (DSDs) and surface chemical compositions. Our results show that fresh emulsions are largely mono-modal with hydrocarbon fragments as the main surface composition. Aged emulsions are also monomodal with slightly larger size. Both SIMS spectral comparison and Principal Component Analysis (PCA) show that some surfactant components appear on the fresh emulsion surface immediately while larger molecular weight components appear at the aged bilge droplet surface. Our results indicate that the oil-water interface evolves after emulsion droplet forms. More importantly, surface evolution not only changes the bilgewater DSD but also surface chemical composition and reactivity.

bilgewater, droplet size distribution, ToF-SIMS, s↗

Synchrotron-based measurement of aluminum agglomerates at motor conditions

Solid rocket propellant combustion is hindered by agglomeration of aluminum particles on its burning surface and determining the particle size has been a problem for half a century. The actual size of the agglomerates at motor pressures is unknown due to the opacity of the combustion plume, particularly at the elevated pressures seen in operational rocket motors. Sampling techniques can provide data at elevated pressure but may be biased due to the sampling method and do not provide information on the dynamics of agglomerate formation. Here this study uses time-resolved synchrotron x-ray imaging (with both absorption and phase contrast) to view aluminum agglomerate formation in situ at relevant rocket pressures. Finally, we have for the first time observed agglomerate formation at motor-relevant pressures in real time with unprecedented fidelity, providing critical data for understanding the combustion of aluminized solid rocket propellants.

36 MATERIALS SCIENCE↗

A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm

For decades, the Vlasov–Darwin model has been recognized to be attractive for particle-in-cell (PIC) kinetic plasma simulations in non-radiative electromagnetic regimes, to avoid radiative noise issues and gain computational efficiency. However, the Darwin model results in an elliptic set of field equations that renders conventional explicit time integration unconditionally unstable. We explore a fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions, which overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. The finite-difference scheme for Darwin field equations and particle equations of motion is space–time-centered, employing particle sub-cycling and orbit-averaging. This algorithm conserves total energy, local charge, canonical-momentum in the ignorable direction, and preserves the Coulomb gauge exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. Finally, we demonstrate the accuracy and efficiency properties of the algorithm with various numerical experiments in 2D–3V.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A scoping study of far-SOL main-wall protection limiters for steady-state operation of compact pilot plant tokamaks

We present a novel method for handling steady-state heat fluxes incident on the main wall of pilot plant-scale magnetic fusion devices, based on the utilization of protection limiters in the far scrape-off layer (SOL). This method helps avoid large plasma-wall gaps, without excessively compromising blanket performance. We present an optimization algorithm for determining the appropriate size and scale of these protection limiters given (1) probability distributions of SOL plasma parameters and (2) assumed risk tolerance. As part of this optimization, we have developed an analytic description of parallel heat fluxes across limiter shadows, and an objective cost function (the ‘Far-SOL Marginal Cost’) to quantify the impact that different main-wall thermal management design choices have on reactor capital cost. Applying the model to a midscale fusion pilot plant concept shows that making use of far-SOL protection limiters can reduce capital costs on the order of $500 M, relative to naively increasing the plasma-wall gap. Our analysis demonstrates that the far-SOL power decay length is the highest-leverage plasma assumption for thermal loading of the first wall, and the primary cost driver for main wall thermal management. The relative cost efficiency of protection limiters increases as assumptions on the far-SOL heat flux become more pessimistic. The concepts described in this paper motivate the further development of far-SOL protection limiters as part of larger efforts to design economical core-edge-wall compatible solutions for a fusion pilot plant.

Design under uncertainty↗

Chronological changes in soil biogeochemical properties of the glacier foreland of Midtre Lovénbreen, Svalbard, attributed to soil-forming factors

We report that Glacier forelands provide an excellent opportunity to investigate vegetation succession and soil development along the chronosequence; however, there are few studies on soil biogeochemical changes from environmental factors, aside from time. This study aimed to investigate soil development and biogeochemical changes in the glacier foreland of Midtre Lovenbreen, Svalbard, by considering various factors, including time. Eighteen vegetation and soil variables were measured at 38 different sampling sites of varying soil age, depth, and glacio-fluvial activity. Soil organic matter (SOM) was quantitatively measured, and the compositional changes in SOM were determined following size-density fractionation. In the topsoil, the soil organic carbon (SOC) and total nitrogen (N) content was found to increase along the soil chronosequence and were highly correlated with vegetation-associated variables. These findings suggest that plant-derived material was the main driver of the light fraction of SOM accumulation in the topsoil. The heavy fractions of SOM were composed of microbially transformed organic compounds, eventually contributing to SOM stabilization within short 90-yr deglaciation under harsh climatic conditions. In addition to time, the soil vertical profiles showed that other environmental parameters, also affected the soil biogeochemical properties. The high total phosphorous (P) content and electrical conductivity in the topsoil were attributed to unweathered subglacial materials and a considerable amount of inorganic ions from subglacial meltwater. The high P and magnesium content in the subsoil were attributed to parent materials, while the high sodium and potassium content in the surface soil were a result of sea-salt deposition. Glacio-fluvial runoff hampered ecosystem development by inhibiting vegetation development and SOM accumulation. This study emphasizes the importance of considering various soil-forming factors, including parent/subglacial materials, aeolian deposition, and glacio-fluvial runoff, as well as soil age, to obtain a comprehensive understanding of the ecosystem development in glacier forelands.

glacio-fluvial runoff↗

Numerical simulation of asteroid geometry variance on airburst threat

For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.

Airburst↗

A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions

Here, we extend a recently proposed fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (Φ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇ • A = 0 exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗