Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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

An automated approach to the alignment of compound refractive lenses

Compound refractive lenses (CRLs) are established X-ray focusing optics, and are used to focus the beam or image the sample in many beamlines at X-ray facilities. While CRLs are quite established, the stack of single lens elements affords a very small numerical aperture because of the thick lens profile, making them far more difficult to align than classical optical lenses that obey the thin-lens approximation. This means that the alignment must be very precise and is highly sensitive to changes to the incident beam, often requiring regular readjustments. Some groups circumvent the full realignment procedure by using engineering controls ( e.g. mounting optics) that sacrifice some of the beam's focusing precision, i.e. spot size, or resolution. While these choices minimize setup time, there are clear disadvantages. This work presents a new automated approach to align CRLs using a simple alignment apparatus that is easy to adapt and install at different types of X-ray experiments or facilities. This approach builds on recent CRL modeling efforts, using an approach based on the Stochastic Nelder–Mead (SNM) simplex method. This method is outlined and its efficacy is demonstrated with numerical simulation that is tested in real experiments conducted at the Advanced Photon Source to confirm its performance with a synchrotron beam. This work provides an opportunity to automate key instrumentation at X-ray facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quantifying the Influence of Size, Shape, and Density of Microplastics on Their Transport Modes: A Modeling Approach

Microplastics (MPs) pose significant risks to marine ecosystems and human health, necessitating accurate predictions of their distributions in aquatic environments for effective risk mitigation. However, understanding MP transport dynamics is challenging because of the inadequate representation of MP characteristics such as size, shape, and density in numerical models. Further, the accuracy of the MP vertical profiles in existing models has not been thoroughly validated. Thus, we developed an MP transport model within the Finite Volume Community Ocean Model framework (FVCOM-MP) by integrating MP characteristics. We validated FVCOM-MP against experimental and analytical data, focusing on various MP transport modes and transitions. FVCOM-MP successfully replicates MP profiles in different transport modes, including the bedload, surface-load, suspended-load, and mixed-load modes. Additionally, we introduce phase diagrams for classifying MP transport modes based on particle characteristics, enhancing our understanding of MP dynamics in aquatic systems. The transport modes for a number of real-world MP particles, including fishing line, plastic bag/bottle fragments, synthetic fibers, tire wear particles, polyvinyl chloride and expanded polystyrene pellets, were analyzed with these phase diagrams.

Microplastic transport, Settling velocity, Rising ↗
DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Effect of toolpath in large-format additive manufacturing with bio-derived composites

There has been growing interest in integrating bio-derived composites into Large-format additive manufacturing (LFAM) feedstocks to reduce the use of petroleum-derived materials and reduce the overall carbon footprint of LFAM. However, these materials present unique challenges during manufacturing due to their variability, which can lead to unintended deformations and failures attributable to suboptimal process conditions. While numerical modelling has been extensively employed to simulate numerous manufacturing processes, its application in LFAM with bio-derived composites remains limited. This study addresses this gap by systematically developing a numerical model to simulate the LFAM process using bio-based materials, specifically wood fiber-reinforced polylactic acid (PLA/WF). Experimental investigations were conducted to characterise the thermal and mechanical properties of additively manufactured PLA/WF specimens. Numerical simulations were performed to predict temperature profiles and deformations during LFAM. The effect of varying infill patterns, internal structures, and tool paths on the temperature distribution and deformation of printed parts was explored using the developed model. This article aims to advance the utilisation of bio-derived composites in LFAM systems and provide a comprehensive understanding of the LFAM process. The findings offer valuable insights for optimising process parameters and enhancing the performance of LFAM with bio-based composites.

Large-format additive manufacturing↗
DOE OSTI2022

Energy and mobility impacts of connected autonomous vehicles with co-optimization of speed and powertrain on mixed vehicle platoons

Intersections are known to be traffic bottlenecks where a significant amount of energy consumption could be caused due to deceleration/acceleration in the presence of red signals. With an increased level of connectivity and automation of intelligent transportation systems, connected autonomous vehicles (CAVs) are expected to be able to proactively adjust their driving strategies subject to constraints imposed by the predicted future traffic. As a result, many potential benefits can be achieved, such as improved energy efficiency, enhanced traffic safety, among many others. Notably, the way CAVs are controlled affects the following legacy vehicles (LVs) due to complex traffic dynamics. Here, we are particularly interested in studying the energy and mobility impact of CAVs with an improved traffic prediction method on mixed vehicle platoons at various market penetration rates. Leveraging traffic prediction, CAVs are controlled with co-optimization of their speed and gear position. Specifically, a traffic prediction framework in a rolling horizon fashion is employed based upon a modified Payne–Whitham (PW) model capable of handling mixed traffic consisting of CAVs and LVs. The prediction error of the modified PW model is reduced by 53.62% compared to that of the standard PW model under test scenarios. According to the predicted traffic conditions, speed and gear position of CAVs are co-optimized with the primary goal of minimizing energy consumption when driving on a signalized arterial. The energy benefits achieved by CAVs and the impact of CAVs on LVs behind are studied comprehensively for mixed vehicle platoons. The lead LV follows a real-world speed profile collected on TH-55 in Minnesota. Numerical results show that energy benefits achieved by the vehicle platoon range from 2% to 16%, and a 1% to 5% reduction in travel time for LVs behind CAVs is also observed, at different penetration rates of CAVs in various traffic scenarios. Furthermore, it is observed that CAVs using the proposed eco-driving approach appear to have a positive impact on the LVs behind in terms of energy consumption, regardless of the driving styles of the LVs ahead.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Response of Atmospheric River Width and Intensity to Aquaplanet Warming: A Detection Algorithm‐ and Background Moisture‐Independent Approach

The width of an atmospheric river (AR) is an important parameter when evaluating its impact. Although previous research suggests ARs will widen with global warming, a precise response has been muddled by the large sensitivity of width to a diverse set of AR detection techniques (ARDTs). Here, we propose a methodology that removes the influence of the ARDT by modeling AR‐integrated vapor transport (IVT) profiles as idealized exponential curves with free parameters given by background IVT, intensity above background IVT, and profile width. Notably, our definition for AR profile width does not include any explicit numerical thresholds, relative or absolute, for IVT. We apply our approach to a series of idealized aquaplanet experiments, first with a baseline sea surface temperature (SST) distribution, and then with +2K, +4K, and +6K uniform warming, so as to determine the contributions of each free parameter to AR width. We also apply our approach to high‐resolution atmosphere‐only models forced with SSTs modified to preserve historical variability but following projected warming over 2016–2050. Our results show that contributions to impacts‐relevant AR widening comes primarily from enhancements in background IVT and AR intensity, as opposed to from dynamic width changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data associated with Quaternary Research manuscript Latest Pleistocene glacial chronology and paleoclimate reconstruction for the East River watershed, Colorado, USA

The data here are associated with Quaternary Research manuscript Latest Pleistocene glacial chronology and paleoclimate reconstruction for the East River watershed, Colorado, USA and include data associated with cosmogenic exposure and depth profile dating as well as glacier-climate numerical modeling. Reconstructing Pleistocene glaciation timing and extent is vital for understanding paleoclimate. While late Pleistocene glaciation has been studied extensively in western North American mountain ranges, the glacial history of the western Elk Range in Colorado remains understudied, particularly in the East River watershed, a site of intense scientific focus. Here we use cosmogenic nuclide exposure and depth–profile dating methods to determine the timing of glaciation in the East River watershed. We use glacier modeling to reconstruct paleoglacier extents and quantify past climate conditions. Our findings indicate that the East River glacier retreated from its maximum position approximately 17–18 ka, moving to recessional positions between 13–15 ka, before experiencing more substantial retreat to high-elevation cirques around 12.9 ka. Glacier modeling suggests that the maximum ice extents at 17–18 ka could have been sustained by temperature depressions of approximately −6.5°C compared to modern conditions, assuming consistent precipitation. Additionally, the ice position at 13–15 ka could have been supported by temperature depressions of around −4.0°C. These results offer insights into the deglaciation timeline in the East River watershed and broader western Elk Range as well as paleoclimate conditions during the late Pleistocene, which may aid future research on critical zone evolution in the East River watershed.The data files include: 1. Table 1 In situ-produced 10Be sample data (.pdf, .csv, .xlsx)2. Table 2 In situ-produced 10Be exposure age results (.pdf, .csv, .xlsx)3. S1_A CRONUS calculator input assuming no erosion (.csv)4. S1_B CRONUS calculator input assuming erosion (.csv)5. S1_C CRONUS calculator results and comparison (.csv)6. S2_A In situ-produced 10Be sample data for depth profile (.csv)7. S2_B In situ-produced 10Be depth profile model input (.csv)8. S2_C In situ-produced 10Be depth profile model results (.csv)9. S3 Monthly cloudiness, rH, and windspeed data used in glacier climate model (.csv)10. S4 Glacier climate model results (.csv)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Preventive Power Outage Estimation Based on A Novel Scenario Clustering Strategy: Preprint

The increasing occurrence of extreme weather events is challenging the power grid operation. In front of the extreme weather, the system operator is responsible for estimating the power outage and scheduling the restoration resources. This paper proposes an outage evaluation framework to identify the possible unserved load profiles, vulnerable areas, and mobile energy adequacy. The predicted vulnerable lines of an outage prediction model tool are utilized to generate numerous faulted line scenarios. Next, each scenario's nodal unserved load profile is obtained by solving a three-phase restoration model that considers the schedule of repair crews and mobile energy resources. Then, a novel scenario clustering strategy is developed to cluster the unserved load profiles into multiple representative ones for straightforward analysis. Finally, case studies on a distribution system evaluate the damage level brought by extreme weather and verify the effectiveness of the proposed scenario clustering strategy.

mobile energy resources↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High Temperature Acclimation of Leaf Gas Exchange, Photochemistry, and Metabolomic Profiles in Populus trichocarpa

High temperatures alter the thermal sensitivities of numerous physiological and biochemical processes that impact tree growth and productivity. Foliar and root applications of methanol have been implicated in plant acclimation to high temperature via the C 1 pathway. Here, we characterized temperature acclimation at 35°C of leaf gas exchange, chlorophyll fluorescence and extractable metabolites of potted Populus trichocarpa saplings and examined potential influences of mM concentrations of methanol added during soil watering over a two-month period. Relative to plants grown under the low growth temperature (LGT), high growth temperature (HGT) plants showed a suppression of leaf water use and carbon cycling including transpiration (E), net photosynthesis (P n ), an estimate of photorespiration (R p ), and dark respiration (R d ), attributed to reductions in stomatal conductance and direct negative effects on gas exchange and photosynthetic machinery. In contrast, HGT plants showed an upregulation of non-photochemical quenching (NPQ t ), the optimum temperature for ETR, and leaf isoprene emissions at 40°C. A large number of metabolites (867) were induced under HGT, many implicated in flavonoid biosynthesis highlighting a potentially protective role for these compounds. Methanol application did not significantly alter leaf gas exchange but slightly reduced the suppression of R d and R p by the high growth temperature while slightly impairing ETR, Fv’/Fm’, and q p . However, we were unable to determine if soil methanol was sufficiently taken up by the plant to have a direct effect on foliar processes. A small number of extracted leaf tissue metabolites (55 out of 10,015) showed significantly altered abundances under LGT and methanol treatments relative to water controls, and this increased in compound number (222) at the HGT. The results demonstrate the large physiological and biochemical impacts of high growth temperature on poplar seedlings and highlights the enhancement of the optimum temperature of ETR as a rapid thermal acclimation mechanism. Finally, although no large effect on leaf physiology was observed, the results are consistent with methanol both impairing photochemistry of the light reactions via formaldehyde toxicity, and stimulating photosynthesis and dark respiration through formate oxidation to CO 2 .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Impurity leakage and radiative cooling in the first nitrogen and neon seeding study in the closed DIII-D SAS configuration

A comparative study of nitrogen versus neon has been carried out to analyze the impact of the two radiative species on power dissipation, SOL impurity distribution, divertor and pedestal characteristics. The experimental results show that N remains compressed in the divertor, thereby providing high radiative losses without affecting the pedestal profiles and displacing carbon as dominant radiator. Neon, instead, radiates more upstream than N thus reducing the power flux through the separatrix leading to a reduced ELM frequency and compression in the divertor. A significant amount of neon is measured in the plasma core leading to a steeper density gradient. The different behavior between the two impurities is confirmed by SOLPS-ITER modeling which for the first time at DIII-D includes multiple impurity species and a treatment of full drifts, currents and neutral–neutral collisions. The impurity transport in the SOL is studied in terms of the parallel momentum balance showing that N is mostly retained in the divertor whereas Ne leaks out consistent with its higher ionization potential and longer mean free path. This is also in agreement with the enrichment factor calculations which indicate lower divertor enrichment for neon. The strong ionization source characterizing the SAS divertor causes a reversal of the main ions and impurity flows. The flow reversal together with plasma drifts and the effect of the thermal force contribute significantly in the shift of the impurity stagnation point affecting impurity leakage. This work provides a demonstration of the impurity leakage mechanism in a closed divertor structure and the consequent impact on pedestal. Since carbon is an intrinsic radiator at DIII-D, in this paper we have also demonstrated the different role of carbon in the N vs Ne seeded cases both in the experiments and in the numerical modeling. Here, carbon contributes more when neon seeding is injected compared to when nitrogen is used. Finally, the results highlight the importance of accompanying experimental studies with numerical modeling of plasma flows, drifts and ionization profile to determine the details of the SOL impurity transport as the latter may vary with changes in divertor regime and geometry. In the cases presented here, plasma drifts and flow reversal caused by high level of closure in the slot upper divertor at DIII-D play an important role in the underlined mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Lidar-Based Evaluation of HRRR Performance in California’s Diablo Range

The performance of the NOAA High-Resolution Rapid Refresh (HRRR) model for capturing low-level winds near a wind energy production site during summer 2019 is evaluated. This study catalogs the ability of HRRR to predict boundary layer dynamics relevant to wind energy interests over complex terrain, which has presented challenges for weather and energy forecasting. Performance is evaluated by comparing HRRR output to wind-profiling Doppler lidars at Lawrence Livermore National Laboratory Site 300. HRRR captured the diurnal profile of horizontal winds in the observed 150-m layer, despite strong underpredictions (∼4 m s −1 ) during evening and nighttime hours. These underpredictions may be a result of local speedup flows observed by the lidars, which were unresolved in HRRR due to their small spatial extent. HRRR bias magnitude relative to observations was found to be minimal during days with synoptic-scale troughs and strong 850-hPa geopotential gradients, while bias magnitude was maximal during days with synoptic ridging and weak 850-hPa geopotential gradients. To translate wind speed predictions to energy forecasting, generic turbine models were used to estimate power generation for turbines characteristic of the nearby Altamont Pass Wind Resource Area. Results show that HRRR-based energy estimates predicted daytime power generation adequately relative to lidar-based estimates with an 18-h lead time (bias magnitude < 0.4 MW from 0900 to 1400 LT) but overpredicted power during the rest of the diurnal cycle (bias > 1 MW). These results demonstrate conditions under which HRRR performs well for wind energy applications in complex terrain, while highlighting biases that require further investigation to support usage of a high-resolution model for wind energy forecasts.

Boundary layer↗
DOE OSTI2022

Comparison of scattering ratio profiles retrieved from ALADIN/Aeolus and CALIOP/CALIPSO observations and preliminary estimates of cloud fraction profiles

The space-borne active sounders have been contributing invaluable vertically resolved information of atmospheric optical properties since the launch of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) in 2006. To build long-term records from space-borne lidars useful for climate studies, one has to understand the differences between successive space lidars operating at different wavelengths, flying on different orbits, and using different viewing geometries, receiving paths, and detectors. In this article, we compare the results of Atmospheric Laser Doppler INstrument (ALADIN) and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) lidars for the period from 28 June to 31 December 2019. First, we build a dataset of ALADIN–CALIOP collocated profiles (Δdist<1°; Δtime<6 h). Then we convert ALADIN's 355 nm particulate backscatter and extinction profiles into the scattering ratio vertical profiles SR(z) at 532 nm using molecular density profiles from Goddard Earth Observing System Data Assimilation System, version 5 (GEOS-5 DAS). And finally, we build the CALIOP and ALADIN globally gridded cloud fraction profiles CF(z) by applying the same cloud detection threshold to the SR(z) profiles of both lidars at the same spatial resolution. Before comparing the SR(z) and CF(z) profiles retrieved from the two analyzed lidar missions, we performed a numerical experiment to estimate the best achievable cloud detection agreement CDA norm (z) considering the differences between the instruments. We define CDA norm (z) in each latitude–altitude bin as the occurrence frequency of cloud layers detected by both lidars, divided by a cloud fraction value for the same latitude–altitude bin. We simulated the SR(z) and CF(z) profiles that would be observed by these two lidars if they were flying over the same atmosphere predicted by a global model. By analyzing these simulations, we show that the theoretical limit for CDA$_{norm}^{theor}$(z) for a combination of ALADIN and CALIOP instruments is equal to 0.81±0.07 at all altitudes. In other words, 19 % of the clouds cannot be detected simultaneously by two instruments due to said differences. The analyses of the actual observed CALIOP–ALADIN collocated dataset containing ~78 000 pairs of nighttime SR(z) profiles revealed the following points: (a) the values of SR(z) agree well up to ~3 km height. (b) The CF(z) profiles show agreement below ~3 km, where ~80 % of the clouds detected by CALIOP are detected by ALADIN as expected from the numerical experiment. (c) Above this height, the CDA$_{norm}^{obs}$(z) reduces to ~50 %. (d) On average, better sensitivity to lower clouds skews ALADIN's cloud peak height in pairs of ALADIN–CALIOP profiles by ~0.5±0.6 km downwards, but this effect does not alter the heights of polar stratospheric clouds and high tropical clouds thanks to their strong backscatter signals. (e) The temporal evolution of the observed CDA$_{norm}^{obs}$(z) does not reveal any statistically significant change during the considered period. This indicates that the instrument-related issues in ALADIN L0/L1 have been mitigated, at least down to the uncertainties of the following CDA$_{norm}^{obs}$(z) values: 68±12 %, 55±14 %, 34±14 %, 39±13 %, and 42±14 % estimated at 0.75, 2.25, 6.75, 8.75, and 10.25 km, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Blade Designs For Improved Multi-Phase Performance In sCO 2 Compressors; Part II - Optical Diagnostics In sCO 2 And Experimental Evaluation With Particle Image Velocimetry

This paper presents the second part of a study in which the leading-edge and suction surface of a compressor blade was modified to delay onset of phase change for sCO 2 compressors operating near the critical point. Using a first-of-its-kind apparatus for the measurement of sCO 2 flow fields, Particle Image Velocimetry (PIV) is used for local flow field measurements of two compressor blade geometries: the modified “biased-wedge,” and a conventional constant thickness blade. Utilizing the developed hardware, the feasibility of a simple, laser-based diagnostic for qualitatively measuring liquid phase regions, is also presented. The design of the optical diagnostics rig, a discussion of numerous challenges, and necessary considerations involved in performing optical-based measurements like PIV, in sCO 2 , are discussed. Velocity field measurements for the modified compressor profile show a much lower suction peak compared to a conventional blade. Furthermore, these results validate numerical results at the tested conditions, where the suction side profile of the biased wedge works to minimize the local pressure gradient.

14 SOLAR ENERGY↗
DOE OSTI2023

Defect‐Free Sound Insulator Using Single Metal‐Based Friction Stir Process Array

Metals are excellent conductors for phonon transportation such as vibration, sound, and heat. Generally, metal sound insulators require multimaterial structure or defects and unimetal sound insulators are challenging. Therefore, a design of a defect-free sound insulator made by single alloys with multiple friction stir processes (FSPs) is proposed. Periodic friction stir processing can induce superlattice-like local mechanical properties’ modifications. By experimental acoustic characterization, it is observed that FSP can introduce clear acoustic–elastic property contrast on an aluminum plate by the presence of stir zone and heat-affected zones. In numerical simulations, the signature FSP-induced property profile is periodically and parallelly arranged on a long aluminum plate. The transmission gap frequencies are present on the frequency spectrum with the sound propagation direction perpendicular to the FSP paths. Disorder offsets on FSP periodicity are further introduced. Anderson localization is found on a resonance frequency, which provides −11 dB sound reduction by an exponential decay. Due to the finite design length, the slight disorder can also enhance sound insulation in the periodic transmission gap frequency. In conclusion, with analysis and comparison with different configurations, the best performance in the models can achieve −30 dB sound insulation in the 350 mm-long aluminum alloy plate with 14 parallel FSPs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2023

Effects of collisional ion orbit loss on tokamak radial electric field and toroidal rotation in an L-mode plasma

Abstract Ion orbit loss has been used to model the formation of a strong negative radial electric field E r in the tokamak edge, as well as edge momentum transport and toroidal rotation. To quantitatively measure ion orbit loss, an orbit-flux formulation has been developed and numerically applied to the gyrokinetic particle-in-cell code XGC. We study collisional ion orbit loss in an axisymmetric DIII-D L-mode plasma using gyrokinetic ions and drift-kinetic electrons. Numerical simulations, where the plasma density and temperature profiles are maintained through neutral ionization and heating, show the formation of a quasisteady negative E r in the edge. We have measured a radially outgoing ion gyrocenter flux due to collisional scattering of ions into the loss orbits, which is balanced by the radially incoming ion gyrocenter flux from confined orbits on the collisional time scale. This suggests that collisional ion orbit loss can shift E r in the negative direction compared to that in plasmas without orbit loss. It is also found that collisional ion orbit loss can contribute to a radially outgoing (counter-current) toroidal-angular-momentum flux, which is not balanced by the toroidal-angular-momentum flux carried by ions on the confined orbits. Therefore, the edge toroidal rotation shifts in the co-current direction on the collisional time scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗