Search NASA⌕ Search

SEARCH · Search NASA

Results for “Overlapping grid”

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 181 records · Page 10

Summertime Coincident Observations of Ice Water Path in the Visible/Near-IR, Radar, and Microwave Frequencies

Accurate representation of the physical and radiative properties of clouds in climate models continues to be a challenge. At present, both remote sensing observations and modeling of microphysical properties of clouds rely heavily on parameterizations or assumptions on particle size distribution (PSD) and cloud phase. In this study, we compare Ice Water Path (IWP), an important physical and radiative property that provides the amount of ice present in a cloud column, using measurements obtained via three different retrieval strategies. The datasets we use in this study include Visible/Near-IR IWP from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument flying aboard the Aqua satellite, Radar-only IWP from the CloudSat instrument operating at 94 GHz, and NOAA/NESDIS operational IWP from the 89 and 157 GHz channels of the Microwave Humidity Sounder (MHS) instrument flying aboard the NOAA-18 satellite. In the Visible/Near-IR, IWP is derived from observations of optical thickness and effective radius. CloudSat IWP is determined from measurements of cloud backscatter and assumed PSD. MHS IWP retrievals depend on scattering measurements at two different, non-water absorbing channels, 89 and 157 GHz. In order to compare IWP obtained from these different techniques and collected at different vertical and horizontal resolutions, we examine summertime cases in the tropics (30S - 30N) when all 3 satellites are within 4 minutes of each other (approximately 1500 km). All measurements are then gridded to a common 15 km x 15 km box determined by MHS. In a grid box comparison, we find CloudSat to report the highest IWP followed by MODIS, followed by MHS. In a statistical comparison, probability density distributions show MHS with the highest frequencies at IWP of 100-1000 g/m(exp 2) and CloudSat with the longest tail reporting IWP of several thousands g/m(exp 2). For IWP greater than 30 g/m(exp 2), MODIS is consistently higher than CloudSat, and it is higher at the lower IWPs but lower at the higher IWPs that overlap with MHS. Some of these differences can be attributed to the limitations of the measuring techniques themselves, but some can result from the assumptions made in the algorithms that generate the IWP product. We investigate this issue by creating categories based on various conditions such as cloud type, precipitation presence, underlying liquid water content, and surface type (land vs. ocean) and by comparing the performance of the IWP products under each condition.

Pittman, Jasna V.↗

A Hybrid Line List for CH4 and Hot Methane Continuum

Molecular line lists (a catalogue of transition frequencies and line strengths) are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here we suggest a new, hybrid line list format to mitigate this problem, based on the idea of temperature-dependent absorption continuum. Methods. The line list is partitioned into a large set of relatively weak lines and a small set of important, stronger lines. The weaker lines are then used either to construct a temperature-dependent (but pressure-independent) set of intensity cross sections or are blended into a greatly reduced set of super-lines. The strong lines are kept in the form of temperature independent Einstein A coefficients. Results. A line list for methane is constructed as a combination of 17 million strong absorption lines relative to the reference absorption spectra and a background methane continuum in two temperature-dependent forms, of cross sections and super-lines. This approach eases the use of large high temperature line lists significantly as the computationally expensive calculation of pressure dependent profiles only need to be performed for a relatively small number of lines. Both the line list and cross sections were generated using a new 34 billion methane line list (34 to10), which extends the 10to10 line list to higher temperatures (up to 2000 K). The new hybrid scheme can be applied to any large line lists containing billions of transitions. We recommend to use super-lines generated on a high resolution grid based on resolving power R = 1,000,000 to model the molecular continuum as a more flexible alternative to the temperature dependent cross sections.

Line Profiles↗

NASCAP/LEO calculations of current collection

NASCAP/LEO is a 3-dimensional computer code for calculating the interaction of a high-voltage spacecraft with the cold dense plasma found in Low Earth Orbit. Although based on a cubic grid structure, NASCAP/LEO accepts object definition input from standard computer aided design (CAD) programs so that a model may be correctly proportioned and important features resolved. The potential around the model is calculated by solving the finite element formulation of Poisson's equation with an analytic space charge function. Five previously published NASCAP/LEO calculations for three ground test experiments and two space flight experiments are presented. The three ground test experiments are a large simulated panel, a simulated pinhole, and a 2-slit experiment with overlapping sheaths. The two space flight experiments are a solar panel biased up to 1000 volts, and a rocket-mounted sphere biased up to 46 kilovolts. In all cases, the authors find good agreement between calculation and measurement.

Mandell, Myron J.↗

Regimes of Cloud Vertical Structure From Active Observations

We introduce a new type of cloud class, which we call “active cloud regime” (ACR), owing to its provenance from active (lidar and cloud radar) spaceborne cloud observations. ACRs intend to provide a climatological description based on cloud vertical structure (CVS) of the most prevalent monthly CVS mixtures encountered at large spatial scales of ~400 km. ACRs are thus a way to create a gridded dataset of a vertically-resolved cloud mask that can facilitate joint analysis with other gridded datasets. The detailed 2D cloud mask comes from the 2B-CLDCLASS-LIDAR CloudSat dataset fusing CALIPSO (lidar) and CloudSat (cloud radar) cloud detections. We show that the global classification of cloudiness under the ACR framework provides valuable insights on how the world’s dominant cloud systems regulate the two major components of atmospheric energetics, precipitation and radiative cooling. NASA’s GEOS model allows us to demonstrate the feasibility of applying the ACR concept in Earth System Models that have the capability to produce subgrid cloudiness obeying pre-specified vertical overlap rules. Comparison of observed and simulated ACRs provides thus another means to assess the realism of modeled clouds.

cloud vertical structure↗

Generating Solid Models from Topographical Data

A method of generating solid models of terrain involves the conversion of topographical data into a form useable by a rapid-prototyping (RP) machine. The method was developed to enable the use of the RP machine to make solid models of Martian terrain from Mars Orbiter laser-altimeter topographical data. The method is equally applicable to the generation of models of the terrains of other astronomical bodies, including other planets, asteroids, and Earth. Topographical data describe a terrain in terms of a set of three-dimensional coordinates [e.g., Cartesian (x,y,z) or polar (latitude, longitude, radius) coordinates] of points or nodes on the terrain surface. The input data for the RP machines are required to provide a three-dimensional description, not of a single surface, but of a volume in this case, a ground volume that underlies the terrain surface. The description is required to be in the form of triangular elements that connect the nodes of all the surfaces and that completely bound the volume, with no open areas, no overlap of triangles, and no extraneous geometric elements. The software used in the present model-generation method was written in IDL - an advanced programming language that affords a number of tools, including subroutines that triangularize surfaces. The software creates a volume from the topographical surface data by adding sides to the edges of the terrain surface and joining the sides with a bottom surface. Each of the sides is triangularized by use of IDL subroutines, and then the software searches for extraneous elements and removes them. Topographical data are usually presented in a grid corresponding to polar coordinates, so that a model generated from such data is equivalent to a topographical map in Mercator projection. However an RP machine is fully capable of including the curvature of a planetary body in a model that it makes. Therefore, the software also offers a capability to transform the topographical data to a projection onto a surface having a curvature corresponding to that of the surface of the modeled planet.

Keller, John W.↗

The Surface Longwave Downward Fluxes of the NASA GEWEX SRB Release 4.0 IP Products: Validation Against the Surface-Based BSRN and PMEL Observed Data

Since the NASA Global Energy Water Exchanges (GEWEX) Surface Radiation Budget (SRB) project released its 3rd version of products in 2010, the GEWEX Data Assessments Panel (GDAP) has been working on integrating various data products to address issues in the closing of the global energy and water cycles. The 4th version of the SRB products, Rel. 4.0-IP, has integrated data products from the cloud, aerosol, atmosphere, ocean surface, and land surface projects, coordinating within GDAP, to produce a long-term time series of TOA and surface radiative estimates. The Rel. 4-IP shortwave products span 34 years continuously from July 1983 to June 2017 on a quasi-equal-area 1degree longitude by 1degree latitude grid system. The longwave products are for land only from 1983-07 to 1987-12; for both land and ocean from 1988-01 to 2009-12; and ocean only from 2010-01 to 2017-06. The data are provided at 3 hourly, 3-hourly-monthly, daily and monthly means. The ISCCP HXS clouds and radiances are the key cloud input of the current GEWEX SRB algorithms. In addition, the longwave algorithm has also made changes in cloud microphysical property, surface skin temperature input, surface emissivity, atmospheric profile, adding longwave aerosol optical properties, revising cloud overlap procedure, and so on. Details of changes in both inputs and algorithms are documented in a NASA Algorithm Theoretical Basis Document (ATBD). We have validated the surface longwave downward fluxes against the surface-based Baseline Surface Radiation Network (BSRN) and the Pacific Marine Environmental Laboratory (PMEL) buoy data. As of 2020, the BSRN archive has 12,116 site-months of observed records from 73 stations on all seven continents, and as of 2017, PMEL archive has 4389 buoy months of observed records from 64 buoys deployed in the tropics of Pacific, Atlantic and Indian Oceans. This paper presents how the SRB Rel. 4.0-IP surface longwave downward fluxes compare with these surface-based measurements and how the comparison statistics differ from that of Rel. 3.0.

GEWEX SRB↗

Creating A Consistent Historical NASA POWER Solar Radiation Dataset to Support Renewable Energy, Building Energy Efficiency and Agro-Climatology Decisions

Prediction of Worldwide Energy Resources (POWER) project provides irradiance dataset to support renewable energy, building energy efficiency and agricultural needs. These datasets are derived from Global Energy and Water Cycle Experiment Surface Radiation Budget (GEWEX SRB) and Clouds and the Earth’s Radiant Energy System (CERES SYN1Deg). A systematic bias has been reported between these two datasets for the years with overlapping observations. For obtaining a consistent climate data record spanning the entire time record of observations, it is crucial to understand and remove the bias in the irradiance dataset. Inconsistency in solar radiation data can lead to inaccurate conclusions about solar energy potential and obscure real trends in solar radiation patterns that would impact energy availability assessments. In this study, we adapt quantile mapping approach to remove the systematic bias and to improve reliability of shortwave and longwave irradiance data. We present a validation of the bias corrected data against ground truth. For each 1° latitude and 1° longitude grid box across the globe, we match the CDFs of the reference dataset (CERES SYN1Deg) to that of the SRB dataset, thereby, adjusting the irradiance values to match the empirical distribution of two different measurements. The performance of quantile mapping is evaluated by using the metrics such as Mean Absolute Deviation (MAD). The results indicate that the quantile mapping significantly improves the accuracy and reliability of solar irradiance dataset especially for the weather conditions associated with high cloud cover and extreme irradiance values. The initial range of MAD for the studied sites for daily data was 4 to 19 Wm-2. After correction these reduced to 3 to 7 Wm-2. The findings from this study have important implications for solar energy system design, agricultural planning, and climate modeling community. Reducing the inconsistency and biases in solar irradiance dataset can enable better planning and operation of solar energy systems, leading to increased efficiency and cost-effectiveness. Additionally, this work also contributes to the statistical post-processing techniques in the renewable energy domain and highlights the potential of historical and near-real-time NASA POWER dataset as a valuable resource for solar energy research and applications.

POWER↗

Sea ice feature and type identification in merged ERS-1 SAR and LANDSAT Thematic Mapper imagery

LANDSAT Thematic Mapper (TM) and ERS-1 SAR (Synthetic Aperture Radar) C band images were acquired for the same area in the Beaufort Sea, 18 Apr. 1992. The two images were co-located to the same grid (25 m resolution) and supervised classification was performed on the TM channel 3 scene in order to classify open water, nilas, grey ice, first year ice, and multiyear ice. Comparison of the LANDSAT classification and the corresponding dB values from the SAR scene showed that, under the given circumstances (high surface winds), open water/nilas/grey ice as defined by a single ice category by the SAR classifier could not be distinguished from first year ice. Surface roughening due to wind appears to be a major problem for the SAR classifier, as the range of the dB values is large enough to overlap all the other ice categories. Additional information, such as surface wind speed is necessary to overcome part of this problem.

Steffen, K.↗

Newton-Krylov-Schwarz: An implicit solver for CFD

Newton-Krylov methods and Krylov-Schwarz (domain decomposition) methods have begun to become established in computational fluid dynamics (CFD) over the past decade. The former employ a Krylov method inside of Newton's method in a Jacobian-free manner, through directional differencing. The latter employ an overlapping Schwarz domain decomposition to derive a preconditioner for the Krylov accelerator that relies primarily on local information, for data-parallel concurrency. They may be composed as Newton-Krylov-Schwarz (NKS) methods, which seem particularly well suited for solving nonlinear elliptic systems in high-latency, distributed-memory environments. We give a brief description of this family of algorithms, with an emphasis on domain decomposition iterative aspects. We then describe numerical simulations with Newton-Krylov-Schwarz methods on aerodynamics applications emphasizing comparisons with a standard defect-correction approach, subdomain preconditioner consistency, subdomain preconditioner quality, and the effect of a coarse grid.

Cai, Xiao-Chuan↗

Parallel Processing of Adaptive Meshes with Load Balancing

Many scientific applications involve grids that lack a uniform underlying structure. These applications are often also dynamic in nature in that the grid structure significantly changes between successive phases of execution. In parallel computing environments, mesh adaptation of unstructured grids through selective refinement/coarsening has proven to be an effective approach. However, achieving load balance while minimizing interprocessor communication and redistribution costs is a difficult problem. Traditional dynamic load balancers are mostly inadequate because they lack a global view of system loads across processors. In this paper, we propose a novel and general-purpose load balancer that utilizes symmetric broadcast networks (SBN) as the underlying communication topology, and compare its performance with a successful global load balancing environment, called PLUM, specifically created to handle adaptive unstructured applications. Our experimental results on an IBM SP2 demonstrate that the SBN-based load balancer achieves lower redistribution costs than that under PLUM by overlapping processing and data migration.

Das, Sajal K.↗

Patched Local Lunar Gravity Solutions using GRAIL Data

We present a method to determine local gravity fields for the Moon using Gravity Recovery and Interior Laboratory (GRAIL) data. We express gravity as gridded gravity anomalies on a sphere, and we estimate adjustments to a background global start model expressed in spherical harmonics. We processed GRAIL Ka-band range-rate data with a short-arc approach, using only data over the area of interest. We determine our gravity solutions using neighbor smoothing constraints. We divided the entire Moon into twelve regions and two polar caps, with a resolution of 0.15˚ x 0.15˚ (which is equivalent to degree and order 1199 in spherical harmonics), and determined the optimal smoothing parameter for each area by comparing localized correlations between gravity and topography for each solution set. Our selected areas share nodes with surrounding areas and they are overlapping. To mitigate boundary effects, we patch the solutions together by symmetrically omitting the boundary parts of overlapping solutions. Our new solution has been iterated, and it has improved correlations with topography when compared to a fully iterated global model. Our method requires fewer resources, and can easily handle regionally varying resolution or constraints. The smooth model describes small-scale features clearly, and can be used in local studies of the structure of the lunar crust.

Sander Goossens↗

NASA Surface Radiation Budget Project: A Look Back and A Look Forward

The NASA Langley Surface Radiation Budget (SRB) project, historically associated with GEWEX, has provided estimates of surface and TOA shortwave and longwave radiative fluxes using International Satellite Cloud Climatology Project (ISCCP) cloud properties and radiances as primary input. The first release was in 1995 (Whitlock et al., 1995) at daily 2.5° resolution. Subsequent releases have improved physics, temporal and spatial resolution, as well as increasing the total length of the time series. Release 4 in 2019 (Kummerow et al., 2019) is on a 1° nested grid, with temporal resolutions from 3-hourly to daily and monthly. It originally spanned from July 1983 to December 2017 and has since been extended through 2019. Release 4 uses the ISCCP H series of data, a complete reprocessing which retains more pixels than the previous D series, used in prior SRB versions. SRB data has been used extensively by a wide range of customers, including but not limited to the solar energy community, climate modelers, architects, oceanographers, and the agricultural community. It is an especial feature of SRB that global radiation budget data are available from the 1983 through 2000 period, which is prior to CERES providing alternate data, and then 2000 through near present to provide overlap coverage. Here we use the 36.5 year record of SRB Release 4 to examine the changes in shortwave and longwave fluxes over the period, including the changes in different regions of the globe. We examine energy budget features such as hemispheric symmetry and implied meridional transport of energy. We look at the variability of these features prior to and during the CERES period and compare to later products. We also look at Release 5, currently under development. This version will provide 0.5° resolution. In addition to surface and TOA radiative fluxes it will include fluxes at several atmospheric layers, and some spectral features beyond the broadband shortwave and longwave fluxes in Release 4 and prior. The shortwave algorithm is undergoing a complete change from the Pinker-Laszlo algorithm. Shortwave fluxes will now be calculated from a forward call to the Fu-Liou radiative transfer model as modified by the CERES team (Rose et al., 2006). Input datasets are being improved, including surface properties, atmospheric profiles of temperature and humidity, aerosols, and total solar irradiance. Preliminary results from Release 5 and comparisons to Release 4 will be shown.

Earth Radiation Budget↗

The Global Impact of Clouds on the Production of MODIS Bidirectional Reflectance Model-based Composites for Terrestrial Monitoring

A global data set of cloud occurrence probability derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua gridded daily data is analyzed to investigate the probability of obtaining at least a minimum number of cloud-free observations within various compositing periods. The probabilities derived from Terra and Aqua, with morning and afternoon overpass times, respectively, are similar and increase with compositing period. Compositing both Terra and Aqua observations results is considerably higher probabilities of obtaining a sufficient number of observations for bidirectional reflectance model-based compositing. Given that the only alternative is obtaining sufficient samples to extend the observation period, which can cause significant problems when the surface state changes, it is concluded that using data from the two MODIS sensors provides the most effective way of generating composited products. Findings with respect to the availability of cloud free composites when n-day composites are generated on a temporally overlapping daily rolling basis, i.e., every day, rather than every n-days, are also discussed for regional and global applications.

Roy, D. P.↗

Retrieval of Aerosol Optical Properties From GOCI-II Observations: Continuation of Long-Term Geostationary Aerosol Monitoring Over East Asia

Since the Geostationary Ocean Color Imager (GOCI) was successfully launched in 2010, the GOCI Yonsei aerosol retrieval (YAER) algorithm has been continuously updated to retrieve hourly aerosol optical properties. GOCI-II has 4 more channels including UV, finer spatial resolution (250 m), and daily full disk coverage as compared to GOCI, and was launched in February 2020, onboard the GEO-KOMPSAT-2B (GK-2B) satellite. In this study, we extended the YAER algorithm to GOCI-II data based on its improved performance in many aspects and present the first results of aerosol optical properties retrieved from GOCI-II data. Utilizing the overlapping period between the GOCI-II and GOCI in geostationary Earth orbit, we present GOCI-II aerosol retrievals for high aerosol-loading cases over East Asia and show that these have a consistent spatial distribution with those from GOCI. Furthermore, GOCI-II provides AOD at an even higher spatial resolution, revealing finer changes in aerosol concentrations. Validation results for one year data show that the GOCI-II AOD has a correlation coefficient of 0.83 and a ratio within the expected error (EE) of 59.4 % when compared with the aerosol robotic network (AERONET) data. We compared statistical metrics for the GOCI and GOCI-II AODs to assess the consistency between the two datasets. In addition, it was found that there is a strong correlation between the two datasets from the comparison of gridded GOCI and GOCI-II AOD products. It is expected that data from GOCI-II will continue long-term aerosol records with high accuracy that can be used to address air-quality issues over East Asia.

AOD↗

Fabrication of a Kilopixel Array of Superconducting Microcalorimeters with Microstripline Wiring

A document describes the fabrication of a two-dimensional microcalorimeter array that uses microstrip wiring and integrated heat sinking to enable use of high-performance pixel designs at kilopixel scales (32 X 32). Each pixel is the high-resolution design employed in small-array test devices, which consist of a Mo/Au TES (transition edge sensor) on a silicon nitride membrane and an electroplated Bi/Au absorber. The pixel pitch within the array is 300 microns, where absorbers 290 microns on a side are cantilevered over a silicon support grid with 100-micron-wide beams. The high-density wiring and heat sinking are both carried by the silicon beams to the edge of the array. All pixels are wired out to the array edge. ECR (electron cyclotron resonance) oxide underlayer is deposited underneath the sensor layer. The sensor (TES) layer consists of a superconducting underlayer and a normal metal top layer. If the sensor is deposited at high temperature, the ECR oxide can be vacuum annealed to improve film smoothness and etch characteristics. This process is designed to recover high-resolution, single-pixel x-ray microcalorimeter performance within arrays of arbitrarily large format. The critical current limiting parts of the circuit are designed to have simple interfaces that can be independently verified. The lead-to-TES interface is entirely determined in a single layer that has multiple points of interface to maximize critical current. The lead rails that overlap the TES sensor element contact both the superconducting underlayer and the TES normal metal

Chervenak, James↗

Radiative Effect of Clouds on Tropospheric Chemistry in a Global Three-Dimensional Chemical Transport Model

Clouds exert an important influence on tropospheric photochemistry through modification of solar radiation that determines photolysis frequencies (J-values). We assess the radiative effect of clouds on photolysis frequencies and key oxidants in the troposphere with a global three-dimensional (3-D) chemical transport model (GEOS-CHEM) driven by assimilated meteorological observations from the Goddard Earth Observing System data assimilation system (GEOS DAS) at the NASA Global Modeling and Assimilation Office (GMAO). We focus on the year of 2001 with the GEOS-3 meteorological observations. Photolysis frequencies are calculated using the Fast-J radiative transfer algorithm. The GEOS-3 global cloud optical depth and cloud fraction are evaluated and generally consistent with the satellite retrieval products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the International Satellite Cloud Climatology Project (ISCCP). Results using the linear assumption, which assumes linear scaling of cloud optical depth with cloud fraction in a grid box, show global mean OH concentrations generally increase by less than 6% because of the radiative effect of clouds. The OH distribution shows much larger changes (with maximum decrease of approx.20% near the surface), reflecting the opposite effects of enhanced (weakened) photochemistry above (below) clouds. The global mean photolysis frequencies for J[O1D] and J[NO2] in the troposphere change by less than 5% because of clouds; global mean O3 concentrations in the troposphere increase by less than 5%. This study shows tropical upper tropospheric O3 to be less sensitive to the radiative effect of clouds than previously reported (approx.5% versus approx.20-30%). These results emphasize that the dominant effect of clouds is to influence the vertical redistribution of the intensity of photochemical activity while global average effects remain modest, again contrasting with previous studies. Differing vertical distributions of clouds may explain part, but not the majority, of these discrepancies between models. Using an approximate random overlap or a maximum-random overlap scheme to take account of the effect of cloud overlap in the vertical reduces the impact of clouds on photochemistry but does not significantly change our results with respect to the modest global average effect.

Liu, Hongyu↗

Multiwavelength Modeling of Nova Atmospheres

We have analyzed the early optically thick ultraviolet spectra of Nova OS And 1986 using a grid of spherically symmetric, non-LTE, line-blanketed, expanding model atmospheres and synthetic spectra with the following set of parameters: 5,000 less than or equal to T(sub model) less than or equal to 60,000K, solar abundances, (rho)(alpha) r(sup -3), v(sub max) = 2000 km/s, L = 6 x 10(exp 4) solar luminosity, and a statistical or microturbulent velocity of 50 km/s. We used the synthetic spectra to estimate the model parameters corresponding to the observed IUE spectra. The fits to the observations were then iteratively improved by changing the parameters of the model atmospheres, in particular T(sub model) and the abundances, to arrive at the best fits to the optically thick pseudo-continuum and the features found in the IUE spectra. The IUE spectra show two different optically thick subphases. The earliest spectra, taken a few days after maximum optical light, show a pseudo-continuum created by overlapping absorption lines. The later observations, taken approximately 3 weeks after maximum light, show the simultaneous presence of allowed, semi-forbidden, and forbidden lines in the observed spectra. Analysis of these phases indicate that OS And 86 had solar metallicities except for Mg which showed evidence of being underabundant by as much as a factor of 10. We determine a distance of 5.1 kpc to OS And 86 and derive a peak bolometric luminosity of approximately 5 x 10(exp 4) solar luminosity. The computed nova parameters provide insights into the physics of the early outburst and explain the spectra seen by IUE. Lastly, we find evidence in the later observations for large non-LTE effects of Fe II which, when included, lead to much better agreement with the observations.

Hauschildt, P. H.↗

Flow Solver for Incompressible Rectangular Domains

This is an extension of the Flow Solver for Incompressible 2-D Drive Cavity software described in the preceding article. It solves the Navier-Stokes equations for incompressible flow using finite differencing on a uniform, staggered grid. There is a runtime choice of either central differencing or modified upwinding for the convective term. The domain must be rectangular, but may have a rectangular walled region within it. Currently, the position of the interior region and exterior boundary conditions are changed by modifying parameters in the code and recompiling. These features make it possible to solve a variety of classical fluid flow problems such as an L-shaped cavity, channel flow, or wake flow past a square cylinder. The code uses fourth-order Runge-Kutta time-stepping and overall second-order spatial accuracy. This software permits the walled region to be positioned such that flow past a square cylinder, an L-shaped cavity, and the flow over a back-facing step can all be solved by reconfiguration. Also, this extension has an automatic detection of periodicity, as well as use of specialized data structure for ease of configuring domain decomposition and computing convergence in overlap regions.

Kalb, Virginia L.↗