Search NASA⌕ Search

SEARCH · Search NASA

Results for “control simulation testing”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Software Development for Data Visualization and Analysis of PN-Diodes & MOSFET Devices

The U.S. Department of Energy has a broad mission to ensure the nation’s security by addressing ongoing environmental challenges, developing novel energy production technologies, and mitigating nuclear security concerns. These efforts benefit from developing more effective semiconductor materials and devices. The research discussed in this report aims to further the development of these new electronics by improving implementations of vertical device structures, which offer theoretically better performance compared to lateral structures. Specifically, it reviews the background and development of a new and easy-to-use program for data processing and analysis from experiments run on vertical gallium nitride power devices. Current vertical structures fall short in critical performance metrics, such as On-Resistance and Breakdown Voltage, due to poor management of the device’s electric field. Junction termination extensions, or JTEs, are a field management technique that increases a device’s resilience to expected failure modes by controlling its surface electric fields. Each JTE design requires significant experimental validation, consisting of hundreds of tested devices, each outputting an enormous set of data points. The new software includes multiple methods of filtering large sets of device testing data to identify and remove flawed devices, automate device analyses, and provide statistical breakdowns for four metrics used to define the effectiveness of the JTE and the current carrying capacity of the device. These results can then be related back to fabrication techniques and device design, and in the future, be combined with processing simulations to further improve our understanding of the JTE mechanisms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Detecting Coastal Wetland Degradation by Combining Remote Sensing and Hydrologic Modeling

Sea-level rise and climate change stresses pose increasing threats to coastal wetlands that are vital to wildlife habitats, carbon sequestration, water supply, and other ecosystem services with global significance. However, existing studies are limited in individual sites, and large-scale mapping of coastal wetland degradation patterns over a long period is rare. Our study developed a new framework to detect spatial and temporal patterns of coastal wetland degradation by analyzing fine-scale, long-term remotely sensed Normalized Difference Vegetation Index (NDVI) data. Then, this framework was tested to track the degradation of coastal wetlands at the Alligator River National Wildlife Refuge (ARNWR) in North Carolina, United States, during the period from 1995 to 2019. We identified six types of coastal wetland degradation in the study area. Most of the detected degradation was located within 2 km from the shoreline and occurred in the past five years. Further, we used a state-of-the-art coastal hydrologic model, PIHM-Wetland, to investigate key hydrologic processes/variables that control the coastal wetland degradation. The temporal and spatial distributions of simulated coastal flooding and saltwater intrusion confirmed the location and timing of wetland degradation detected by remote sensing. The combined method also quantified the possible critical thresholds of water tables for wetland degradation. The remote sensing–hydrologic model integrated scheme proposed in this study provides a new tool for detecting and understanding coastal wetland degradation mechanisms. Our study approach can also be extended to other coastal wetland regions to understand how climate change and sea-level rise impact wetland transformations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

All-Digital Plug and Play Passive RFID Sensors for Energy Efficient Building Control

This is the final report for the project “All-Digital Plug and Play Passive RFID Sensors for Energy Efficient Building Control”, funded by DOE, and performed by Clemson University, Phase IV Engineering and Harvard University from October 1, 2016 to December 31, 2020. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, plug & play, battery-free, wireless sensor technology to enable low-cost (<$10 per node) indoor and outdoor temperature and humidity measurement for energy efficient building controls and operations. The proposed technology is based on the novel concept of all-digital sensing and its seamless integration with the passive RFID technology. This project focuses on the design, fabrication, material optimization, interrogation electronics, validation, and demonstration of the novel sensor nodes for building applications. The specific objectives of this research program include: (1) Design, fabricate and optimize a compact, robust, and high-resolution digitizer, which could convert rotation angle into digital numbers. (2) Develop the multi-physics-based modeling and simulation of the temperature/humidity transducer to derive a rational design of the architecture, dimension, structure, materials (i.e., mechanical, electrical, thermal and hygroscopic) properties and functions of the sensor node. (3) Design, fabricate and optimize the bi-material based humidity sensitive coil which could linearly transduce the environmental relative humidity variations to rotation angles. (4) Design, fabricate and optimize an UHF RFID platform which could support long-range passive wireless communications of 8-bit digital numbers. (5) Design, fabricate and optimize the miniaturized all-digital sensor using MEMS technology. (6) Validate the integrated all-digital sensing system in a real building environment to test the system’s performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

Model-based analysis of environmental controls over ecosystem primary production in an alpine tundra dry meadow

We investigated several key limiting factors that control alpine tundra productivity by developing an ecosystem biogeochemistry model. The model simulates the coupled cycling of carbon (C), nitrogen (N), and phosphorus (P) and their interactions with gross primary production (GPP). It was parameterized with field observations from an alpine dry meadow ecosystem using a global optimization strategy to estimate the unknown parameters. The model, along with the estimated parameters, was first validated against independent data and then used to examine the environmental controls over plant productivity. Our results show that air temperature is the strongest limiting factor to GPP in the early growing season, N availability becomes important during the middle portion of the growing season, and soil moisture is the strongest limiting factors by late in the growing season. Overall, the controls over GPP during the growing season, from strongest to weakest, are soil moisture content, air temperature, N availability, and P availability. This simulation provides testable predictions of the shifting nature of physical and nutrient limitations on plant growth. The model also indicates that changing environmental conditions in the alpine will likely lead to changes in productivity. For example, warming eliminates the control of P availability on GPP and makes N availability surpass air temperature to become the second strongest limiting factor. In contrast, an increase in atmospheric nutrient deposition eliminates the control of N availability and enhances the importance of P availability. Furthermore, these analyses provide a quantitative and conceptual framework that can be used to test predictions and refine ecological analyses at this long-term ecological research site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Optimization of operational strategies for rich gas enhanced oil recovery based on a pilot test in the Bakken tight oil reservoir

Horizontal well drilling and multistage hydraulic fracturing have been demonstrated as effective approaches for stimulating oil production in the Bakken tight oil reservoir. However, after multiple years of production, primary oil recovery in the Bakken is generally less than 10% of the estimated original oil in place. Gas huff ‘n’ puff (HnP) has been tested in the Bakken Formation as an enhanced oil recovery (EOR) method; however, most field pilot test results showed no significant incremental oil production. One of the factors affecting HnP EOR performance is premature gas breakthrough, which is one of the most critical issues observed in the field because of the presence of interwell fractures. Consequently, injected gas rapidly reaches adjacent production wells without contacting reservoir rock and increasing oil recovery. Proper conformance control is therefore needed to avoid early gas breakthrough and improve EOR performance. In this study, a rich gas EOR pilot in the Bakken was carefully analyzed to collect the essential reservoir and operational data. A simulation model with 16 wells was then developed to reproduce the production history and predict the EOR performance with and without conformance control. EOR operational strategies, including single- and multiple-well HnP, with different gas injection constraints were investigated. The simulation results of single-well HnP without conformance control showed that a rich gas injection rate of at least 10 MMscfd was needed to yield meaningful incremental oil production. The strategy of conformance control via water injection could significantly improve oil production in the HnP well, but injecting an excessive amount of water also leads to water breakthrough and loss of oil production in the offset wells. By analyzing the production performance of the wells individually, the arrangement of wells was optimized for multiple-well HnP EOR. The multiwell results showed that rich gas EOR could improve oil production up to 7.4% by employing conformance control strategies. Furthermore, replacing rich gas with propane as the injection gas could result in 14% of incremental oil production.

02 PETROLEUM↗
DOE OSTI2025

Computational Fluid Dynamics Analysis of the Molten Salt Tritium Transport Experiment Test Section

Tritium, a radionuclide produced through neutron capture by lithium and other elements (beryllium and fluoride) in molten salts, presents unique challenges to radionuclide release. This is true for both fusion energy breeder blankets and molten salt fission reactors. The fundamental understanding of tritium transport is crucial to the safe design and operation of these reactors. The Molten Salt Tritium Transport Experiment (MSTTE), currently under construction at Idaho National Laboratory, aims to investigate tritium transport phenomena using a forced-convection fluoride salt loop. This loop is designed to study various transport mechanisms, such as permeation through metals and gas-liquid interactions, and is intended to support future research on tritium extraction units. A critical aspect of the MSTTE loop design is ensuring a fully developed velocity profile before the fluid reaches the permeation test section where measurements are made. This study employs computational fluid dynamics to model the salt flow behavior within the MSTTE permeation test section. A realizable k-ε turbulent model with enhanced wall treatment is used to simulate the single-phase, vertical upward flow of molten salt FLiNaK under isothermal conditions. The simulation results indicated flow distortion and underdeveloped profiles at all planned flow rates within the test section due to the 85-deg sharp bend. To address this issue, a reduced diameter with a reducer and expander and a flow conditioner are investigated to achieve fully developed flow. The analysis showed that the flow conditioner successfully corrected the flow profile, achieving fully developed behavior at a flow rate of 50 liters per minute (LPM). In conclusion, this research enhances our understanding of flow dynamics in molten salt systems and contributes to optimizing tritium transport control technologies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Optimal Submodule Capacitor Sizing for Modular Multilevel Converters with Common Mode Voltage Injection and Circulating Current Control

Abstract - Modular multilevel converters (MMCs), with the filterless and transformer-less features, are considered as a promising topology for high-efficiency medium voltage (MV) high voltage (HV) drive systems. However, at low frequencies (< 30 Hz), a large current ripple flowing through the capacitors of each submodule will cause significant capacitor voltage fluctuation. This paper proposes an optimal submodule capacitor selection based on the common mode voltage injection. The proposed sizing model includes two functions: a) accurately estimating the capacitor voltage ripple based on the drive system parameters; and b) identifying the minimum capacitance of the submodule capacitors based on the allowed maximum capacitor voltage fluctuations. Simulation and hardware-in-the-loop (HIL) test results show that the error between the actual and estimated values of the voltage fluctuation is within 3% using the proposed model.

Capacitor Sizing↗
DOE OSTI2022

A Fairness-based Distributed Energy Coordination for Voltage Regulation in Distribution Systems

The rapid adoption of photovoltaic (PV) systems combined with falling prices of energy storage is paving the way for a future in which customers could locally supply their energy needs and export surplus power to the grid. However, reverse power flow from multiple sites in a network can cause overvoltage that decreases system reliability and the utilization level of PV system. In this paper, we present a two-stage approach to solve the energy coordination problem in distribution networks. In the first stage, a household with PV and energy storage is optimized to allow customers to adopt different control strategies such as load following and tariff arbitrage. We show that such customer-led control strategies can cause reverse power flows during peak PV generation periods, thus resulting in overvoltage. In such cases, the second stage of the proposed approach is called upon to solve a fairness-based energy coordination problem that can maintain voltages within limits while maximizing the PV generation while fairly distributing the actions among all the PV systems. As a result, each of the resources in the network has equal opportunity to export power and shares the responsibility of voltage regulation. Simulation studies have been carried out on the IEEE 13-bus test feeder to verify the effectiveness of the proposed approach.

Distributed energy resources, fairness, overvoltag↗
DOE OSTI2016

Predicting the Responses of Soil Nitrite-Oxidizers to Multi-Factorial Global Change: A Trait-Based Approach

Soil microbial diversity is huge and a few grams of soil contain more bacterial taxa than there are bird species on Earth. This high diversity often makes predicting the responses of soil bacteria to environmental change intractable and restricts our capacity to predict the responses of soil functions to global change. Here, using a long-term field experiment in a California grassland, we studied the main and interactive effects of three global change factors (increased atmospheric CO 2 concentration, precipitation and nitrogen addition, and all their factorial combinations, based on global change scenarios for central California) on the potential activity, abundance and dominant taxa of soil nitrite-oxidizing bacteria (NOB). Using a trait-based model, we then tested whether categorizing NOB into a few functional groups unified by physiological traits enables understanding and predicting how soil NOB respond to global environmental change. Contrasted responses to global change treatments were observed between three main NOB functional types. In particular, putatively mixotrophic Nitrobacter, rare under most treatments, became dominant under the 'High CO 2 +Nitrogen+Precipitation' treatment. The mechanistic trait-based model, which simulated ecological niches of NOB types consistent with previous ecophysiological reports, helped predicting the observed effects of global change on NOB and elucidating the underlying biotic and abiotic controls. Our results are a starting point for representing the overwhelming diversity of soil bacteria by a few functional types that can be incorporated into models of terrestrial ecosystems and biogeochemical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Development of a hardware-in-the-loop testbed for laboratory performance verification of flexible building equipment in typical commercial buildings

The goals of reducing energy costs, shifting electricity peaks, increasing the use of renewable energy, and enhancing the stability of the electric grid can be met in part by fully exploiting the energy flexibility potential of buildings and building equipment. The development of strategies that exploit these flexibilities could be facilitated by publicly available high-resolution datasets illustrating how control of HVAC systems in commercial buildings can be used in different climate zones to shape the energy use profile of a building for grid needs. This article presents the development and integration of a Hardware-In-the-Loop Flexible load Testbed (HILFT) that integrates physical HVAC systems with a simulated building model and simulated occupants with the goal of generating datasets to verify load flexibility of typical commercial buildings. Compared to simulation-only experiments, the hardware-in-the-loop approach captures the dynamics of the physical systems while also allowing efficient testing of various boundary conditions. The HILFT integration in this article is achieved through the co-simulation among various software environments including LabVIEW, MATLAB, and EnergyPlus. Although theoretically viable, such integration has encountered many real-world challenges, such as: 1) how to design the overall data infrastructure to ensure effective, robust, and efficient integration; 2) how to avoid closed-loop hunting between simulated and emulated variables; 3) how to quantify system response times and minimize system delays; and 4) how to assess the overall integration quality. Lessons-learned using the examples of an AHU-VAV system, an air-source heat pump system, and a water-source heat pump system are presented.

Chen, zhelun↗
DOE OSTI2016

Stripline kicker for integrable optics test accelerator

We present a design of a stripline kicker for Integrable Optics Test Accelerator (IOTA). For its experimental program IOTA needs two full-aperture kickers, capable to create an arbitrary controllable kick in 2D. For that reason their strengths are variable in a wide range of amplitudes up to 16 mrad, and the pulse length 100 ns is less than a revolution period for electrons. In addition, the kicker should have a physical aperture of 40 mm for a proposed operation with proton beam, and an outer size of 70 mm to fit inside existing quadrupole magnets to save space in the ring. Computer simulations using CST Microwave Studio show high field uniformity and wave impedance close to 50 {\Omega}.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Modernizing the Legacy Fission Wire Measurement System for the Advanced Test Reactor-Critical Facility

Operational lifetime extensions of existing research reactors have emphasized the need for refurbishment, replacements, and upgrades to supporting equipment and instrumentation. The Advanced Test Reactor (ATR) at Idaho National Laboratory (INL), which entered service in 1967, has recently completed the sixth core internals change-out and has scheduled operations until at least 2040. Reactor maintenance and operational risk management is critically important in the research reactor community, however supporting measurement systems sometimes get overlooked when maintenance is planned. The Fission Wire Measurement System (FWMS) is a custom measurement system designed in the 1960s to measure the beta-particle activity of irradiated uranium-aluminum fission wires. This measurement is conducted to determine the fission rate profile of the Advanced Reactor Test Critical (ATR-C) facility. The ATR-C is an open-pool, low-power test reactor that was purpose driven to resemble ATR and is used to qualify experiment configurations and verify core models prior to full-power experiment irradiations in ATR. A power distribution measurement in ATR-C uses uranium-aluminum wires that are distributed throughout the ATR-C core to validate simulation and modeling results. These measurements require 340 to 1500 wires to be irradiated and measured within a 12-hour window. The activity of the wires is measured in the required time with the FWMS, which was put into service in 1965 at the Radiation Measurements Laboratory (RML). The system consists of 4 measurement channels and one reference channel, each with a 2-pi proportional gas flow detector and the measurement channels each have an automated sample changer. This legacy system is crucial to the continued operations of ATR and has undergone some minor hardware upgrades since 1965, however the system presently relies on custom control boards, custom gas ion chambers, analog amplifiers/discriminators, and a user interface (UI) for the system written in outdated code. Much of the equipment and software is custom with no commercial replacements or support and limited documentation. The existing control software requires an operating system that is no longer supported, creating more vulnerabilities to continued operations. A project is underway with a third-party vendor to design, build, and document a new control and data acquisition system (CDAS) for the FWMS. The new upgrade will replace the control system, computer, UI, sample changer motors, and main power supply while maintaining the interface with existing detector hardware. The upgraded system will be operated in parallel with the current hardware and software to conduct validation testing. This equipment upgrade demonstrates the commitment at ATR to ensuring successful operations and potential future research reactors at INL.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2021

In-pixel automatic threshold calibration for the CMS Endcap Timing Layer readout chip

We present the implementation and verification of an in-pixel automatic threshold calibration circuit for the CMS Endcap Timing Layer (ETL) in the High-Luminosity LHC upgrade. The discriminator threshold of the ETL readout chip (ETROC) needs to be calibrated regularly to mitigate the circuit baseline change. Traditional methods need a lot of communication through a slow control system hence are time-consuming. This paper describes an in-pixel automatic scheme with improvements in operating time and usability. In this scheme, a sample-accumulation circuit is used to measure the average discriminator output. A binary successive approximation and linear combination scan are applied to find the equivalent baseline. The actual calibration procedure has been first implemented in FPGA firmware and tested with the ETROC front-end prototype chip (ETROC0). The calibration circuit has been implemented with Triple Modular Redundancy (TMR) and verified with Single Event Effects (SEEs) simulation. A complete calibration process lasts 35 ms with a 40 MHz clock. In the worst case, the dynamic and static power consumption are estimated to be 300 uW and 10.4 uW, respectively. The circuit design, implemented in a 65 CMOS technology, will be integrated into ETROC2, the next iteration of the ETROC with a 16 x16 pixel matrix.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

An automated approach to the alignment of compound refractive lenses

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Virtualizing Industrial Control Networks for Cyber Resilience Experiments

Industrial control systems (ICS) networks are undergoing constant shifts to accommodate new security measures. It is challenging to test varying network configurations and security tools with physical systems as they typically include large, expensive equipment. Not only this, but researchers often do not have access to this type of equipment for development of new security tools and techniques. As a solution to these issues, this work presents a set of tools for utilizing GNS3 and Docker as a virtual ICS network. Additionally, the virtual network can be attached to physical devices including network switches, hardware simulations, and intelligent electronic devices (IEDs). Two case studies showcase a relatively complex automatically generated network and an attack on a simple ICS network with an example mitigation.

42 ENGINEERING↗
DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗
DOE OSTI2024

Progress of Gas Injection EOR Surveillance in the Bakken Unconventional Play—Technical Review and Machine Learning Study

Although considerable laboratory and modeling activities were performed to investigate the enhanced oil recovery (EOR) mechanisms and potential in unconventional reservoirs, only limited research has been reported to investigate actual EOR implementations and their surveillance in fields. Eleven EOR pilot tests that used CO2, rich gas, surfactant, water, etc., have been conducted in the Bakken unconventional play since 2008. Gas injection was involved in eight of these pilots with huff ‘n’ puff, flooding, and injectivity operations. Surveillance data, including daily production/injection rates, bottomhole injection pressure, gas composition, well logs, and tracer testing, were collected from these tests to generate time-series plots or analytics that can inform operators of downhole conditions. A technical review showed that pressure buildup, conformance issues, and timely gas breakthrough detection were some of the main challenges because of the interconnected fractures between injection and offset wells. The latest operation of co-injecting gas, water, and surfactant through the same injection well showed that these challenges could be mitigated by careful EOR design and continuous reservoir monitoring. Reservoir simulation and machine learning were then conducted for operators to rapidly predict EOR performance and take control actions to improve EOR outcomes in unconventional reservoirs.

Energy & Fuels↗
DOE OSTI2017

Streamline integration as a method for two-dimensional elliptic grid generation

We propose a new numerical algorithm to construct a structured numerical elliptic grid of a doubly connected domain. Our method is applicable to domains with boundaries defined by two contour lines of a two-dimensional function. Furthermore, we can adapt any analytically given boundary aligned structured grid, which specifically includes polar and Cartesian grids. The resulting coordinate lines are orthogonal to the boundary. Grid points as well as the elements of the Jacobian matrix can be computed efficiently and up to machine precision. In the simplest case we construct conformal grids, yet with the help of weight functions and monitor metrics we can control the distribution of cells across the domain. Our algorithm is parallelizable and easy to implement with elementary numerical methods. We assess the quality of grids by considering both the distribution of cell sizes and the accuracy of the solution to elliptic problems. Among the tested grids these key properties are best fulfilled by the grid constructed with the monitor metric approach. - Graphical abstract: - Highlights: • Construct structured, elliptic numerical grids with elementary numerical methods. • Align coordinate lines with or make them orthogonal to the domain boundary. • Compute grid points and metric elements up to machine precision. • Control cell distribution by adaption functions or monitor metrics.

97 MATHEMATICS AND COMPUTING↗