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 163 records · Page 9

SOWFA Large-Eddy Simulation of a Diurnal Cycle at the SWiFT Site

These are large-eddy simulation (LES) data output from the Simulator for Wind Farm Applications (SOWFA) with internal mesoscale-to-microscale coupling. A diurnal cycle was simulated over the Scaled Wind Farm Technology (SWIFT) facility near Lubbock, Texas, from 2013-11-08 12:00:00 UTC to 2013-11-09 12:00:00 UTC. Indirect profile assimilation was used to enforce nonstationary mesoscale conditions. Mesoscale forcing data were extracted from WRF numerical weather prediction outputs or derived from observations at SWiFT.

17 WIND ENERGY↗

Turbulence dissipation rate estimated from lidar observations during the LAPSE-RATE field campaign

The International Society for Atmospheric Research using Remotely-piloted Aircraft (ISARRA) hosted a flight week in July 2018 to demonstrate unmanned aircraft systems' (UASs) capabilities in sampling the atmospheric boundary layer. This week-long experiment was called the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE) field campaign. Numerous remotely piloted aircraft and ground-based instruments were deployed with the objective of capturing meso- and microscale phenomena in the atmospheric boundary layer. The University of Oklahoma deployed one Halo Streamline lidar, and the University of Colorado Boulder deployed two WindCube lidars. In this paper, we use data collected from these Doppler lidars to estimate turbulence dissipation rate throughout the campaign. We observe large temporal variability of turbulence dissipation close to the surface with the WindCube lidars that is not detected by the Halo Streamline. However, the Halo lidar enables estimating dissipation rate within the whole boundary layer, where a diurnal variability emerges. We also find a higher correspondence in turbulence dissipation between the WindCube lidars, which are not co-located, compared to the Halo and WindCube lidar that are co-located, suggesting a significant influence of measurement volume on the retrieved values of dissipation rate.

WIND ENERGY↗

Ionization wave propagation in a He plasma jet in a controlled gas environment

Characterizing ionization wave propagation in low temperature plasma jets is critical to predicting production of reactive species and plasma–surface interactions for biomedical applications and surface functionalization. In this work, results from optical emission and laser induced fluorescence measurements of the ionization wave in a He plasma jet operating in a controlled gas environment are discussed and used for comparison with numerical modeling. Additionally, the ionization wave was observed using ICCD (Intensified Charge Coupled Device) imaging and characterized by time and spatially resolved electron density measurements using laser-collision-induced fluorescence. The plasma jet was initially characterized using pure He (nominally at 200 Torr), while varying pressure and voltage. When operating in pure He, the ionization wave broadly expands exiting the plasma tube. Increasing the operating pressure reduces the speed and isotropic expansion of the ionization wave. Furthermore, the jet operated with a humid He shroud was also studied. The humid He shroud results in the electron density increasing and having an annular profile due to the lower ionization potential of H2O compared to He and localized photoionization in the mixing region. Numerical modeling highlighted the importance of resonance radiation emitted by excited states of He, photoelectron emission from the quartz tube, and the kinetic behavior of the electrons produced by photoionization ahead of the ionization front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Construction of wave dark matter halos: Numerical algorithm and analytical constraints

Here we present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos—such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes and show how it simplifies in the high-energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schrödinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Inferring the temperature profile of the radiative shock in the COAX experiment with shock radiography, Dante, and spectral temperature diagnostics

Predicting and modeling the behavior of experiments with radiation waves propagating through low-density foams require a detailed quantification of the numerous uncertainties present. In regimes where a prominent radiative shock is produced, key dynamical features include the shock position, temperature, and curvature and the spatial distribution and temperature of the corresponding supersonic radiation wave. The COAX experimental platform is designed to constrain numerical models of such a radiative shock propagating through a low-density foam by employing radiography for spatial and shock information, Dante for characterizing the x-ray flux from the indirectly driven target, and a novel spectral diagnostic designed to probe the temperature profile of the wave. In this work, we model COAX with parameterized 2D simulations and a Hohlraum-laser modeling package to study uncertainties in diagnosing the experiment. The inferred temperature profile of the COAX radiation transport experiments has been shown to differ from simulations more than expected from drive uncertainties that have been constrained by simultaneous soft x-ray flux and radiography measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of high-fidelity ion cyclotron range of frequencies antenna coupling simulations in full 3D geometry against experiments in the ASDEX Upgrade tokamak

We address the validation of finite elements ion cyclotron range of frequencies (ICRF) antenna coupling simulation against experiments performed in the ASDEX Upgrade tokamak. Measurements of the loading resistance in ICRF-heated, magnetically-perturbed 3D plasma discharges are compared against numerical predictions of the RAPLICASOL code. To this end, the 3D induction field and the 3D density profile are modeled by concatenating the PARVMEC, BMW and EMC3-EIRENE codes. The 3D density is input to RAPLICASOL, where full-wave simulations are performed on a finite element mesh retaining full 3D geometry in the ICRF antenna model and the plasma description. The results are further compared with RAPLICASOL simulations using a 1D density profile as measured at the outboard midplane in the same experiments. We find that simulations using a 1D density profile overestimate the change in loading resistance by a factor of ~197 -248%, while simulations using the full 3D density profile are in agreement with experiments within a factor ~16%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resilient Operating Constraints for Power Distribution Systems under Setpoint Attacks

Integration and operation of distributed generation (DG) and energy storage (ES) in power distribution systems are enabled by communication networks and embedded sensor and control devices that increase the vulnerability of the systems to cyber-threats, broadening the attack surface and making adversary actions more unpredictable. This paper proposes a methodology that uses ellipsoidal approximations to quantify the potential damage caused by successful attacks that affect, directly or indirectly, the desired operation setpoints and may drive the power distribution operation to unsafe states by violating the limits of voltage or line flows. More specifically, a new methodology is introduced to find the optimal non-symmetric operating constraints that can be imposed to each DG and ES in order to guarantee that the power distribution system is resilient to any malicious setpoints. The proposed method takes as inputs the system topology, DG and ES capabilities, and load limits to solve a convex optimization problem formulated using linear matrix inequalities (LMIs) and the power flow equations. The proposed solution is agnostic to the attacker's action or load profile and it does not require any assumption about the location or means of the attack. The numerical results on a test distribution feeder with several DG and ES illustrate how the proposed resilient operating constraints guarantee the security of the power distribution system under setpoint attacks.

Giraldo, Jairo↗

Thermal chains and entrainment in cumulus updrafts, Part 2: Analysis of idealized simulations

Research has suggested that the structure of deep convection often consists of a series of rising thermals, or “thermal chain”, which contrasts with existing conceptual models that are used to construct cumulus parameterizations. In this work, simplified theoretical expressions for updraft properties obtained in Part 1 of this study are used to develop a hypothesis explaining why this structure occurs. In this hypothesis, cumulus updraft structure is strongly influenced by organized entrainment below the updraft’s vertical velocity maximum. In a dry environment, this enhanced entrainment can locally reduce condensation rates and increase evaporation, thus eroding buoyancy. For moderate-to-large initial cloud radius R, this breaks up the updraft into a succession of discrete pulses of rising motion (i.e., a thermal chain). For small R, this leads to the structure of a single, isolated rising thermal. In contrast, moist environments are hypothesized to favor plume-like updrafts for moderate-to-large R. In a series of axisymmetric numerical cloud simulations, R and environmental relative humidity (RH) are systematically varied to test this hypothesis. Vertical profiles of fractional entrainment rate, passive tracer concentration, buoyancy, and vertical velocity from these runs agree well with vertical profiles calculated from the theoretical expressions in Part 1. Analysis of the simulations supports the hypothesized dependency of updraft structure on R and RH, that is, whether it consists of an isolated thermal, a thermal chain, or a plume, and the role of organized entrainment in driving this dependency. Additional 3-dimensional (3-D) turbulent cloud simulations are analyzed, and the behavior of these 3-D runs is qualitatively consistent with the theoretical expressions and axisymmetric simulations.

54 ENVIRONMENTAL SCIENCES↗

Finite-difference multiple fluid solution for source-driven rotation in highly magnetized linear plasma device

The rotation profile of a magnetized plasma cylinder composed of multiple fluids is investigated analytically, expanding on previous results. The analytic steady-state solution is used as a benchmark for a time-dependent multiple-fluid finite-difference code, MITNS: Multiple-Ion Transport Numerical Solver. Magnetic field evolution is taken into account, both analytically and numerically. Its details are shown to be of importance when particles are allowed out of the domain. MITNS reproduces the asymptotic expansion results for a small parameter $δ<<<1$. For $\sqrt{m_elm_i} ~ δ << 1$ , a slightly different regime, dominated by viscosity-induced transport of ions, is found numerically and analytically. Finally, this verification supports the use of this code for more complex time-dependent calculations in the future. Additionally, we derive the angular velocity profile of each species due to radial particle and charge fluxes of various strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fine Root and Soil Organic Carbon Depth Distributions are Inversely Related Across Fertility and Rainfall Gradients in Lowland Tropical Forests

Humid tropical forests contain some of the largest soil organic carbon (SOC) stocks on Earth. Much of this SOC is in subsoil, yet variation in the distribution of SOC through the soil profile remains poorly characterized across tropical forests. We used a correlative approach to quantify relationships among depth distributions of SOC, fine root biomass, nutrients and texture to 1 m depths across 43 lowland tropical forests in Panama. The sites span rainfall and soil fertility gradients, and these are largely uncorrelated for these sites. We used fitted b parameters to characterize depth distributions, where b is a numerical index based on an asymptotic relationship, such that larger b values indicate greater concentrations of root biomass or SOC at depth in the profile. Root b values ranged from 0.82 to 0.95 and were best predicted by soil pH and extractable potassium (K) stocks. For example, the three most acidic (pH < 4) and Kpoor (< 20 g K m -2 ) soils contained 76 ± 5% of fine root biomass from 0 to 10 cm depth, while the three least acidic (pH > 6.0) and most K-rich (> 50 g K m -2 ) soils contained only 41 ± 9% of fine root biomass at this depth. Root b and SOC b values were inversely related, such that a large fine root biomass in surface soils corresponded to large SOC stocks in subsoils (50–100 cm). SOC b values were best predicted by soil pH and base cation stocks, with the three most base-poor soils containing 34 ± 8% of SOC from 50 to 100 cm depth, and the three most base-rich soils containing just 9 ± 2% of SOC at this depth. Nutrient depth distributions were not related to Root b or SOC b values. Here, these data show that large surface root biomass stocks are associated with large subsoil C stocks in strongly weathered tropical soils. Further studies are required to evaluate why this occurs, and whether changes in surface root biomass, as may occur with global change, could in turn influence SOC storage in tropical forest subsoils.

54 ENVIRONMENTAL SCIENCES↗

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗

Radiative heating rate profiles over the southeast Atlantic Ocean during the 2016 and 2017 biomass burning seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96° S, 14.35° W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d -1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

54 ENVIRONMENTAL SCIENCES↗

Effect of negative triangularity on tearing mode stability in tokamak plasmas

Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beam-Based Diagnostics of Electric Guide Fields and Lattice Parameters for Run-1 of the Muon g-2 Storage Ring at Fermilab

A portion of the Muon g-2 Storage Ring electric system, which provides vertical beam focusing, exhibited an unexpected time dependence that produced a characteristic evolution of the stored beam during Run-1 of the Muon g-2 Experiment at Fermilab (E989). A method to reconstruct the Run-1 electric guide fields has been developed, which is based on a numerical model of the muon storage ring and optimization algorithms supported by COSY INFINITY. This method takes beam profile measurements from the Muon g-2 straw tracking detectors as input, and it produces a full reconstruction of the time-dependent fields. The fields can then be used for the reproduction of detailed beam tracking simulations and the calculation of ring lattice parameters for acceptance studies and systematic error evaluations.

43 PARTICLE ACCELERATORS↗

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗

A2Cloud‐RF : A random forest based statistical framework to guide resource selection for high‐performance scientific computing on the cloud

Summary This article proposes a random‐forest based A2Cloud framework to match scientific applications with Cloud providers and their instances for high performance. The framework leverages four engines for this task: PERF engine, Cloud trace engine, A2Cloud‐ext engine, and the random forest classifier (RFC) engine. The PERF engine profiles the application to obtain performance characteristics, including the number of single‐precision (SP) floating‐point operations (FLOPs), double‐precision (DP) FLOPs, x87 operations, memory accesses, and disk accesses. The Cloud trace engine obtains the corresponding performance characteristics of the selected Cloud instances including: SP floating point operations per second (FLOPS), DP FLOPS, x87 operations per second, memory bandwidth, and disk bandwidth. The A2Cloud‐ext engine uses the application and Cloud instance characteristics to generate objective scores that represent the application‐to‐Cloud match. The RFC engine uses these objective scores to generate two types of random forests to assist users with rapid analysis: application‐specific random forests (ARF) and application‐class based random forests. The ARF consider only the input application's characteristics to generate a random forest and provide numerical ratings to the selected Cloud instances. To generate the application‐class based random forests, the RFC engine downloads the application profiles and scores of previously tested applications that perform similar to the input application. Using these data, the RFC engine creates a random forest for instance recommendation. We exhaustively test this framework using eight real‐world applications across 12 instances from different Cloud providers. Our tests show significant statistical agreement between the instance ratings given by the framework and the ratings obtained via actual Cloud executions.

Samuel, David↗

Integrated modeling of alpha particle losses in JET DT plasmas

Alpha particle confinement is crucial for sustaining burning plasmas and designing future reactor concepts. Along with classical/prompt losses, various magnetohydrodynamic instabilities can lead to wave-particle interactions which can transport alpha particles outward from the plasma. This can result in a reduction in plasma heating/performance, and, at worst, damage in-vessel components. Joint European Torus’s recent deuterium–tritium campaigns in 2021–2023 have produced numerous alpha particle loss measurements with its scintillator probe and Faraday cup array fast ion loss detectors as discussed in Bonofiglo et al (2024 Nucl. Fusion 64 096038). This paper will report on integrated energetic particle transport modeling in support of those measurements. The modeling is accomplished with the TRANSP and ORBIT-kick codes with the use of recently developed reduced models which calculate mode structure, amplitude, and the evolving dynamics. When possible, constraints and comparisons to experiment are conducted. Case studies are performed on a variety of magnetohydrodynamic activity, including: fishbones, tearing modes (TMs), and sawtooth crashes. Additionally, a special case of an alpha-driven toroidal Alfvén eigenmode is briefly discussed, where modeling showed marginally weak alpha losses and was unable to support experimental observations. Coupled effects between a TM and toroidal field ripple are presented and were unable to replicate the observations in lost particle pitch but did duplicate the localized flattening of the measured neutron profile. Additional modeling results compare the magnitude of losses and energy/velocity-space sensitivities against experimental observations/measurements for each scenario. This work corroborates numeric alpha transport modeling while also identifying model deficiencies. While this report details alpha transport, it also presents open issues for discussion in assessing the validity of our numerical models towards burning plasmas.

Bonofiglo, Phillip J. [Princeton Plasma Physics La↗

Scaling Turbulent Combustion Fields in Explosions

We considered the topic of explosions from spherical high-explosive (HE) charges. We studied how the turbulent combustion fields scale. On the basis of theories of dimensional analysis by Bridgman and similarity theories of Sedov and Barenblatt, we found that all fields scaled with the explosion length scale r0. This included the blast wave, the mean and root mean squared (RMS) profiles of thermodynamic variables, combustion variables, velocities, vorticity, and turbulent Reynolds stresses. This was a consequence of the formulation of the problem and our numerical method, which both satisfied the similarity conditions of Sedov. We performed numerical simulations of 1 g charges and 1 kg charges; the solutions were identical (within roundoff error) when plotted in scaled variables. We also explored scaling laws related to three-phase pyrotechnic explosions. We show that although the scaling formally broke down, the fireball still essentially scaled with the explosion length scale r0. However, the discrete Lagrange particles (DLP) (phase 2) and the heterogeneous continuum model (HCM) of the DLP wakes (phase 3) did not scale with r0, and mean and RMS profiles could differ by a factor of 10 in some regions. This was because the DLP particles and wakes introduced an additional scale that broke the similarity conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗