Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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

Full time-dependent SOLPS-ITER simulation of the SPARC tokamak: actuator design for particle and divertor condition control *

Abstract This paper presents the application of full time-dependent SOLPS-ITER simulations for actuator design in the SPARC tokamak. This study employs both the EIRENE module, a neutral solver, and the B2.5 plasma module in a time-dependent mode. This is in contrast to most SOLPS simulations, which focus on steady-state solutions, where the neutral distribution is evolved without any time limit or for a time step of 1 ⋅ 10 − 3 second, which is several orders of magnitude larger than the fluid plasma time step. The time-dependent EIRENE was tested with a fixed B2.5 background and compared with a simple conductance based model in a simplified pump chamber geometry. This comparison aimed to verify the reliability of the neutral relaxation timescale derived from the time-dependent EIRENE. Subsequently, a full time-dependent simulation was performed in a realistic geometry, with the Monte-Carlo neutral time step synchronized with the plasma fluid time step. The numerical setup of the code, including relative time steps and the size of the census data used to store Monte-Carlo particle information is considered. The full-time dependent simulations are then applied to inform the design of the SPARC louver structure, which affects divertor plasma parameters by regulating the neutral conductance from the divertor to the pump. The response of the plasma and neutral parameters was captured on a timescale that enables the design of the actuator to consider time-dependent control capability. It was found that changing the louver opacity has an equivalent effect as varying the gas throughput via puff actuation. Therefore, equivalent divertor plasma conditions can be obtained from both actuators, while the neutral pressure distribution in the pump and divertor differs for each actuator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Gravitational waves from asymmetric oscillon dynamics?

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of e-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This new finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multiple Magnetorotons and Spectral Sum Rules in Fractional Quantum Hall Systems

We study, numerically, the charge neutral excitations (magnetorotons) in fractional quantum Hall systems, concentrating on the two Jain states near quarter filling, ν = 2/7 and ν = 2/9, and the ν = 1/4 Fermi-liquid state itself. In contrast to the ν = 1/3 states and the Jain states near half filling, on each of the two Jain states ν = 2/7 and ν = 2/9 the graviton spectral densities show two, instead of one, magnetoroton peaks. The magnetorotons have spin 2 and have opposite chiralities in the ν = 2/7 state and the same chirality in the ν = 2/9 state. Here, we also provide a numerical verification of a sum rule relating the guiding center spin s¯ with the spectral densities of the stress tensor.

2-dimensional systems↗
DOE OSTI2024

Electrical Control of Spin and Valley in Spin-Orbit Coupled Graphene Multilayers

Electrical control of magnetism has been a major technological pursuit of the spintronics community, owing to its far-reaching implications for data storage and transmission. Here, in this study, we propose and analyze a new mechanism for electrical switching of isospin, using chiral-stacked graphene multilayers, such as Bernal bilayer graphene or rhombohedral trilayer graphene, encapsulated by transition metal dichalcogenide (TMD) substrates. Leveraging the proximity-induced spin-orbit coupling from the TMD, we demonstrate electrical switching of correlation-induced spin and/or valley polarization, by reversing a perpendicular displacement field or the chemical potential. We substantiate our proposal with both analytical arguments and self-consistent Hartree-Fock numerics. Finally, we illustrate how the relative alignment of the TMDs, together with the top and bottom gate voltages, can be used to selectively switch distinct isospin flavors, putting forward correlated Van der Waals heterostructures as a promising platform for spintronics and valleytronics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Glass formulation and composition optimization with property models: A review

Abstract Glass is a versatile material with a remarkable history and many practical applications. It plays a critical role in our everyday lives, the advancement of science, and the development of many technologies. The Edisonian type trial‐and‐error method was commonly used for conventional design of glass compositions, which was time‐consuming and costly. With the urgent need to develop new glass compositions for technology applications rapidly, it has become necessary to develop precise property models with predictive powers using large databases and efficient formulation approaches. This paper reviews the design of glass compositions using these analytical and numerical models of composition–structure–property relations of glasses, some based on large databases and machine learning approaches. Aspects of data collection, model fitting, feature extraction, model evaluation, and uncertainty quantification will be covered. Furthermore, advances in the glass optimization framework and available tools are summarized with examples. The outlook and perspective for further glass property model development and formulation approaches are discussed.

Lu, Xiaonan↗
DOE OSTI2025

Local adaptation of both plant and pathogen: an arms‐race compromise in switchgrass rust

In coevolving species, parasites locally adapt to host populations as hosts locally adapt to resist parasites. Parasites often outpace host local adaptation since they have rapid life cycles, but host diversity, the strength of selection, and external environmental influence can result in complex outcomes. To better understand local adaptation in host–parasite systems, we examined locally adapted switchgrass (Panicum virgatum), and its leaf rust pathogen (Puccinia novopanici) across a latitudinal range in North America. We grew switchgrass genotypes in 10 replicated multiyear common gardens, measuring rust severity from natural infection in a ‘host reciprocal transplant’ framework for testing local adaptation. We conducted genome-wide association mapping to identify genetic loci associated with rust severity. Genetically differentiated rust populations were locally adapted to northern and southern switchgrass, despite host local adaptation to environmental conditions in the same regions. Rust resistance was polygenic, and distinct loci were associated with rust severity in the north and south. We narrowed a previously identified large-effect quantitative trait locus for rust severity to a candidate YELLOW STRIPE-LIKE gene and linked numerous other loci to defense-related genes. Overall, our results suggest that both hosts and parasites can be simultaneously locally adapted, especially when parasites impose less selection than other environmental factors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Local and nonlocal phase-field models of tumor growth and invasion due to ECM degradation

We present and analyze new multi-species phase-field mathematical models of tumor growth and ECM invasion. The local and nonlocal mathematical models describe the evolution of volume fractions of tumor cells, viable cells (proliferative and hypoxic cells), necrotic cells, and the evolution of matrix-degenerative enzyme (MDE) and extracellular matrix (ECM), together with chemotaxis, haptotaxis, apoptosis, nutrient distribution, and cell-to-matrix adhesion. We provide a rigorous proof of the existence of solutions of the coupled system with gradient-based and adhesion-based haptotaxis effects. In addition, we discuss finite element discretizations of the model, and we present the results of numerical experiments designed to show the relative importance and roles of various effects, including cell mobility, proliferation, necrosis, hypoxia, and nutrient concentration on the generation of MDEs and the degradation of the ECM.

Mathematics↗
DOE OSTI2020

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Casing-integrity mapping using top-casing electrodes and surface based EM fields

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage and hazardous waste disposal. We introduce a simple non-invasive technology for casing integrity screening, based on the continuity of electrical current flow. Applying low frequency current to a wellhead, with a distant return electrode, produces a casing current dependent on the properties and depth extent of the well casing as well as the background formation. These currents in-turn generate surface electrical fields that can be captured in a radial profile and be used to analyze properties of the well casing. Numerical modeling results reveal a strong relation of the electric field to the casing properties and depth extent of the well. A small breakage in the casing produces a profile coincident to a cased well with a completion depth above the break. A corroded patch, where the casing conductivity is reduced, also alters the field profiles and its depth may be estimated by comparing to the profile expected from the well completion diagrams. The electric field profiles are also strongly dependent on background resistivity distributions and on whether the well was drilled using water or oil-based drilling fluids. We validate the proof of concept in a field experiment, where we applied currents at the wellheads of two wells with different casing lengths. The two profiles were similar in appearance but offset in amplitude by more than a factor of 5, consistent with the theoretical analysis as well as the 3D modeling results. These results demonstrate that our proposed approach has promise for mapping the general casing condition without well intervention. This approach can be a practical and effective tool for rapidly screening a number of wells before expensive logging-based technologies are employed for casing inspection in detail.

58 GEOSCIENCES↗
DOE OSTI2021

Paying it Forward: How an Emerging Unconventional Play can Hit the Ground Running

While every tight oil play is unique, there are lessons that can be transferred from one play to another to improve the efficiency and pace of production operations and development. These improvements may not fit precisely in every basin or play but generally hold to themes that can be tested against and built upon. Themes such as the quantity of proppant, longer lateral length, or the number of stages can be directly tied to increased productivity. However, there are diminishing returns on these investment activities for which each operator, within a given play, will be required to identify and mitigate against. This is especially true as the industry steps in and begins developing new tight oil plays. In their nascent stages, these plays may have limited well penetrations and, as a result, limited geological and performance data from which to extrapolate. Pulling together an understanding of where the industry currently resides in terms of how to exploit these resources can provide a boost in terms of working towards greatly improved well completions. In 2019, the US EIA estimated that nearly 8 million barrels of oil per day were produced from tight oil reservoirs in the United States (US EIA, 2020). This represents over 60% of the domestic crude production, originating from multiple reservoirs in the Permian Basin (TX) as well as the Bakken (MT, ND), Eagle Ford (TX), Niobrara (CO, WY), and Anadarko Basin (OK) formations, among others. As such, there are 1,000s of wells across these numerous tight oil plays that can relate an informative story. To build this story, the interplay of geology, well spacing, lateral length, and stimulation, all critical to economic success, will be explored. This paper proposes to look back at these mature tight oil (and gas) basins and bring forth an understanding of what lessons can be applied to the emerging Powder River Basin tight oil reservoirs (Mowry and the Turner/Frontier). The authors will delve into the four broad topics of geology, well spacing, lateral length, and stimulation, highlighting case studies to demonstrate those lessons from established tight oil plays that will underpin planned activities at a Field Laboratory Test Site in the southern Powder River Basin.

02 PETROLEUM↗
DOE OSTI2023

Oblique propagation and temperature effects on the resonant right-hand ion beam instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A New Approach to Quantify Shallow Water Hydrologic Exchanges in a Large Regulated River Reach

Hyporheic exchange is a crucial component in the water cycle. The strength of the exchange directly affects the biogeochemical and ecological processes occurred in the hyporheic zone from micro to reach scale. Hyporheic fluxes can be quantified using many direct and indirect measurements as well as analytical and numerical modeling tools. However, in a relatively large river, these methods are limited by accessibility, the difficulty of performing representative sampling, and complexity of geomorphologic features and subsurface properties. In rivers regulated by hydroelectric dams, quantifying hyporheic fluxes becomes more challenging due to frequent hydropeaking events, featured by hourly to daily variations in flow and river stages created by dam operations(Hancock 2002). In this study, we developed and validated methods that based on field measurements to estimate shallow water hyporheic fluxes across the river bed at five locations along the shoreline of the Columbia River. Vertical thermal profiles measured by self-recording thermistors were combined with time series of hydraulic gradients derived from river stage and water level at in-land wells to estimate the hyporheic flux rate. The results suggested that the hyporheic exchange rate had high spatial and temporal heterogeneities over the riverbed, with predicted flux rate varies from +1×10 -6 m s-1 to -1.5×10 -6 m s -1 under various flow conditions at the some locations, and with a magnitude of fluxes 6-9 times higher in the primary channel than that in the secondary channel. Furthermore, the variations on shallow water hyporheic flow dynamics may further lead to different biogeochemical and ecological consequences at different river segments.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1971

An analysis of the flow field near the fuel injection location in a gas core reactor.

An analytical study is presented which shows the effects of large energy release and the concurrent high acceleration of inner stream fluid on the coaxial flow field in a gas core reactor. The governing equations include the assumptions of only radial radiative transport of energy represented as an energy diffusion term in the Euler equations. The method of integral relations is used to obtain the numerical solution. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Weinstein, H.↗
NASA NTRS1973

Outline of a new approach to the analysis of complex systems and decision processes.

Development of a conceptual framework for dealing with systems which are too complex or too ill-defined to admit of precise quantitative analysis. The approach outlined is based on the premise that the key elements in human thinking are not numbers, but labels of fuzzy sets - i.e., classes of objects in which the transition from membership to nonmembership is gradual rather than abrupt. The approach in question has three main distinguishing features - namely, the use of so-called 'linguistic' variables in place of or in addition to numerical variables, the characterization of simple relations between variables by conditional fuzzy statements, and the characterization of complex relations by fuzzy algorithms.

Zadeh, L. A.↗
NASA NTRS1974

Evaluation of analytical procedures for prediction of turbulent boundary layers on a porous wall

An analytical study has been made to determine how well current boundary layer prediction techniques work when there is mass transfer normal to the wall. The data that were considered in this investigation were for two-dimensional, incompressible, turbulent boundary layers with suction and blowing. Some of the bleed data were taken in an adverse pressure gradient. An integral prediction method was used three different porous wall skin friction relations, in addition to a solid-surface relation for the suction cases. A numerical prediction method was also used. Comparisons were made between theoretical and experimental skin friction coefficients, displacement and momentum thicknesses, and velocity profiles. The integral method with one of the porous wall skin friction laws gave very good agreement with data for most of the cases considered. The use of the solid-surface skin friction law caused the integral to overpredict the effectiveness of the bleed. The numerical techniques also worked well for most of the cases.

Towne, C. E.↗