Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical stability”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Twenty-Year Analysis of Winds in California for Offshore Wind Energy Production Using WRF v4.1.2

Offshore wind resource characterization in the United States relies heavily on simulated winds from numerical weather prediction (NWP) models, given the lack of hub-height observations offshore. One such NWP data set used extensively by U.S. stakeholders is the Wind Integration National Dataset (WIND) Toolkit, a 7-year time-series data set produced in 2013 by the National Renewable Energy Laboratory. In this study, we present an update to that data set for offshore California that leverages recent advancements in NWP modeling capabilities and extends the period of record to a full 20 years. The data set predicts a significantly larger wind resource (0.25–1.75 m s-1 stronger), including in three Call Areas that the Bureau of Ocean Energy Management is considering for commercial activity. We conduct a set of yearlong simulations to study factors that contribute to this increase in the modeled wind resource. The largest impact arises from a change in the planetary boundary layer parameterization from the Yonsei University scheme to the Mellor-Yamada-Nakanishi-Niino scheme and their diverging wind profiles under stable stratification. Additionally, we conduct a refined wind resource assessment at the three Call Areas, characterizing distributions of wind speed, shear, veer, stability, frequency of wind droughts, and power production. We find that, depending on the attribute, the new data set can show substantial disagreement with the WIND Toolkit, thereby driving important changes in predicted power.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2023

Microphysical Modeling of Fault Slip and Stability Transition in Hydrothermal Conditions

Abstract Field and laboratory observations indicate that the frictional behaviors of faults depend on hydrothermal conditions. We extend the microphysical Chen‐Niemeijer‐Spiers (CNS) model to hydrothermal conditions by using the observed temperature variation of indentation hardness to infer the temperature dependence of a microphysical parameter . This parameter is assumed constant in previous versions of the CNS model. A simple spring‐slider system is used to simulate the fault system and investigate the steady‐state frictional behaviors of wet granite gouges. Our numerical results quantitatively reproduce experimental data showing the frictional‐plastic transition. The results also describe the transition from velocity‐strengthening at low temperatures (<160°C), to velocity‐weakening at intermediate temperatures (160°C–370°C), then back to velocity‐strengthening at high temperatures (>370°C). In our extended CNS model, these results suggest that the dominant shear deformation mechanism does transition from frictional granular flow to fully plastic creep with increasing temperature.

58 GEOSCIENCES↗
DOE OSTI2022

Towards fast-charging high-energy lithium-ion batteries: From nano- to micro-structuring perspectives

Electric vehicles (EVs) have been playing an indispensable role in reducing greenhouse gas emissions for our modern society. However, current EVs are difficult to meet people’s diverse travel needs, especially in long endurance and fast-charging capacities. At the heart of this issue is the physicochemical limit of current lithium-ion batteries (LIBs), which are the core parts for powering the vehicles. Hence, LIBs with simultaneous high energy and power are critically required to further promote the development of EVs. Here, in this review, we first summarize the key electrochemical processes in electrochemical reactions which lead to the corresponding overpotentials in or between multiple battery components. Furthermore, numerical simulations are employed to quantitatively analyze the effects of versatile electrode parameters on electrochemical properties in high-energy NMC811//graphite systems. On the basis of the in-depth understandings from simulation, recent experimental efforts on designing electroactive materials and electrode architectures across multiple length scales are discussed. Among them, nano-structuring can promote local mass transport and stabilize the interfaces at the particle level, while micro-structuring can establish efficient pathways for charge carriers at the electrode level. Finally, we conclude that a tight feedback loop among structure engineering, characterization and simulation should be followed to speed up the understanding of the deficiencies existing in current electrode designs as well as point out the possible electrode optimization routes for next-generation fast-charging LIBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators

Abstract The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using numerical optimization methods. It is found that while the near-axis expansion is able to replicate the growth rates, real frequencies and perpendicular wave-vector at the inner core (both using simplified dispersion relations and first-principle gyrokinetic simulations), it leads to an overestimation of the growth rate at larger radii. An approximate analytic solution of the ITG dispersion relation for the non-resonant limit suggests growth rates could be systematically higher in quasi-axisymmetric (QA) configurations compared to quasi-helically (QH) symmetric ones. However except for very close to the axis, linear gyrokinetic simulations do not show systematic differences between QA and QH configurations.

Physics↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2022

Implications of Shear and Thermal Stratification on Wind Turbine Tip-Vortex Stability

The interaction between wind turbines in a wind farm through their wakes is a phenomenon that has been studied for decades and is still relevant today. Turbines clustered together in arrays will often operate in the wake of other upstream turbines which may lead to significant power losses and fatigue loads. For modern large-scale wind turbines, the mean shear velocity profile and thermal stratification are major components of the atmospheric boundary layer so it is important to understand their impact on near-wake development. Additionally, veer is present due to the rotation of the Earth. The impact of shear, thermal stratification and veer on the stable wake length of turbines with a dynamic control strategy is studied numerically in this work using a suite of highly resolved large-eddy simulations. Instantaneous flow fields are extracted from the simulations and used to conduct proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip-vortex instability mechanisms. Our findings show that the dynamic pitch control scheme is able to shorten the stable wake length to about 1.5R in uniform flow. Shear can significantly affect the break up of wind turbine tip-vortices as well as the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near-wake field. Veer causes the wake boundary to skew but has a limited impact on the wake length.

length thermal stratification↗
DOE OSTI2024

Operational Analytics Studies for ATLAS Distributed Computing: Data Popularity Forecast and Utilization of the WLCG Centers

Operational analytics is the direction of research related to the analysis of the current state of computing processes and the prediction of future states in order to anticipate imbalances and take timely measures to stabilize a complex system. There are two relevant areas in ATLAS Distributed Computing that are currently the focus of studies: user physics analysis including the forecast of popularity of data samples among users, and evaluating WLCG centers for their readiness to process user analysis payloads. Studying these areas is challenging due to the complexity involved, as it requires a comprehensive understanding of numerous boundary conditions typically found in large-scale distributed computing infrastructures. Forecasts of data popularity are problematic without the categorization of user tasks by their types (data transformation or physics analysis), which do not always appear on the surface but may induce noise, which introduces significant distortions for predictive analysis. Evaluating the WLCG resources by their analysis workloads is also a challenging task as it is necessary to find a balance between the workload of the resource, its performance, the waiting time for jobs on it, as well as the volume of jobs that it processes. This is especially difficult in a heterogeneous computing environment, where legacy resources are used along with modern high-performance machines. We will look at these areas of research in detail and discuss what tools and methods are used in our work, demonstrating results already obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Neural network based fast prediction of β N limits in HL-2M

Artificial neural networks (NNs) are trained, based on the numerical database, to predict the no-wall and ideal-wall β N limits, due to onset of the n = 1 ( n is the toroidal mode number) ideal external kink instability, for the HL-2M tokamak. The database is constructed by toroidal computations utilizing both the equilibrium code CHEASE (Lütjens et al 1992 Comput. Phys. Commun. 69 287) and the stability code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The stability results show that (1) the plasma elongation generally enhances both β N limits, for either positive or negative triangularity plasmas; (2) the effect is more pronounced for positive triangularity plasmas; (3) the computed no-wall β N limit linearly scales with the plasma internal inductance, with the proportionality coefficient ranging between 1 and 5 for HL-2M; (4) the no-wall limit substantially decreases with increasing pressure peaking factor. Furthermore, both the NN model and the convolutional neural network (CNN) model are trained and tested, producing consistent results. The trained NNs predict both the no-wall and ideal-wall limits with as high as 95% accuracy, compared to those directly computed by the stability code. Additional test cases, produced by the Tokamak Simulation Code (Jardin et al 1993 Nucl. Fusion 33 371), also show reasonable performance of the trained NNs, with the relative error being within 10%. The constructed database provides effective references for the future HL-2M operations. The trained NNs can be used as a real-time monitor for disruption prevention in the HL-2M experiments, or serve as part of the integrated modeling tools for ideal kink stability analysis.

Physics↗
DOE OSTI2025

A Predictive Deep-Reinforcement-Learning-Based Connected Automated Vehicle Anticipatory Longitudinal Control in a Mixed Traffic Lane Change Condition

Maintaining safety and efficiency for mixed traffic consisting of connected automated vehicles (CAVs) and human-driven vehicles (HDVs) is an arduous task due to the inherent HDVs’ stochasticity. Especially for longitudinal control, which is the basic function of vehicle automation, prevailing research primarily considers CAV’s car-following control merely the acceleration and deceleration of leading vehicles. However, this approach overlooks the potential disruptions caused by surrounding vehicles executing lane changes, which can significantly impact the control vehicle’s stability and overall safety. Hence, our study introduces a predictive deep reinforcement learning (DRL) longitudinal CAV controller. This innovative approach leverages prediction from a physics-informed neural network as well as the control capability of DRL to better anticipate and mitigate issues arising from lane-changing, enhancing the safety and efficiency of CAVs in such scenarios. Finally, validated by the numerical simulations embedded with the real-world data, the results indicate that the proposed controller significantly enhances the safety and efficiency of CAVs in situations involving lane changes by other vehicles, showcasing its potential as a valuable tool in advancing CAV technology in mixed traffic.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Passive Ultrasonic Deterrents to Reduce Bat Mortality in Wind Farms (Final Report)

The overall objective of this project was to develop active and passive ultrasonic deterrent concepts to mitigate bat mortality at wind turbines. The research was supported by the US Department of Energy under the contract DE-EE0008731. Novel ultrasonic bat deterrents are investigated in this study. The deterrents are based on the idea of aerodynamic whistles wherein flow-acoustic resonance is used to produce high-amplitude tonal sound at desired ultrasonic frequencies. The deterrents can be grouped into active and passive. The active deterrents/whistles are driven by pressurized air supplied by an external source (e.g., an air compressor). The passive deterrents/whistles are designed to be mounted on wind turbine blades and are “powered” by the energy in the air moving past the rotor blades. Multiple active and passive deterrent ideas are proposed and investigated numerically and experimentally. The computational analysis involves solving the compressible unsteady Reynolds-averaged Navier-Stokes (URANS) equations and coupling the near-field aerodynamic solution with an integral method based on the Ffowcs Williams-Hawkings acoustic analogy to predict ultrasound in the farfield. Experiments are performed in the anechoic chamber at ISU (for the active whistles) and in the Stability wind tunnel at Virginia Tech (for the passive whistles).

17 WIND ENERGY↗
DOE OSTI2023

Theoretical Investigation of Delafossite-Cu 2 ZnSnO 4 as a Promising Photovoltaic Absorber

In the quest for efficient and cost-effective photovoltaic absorber materials beyond silicon, considerable attention has been directed toward exploring alternatives. One such material, zincblende-derived Cu2ZnSnS4 (CZTS), has shown promise due to its ideal band gap size and high absorption coefficient. However, challenges such as structural defects and secondary phase formation have hindered its development. In this study, we examine the potential of another compound, Cu2ZnSnO4 (CZTO), with a similar composition to CZTS as a promising alternative. Employing ab initio density function theory (DFT) calculations in combination with an evolutionary structure prediction algorithm, we identify that the crystalline phase of delafossite structure is the most stable among the 900 (meta)stable CZTO. Its thermodynamic stability at room temperature is also confirmed by the molecular dynamics study. Excitingly, this new phase of CZTO displays a direct band gap where the dipole-allowed transition occurs, making it a strong candidate for efficient light absorptions. Furthermore, the estimation of spectroscopic limited maximum efficiency (SLME) directly demonstrates the high potential of delafossite-CZTO as a photovoltaic absorber. Our numerical results suggest that delafossite-CZTO holds promise for future photovoltaic applications.

14 SOLAR ENERGY↗
DOE OSTI2024

Flashback characterization of additively manufactured swirl-stabilized fuel injector with varying surface roughness

This work investigates the effects of as-built surface roughness on the flashback propensity in additively manufactured (AM) swirl-stabilized lean premixed (LPM) fuel injectors. Adoption of AM for rapid prototyping and fabrication of complex fuel-flexible injector designs requires investigating surface roughness effects on the flow and flame stability characteristics of the combustor. Wall roughness increases the near-wall shear, which could alter the boundary layer structure and change the propensity for flame flashback. Accounting for the realistic as-built surface roughness is crucial in carrying out computational modeling and experimental analysis to establish a feedback loop for the precision designing of fuel-flexible injectors. The presented numerical analysis of as-built AM injector is an essential consideration to optimize injector design for aerodynamics and additive manufacturability.

Mohanty, Pratikshya↗
DOE OSTI2021

Unraveling the Combined Effects of Ice and Permafrost on Arctic Delta Morphodynamics

Arctic delta channel dynamics and coastal deposition patterns are uniquely affected by ice and permafrost. Building on previous experiments that documented the individual effects of ice and permafrost, we conducted new numerical modeling experiments to examine the combined influences of ice thickness and permafrost erodibility. We quantified the effects of ice, permafrost, and their combination on channel migration, avulsion, and persistence; onshore and nearshore deposition; and large-scale delta morphology. Both ice and permafrost limited channel mobility, enhanced overbank deposition, and roughened shorelines. Thick ice increased offshore channel incision and sediment delivery. Permafrost stabilized channels and focused deposition that led to enhanced channel levee development. The increased levee deposition limited widespread channel infilling leading to the preservation of abandoned channels. Levee growth was enhanced by overbank flooding due to ice-driven changes in water surface elevations. Furthermore, this flooding led to more extensive levees on deltas with both ice and permafrost compared to deltas with only permafrost. Ice also preserved small channels, leading to size-dependent channel persistence: the presence of ice in the channels reduced the fraction of time small channels remained active relative to larger channels. Resistant permafrost damped this effect, as smaller and more distal channels were active earlier in the flood on deltas with resistant permafrost likely due to the smaller number of channels on permafrost-laden deltas. Finally, our results suggest that the loss of both ice and permafrost will increase channel mobility, decrease overbank flooding, increase shoreline progradation rates, and decrease offshore transport of riverine fluxes.

58 GEOSCIENCES↗
DOE OSTI2024

Bayesian stability and force modeling for uncertain machining processes

Accurately simulating machining operations requires knowledge of the cutting force model and system frequency response. However, this data is collected using specialized instruments in an ex-situ manner. Bayesian statistical methods instead learn the system parameters using cutting test data, but to date, these approaches have only considered milling stability. This paper presents a physics-based Bayesian framework which incorporates both spindle power and milling stability. Initial probabilistic descriptions of the system parameters are propagated through a set of physics functions to form probabilistic predictions about the milling process. The system parameters are then updated using automatically selected cutting tests to reduce parameter uncertainty and identify more productive cutting conditions, where spindle power measurements are used to learn the cutting force model. The framework is demonstrated through both numerical and experimental case studies. Results show that the approach accurately identifies both the system natural frequency and cutting force model.

42 ENGINEERING↗
DOE OSTI2020

Widely Tunable Quantum Phase Transition from Moore-Read to Composite Fermi Liquid in Bilayer Graphene

Here, we develop a proposal to realize a widely tunable and clean quantum phase transition in bilayer graphene between two paradigmatic fractionalized phases of matter: the Moore-Read fractional quantum Hall state and the composite Fermi liquid metal. This transition can be realized at total fillings $\textit{ν}$ = ±3 + 1/2 and the critical point can be controllably accessed by tuning either the interlayer electric bias or the perpendicular magnetic field values over a wide range of parameters. We study the transition numerically within a model that contains all leading single particle corrections to the band structure of bilayer graphene and includes the fluctuations between the $\textit{n}$ = 0 and $\textit{n}$ = 1 cyclotron orbitals of its zeroth Landau level to delineate the most favorable region of parameters to experimentally access this unconventional critical point. We also find evidence for a new anisotropic gapless phase stabilized near the level crossing of $\textit{n}$ = 0/1 orbits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗