Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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

Enhanced diffusion bonding of alloy 617 using electric field-assisted sintering

The development of compact heat exchangers (CHXs) has gained increasing interest in many industries owing to their high thermal efficiency and reduced size. Diffusion bonding (DB) is an advantageous technique for fabricating CHXs. Alloy 617 is a candidate for manufacturing CHXs for high-temperature advanced nuclear reactors due to its elevated-temperature properties. Previous endeavors in DB of Alloy 617 were conducted by hot pressing (HP), which reported precipitates at the diffusion-bond interface, limited grain boundary (GB) migration, and significantly reduced high-temperature mechanical properties. To overcome these challenges, this study investigated DB of Alloy 617 using electric field-assisted sintering (EFAS). Stacks composed of three sheets were bonded with EFAS using different temperatures, pressures, and hold times. DB using HP as the zero-current analog of EFAS was also performed for comparison. The result shows that Cr- and Mo-rich precipitates were formed at the interface of the hot-pressed samples. The electric current and temperature in EFAS play a significant role in precipitation and GB migration. The electric current coupled with correct temperatures can effectively prevent precipitate formation at the interface and achieve excellent GB migration. Nanoscale Al-rich oxide was formed at the interface of the samples made by both HP and EFAS, but grain boundaries can ignore the nanoscale Al-oxide and migrate across the interface. The temperature, pressure, and hold time also affected diffusion. The temperature is a prerequisite for a successful GB migration, and GB migration can be enhanced by increasing pressure and hold time.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Moisture ingress in commercial steel drums: Water content determination, diffusion modelling and predicted permeation rates

Commercial steel drums underpin the global economy, playing a pivotal role in the storage and transportation of critical materials. Transported and stored materials, such as food, chemical and nuclear waste, can be sensitive to ambient conditions, particularly moisture that can enhance negative effects such as corrosion and material degradation. Although international standards and regulations are in place for the qualification of steel drums, there are no current testing requirements, established limits or boundaries for the permeation of moisture into the drums during transportation or storage. This work aims to provide insights into the moisture ingress over time into properly sealed steel drums and provides estimated moisture ingress rates over time through extrapolation. Water vapour transmission rate (WVTR) measurements through the gasket material at 10–40°C were 0.11–2.1 g/m 2 /day resulting in a permeation activation energy of 30.2 kJ/mol. Water sorption measurements and Karl Fischer titration (KFT) on ethylene propylene diene monomer (EPDM) gasket material revealed a decrease in equilibrium moisture saturation with increasing temperature. KFT measurements also revealed the presence of moisture within the adhesive and drum wall after exposure to ambient conditions. KFT and Fourier transform infrared spectroscopy (FTIR) show that moisture will desorb from the EPDM and drum wall after exposure to desiccating conditions, although a minimal amount of moisture will remain present. When sealed to the manufacturer's recommendations, the steel drums are effective in minimizing moisture ingress. In conclusion, in sealed empty drums, moisture ingress rates for 19-L drums were 0.4–1.5 mg/day at 25°C 15% relative humidity (RH) and increased to 7.1–8.8 mg/day at 40°C 90% RH, and moisture ingress rates for 210-L drums were 2.5 and 3.5 mg/day at field deployment conditions of 15.5°C 51.5% RH and 23°C 40% RH, respectively.

42 ENGINEERING↗
DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗
DOE OSTI2023

Potential Benefits in Remapping the Special Flood Hazard Area: Evidence from the U.S. Housing Market

A typical U.S. homebuyer's understanding of whether a property faces flood risk is based on whether the property is located inside the National Flood Insurance Program's (NFIP) Special Flood Hazard Area (SFHA). The SFHA boundary, however, may bias homebuyers’ perceptions of flood risk relative to unobserved true risk because the SFHA is an incomplete, and sometimes inaccurate, representation of flood hazards. Using a national dataset of property transactions, flood hazard data from the First Street Foundation, state-level measures of flood disclosure laws, and a spatially restricted triple-difference design, we distinguish capitalization effects of policy-driven (SFHA) versus quasi-objective (First Street) indicators of flood hazard. We identify these effects by assessing thousands of local spatial interactions between property SFHA designations, measures of quasi-objective flood hazard, and being in a state that mandates flood history disclosures. Being inside the SFHA generates a risk discount, but the signal is muted relative to underlying hazard exposure. Further, the SFHA signal can result in inefficient discounts for properties erroneously mapped in the SFHA. Disclosure requirements about flood history accentuate price effects for hazardous properties and result in more adequate risk internalization. In the absence of disclosures, however, we find either no or a weak risk signal for houses outside the SFHA that face flood hazard. Our results highlight potential benefits of updating SFHA boundaries to include all houses that may experience flooding, and argue in favor of requiring flood-related disclosures from sellers to improve the market's ability to internalize flood risk.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Giant Magnetoelectric Coupling and Magnetic-Field-Induced Permanent Switching in a Spin Crossover Mn(III) Complex

We investigate giant magnetoelectric coupling at a Mn 3+ spin crossover in [Mn III L]BPh 4 (L = (3,5-diBr-sal) 2 323) with a field-induced permanent switching of the structural, electric, and magnetic properties. An applied magnetic field induces a first-order phase transition from a high spin/low spin (HS-LS) ordered phase to a HS-only phase at 87.5 K that remains after the field is removed. We observe this unusual effect for DC magnetic fields as low as 8.7 T. The spin-state switching driven by the magnetic field in the bistable molecular material is accompanied by a change in electric polarization amplitude and direction due to a symmetry-breaking phase transition between polar space groups. The magnetoelectric coupling occurs due to a γη 2 coupling between the order parameter γ related to the spin-state bistability and the symmetry-breaking order parameter η responsible for the change of symmetry between polar structural phases. We also observe conductivity occurring during the spin crossover and evaluate the possibility that it results from conducting phase boundaries. We perform ab initio calculations to understand the origin of the electric polarization change as well as the conductivity during the spin crossover. Thus, in this work, we demonstrate a giant magnetoelectric effect with a field-induced electric polarization change that is 1/10 of the record for any material.

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

Evaluating wind profiles in a numerical weather prediction model with Doppler lidar

We use Doppler lidar wind profiles from six locations around the globe to evaluate the wind profile forecasts in the boundary layer generated by the operational global Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts (ECMWF). The six locations selected cover a variety of surfaces with different characteristics (rural, marine, mountainous urban, and coastal urban). We first validated the Doppler lidar observations at four locations by comparison with co-located radiosonde profiles to ensure that the Doppler lidar observations were of sufficient quality. The two observation types agree well, with the mean absolute error (MAE) in wind speed almost always less than 1 m s –1 . Large deviations in the wind direction were usually only seen for low wind speeds and are due to the wind direction uncertainty increasing rapidly as the wind speed tends to zero. Time–height composites of the wind evaluation with 1 h resolution were generated, and evaluation of the model winds showed that the IFS model performs best over marine and coastal locations, where the mean absolute wind vector error was usually less than 3 m s –1 at all heights within the boundary layer. Larger errors were seen in locations where the surface was more complex, especially in the wind direction. For example, in Granada, which is near a high mountain range, the IFS model failed to capture a commonly occurring mountain breeze, which is highly dependent on the sub-grid-size terrain features that are not resolved by the model. The uncertainty in the wind forecasts increased with forecast lead time, but no increase in the bias was seen. At one location, we conditionally performed the wind evaluation based on the presence or absence of a low-level jet diagnosed from the Doppler lidar observations. The model was able to reproduce the presence of the low-level jet, but the wind speed maximum was about 2 m s –1 lower than observed. This is attributed to the effective vertical resolution of the model being too coarse to create the strong gradients in wind speed observed. Our results show that Doppler lidar is a suitable instrument for evaluating the boundary layer wind profiles in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Assessment of Contaminant Bioaccumulation in Aquatic Biota on and Adjacent to the Oak Ridge Reservation—2020

This report provides information on contaminant concentrations in multiple wildlife prey species inhabiting or associated with water bodies on and downstream from the Oak Ridge Reservation (ORR), including regional reference sites. This information can be used to understand the nature and extent of contaminant exposure and transfer through the food web and to provide a baseline for evaluating temporal trends in contaminant exposure and accumulation. This information was gathered as part of the Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA) Five-Year Review (FYR) process, and it also addresses regulator comments on the annual Remediation Effectiveness Report (RER). This report summarizes the results of biological sampling completed in the spring and summer of 2020, prior to the 2021 FYR. Sampled locations included key sites upstream from the point sources of contaminants, along with watershed integration points on and off the ORR. Integration points are spatially the most important sampling sites because they capture biological responses to upstream contaminant inputs from both aqueous and groundwater sources, as well as responses to various remedial actions. These sites are most often located at a watershed administration unit boundary and may have regulatory significance for performance measures. In many cases, there are long-term biological and other monitoring data associated with these sites; these data generally fall within the areas of exposure and effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling Ti–6Al–4V using crystal plasticity, calibrated with multi-scale experiments, to understand the effect of the orientation and morphology of the α and β phases on time dependent cyclic loading

Classically, crystal plasticity modeling has used a range of constitutive equations, in which the incorporation of additional physics-based relationships typically results in additional model parameters. These additional parameters need to be reliably calibrated, which often necessitates the use of a range of experimental data acquired at multiple length scales. In this work, a crystal plasticity based finite element (CPFE) model for a dual-phase Titanium alloy, Ti–6Al–4V, is developed. The α and β phases of the microstructure are explicitly modeled. The model is calibrated using a systematic optimization routine and experimental data that consist of macroscopic stress-strain curves coupled with lattice strains on different crystallographic planes for the two phases. These experimental data were obtained from in situ high energy X-ray diffraction experiments for multiple material pedigrees, with varying crystallographic orientation distribution and β volume fractions. Depending on the thermomechanical-processing route and the heat treatment used to manufacture the alloy, Ti–6Al–4V can exist in a wide number of microstructural forms, which often results in the α and β phases either having well aligned slip systems (following the Burgers orientation relationship (BOR)) or possessing no alignment of the slip systems across the interphase boundary (not following the BOR). In this study, the fully-calibrated CPFE model is used to gain a comprehensive understanding of the deformation behavior of Ti–6Al–4V, specifically, the effect of microstructures that follow the BOR (or not) on time-dependent cyclic loading (including the effects of dwell hold times).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase-field modeling for restructuring in the dark zone of high burnup UO 2

This report summarizes the mesoscale modeling work performed in fiscal year 2024 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstructural evolution and restructuring observed in the dark regions of high burnup UO 2 nuclear fuel. This is the first attempt to realistically simulate the restructuring behavior observed in different region of a high burnup fuel. We employ a grand-potential based phase-field model to concurrently evaluate the formation of subgrains and growth of fission bubbles within the fuel. A energy-based subgrain formation criteria is introduced to simulate the restructuring process. Effect of different initial conditions and different modeling parameters are studies systematically to capture how each of these parameters influence the characteristics of the restructured fuel. It is observed that the subgrain formation begins around existing fission gas bubbles and then proceeds towards triple junctions, grain boundaries and grain interiors. It is demonstrated that restructuring is influenced by a combination of initial dislocation densities, subgrain formation rate, and temperature. Rate of restructuring increases with increase in fuel temperature. A restructuring bias is observed within the microstructure due to variation in defect accumulation among different grains. Furthermore, bubble sizes and distribution does not have a significant effect on rate of restructuring. The predicted microstructures resembles the characteristics of the restructured regions as observed in experiments. Finally, a correlation is presented that demonstrates the evolution of the restructuring volume fraction as a function of local effective burnup. This work provides a first of its kind restructuring model for darkzone that can be used by BISON for performance prediction of high burnup UO 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Effect of annealing and supercritical CO 2 exposure at 750°C on the tensile properties of stainless steel and Ni-based structural alloys

Subsize dogbone tensile specimens of alloys 304H, 310HCbN, and 740H were exposed to 30 MPa supercritical CO 2 (sCO 2 ) for up to 4000 h at 750°C. Similar exposures in laboratory air were used to assess the specific role of sCO 2 exposure on the alloy tensile properties. Type 304H specimens were only embrittled by exposure in sCO 2 due to the formation of a thick Fe-rich oxide and C ingress in the bulk of the alloy. For 310HCbN and 740H, no effect of the sCO 2 environment was observed and a decrease of the alloy's ductility was associated with previously observed microstructure evolutions, for example, grain boundary carbide formation. The 740H was tested in the solution annealed condition and self-aging took place during 750°C exposures in sCO 2 and air with the formation of fine gamma prime precipitates leading to significant tensile strength increase. It was found that both 740H and 310HCbN have adequate sCO 2 compatibility for high-temperature commercial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The R-matrix bootstrap for the 2d O( N ) bosonic model with a boundary

The S-matrix bootstrap is extended to a 1+1d theory with O(N) symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Validation and moisture content sensitivity analysis of cross-laminated timber wall assemblies in EnergyPlus

Cross-laminated timber buildings are becoming more common in North America, with many numerical studies showing potential energy savings. However, no studies have validated any EnergyPlus heat transfer algorithms or quantified their accuracy in simulating CLT in building envelopes. This study empirically validates the heat flux predictions for each of EnergyPlus's heat transfer algorithms (Conduction Transfer Functions (CTF), Effective Moisture Penetration Depth (EMPD), Conduction Finite Difference (CondFD), and Heat and Moisture Transfer (HAMT)) for two different CLT ply thicknesses with both summer and winter boundary conditions measured in controlled lab experiments. It also evaluates the model sensitivity of heat flux and heating and cooling loads to moisture content. The 1D validation shows that the HAMT model is the most accurate among all algorithms. All EnergyPlus's heat flux predictions are accurate independent of CLT plate thickness for summer conditions. However, the three constant property algorithms (CTF, EMPD, and CondFD) underpredict heat flux throughout the whole day during winter conditions. The 1D sensitivity analysis indicates that elevated moisture content can increase peak heat fluxes through the material by up to 20 %. Finally, the whole building model sensitivity analysis shows increased heating load and slight cooling load variation due to increased moisture content when using constant property models. The analysis shows significantly lower peak thermal demand (7 % lower heating and 6 % lower cooling) and monthly thermal load (8 % less cooling and 6 % less heating) predictions when using HAMT vs a constant property model.

42 ENGINEERING↗
DOE OSTI2022

Dirac Nodal Line in Hourglass Semimetal Nb 3 SiTe 6

Glide-mirror symmetry in nonsymmorphic crystals can foster the emergence of novel hourglass nodal loop states. Here, we present spectroscopic signatures from angle-resolved photoemission of a predicted topological hourglass semimetal phase in Nb 3 SiTe 6 . Linear band crossings are observed at the zone boundary of Nb 3 SiTe 6 , which could be the origin of the nontrivial Berry phase and are consistent with a predicted glide quantum spin Hall effect; such linear band crossings connect to form a nodal loop. Furthermore, the saddle-like Fermi surface of Nb 3 SiTe 6 observed in our results helps unveil linear band crossings that could be missed. Additionally, in situ alkali-metal doping of Nb 3 SiTe 6 also facilitated the observation of other band crossings and parabolic bands at the zone center correlated with accidental nodal loop states. Overall, our results complete the system's band structure, help explain prior Hall measurements, and suggest the existence of a nodal loop at the zone center of Nb 3 SiTe 6 .

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Charged Biological Membranes Repel Large Neutral Molecules by Surface Dielectrophoresis and Counterion Pressure

Macromolecular crowding is the usual condition of cells. The implications of the crowded cellular environment for protein stability and folding, protein–protein interactions, and intracellular transport drive a growing interest in quantifying the effects of crowding. While the properties of crowded solutions have been extensively studied, less attention has been paid to the interaction of crowders with the cellular boundaries, i.e., membranes. However, membranes are key components of cells and most subcellular organelles, playing a central role in regulating protein channel and receptor functions by recruiting and binding charged and neutral solutes. While membrane interactions with charged solutes are dominated by electrostatic forces, here we show that significant charge-induced forces also exist between membranes and neutral solutes. Using neutron reflectometry measurements and molecular dynamics simulations of poly(ethylene glycol) (PEG) polymers of different molecular weights near charged and neutral membranes, we demonstrate the roles of surface dielectrophoresis and counterion pressure in repelling PEG from charged membrane surfaces. The resulting depletion zone is expected to have consequences for drug design and delivery, the activity of proteins near membrane surfaces, and the transport of small molecules along the membrane surface.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Substantial contribution of iodine to Arctic ozone destruction

Unlike bromine, the effect of iodine chemistry on the Arctic surface ozone budget is poorly constrained. We present ship-based measurements of halogen oxides in the high Arctic boundary layer from the sunlit period of March to October 2020 and show that iodine enhances springtime tropospheric ozone depletion. We find that chemical reactions between iodine and ozone are the second highest contributor to ozone loss over the study period, after ozone photolysis-initiated loss and ahead of bromine.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗
DOE OSTI2022

Numerical simulation of a relativistic magnetron using a fluid electron model

An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. Here we show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child–Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗