Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Role of structural anisotropy in geometric birefringence of high-birefringence index-guiding PCFs: analysis by moment of inertia

We have studied the role of structural anisotropy on the geometric birefringence of four typical high-birefringence (HB) index-guiding photonic crystal fibers (PCFs) with the same air-silica structure and circular, ellipse, and rectangle air holes, respectively. The normalized difference between the moment of inertia ΔI in x and y direction of the PCF structure is utilized to quantificationally describe the structural anisotropy. We demonstrate that the phase birefringence B increases monotonously when the normalized ΔI rises in different cases of PCFs structures. In order to obtain high birefringence B, it can be an effective method to optimize the parameters of the structure to increase the normalized ΔI. This work can be not only a reference of founding the essence of geometric birefringence and proposing the optimal design of HB-IG-PCFs, but also an enlightenment for bringing the different method (compared to the finite element) into the investigation of PCFs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Designing a high sensitivity hexagonal nano-cavity photonic crystal resonator for the purpose of seawater salinity sensing

In this paper a sensing platform based on nano-cavity photonic crystal resonator (NC-PCR) is proposed for seawater salinity concentration detection. It has been demonstrated that the proposed NC-PCR has high sensitivity and selectivity for various concentrations of seawater salinity. The sensitivity is computed to be as high as 544 nm/RIU (nm per refractive index unit) and the quality factor (Q) is 3860. Also it has been shown that the figure of merit for the proposed structure calculated to be as high as ~ 1.234 × 10{sup 3} RIU{sup −1} which is very appropriate for sensing applications. Moreover, the minimum detection limit for the presented structure is obtained as 8.1 × 10{sup −5} RIU. The finite-difference time-domain method is employed to investigate the transmission characteristics of the device for a verity of refractive indices. It should also be mentioned that our design is comparatively easy to fabricate and certainly the advantages of our design justify its usage in sensing applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Effect of electron–electron interactions on optical properties of GaN/AlN quantum wells: a nonlinear Schrödinger equation approach

Researchers usually neglect the electron–electron interaction effect when they study the optical properties of semiconducting nanostructures through the compact density matrix approach. In the existing papers, this work has also been done through self-consistent solution of the Schrödinger and Poisson equations. For the first time, we have investigated the electronic and optical properties of wurtzite GaN/AlN quantum wells in the presence of electron–electron interactions assuming weak electron–electron scattering through a numerical solution of a nonlinear Schrödinger equation. It is shown that the electron–electron interaction effect leads to the removal of the N-fold degeneracy in the energy levels of N-well multiple quantum well systems. Effects of different parameters such as electric field, the strength of the electron–electron interaction, a total length of the system and also two different numbers of wells (1 and 2) on the mentioned optical properties have been investigated. Also, we observed redshifts in the optical rectification coefficient peak positions when the nonlinear parameter increased. Finite difference self-consistent method was employed to solve the nonlinear Schrödinger equation and we demonstrated the method convergence. Finally, we illustrated the accuracy of the method by comparing our energy values with those of the Sinc method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A new proposal for a high-performance 4-channel demultiplexer based on 2D photonic crystal using three cascaded ring resonators for applications in advanced optical systems

In this paper, a high-performance 4-channel demultiplexer (DEMUX) based on 2D photonic crystal using three cascaded ring resonators is proposed. The proposed DEMUX is composed of square lattice of silicon dielectric rods with a refractive index of 3.46 suspended in air. For the proposed device, the lattice constant a = 590 nm and the filling factor r/a = 0.2. In this study, the plane wave expansion and finite-difference time-domain methods are used respectively to characterize the photonic band gap and to investigate the optical behaviour of the structure. It is demonstrated that by choosing appropriate values for the radius of the inner rods, scattering rods, coupling rods and the coupling length between the bus and the drop waveguides, the desired wavelengths can effectively be separated. The resonant wavelengths of the proposed DEMUX are 1528.7 nm, 1535.3 nm, 1531 nm, and 1533 nm. The mean quality factor, spectral bandwidth and the transmission efficiency of our device are computed to be respectively 7358.5, 0.217 nm and 99.25%. The proposed DEMUX has a mean channel spacing, minimum and maximum values of crosstalk of about 2.2 nm, − 46.68 dB and − 9.79 dB, respectively. Furthermore, the footprint of the device is 294.4 μm{sup 2}, which is suitable for photonic integrated circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A 2 × 1 all-optical multiplexer using Kerr nonlinear nano-plasmonic switch

In this paper, a novel surface plasmon polariton based 2 × 1 multiplexer, which is based on an adder module and an all-optical controllable nano-plasmonic switching module is presented. It is shown that by applying a pump light, the transmissions of the input signals into the output waveguide can be controlled. The operation of the proposed device is investigated through the use of coupled-mode theory and finite-difference time-domain method. In the proposed adder module, the transmissions of the input signals at the working wavelengths are about 0.7. Based on the simulation results, by applying a pump light intensity of 150 KW/cm{sup 2} the transmissions of input signals alter with the extinction ratio of 22.04 dB, which shows a significant switching effect. For a pump light intensity of 150 KW/cm{sup 2} the switching times of input signals are measured about 560 fs. It is expected that the proposed 2 × 1 multiplexer can be used in all-optical signal processing, nano-plasmonic circuits, communication, etc.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Diffraction of radio frequency waves by spatially modulated interfaces in the plasma edge in tokamaks

We report the use of radio frequency (RF) waves in fusion plasmas for heating, for non-inductive current generation, for profile control and for diagnostics has been well established. The RF waves, excited by antenna structures placed near the wall of a fusion device, have to propagate through density fluctuations at the plasma edge. These fluctuations can modify the properties of the RF waves that propagate towards the core of the plasma. A full-wave electromagnetic computational code ScaRF based on the finite difference frequency domain (FDFD) method has been developed to study the effect of density turbulence on RF waves. The anisotropic plasma permittivity used in the scattering studies is that for a magnetized, cold plasma. The code is used to study the propagation of an RF plane wave through a modulated, spatially periodic density interface. Such an interface could arise in the edge region due to magnetohydrodynamic instability or drift waves. The frequency of the plane wave is taken to be in the range of the electron cyclotron frequency. The scattering analysis is applicable to ITER-like plasmas, as well as to plasmas in medium sized tokamaks such as TCV, ASDEX-U and DIII-D. The effect of different density contrasts across the interface and of different spatial modulations are discussed. While ScaRF is used to study a periodic density fluctuation, the code is general enough to include different varieties of density fluctuations in the edge region – such as blobs and filaments, and spatially random fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling Structural Colors from Disordered One-Component Colloidal Nanoparticle-Based Supraballs Using Combined Experimental and Simulation Techniques

Bright, saturated structural colors in birds have inspired synthesis of self-assembled, disordered arrays of assembled nanoparticles with varied particle spacings and refractive indices. However, predicting colors of assembled nanoparticles, and thereby guiding their synthesis, remains challenging due to the effects of multiple scattering and strong absorption. In this work, we use a computational approach to first reconstruct the nanoparticles’ assembled structures from small-angle scattering measurements and then input the reconstructed structures to a finite-difference time-domain method to predict their color and reflectance. This computational approach is successfully validated by comparing its predictions against experimentally measured reflectance and provides a pathway for reverse engineering colloidal assemblies with desired optical and photothermal properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

A practical approach to the sensitivity analysis for kinetic Monte Carlo simulation of heterogeneous catalysis

Lattice kinetic Monte Carlo simulations have become a vital tool for predictive quality atomistic understanding of complex surface chemical reaction kinetics over a wide range of reaction conditions. In order to expand their practical value in terms of giving guidelines for atomic level design of catalytic systems, it is very desirable to readily evaluate a sensitivity analysis for a given model. The result of such a sensitivity analysis quantitatively expresses the dependency of the turnover frequency, being the main output variable, on the rate constants entering the model. In the past the application of sensitivity analysis, such as Degree of Rate Control, has been hampered by its exuberant computational effort required to accurately sample numerical derivatives of a property that is obtained from a stochastic simulation method. Here in this study we present an efficient and robust three stage approach that is capable of reliably evaluating the sensitivity measures for stiff microkinetic models as we demonstrate using CO oxidation on RuO 2 (110) as a prototypical reaction. In a first step, we utilize the Fisher Information Matrix for filtering out elementary processes which only yield negligible sensitivity. Then we employ an estimator based on linear response theory for calculating the sensitivity measure for non-critical conditions which covers the majority of cases. Finally we adopt a method for sampling coupled finite differences for evaluating the sensitivity measure of lattice based models. This allows efficient evaluation even in critical regions near a second order phase transition that are hitherto difficult to control. The combined approach leads to significant computational savings over straightforward numerical derivatives and should aid in accelerating the nano scale design of heterogeneous catalysts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Development of an Optical Library for Coevaporated CdSe x Te 1– x

The conversion efficiency of CdTe solar cells may be improved by bandgap engineering, i.e., changing the bandgap value through the addition of Se in the absorber. The Se alloying enables a short-circuit current density improvement, as it leads to a bandgap energy value decrease. Furthermore, it has been associated with increased minority carrier lifetimes, assuring high open-circuit voltage values. An Se gradient profile control can further optimize the solar cell performance. Thus, an optical model baseline of the CdSe x Te 1–x (CST) compound was developed. Spectroscopic ellipsometry measurements were conducted to accurately extract the optical constants of ten CST layers deposited through coevaporation with x varying from 0 to 1. Using the measured dielectric function spectra from the discrete CST layers with varying x, and considering the composition-induced shift in the critical point energies, an energy-shift model was employed to develop the accurate optical library for the CST compound for any x value to provide data for future modeling and optimization. Furthermore, the library accuracy was validated through optical simulations of the quantum efficiency of a graded CST solar cell using the finite-difference time-domain method by replicating the Se profile in the absorber layer measured through secondary ion mass spectrometry.

14 SOLAR ENERGY↗
DOE OSTI2017

Multistep inversion workflow for 3D long-offset damped elastic waves in the Fourier domain

We present a new workflow for imaging damped three-dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short-offset and low-frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine-scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite-difference frequency-domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full-waveform inversion workflow is applied to the 3D Marmousi-2 and SEG/EAGE Salt models with long-offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one-third of the maximum offset length at a resolution consistent with the inverted frequencies.

3D full waveform inversion↗
DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗
NASA NTRS1972

Convection in the tanks of a rotating spacecraft

A general study of convection and mixing of a stratified fluid in a rotating container is described, with application to the special problem of fluid heating and convection in a spacecraft tank. The analysis was based on a set of approximate equations for the Navier-Stokes description of fluid convection with small density variations in a rotating system, including effects of body forces due to temperature stratification (caused by a heater) and arbitrary time-dependent rotation of the tank about a noncentral axis. An efficient numerical finite difference scheme and computational method are described for the convection of vorticity and energy in a two-dimensional tank. Special procedures were developed for analysis of the thermodynamic states resulting from the approximate flow equations derived for small density variations. Results of the numerical simulation are presented for studying the effectiveness of rotation maneuvers in mixing stratified oxygen in the tanks of an Apollo spacecraft. Significant effects of the rotation maneuvers are discussed, including the reduction of the potential for pressure collapse.

Source record↗
NASA NTRS1972

Effects of circular geometry in simulation of convection in rotating spacecraft tanks

For computational simulation of the convection and mixing of stratified fluid in a rotating tank (such as used in Apollo flights) with time-dependent rotation, the Navier-Stokes convection problem was formulated for a circular tank configuration. The final equations results from a general approximate theory for combined forced and contained natural convection in a time-dependent rotating system. The equations are cast in terms of vorticity and stream function in a form convenient for computation, with a transformed coordinate system, and appropriate boundary conditions are derived. Accurate representations for the cryogenic supercritical oxygen thermodynamic properties are used in the computations, and an efficient numerical finite difference scheme and computational method are employed.

Martin, E. D.↗
NASA NTRS1975

Computation of three-dimensional, inviscid supersonic flows

The paper sets forth in detail a method for the finite-difference computation of three-dimensional supersonic fields in an Eulerian mesh. First-, second-, and third-order finite difference schemes are examined. Attention is given to proper treatment of the impermeable and permeable boundaries encompassing the computational plane. Numerical results are presented for certain specific configurations: a conical wing-body combination, internal corner flow, a two-dimensional blunt body, an interfering shock problem, and three-dimensional inviscid supersonic flow past a shuttle-orbiter type vehicle.

Kutler, P.↗
NASA NTRS1978

Theoretical and experimental study of the drag of multielement airfoils

The viscous/potential flow past single-element and multielement airfoils is studied theoretically and experimentally. A computerized analysis, based on iteratively coupled potential-flow and boundary-layer analysis, is used to predict the flow field of the airfoil. The method yields detailed characteristics of conventional laminar and turbulent boundary layers, turbulent wakes, and confluent boundary layers. The viscous flows are analyzed with a method that uses finite-difference solutions of the boundary-layer equations. Reynolds stress in the boundary layers and wakes is simulated with eddy viscosity models for the various flow zones. The viscous calculations are carried into the wake of the airfoil where the drag is found from the defect in the wake momentum.

Olson, L. E.↗