Search NASA⌕ Search

SEARCH · Search NASA

Results for “Utility and Grid Operations”

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 451 records · Page 25

Feasibility Analyses for the Microgrid of the Mountain in the Cordillera Central Region of Puerto Rico

Distributed energy resource (DER) development can benefit communities through improvement of access to reliable electricity and a resilient energy system. Many place-specific considerations must be accounted for when designing electrical systems with co-located generation and loads. This paper presents an overview of four studies focusing on different aspects for regional DER development in the Cordillera Central region of Puerto Rico, including a grid and load stability study with the introduction of solar energy generation, microgrid design based on the solar resource in the area, resilience considerations during power outage scenarios, and qualitative risk aspects for components in the microgrid. Additionally, a conceptual substation microgrid analysis is conducted using the Microgrid Design Toolkit (MDT) to evaluate trade-offs between cost and energy availability, providing insights into optimal design strategies. A microgrid resilience assessment is performed using the Resilient Nodal Cluster Analysis Tool (ReNCAT) to identify critical loads and evaluates their accessibility during outages, emphasizing the importance of community needs. Finally, a qualitative risk assessment utilizing Hazard and Operability Analysis (HAZOP) methods highlights potential hazards and mitigations for microgrid components, ensuring robust system design and operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

WE-Validate: An Open-Source Framework For Wind Power Validation

Grid operators rely on historical weather time series at existing and planned wind power plants to make informed decisions when planning for a future power grid with very high penetration of renewable power. While synthetic wind power time series have been developed based on historical weather models, their validation with actual power production data remains complex due to variations in modeling practices and methodologies. This paper introduces the WE-Validate framework, originally designed for wind speed validation and now enhanced for wind power validation with a graphical user interface to support users with minimal programming experience. Validation of wind power with WE-Validate is based on robust metrics consisting of RMSE, centered RMSE, average bias, average percent bias, mean absolute error, mean absolute percent error, cross correlation, and calculation of ramping magnitude, rate, and duration. This paper showcases WE-Validate with validation of synthetically derived power for a wind plant in Washington state for one month in 2018. Validation of the synthetic power from two comparison data sets compared with observations shows both comparison series have strong correlation with observed across weekly and monthly aggregations while suffering from persistent negative bias. The suite of metrics within WE-Validate facilitates immediate insight into the utility of the comparison data sets through compression across multiple axes. This user-friendly, open-source tool can be extended beyond wind power, making it a valuable resource for system planners and operators in different domains.

Moncheur de Rieudotte, Malcolm P.↗

Virtual Power Plant Profiles and Inventory

Increased adoption of distributed energy resources (DERs), growing quantities of renewable energy on the grid, and load growth from buildings, industry, and data centers are driving the need for increased demand flexibility and solutions to reliably and cost-effectively balance electricity supply and demand. In order to deploy DERs on a similar scale as conventional supply-side resources, utilities and state decision-makers are exploring a range of programmatic approaches and market designs, including Virtual Power Plants (VPP). Based on interviews with more than twenty subject matter experts on VPPs, the Insights into Scaling Virtual Power Plants report and appendix outlines actions that utilities and regulators can take to expand VPPs, organized around successful VPP designs, deployment, and operation. The companion report, Virtual Power Plant Profiles and Inventory describes three conceptual profiles that can be used as a starting point to facilitate clear discussions among regulators, utilities, solution providers, and stakeholders about various VPP strategies. It also includes details on more than 790 demand response programs and net energy metering solar plus battery programs, and 180 VPP programs. The Inventory is provided below as a separate file.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Responsible Adoption of Artificial Intelligence (AI) in Electric Grid Operations

The future of the grid will be powered by AI—or undermined by it. Artificial intelligence is rapidly reshaping grid operations, improving fault detection, forecasting accuracy, and real-time optimization. As AI systems move closer to operational decision loops, however, they introduce new consequence pathways: expanded attack surfaces, model integrity risks, regulatory exposure, and human-automation challenges. This talk presents a consequence-driven framework for deploying AI responsibly in the electric grid. Attendees will gain practical strategies to strengthen resilience, boost reliability, and deploy AI securely — ensuring the grid of the future is not only smarter but safer.

25 - ENERGY STORAGE↗

A Review on Protection Challenges on Transmission Lines Connected to Inverter Based Resources

The increased penetration of inverter-based resources (IBRs) connected to transmission grids is challenging how traditional electric power systems operate. Such systems were originally designed based on the dominant presence of synchronous generators. Some main challenges include the increased concerns about system stability due to the reduced system inertia caused by IBRs and the reduced protection reliability caused by the non-deterministic fault current response of IBRs. The protection relays’ algorithms have been designed and validated based on well-known fault current characteristics of synchronous generations, while the short-circuit responses of IBRs heavily depend on their vendor-specific control schemes, which are mainly set to protect inverter semiconductor switches and DC components. The magnitude of short current that IBRs contribute is typically between 1.1–1.5 times the rated current that may behave unstably and incoherently with the voltages. The control scheme of the IBRs also controls the active (P) and reactive (Q) support during system voltage sag and swell conditions. These characteristics challenge not only protective relay schemes but also existing short-circuit software tools to correctly model IBRs for fault studies and protection coordination analyses. Utilities, university researchers, research institutes, and relay manufacturers have performed studies to investigate different aspects of high IBR penetration. Almost all studies support 1) creating international and national standards to regulate the IBRs, 2) revising engineering tools to address shortcomings, and 3) adjusting the protective relay settings to resolve protection problems. This paper summarizes the studies and research conducted in recent years to quantify the impact of IBR penetration on transmission line protection and to identify potential solutions to improve the dependability and reliability of grid protection systems.

14 SOLAR ENERGY↗

Modeling Framework for the Assessment of a Sustainable Hydrogen Production and Supply Chain Network in California

The cost-effective and sustainable deployment of hydrogen supply and demand networks, especially in large economic regions like California, can be challenging considering the spatial-temporal availability and variability of the different actors across the network such as production processes, distribution modes, and end-users. In this presentation, we will provide an overview and demonstration of a modeling framework used to assess the environmental, economic, and human health impacts of plausible hydrogen production and supply chain networks in California. Scenarios focus on green hydrogen production pathways using water electrolysis and biomass gasification. End-use applications included in the model are transit, medium and heavy-duty trucking, port authorities, and power and aviation companies that currently consume natural gas, diesel, and aviation fuel for their day-to-day operation. Representative locations for hydrogen production and end-use are based on recent projections of the hydrogen economy in California. All mass and energy flows, as well as estimated emissions, are based on H2A process model designs and projections of technology performance, literature review, and LBNL process, economic and life cycle modeling, and not on company data for the sake of this presentation. Human health impacts are included following methodologies developed for the University of California Irvine HyDeal project. Life cycle phases associated with hydrogen production include feedstock preparation (water and biomass), energy production and consumption (renewable, grid, and combination of renewable and grid electricity), maintenance (chemical utilization in electrolysis and natural gas combustion in gasification), carbon sequestration, hydrogen storage (compression and liquefaction), and distribution (truck and pipeline). We apply the framework utilizing California specific emission factors, financial data, and human health damages and explore the impact of network characteristics on results. Example variations include: the inclusion of policy incentives or not, different representations of the electricity grid and source, electrolysis versus gasification versus combinations of both for production, liquefaction versus compression based on producer capacity cutoffs, transportation truck versus pipeline based on existing infrastructure, and ultimate end use. Comparison of these different scenarios can help inform future projects by demonstrating the trade-offs among environmental, economic, and human health impacts. This model, automated in R, is a starting platform upon which new analysis, modeling capabilities, locations, and emission factors can be rapidly tested and integrated.

Zaki, Mohammed Tamim↗

Boundary-Layer Separation Control under Low-Pressure Turbine Airfoil Conditions using Glow-Discharge Plasma Actuators

Modem low-pressure turbines, in general, utilize highly loaded airfoils in an effort to improve efficiency and to lower the number of airfoils needed. Typically, the airfoil boundary layers are turbulent and fully attached at takeoff conditions, whereas a substantial fraction of the boundary layers on the airfoils may be transitional at cruise conditions due to the change of density with altitude. The strong adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation at the latter low Reynolds number conditions. Large separation bubbles, particularly those which fail to reattach, cause a significant degradation of engine efficiency. A component efficiency drop of the order 2% may occur between takeoff and cruise conditions for large commercial transport engines and could be as large as 7% for smaller engines at higher altitude. An efficient means of of separation elimination/reduction is, therefore, crucial to improved turbine design. Because the large change in the Reynolds number from takeoff to cruise leads to a distinct change in the airfoil flow physics, a separation control strategy intended for cruise conditions will need to be carefully constructed so as to incur minimum impact/penalty at takeoff. A complicating factor, but also a potential advantage in the quest for an efficient strategy, is the intricate interplay between separation and transition for the situation at hand. Volino gives a comprehensive discussion of several recent studies on transition and separation under low-pressure-turbine conditions, among them one in the present facility. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions. If the transition occurs early in the boundary layer then separation may be reduced or completely eliminated. Transition in the shear layer of a separation bubble can lead to rapid reattachment. This suggests using control mechanisms to trigger and enhance early transition. Gad-el-Hak provides a review of various techniques for flow control in general and Volino discusses recent studies on separation control under low-pressure-turbine conditions utilizing passive as well as active devices. As pointed out by Volino, passive devices optimized for separation control at low Reynolds numbers tend to increase losses at high Reynolds numbers, Active devices have the attractive feature that they can be utilized only in operational regimes where they are needed and when turned off would not affect the flow. The focus in the present paper is an experimental Separation is induced on a flat plate installed in a closed-circuit wind tunnel by a shaped insert on the opposite wall. The flow conditions represent flow over the suction surface of a modem low-pressure-turbine airfoil ('Pak-B'). The Reynolds number, based on wetted plate length and nominal exit velocity, is varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low (0.2%) and high (2.5%) Gee-stream turbulence intensities are set using passive grids. A spanwise-oriented phased-plasma-array actuator, fabricated on a printed circuit board, is surface- flush-mounted upstream of the separation point and can provide forcing in a wide frequency range. Static surface pressure measurements and hot-wire anemometry of the base and controlled flows are performed and indicate that the glow-discharge plasma actuator is an effective device for separation control. of active separation control using glow discharge plasma actuators.

Hultgren, Lennart S.↗

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗

ELECTRIFICATION OF A HEAVY-DUTY OFF-ROAD MATERIAL HANDLER: ENERGY SAVINGS AND EMISSION REDUCTIONS

Federal regulations are driving the adoption of electrification technologies to reduce carbon dioxide equivalent (CO2e) emissions, a metric that quantifies the global warming potential of various greenhouse gases in terms of carbon dioxide (CO2). Although no specific CO2 regulations exist for heavy-duty off-road machines, future reductions are likely, given stricter emissions standards for on-road vehicles. The heavy-duty off-road sector offers significant fuel-saving potential, as its focus has traditionally been on reliability and performance rather than fuel efficiency. This dissertation examines fuel and CO2e savings opportunities on a heavy-duty off-road material handler, the Pettibone Cary-Lift 204i, from stock configuration to simple modifications to a complete teardown and reconfiguration of the machine with a plug-in series hybrid architecture using electrified hydraulics. The study begins by modeling the baseline machine’s fuel and energy consumption, calibrating with experimental data from custom operating cycles. An energy analysis identifies key areas for fuel savings. Two simple powertrain modifications result in a combined 16.2% fuel savings. Next, a Pugh-style analysis narrows a list of electrified architectures, leading to high-fidelity models that evaluate total lifetime CO2e and costs. Higher electrification levels reduce CO2e emissions but increase costs, and electricity grid emissions significantly impact CO2e for plug-in architectures. A plug-in series hybrid is chosen for the project. In its base control form, 49% fuel and 29% CO2e savings are expected from the plug-in series hybrid compared to the baseline machine. Further savings are pursued through regenerative braking (6.3%) and load-following hydraulic control (17.8%), totaling 24.1% fuel savings, and leading to a total of 61% fuel and 41% CO2e savings compared to the baseline. Battery chemistries and charging strategies are also analyzed for cost and CO2e impacts, finding LFP batteries as superior due to longevity, and overnight level 2 charging usually at a lower cost but resulting in higher emissions than opportunity DC fast-charging (DCFC). DCFC emissions are highly dependent on grid emissions, and DCFC cost is highly dependent on grid demand charges. Finally, artificial intelligence is applied to operating cycle recognition. Neural network accuracy ranges from 81% to 99%, with applications to worksite efficiency and safety improvements.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Practical challenges of model predictive control (MPC) for grid interactive small and medium commercial buildings

To the urgent call for mitigating climate change, substantial initiatives have been undertaken to deploy grid-interactive heating, ventilation, and air-conditioning (HVAC) controls, such as model predictive control (MPC) for buildings. These efforts typically aim to curtail peak energy demand, shift load and enhance overall energy efficiency. With the recent development of low-cost MPC technologies that don’t require extensive instrumentation or manual modeling, small and medium commercial buildings (SMCBs), which rarely utilize advanced HVAC control systems, have become candidates for grid-interactive efficient buildings (GEBs). However, despite the potential benefits and maturity of the technology itself, several practical challenges remain in real-world implementation. In this paper, we share the practical challenges that we have encountered in implementing and testing three types of MPC solutions (ON/OFF unit, dualfuel, and VRF systems) on multiple SMCB sites. We describe the MPC deployment process and discuss the lessons learned. The site selection, eligibility, and retrofit availability (e.g., utility price structure, thermostat communications, etc.) are the main discussion points at the beginning of the project. Also, the modeling automation and the best practices for interacting with endusers and handling erroneous situations are presented for successful operations.

woo Ham, Sang↗

Titanium Optics for Ion Thrusters

Ion thruster total impulse capability is limited, in part, by accelerator grid sputter erosion. A development effort was initiated to identify a material with a lower accelerator grid volumetric sputter erosion rate than molybdenum, but that could utilize the present NSTAR thruster grid design and fabrication techniques to keep development costs low, and perform as well as molybdenum optics. After comparing the sputter erosion rates of several atomic materials to that of molybdenum at accelerator voltages, titanium was found to offer a 45% reduction in volumetric erosion rates. To ensure that screen grid sputter erosion rates are not higher at discharge chamber potentials, titanium and molybdenum sputter erosion rates were measured at these potentials. Preliminary results showed only a slightly higher volumetric erosion rate for titanium, so that screen grid erosion is insignificant. A number of material, thermal, and mechanical properties were also examined to identify any fabrication, launch environment, and thruster operation issues. Several titanium grid sets were successfully fabricated. A titanium grid set was mounted onto an NSTAR 30 cm engineering model ion thruster and tested to determine optics performance. The titanium optics operated successfully over the entire NSTAR power range of 0.5 to 2.3 kW. Differences in impingement-limited perveances and electron backstreaming limits were found to be due to a larger cold gap for the titanium optics. Discharge losses for titanium grids were lower than those for molybdenum, likely due to a slightly larger titanium screen grid open area fraction. Radial distributions of beam current density with titanium optics were very similar to those with molybdenum optics at all power levels. Temporal electron backstreaming limit measurements showed that titanium optics achieved thermal equilibrium faster than molybdenum optics.

Soulas, George C.↗

Numerical Simulation of the Combustion of Fuel Droplets: Finite Rate Kinetics and Flame Zone Grid Adaptation (CEFD)

The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to support outstanding research endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers has developed a transient numerical model to study the combustion of suspended and moving droplets. The engines that propel missiles, jets, and many other devices are dependent upon combustion. Therefore, data concerning the combustion of fuel droplets is of immediate relevance to aviation and aeronautical personnel, especially those involved in flight operations. The experiments being conducted by Dr. Gogos and Dr. Nayagam s research teams, allow investigators to gather data for comparison with theoretical predictions of burning rates, flame structures, and extinction conditions. The consequent improved hndamental understanding droplet combustion may contribute to the clean and safe utilization of fossil hels (Williams, Dryer, Haggard & Nayagam, 1997, 72). The present state of knowledge on convective extinction of he1 droplets derives fiom experiments conducted under normal gravity conditions. However, any data obtained with suspended droplets under normal gravity are grossly affected by gravity. The need to obtain experimental data under microgravity conditions is therefore well justified and addresses one of the goals of NASA s Human Exploration and Development of Space (HEDS) microgravity combustion experiment.

Gogos, George↗

Analysis of a Cylindrical Specimen Heated by an Impinging Hot Hydrogen Jet

A computational conjugate heat transfer methodology was developed, as a first step towards an efficient and accurate multiphysics, thermo-fluid computational methodology to predict environments for hypothetical solid-core, nuclear thermal engine thrust chamber and components. A solid conduction heat transfer procedure was implemented onto a pressure-based, multidimensional, finite-volume, turbulent, chemically reacting, thermally radiating, and unstructured grid computational fluid dynamics formulation. The conjugate heat transfer of a cylindrical material specimen heated by an impinging hot hydrogen jet inside an enclosed test fixture was simulated and analyzed. The solid conduction heat transfer procedure was anchored with a standard solid heat transfer code. Transient analyses were then performed with ,variable thermal conductivities representing three composites of a material utilized as flow element in a legacy engine test. It was found that material thermal conductivity strongly influences the transient heat conduction characteristics. In addition, it was observed that high thermal gradient occur inside the cylindrical specimen during an impulsive or a 10 s ramp start sequence, but not during steady-state operations.

Wang, Ten-See↗

Parallel computing for power system climate resiliency: Solving a large-scale stochastic capacity expansion problem with mpi-sppy

Here we propose a nodal stochastic generation and transmission expansion planning model that incorporates the output from high-resolution global climate models through load and generation availability scenarios. We implement our model in Pyomo and perform computational studies on a realistically-sized test case of the California electric grid in a high performance computing environment. We propose model reformulations and algorithm tuning to efficiently solve this large problem using a variant of the Progressive Hedging Algorithm. We utilize the parallelization capabilities and overall versatility of mpi-sppy, exploiting its hub-and-spoke architecture to concurrently obtain inner and outer bounds on an optimal expansion plan. Initial results show that instances with 360 representative days on a system with over 8,000 buses can be solved to within 5% of optimality in under 4 h of wall clock time, a first step towards solving a large-scale power system expansion planning problem across a wide range of climate-informed operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analog In-Memory Computing for the Synthetic Aperture Radar Polar Format Algorithm

As the utility of synthetic aperture radar (SAR) systems increases in autonomous vehicles, satellites, and other power- and space-constrained edge applications, there is a growing need for processors that can form SAR images at low power. In recent years, analog in-memory compute (AIMC) has shown immense promise for accelerating neural networks and other matrix-vector multiplication (MVM) heavy workloads at the edge. Here, in this work, we examine how the polar format algorithm (PFA), a popular SAR image formation algorithm, can be mapped to these AIMC systems. The PFA maps readily onto analog MVMs because it primarily consists of two linear operations: interpolation of frequency-domain data to a Cartesian grid, followed by a 2-D Fourier transform. This work presents two approaches to map the interpolation operation onto MVMs in analog hardware: a chirp transform and a modified form of sinc interpolation. These mappings introduce algorithmic errors, and their effect on the quality of SAR image formation is examined, both quantitatively and qualitatively. In addition, the impact of errors introduced by the analog hardware is explored to determine which approach is optimal under varying assumptions about the underlying analog memory devices and circuits.

Analog computing↗

Post-Event Fault Identification with Machine Learning for Protection System Validation

Power system protection devices have transitioned over the past few decades from mechanical to analog devices, then to solid state and finally digital. Relays and their associated critical network of equipment have significantly increased in complexity. Even internally, relays have gained significant intricacy, with relatively simple overcurrent or differential functions now being assisted by a myriad of other functions. This is necessary as the grid becomes more complex, but it brings increased difficulty in monitoring and upkeep. Misoperation caused by accidental improper relay settings or deliberate malicious actions is a constant challenge faced by all utilities. These improper settings can be difficult to identify and may require exhaustive post-mortem analysis, typically after a major outage event has already occurred. A mechanism is needed for monitoring the behavior of protection systems to validate that their performance falls within expectations. Relays that fail to isolate a fault or trip when there is no system disturbance can be flagged for settings review in situations where this behavior may not have been noticed due to manual restoration or backup protection operations. This work presents a concept for a machine learning (ML) system capable of validating the performance of protection systems by identifying fault events and characterizing protection system responses based solely on available current and voltage measurements. As a first step in its development, an experimental dataset is generated, and a random forest model is implemented with high accuracy in distinguishing four power system scenarios.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Toward Immersed Boundary Simulation of High Reynolds Number Flows

In the immersed boundary (IB) method, the surface of an object is reconstructed with forcing terms in the underlying flow field equations. The surface may split a computational cell removing the constraint of the near wall gridlines to be aligned with the surface. This feature greatly simplifies the grid generation process which is cumbersome and expensive in particular for structured grids and complex geometries. The IB method is ideally suited for Cartesian flow solvers. The flow equations written in Cartesian coordinates appear in a very simple form and several numerical algorithms can be used for an efficient solution of the equations. In addition, the accuracy of numerical algorithms is dependent on the underlying grid and it usually deteriorates when the grid deviates from a Cartesian mesh. The challenge for the IB method lies in the representation of the wall boundaries and in providing an adequate near wall flow field resolution. The issue of enforcing no-slip boundary conditions at the immersed surface has been addressed by several authors by imposing a local reconstruction of the solution. Initial work by Verzicco et al. was based on a simple linear, one-dimensional operator and this approach proved to be accurate for boundaries largely aligned with the grid lines. Majumdar et al. used various multidimensional and high order polynomial interpolations schemes. These high order schemes, however, are keen to introduce wiggles and spurious extrema. Iaccarino & Verzicco and Kalitzin & Iaccarino proposed a tri-linear reconstruction for the velocity components and the turbulent scalars. A modified implementation that has proven to be more robust is reported in this paper. The issue of adequate near wall resolution in a Cartesian framework can initially be addressed by using a non-uniform mesh which is stretched near the surface. In this paper, we investigate an unstructured approach for local grid refinement that utilizes Cartesian mesh features. The computation of high Reynolds number wall bounded flows is particularly challenging as it requires the consideration of thin turbulent boundary layers, i.e. near wall regions with large gradients of the flow field variables. For such flows, the representation of the wall boundary has a large impact on the accuracy of the computation. It is also critical for the robustness and convergence of the flow solver.

Kalitzin, Georgi↗