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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Runaway electron plateau current profile reconstruction from synchrotron imaging and Ar-II line polarization angle measurements in DIII-D

Abstract Current profile reconstructions are obtained for high current ( I p ≃ 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Simplified Numerical Approach to Characterize the Thermal Response of a Moving Bed Solar Reactor

Abstract Concentrated solar thermochemical storage in the form of a zero-emission fuel is a promising option to produce long-duration energy storage. Solar fuel is produced using a cavity reactor that captures concentrated solar radiation from a solar field of heliostats. In this paper, heat transfer model of a tubular plug-flow reactor designed and manufactured for a solar fuel production is presented. Experimental data collected from a fixed bed tubular reactor testing are used for model comparison. The system consists of an externally heated tube with counter-current flowing gas and moving solid particles as the heated media. The proposed model simulates the dynamic behavior of temperature profiles of the tube wall, gas, and particles under various gas flowrates and residence times. The heat transfer between gas–wall, solid particle–wall, and gas–solid particle is numerically studied. The model results are compared with the results of experiments done using a 4 kW furnace with a 150 mm heating zone surrounding a horizontal alumina tube (reactor) with 50.8 mm outer diameter and thickness of 3.175 mm. Solid fixed particles of magnesium manganese oxide (MgMn2O4) with the size of 1 mm are packed within the length of 250 mm at the center of the tube length. Simulation results are assessed with respect to fixed bed experimental data for four different gas flowrates, namely, 5, 10, 15, and 20 standard liters per minute of air, and furnace temperatures in the range of 200–1200 °C. The simulation results showed good agreement with maximum steady state error that is less than 6% of those obtained from the experiments for all runs. The proposed model can be implemented as a low-order physical model for the control of temperature inside plug-flow reactors for thermochemical energy storage applications.

Engineering↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2026

photoD with Rubin ’s Data Preview 1: First stellar photometric distances and faint blue star deficits

Aims. We investigate the utility of Rubin’s Data Preview 1 (DP1) for estimating stellar number density profiles across the Milky Way halo. Methods. We used stellar broad-band near-UV to near-IR ugrizy photometry released in Rubin’s DP1 to estimate distance and metallicity for blue main sequence stars brighter than r = 24 in three ~1.1 sq. deg. fields at southern Galactic latitudes. Results. Compared to TRILEGAL simulations of the Galaxy’s stellar content, we found a likely deficit of blue main sequence turn-off stars with 22 < r < 24. We interpreted this discrepancy as a signature of a steeper halo number density profile at galactocentric distances 10–50 kpc than the canonical ~1/r 3 profile assumed in TRILEGAL simulations. Conclusions. This interpretation is consistent with earlier suggestions based on observations of more luminous, but much less numerous, evolved stellar populations, along with a few pencil beam surveys of blue main sequence stars in the northern sky. These results bode well for the future Galactic halo exploration with Rubin’s Legacy Survey of Space and Time (LSST).

Galaxy: fundamental parameters↗
DOE OSTI2024

On the subject of large-scale pool fires and turbulent boundary layer interactions

The role to which a realistic inflow turbulent boundary layer (TBL) influences transient and mean large-scale pool fire quantities of interest (QoIs) is numerically investigated. High-fidelity, low-Mach large-eddy simulations that activate low-dissipation, unstructured numerics are conducted using an unsteady flamelet combustion modeling approach with mutiphysics coupling to soot and participating media radiation transport. Three inlet profile configurations are exercised for a large-scale, high-aspect rectangular pool that is oriented perpendicular to the flow direction: a time-varying, TBL inflow profile obtained from a periodic precursor simulation, the time-mean of the transient TBL, and a steady power-law inflow profile that replicates the mean TBL crosswind velocity of 10.0 m/s at a vertical height of 10 m. Results include both qualitative transient flame evolution and quantitative flame shape with ground-level temperature and convective/radiative heat flux profiles. While transient fire events, which are driven by burst-sweep TBL coupling, such as blow-off and reattachment are vastly different in the TBL case (contributing to increased root mean square QoI fluctuation prediction and disparate flame lengths), mean surface QoI magnitudes are similar. Quadrant analysis demonstrates that the TBL configuration modifies burst-sweep phenomena at windward pool locations, while leeward recovery is found. In conclusion, positive fluctuations of convective heat flux correlate with fast moving fluid away from the pool surface due to intermittent combustion events.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗
DOE OSTI2025

Exceptional points in a passive strip waveguide

Abstract Exceptional points (EPs) in non‐Hermitian systems have attracted significant interest due to their unique behaviors, including novel wave propagation and radiation. While EPs have been explored in various photonic systems, their integration into standard photonic platforms can expand their applicability to broader technological domains. In this work, we propose and experimentally demonstrate EPs in an integrated photonic strip waveguide configuration, exhibiting unique deep wave penetration and uniform‐intensity radiation profiles. By introducing the second‐order grating on one side of the waveguide, forward and backward propagating modes are coupled both directly through second‐order coupling and indirectly through first‐order coupling via a radiative intermediate mode. To describe the EP behavior in a strip configuration, we introduce modified coupled‐mode equations that account for both transverse and longitudinal components. These coupled‐mode formulas reveal the formation of EPs in bandgap closure, achieved by numerically optimizing the grating’s duty cycle to manipulate the first‐ and second‐order couplings simultaneously. Experimental observations, consistent with simulations, confirm the EP behavior, with symmetric transmission spectra and constant radiation profiles at the EP wavelength, in contrast to conventional exponential decay observed at detuned wavelengths. These results demonstrate the realization of EPs in a widely applicable strip waveguide configuration, paving the way for advanced EP applications in nonlinear and ultrafast photonics, as well as advanced sensing technologies.

Materials Science↗
DOE OSTI2021

Radiosonde Observations of Environments Supporting Deep Moist Convection Initiation during RELAMPAGO-CACTI

The Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) and Cloud, Aerosol, and Complex Terrain Interactions (CACTI) projects deployed a high-spatiotemporal-resolution radiosonde network to examine environments supporting deep convection in the complex terrain of central Argentina. This study aims to characterize atmospheric profiles most representative of the near-cloud environment (in time and space) to identify the mesoscale ingredients affecting storm initiation and growth. Spatiotemporal autocorrelation analysis of the soundings reveals that there is considerable environmental heterogeneity, with boundary layer thermodynamic and kinematic fields becoming statistically uncorrelated on scales of 1–2 h and 30 km. Using this as guidance, we examine a variety of environmental parameters derived from soundings collected within close proximity (30 km in space and 30 min in time) of 44 events over 9 days where the atmosphere either: 1) supported the initiation of sustained precipitating convection, 2) yielded weak and short-lived precipitating convection, or 3) produced no precipitating convection in disagreement with numerical forecasts from convection-allowing models (i.e., Null events). There are large statistical differences between the Null event environments and those supporting any convective precipitation. Null event profiles contained larger convective available potential energy, but had low free-tropospheric relative humidity, higher freezing levels, and evidence of limited horizontal convergence near the terrain at low levels that likely suppressed deep convective growth. Here, we also present evidence from the radiosonde and satellite measurements that flow–terrain interactions may yield gravity wave activity that affects CI outcome.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multi-modal Energy-optimal Trip Scheduling in Real-time (METS-R) for Transportation Hubs (Final Report)

This report summarizes the work performed under the award number EE0008524. The project develops the Multi-modal Energy-optimal Trip Scheduling in Real-time (METS-R) platform as the next-generation transportation solution based on autonomous electric vehicles (AEV) serving passenger trips from and to urban transportation hubs, to substantially reduce transportation energy consumption. Extensive data collection and analyses were first conducted to understand the demand patterns and energy consumption of hub-based on-road trips. Then, a data-driven framework that consists of an analytical module and a simulation module was proposed. For the analytical module, five planning + operation tools were developed to support the planning and energy-efficient operations of urban AEV services: the charging station planning that robotically allocates charging supplies based on the stationary charging demand distribution; the transit planning and demand adaptive scheduling model that efficiently generates\ candidate transit routes from hubs to other places and dynamically adjusts the transit time table to fit the current demand; the online energy-efficient routing that learns the energy-optimal paths from observations of link-level energy consumption in real-time; the hub-based ridesharing that matches trip requests together with account for the uncertainty of future trip demand and vehicle supply; and finally, the integrated demand prediction and anomaly detection pipeline that leverages the flight/train time table and support other planning/operation tools. To demonstrate the performance of these tools, a scalable high-performance agent-based simulator was built. We divided the urban space into multiple service zones where each zone was considered as an agent for passenger generation and vehicle charging. Two types of AEV agents were coded to model two types of mobility services: AEV taxi and AEV transit. For the AEV taxi, the team implemented the functions of pickup/drop-off passengers, energy-efficient routing, ridesharing, fleet rebalancing, and recharging. For the AEV bus, the team implemented the functions of demand-adaptive route scheduling, passenger boarding, and recharging. A high-performance computing framework was introduced to receive various profiling information (such as link energy updates, vehicle speed) from the simulator instances and communicate the operational commands back to the instances. The numerical experiments show that each of the proposed operational algorithms can reduce energy consumption and improve system efficiency. Furthermore, there exists the need to collectively consider multiple planning + operational strategies as multiple strategies can influence each other in terms of performance impacts. Recommendations for future work related to AEV planning and simulation are discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Additive Manufacture of Multi-Functional, Large-Scale Components

The poster summarizes the main results for a LDRD focused on developing a arc-based additive manufacturing (AM) system capable of fabricating large components with a composition gradient otherwise known as functionally graded materials (FGMs). FGMs are capable of addressing local extreme conditions in an optimal and cost-effective manner. The AM system (referred to as multi-Wire Arc Additive Manufacturing or m-WAAM) consisted of modifying a standard automated gas tungsten arc welding system to deliver multiple welding wires to the weld pool at controlled rates to deposit desired alloy compositions at critical locations. Numerous shapes and composition gradients were fabricated with the system and the effect of various deposition parameters on the nature of the composition profile were evaluated. Additionally, heat treatment of the FGM was performed to develop graded microstructures which should enhance dynamic impact performance of inexpensive steels. The project highlights the potential of m-WAAM to efficiently fabricate FGMs. However, this project only explored the science behind the technique and additional development work will be required to develop a fully-functional system that can build more complex components.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine Straightness Error Measurement Based on Optical Fiber Fabry–Pérot Interferometer Monitoring Technique

Abstract In this research, we propose an Error Separation Technique (EST) based on optical fiber sensors for on-machine straightness error measurement. Two fiber optic Fabry–Pérot interferometers have been developed serving as two displacement sensors. The displacement distance is computed according to the reflected spectrum from interferometers, which can achieve a sub-micrometer resolution. The two-point method has been employed to separate the straightness error of the slides and the profile error of a fine-polished standard block. The spacing distance between two interferometers is determined by the diameter of optical fibers so that the EST's resolution has the potential to reach the sub-millimeter scale. In the experiment, the straightness error has been measured on a commercially available computer numerical control machine tool, and the measurement has been conducted on its x-axis. The spacing distance between two optical fiber sensors is 1.5 mm which equals the EST's resolution along the machine tool's x-axis. The separated profile error of the measured standard block is around 30 µm which has been verified by a high precision Coordinate Measurement Machine (CMM). The magnitude of the separated straightness error is around 40 µm. This technique is flexible and simple to be conducted, which can contribute to the micro-machine tool calibration and other straightness error applications.

Engineering↗
DOE OSTI2026

Nonlinear saturation of ballooning modes in stellarators

Ballooning mode saturation is investigated in realistic stellarator configurations using the flux tube approach of Ham et al (2018 Plasma Phys. Control. Fusion 60 075017), Ham et al (2016 Phys. Rev. Lett. 116 235001). The method is adapted to account for the lack of exact force balance in stellarator equilibrium solvers that assume existence of nested flux surfaces. A variational approach for calculating flux tube energy is developed to overcome this force error problem in stellarator numerical equilibria. Saturated (equilibrium) flux tube states that cross 10%–20% of the plasma minor radius are shown to exist for linearly ballooning unstable profiles. It is shown that several features of the displaced flux tube structure in a full nonlinear MHD simulation of Wendelstein 7X are reproduced by our model. Saturated states are found in a compact stellarator equilibrium close but below the marginal ballooning linear instability, i.e. the unperturbed equilibrium is metastable. This suggests that edge-localized-mode-like explosive MHD behavior may be possible in stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Distributed Data-Driven Optimization for Voltage Regulation in Distribution Systems

Here, this paper proposes a distributed data-driven optimization framework for voltage regulation in distribution systems. The recursive kernel regression and alternating direction method of multipliers (ADMM) are selected to cover the system learning and distributed optimization tasks. The proposed distributed data-driven framework is capable of having a rapid response to system or load changes while considering the operation optimality. Besides, the distributed algorithm parallels the computation tasks and reduces the computational expense of a single agent. To validate the performance of the proposed method, a hypothetical 7-Bus system and the IEEE 123-Bus system are selected to show the effectiveness of the proposed data-driven framework. According to the numerical study results, the proposed method offers great flexibility for selecting customized kernel models for different regions and can effectively improve the system voltage profile in a distributed manner.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Techno-economic feasibility of borehole thermal energy storage system connected to geothermal heat pumps for seasonal heating load of two buildings in Fairbanks, Alaska

Borehole thermal energy storage (BTES) is an effective solution for managing imbalanced heating and cooling loads in cold regions. This study evaluated the long-term feasibility of a BTES system in the Fairbanks area, Alaska, through building energy modeling, resource characterization, and numerical modeling. The system was designed to store waste heat from a nearby coal power plant during summer and provide thermal energy during winter to geothermal heat pumps supplying heating loads in two buildings. Heating load profiles were modeled for the buildings using EnergyPlus, and the results indicated the annual heating load was 5.6 times greater than the cooling load. 40 borehole heat exchangers were pre-designed approximately 100 m away from the two buildings in terms of land availability and regulatorily optimized depth. The 20-year performance of the designed BTES system under two operational scenarios—one with a 5-year preheating period and one without—was numerically modeled using subsurface temperature and properties characterized through the literature review and thermal response tests. Both scenarios demonstrated that the BTES has the capacity to fully cover the heating loads in the two buildings throughout the 20-year lifetime. Production temperatures at central wells were 33 % higher on average than at outer wells in both scenarios. The 5-year preheating period increased subsurface and extraction temperatures, and correspondingly annual average and total thermal energy production was higher for 8 years than in the scenario without the preheating period. These results highlight the long-term reliability and sustainability of the BTES system in meeting heating demands over its lifetime, with the preheating period offering potential performance improvements. Implementing the BTES system in cold regions with high heating demand, such as Fairbanks, Alaska, could provide a long-term, sustainable energy solution for managing imbalanced heating and cooling loads.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗