Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “1D”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

High-Resolution Genuinely Multidimensional Solution of Conservation Laws by the Space-Time Conservation Element and Solution Element Method

In this overview paper, we review the basic principles of the method of space-time conservation element and solution element for solving the conservation laws in one and two spatial dimensions. The present method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. In contrast to the modern upwind schemes, the approach here does not use the Riemann solver and the reconstruction procedure as the building blocks. The drawbacks of the upwind approach, such as the difficulty of rationally extending the 1D scalar approach to systems of equations and particularly to multiple dimensions is here contrasted with the uniformity and ease of generalization of the Conservation Element and Solution Element (CE/SE) 1D scalar schemes to systems of equations and to multiple spatial dimensions. The assured compatibility with the simplest type of unstructured meshes, and the uniquely simple nonreflecting boundary conditions of the present method are also discussed. The present approach has yielded high-resolution shocks, rarefaction waves, acoustic waves, vortices, ZND detonation waves, and shock/acoustic waves/vortices interactions. Moreover, since no directional splitting is employed, numerical resolution of two-dimensional calculations is comparable to that of the one-dimensional calculations. Some sample applications displaying the strengths and broad applicability of the CE/SE method are reviewed.

Himansu, Ananda↗
NASA NTRS2002

Investigation of N2O Production from 266 and 532 nm Laser Flash Photolysis of O3/N2/O2 Mixtures

Tunable diode laser absorption spectroscopy has been employed to measure the amount of N2O produced from laser flash photolysis of O3/N2/O2 mixtures at 266 and 532 nm. In the 532 nm photolysis experiments very little N2O is observed, thus allowing an upper limit yield of 7 x 10(exp -8) to be established for the process O3 + N2 yield N2O + O2, where O3 is nascent O3 that is newly formed via O(3P(sub J)) + O2 recombination (with vibrational excitation near the dissociation energy of O3). The measured upper limit yield is a factor of approx. 600 smaller than a previous literature value and is approximately a factor of 10 below the threshold for atmospheric importance. In the 266 nm photolysis experiments, significant N2O production is observed and the N2O quantum yield is found to increase linearly with pressure over the range 100 - 900 Torr in air bath gas. The source of N2O in the 266 nm photolysis experiments is believed to be the addition reaction O(1D(sub 2)) + N2 + M yields (k(sub sigma)) N2O + M, although reaction of (very short-lived) electronically excited O3 with N2 cannot be ruled out by the available data. Assuming that all observed N2O comes from the O(1D(sub 2)) + N2 + M reaction, the following expression describes the temperature dependence of k(sub sigma) (in its third-order low-pressure limit) that is consistent with the N2O yield data: k(sub sigma) = (2.8 +/- 0.1) x 10(exp -36)(T/300)(sup -(0-88+0.36)) cm(sup 6) molecule(sup -2)/s, where the uncertainties are 2(sigma) and represent precision only. The accuracy of the reported rate coefficients at the 95% confidence level is estimated to be 30 - 40% depending on the temperature. Model calculations suggest that gas phase processes initiated by ozone absorption of a UV photon represent about 1.4% of the currently estimated global source strength of atmospheric N2O. However, these processes could account for a significant fraction of the oxygen mass-independent enrichment observed in atmospheric N2O, and they appear to be the first suggested photochemical mechanism that is capable of explaining the altitude dependence of the observed mass -independent isotopic signature.

Estupinan, E. G.↗
NASA NTRS2003

A One Billion Year Martian Climate Model: The Importance of Seasonally Resolved Polar Caps and the Role of Wind

Wind deflation and deposition are powerful agents of surface change in the present Mars climate regime. Recent studies indicate that, while the distribution of regions of potential deflation (or erosion) and deposition is remarkably insensitive to changes in orbital parameters (obliquity, timing of perihelion passage, etc.), rates of aeolian surface modification may be highly sensitive to these parameters even if the atmospheric mass remains constant. But previous work suggested the atmospheric mass is likely to be sensitive to obliquity, especially if a significant mass of carbon dioxide can be stored in the regolith or deposited in the form of massive polar caps. Deflation and erosion are highly sensitive to surface pressure, so feedback between orbit variations and surface pressure can greatly enhance the sensitivity of aeolian modification rates to orbital parameters. We used statistics derived from a 1 Gyr orbital integration of the spin axis of Mars, coupled with 3D general circulation models (GCMs) at a variety of orbital conditions and pressures, to explore this feedback. We also employed a seasonally resolved 1D energy balance model to illuminate the gross characteristics of the longterm atmospheric evolution, wind erosion and deposition over one billion years. We find that seasonal polar cycles have a critical influence on the ability for the regolith to release CO2 at high obliquities, and find that the atmospheric CO2 actually decreases at high obliquities due to the cooling effect of polar deposits at latitudes where seasonal caps form. At low obliquity, the formation of massive, permanent polar caps depends critically on the values of the frost albedo, A(sub frost), and frost emissivity, E(sub frost). Using our 1D model with values of A(sub frost) = 0.67 and E(sub frost) = 0.55, matched to the NASA Ames GCM results, we find that permanent caps only form at low obliquities (< 10 degrees). Thus, contrary to expectations, the Martian atmospheric pressure is remarkable static over time, and decreases both at high and low obliquity. Also, from our one billion year orbital model, we present new results on the fraction of time Mars is expected to experience periods of high and low obliquity. Finally, using GCM runs at a variety of pressures, we examine the likely role of wind erosion under an early more massive Martian atmosphere.

J C Armstrong↗
NASA NTRS2003

Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations: A New Approach Based on the Space-Time CE/SE Method

This paper reports on a significant advance in the area of non-reflecting boundary conditions (NRBCs) for unsteady flow computations. As a part of the development of the space-time conservation element and solution element (CE/SE) method, sets of NRBCs for 1D Euler problems are developed without using any characteristics-based techniques. These conditions are much simpler than those commonly reported in the literature, yet so robust that they are applicable to subsonic, transonic and supersonic flows even in the presence of discontinuities. In addition, the straightforward multidimensional extensions of the present 1D NRBCs have been shown numerically to be equally simple and robust. The paper details the theoretical underpinning of these NRBCs, and explains their unique robustness and accuracy in terms of the conservation of space-time fluxes. Some numerical results for an extended Sod's shock-tube problem, illustrating the effectiveness of the present NRBCs are included, together with an associated simple Fortran computer program. As a preliminary to the present development, a review of the basic CE/SE schemes is also included.

Chang, Sin-Chung↗
NASA NTRS2003

Evidence for Three-Dimensional Radiative Effects in MODIS Cloud Optical Depths Retrieved at Back Scattering View Angles

This study addresses the question whether 1D radiative transfer theory describes well the angular distribution of shortwave cloud reflection. The statistical analysis of a large set of MODIS observations indicates that in oblique backward scattering directions, cloud reflection is stronger than 1D theory would predict. After considering a variety of possible causes, the paper concludes that the most likely reason for the increase lies in 3D radiative interactions. The results' main implication is that cloud optical depths retrieved at back scattering view angles larger than about 50 degrees tend to be overestimated and should be used only with great caution.

Varnai, Tamas↗
NASA NTRS2004

Advanced Applications of Adifor 3.0 for Efficient Calculation of First-and Second-Order CFD Sensitivity Derivatives

This final report will document the accomplishments of the work of this project. 1) The incremental-iterative (II) form of the reverse-mode (adjoint) method for computing first-order (FO) aerodynamic sensitivity derivatives (SDs) has been successfully implemented and tested in a 2D CFD code (called ANSERS) using the reverse-mode capability of ADIFOR 3.0. These preceding results compared very well with similar SDS computed via a black-box (BB) application of the reverse-mode capability of ADIFOR 3.0, and also with similar SDs calculated via the method of finite differences. 2) Second-order (SO) SDs have been implemented in the 2D ASNWERS code using the very efficient strategy that was originally proposed (but not previously tested) of Reference 3, Appendix A. Furthermore, these SO SOs have been validated for accuracy and computational efficiency. 3) Studies were conducted in Quasi-1D and 2D concerning the smoothness (or lack of smoothness) of the FO and SO SD's for flows with shock waves. The phenomenon is documented in the publications of this study (listed subsequently), however, the specific numerical mechanism which is responsible for this unsmoothness phenomenon was not discovered. 4) The FO and SO derivatives for Quasi-1D and 2D flows were applied to predict aerodynamic design uncertainties, and were also applied in robust design optimization studies.

Taylor, Arthur C., III↗
NASA NTRS1995

A Portable 3D FFT Package for Distributed-Memory Parallel Architectures

A parallel algorithm for 3D FFTs is implemented as a series of local 1D FFTs combined with data transposes. This allows the use of vendor supplied (often fully optimized) sequential 1D FFTs. The FFTs are carried out in-place by using an in-place data transpose across the processors.

fast Fourrier transform FFT parallel processing↗
NASA NTRS2007

3D Aerosol-Cloud Radiative Interaction Observed in Collocated MODIS and ASTER Images of Cumulus Cloud Fields

3D aerosol-cloud interaction is examined by analyzing two images containing cumulus clouds in biomass burning regions in Brazil. The research consists of two parts. The first part focuses on identifying 3D clo ud impacts on the reflectance of pixel selected for the MODIS aerosol retrieval based purely on observations. The second part of the resea rch combines the observations with radiative transfer computations to identify key parameters in 3D aerosol-cloud interaction. We found that 3D cloud-induced enhancement depends on optical properties of nearb y clouds as well as wavelength. The enhancement is too large to be ig nored. Associated biased error in 1D aerosol optical thickness retrie val ranges from 50% to 140% depending on wavelength and optical prope rties of nearby clouds as well as aerosol optical thickness. We caution the community to be prudent when applying 1D approximations in comp uting solar radiation in dear regions adjacent to clouds or when usin g traditional retrieved aerosol optical thickness in aerosol indirect effect research.

Wen, Guoyong↗
NASA NTRS2011

Study of Unsteady, Sphere-Driven, Shock-Induced Combustion for Application to Hypervelocity Airbreathing Propulsion

A premixed, shock-induced combustion engine has been proposed in the past as a viable option for operating in the Mach 10 to 15 range in a single stage to orbit vehicle. In this approach, a shock is used to initiate combustion in a premixed fuel/air mixture. Apparent advantages over a conventional scramjet engine include a shorter combustor that, in turn, results in reduced weight and heating loads. There are a number of technical challenges that must be understood and resolved for a practical system: premixing of fuel and air upstream of the combustor without premature combustion, understanding and control of instabilities of the shock-induced combustion front, ability to produce sufficient thrust, and the ability to operate over a range of Mach numbers. This study evaluated the stability of the shock-induced combustion front in a model problem of a sphere traveling in a fuel/air mixture at high Mach numbers. A new, rapid analysis method was developed and applied to study such flows. In this method the axisymmetric, body-centric Navier-Stokes equations were expanded about the stagnation streamline of a sphere using the local similarity hypothesis in order to reduce the axisymmetric equations to a quasi-1D set of equations. These reduced sets of equations were solved in the stagnation region for a number of flow conditions in a premixed, hydrogen/air mixture. Predictions from the quasi-1D analysis showed very similar stable or unstable behavior of the shock-induced combustion front as compared to experimental studies and higher-fidelity computational results. This rapid analysis tool could be used in parametric studies to investigate effects of fuel rich/lean mixtures, non-uniformity in mixing, contaminants in the mixture, and different chemistry models.

Axdahl, Erik↗
NASA NTRS2012

2D Radiative Processes Near Cloud Edges

Because of the importance and complexity of dynamical, microphysical, and radiative processes taking place near cloud edges, the transition zone between clouds and cloud free air has been the subject of intense research both in the ASR program and in the wider community. One challenge in this research is that the one-dimensional (1D) radiative models widely used in both remote sensing and dynamical simulations become less accurate near cloud edges: The large horizontal gradients in particle concentrations imply that accurate radiative calculations need to consider multi-dimensional radiative interactions among areas that have widely different optical properties. This study examines the way the importance of multidimensional shortwave radiative interactions changes as we approach cloud edges. For this, the study relies on radiative simulations performed for a multiyear dataset of clouds observed over the NSA, SGP, and TWP sites. This dataset is based on Microbase cloud profiles as well as wind measurements and ARM cloud classification products. The study analyzes the way the difference between 1D and 2D simulation results increases near cloud edges. It considers both monochromatic radiances and broadband radiative heating, and it also examines the influence of factors such as cloud type and height, and solar elevation. The results provide insights into the workings of radiative processes and may help better interpret radiance measurements and better estimate the radiative impacts of this critical region.

Varnai, T.↗
NASA NTRS2012

A Sounding Rocket Mission Concept to Acquire High-Resolution Radiometric Spectra Spanning the 9 nm - 31 nm Wavelength Range

When studying Solar Extreme Ultraviolet (EUV) emissions, both single-wavelength, two- dimensional (2D) spectroheliograms and multi-wavelength, one-dimensional (1D) line spectra are important, especially for a thorough understanding of the complex processes in the solar magnetized plasma from the base of the chromosphere through the corona. 2D image data are required for a detailed study of spatial structures, whereas radiometric (i.e., spectral) data provide information on relevant atomic excitation/ionization state densities (and thus temperature). Using both imaging and radiometric techniques, several satellite missions presently study solar dynamics in the EUV, including the Solar Dynamics Observatory (SDO), Hinode, and the Solar-Terrestrial Relations Observatory (STEREO). The EUV wavelengths of interest typically span 9 nm to 31 nm, with the shorter wavelengths being associated with the hottest features (e.g., intense flares and bright points) and the longer wavelengths associated with cooler features (e.g., coronal holes and filaments). Because the optical components of satellite instruments degrade over time, it is not uncommon to conduct sounding rocket underflights for calibration purposes. The authors have designed a radiometric sounding rocket payload that could serve as both a calibration underflight for and a complementary scientific mission to the upcoming Solar Ultraviolet Imager (SUVI) mission aboard the GOES-R satellite (scheduled for a 2015 launch). The challenge to provide quality radiometric line spectra over the 9-31 nm range covered by SUVI was driven by the multilayer coatings required to make the optical components, including mirrors and gratings, reflective over the entire range. Typically, these multilayers provide useful EUV reflectances over bandwidths of a few nm. Our solution to this problem was to employ a three-telescope system in which the optical components were coated with multilayers that spanned three wavelength ranges to cover the three pairs of SUVI bands. The complete system was designed to fit within the Black Brandt-IX 22.-diameter payload skin envelope. The basic optical path is that of a simple parabolic telescope in which EUV light is focused onto a slit and shutter assembly and imaged onto a normal-incidence diffraction grating, which then disperses the light onto a 2048 2048 CCD sensor. The CCD thus records 1D spatial information along one axis and spectral information along the other. The slit spans 40 arc-minutes in length, thus covering a solar diameter out to +/- 1.3 solar radii. Our operations concept includes imaging at three distinct positions: the north-south meridian, the northeast-southwest diagonal, and real-time pointing at an active region. Six 10-second images will be obtained at each position. Fine pointing is provided by the SPARCS-VII attitude control system typically employed on Black Brandt solar missions. Both before and after launch, all three telescopes will be calibrated with the EUV line emission source and monochromater system at NASA's Stray Light Facility at Marshall Spaceflight Center. Details of the payload design, operations concept, and data application will be presented.

Krause, L. Habash↗
NASA NTRS2013

Understanding the Flow Physics of Shock Boundary-Layer Interactions Using CFD and Numerical Analyses

Computational fluid dynamic (CFD) analyses of the University of Michigan (UM) Shock/Boundary-Layer Interaction (SBLI) experiments were performed as an extension of the CFD SBLI Workshop held at the 48th AIAA Aerospace Sciences Meeting in 2010. In particular, the UM Mach 2.75 Glass Tunnel with a semi-spanning 7.75deg wedge was analyzed in attempts to explore key physics pertinent to SBLI's, including thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses were 3D CFD simulations using the OVERFLOW flow solver, with additional quasi-1D simulations performed with an in house MATLAB code interfacing with the NIST REFPROP code to explore perfect verses non-ideal air. A fundamental exploration pertaining to the effects of particle image velocimetry (PIV) on post-processing data is also shown. Results from the CFD simulations showed an improvement in agreement with experimental data with key contributions including adding a laminar zone upstream of the wedge and the necessity of mimicking PIV particle lag for comparisons. Results from the quasi-1D simulation showed that there was little difference between perfect and non-ideal air for the configuration presented.

Friedlander, David J.↗
NASA NTRS2014

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗
NASA NTRS2013

Improving NASA's Multiscale Modeling Framework for Tropical Cyclone Climate Study

One of the current challenges in tropical cyclone (TC) research is how to improve our understanding of TC interannual variability and the impact of climate change on TCs. Recent advances in global modeling, visualization, and supercomputing technologies at NASA show potential for such studies. In this article, the authors discuss recent scalability improvement to the multiscale modeling framework (MMF) that makes it feasible to perform long-term TC-resolving simulations. The MMF consists of the finite-volume general circulation model (fvGCM), supplemented by a copy of the Goddard cumulus ensemble model (GCE) at each of the fvGCM grid points, giving 13,104 GCE copies. The original fvGCM implementation has a 1D data decomposition; the revised MMF implementation retains the 1D decomposition for most of the code, but uses a 2D decomposition for the massive copies of GCEs. Because the vast majority of computation time in the MMF is spent computing the GCEs, this approach can achieve excellent speedup without incurring the cost of modifying the entire code. Intelligent process mapping allows differing numbers of processes to be assigned to each domain for load balancing. The revised parallel implementation shows highly promising scalability, obtaining a nearly 80-fold speedup by increasing the number of cores from 30 to 3,335.

tropical cyclone interannual variability↗
NASA NTRS2014

Evaluation of Dielectric-Barrier-Discharge Actuator Substrate Materials

A key, enabling element of a dielectric barrier discharge (DBD) actuator is the dielectric substrate material. While various investigators have studied the performance of different homogeneous materials, most often in the context of related DBD experiments, fundamental studies focused solely on the dielectric materials have received less attention. The purpose of this study was to conduct an experimental assessment of the body-force-generating performance of a wide range of dielectric materials in search of opportunities to improve DBD actuator performance. Materials studied included commonly available plastics and glasses as well as a custom-fabricated polyimide aerogel. Diagnostics included static induced thrust, electrical circuit parameters for 2D surface discharges and 1D volume discharges, and dielectric material properties. Lumped-parameter circuit simulations for the 1D case were conducted showing good correspondence to experimental data provided that stray capacitances are included. The effect of atmospheric humidity on DBD performance was studied showing a large influence on thrust. The main conclusion is that for homogeneous, dielectric materials at forcing voltages less than that required for streamer formation, the material chemical composition appears to have no effect on body force generation when actuator impedance is properly accounted for.

Wilkinson, Stephen P.↗
NASA NTRS2009

Power Spectral Density Specification and Analysis of Large Optical Surfaces

The 2-dimensional Power Spectral Density (PSD) can be used to characterize the mid- and the high-spatial frequency components of the surface height errors of an optical surface. We found it necessary to have a complete, easy-to-use approach for specifying and evaluating the PSD characteristics of large optical surfaces, an approach that allows one to specify the surface quality of a large optical surface based on simulated results using a PSD function and to evaluate the measured surface profile data of the same optic in comparison with those predicted by the simulations during the specification-derivation process. This paper provides a complete mathematical description of PSD error, and proposes a new approach in which a 2-dimentional (2D) PSD is converted into a 1-dimentional (1D) one by azimuthally averaging the 2D-PSD. The 1D-PSD calculated this way has the same unit and the same profile as the original PSD function, thus allows one to compare the two with each other directly.

wavefront↗
NASA NTRS2015

On Hilbert-Huang Transform Based Synthesis of a Signal Contaminated by Radio Frequency Interference or Fringes

Norden E. Huang et al. had proposed and published the Hilbert-Huang Transform (HHT) concept correspondently in 1996, 1998. The HHT is a novel method for adaptive spectral analysis of non-linear and non-stationary signals. The HHT comprises two components: - the Huang Empirical Mode Decomposition (EMD), resulting in an adaptive data-derived basis of Intrinsic Mode functions (IMFs), and the Hilbert Spectral Analysis (HSA1) based on the Hilbert Transform for 1-dimension (1D) applied to the EMD IMF's outcome. Although paper describes the HHT concept in great depth, it does not contain all needed methodology to implement the HHT computer code. In 2004, Semion Kizhner and Karin Blank implemented the reference digital HHT real-time data processing system for 1D (HHT-DPS Version 1.4). The case for 2-Dimension (2D) (HHT2) proved to be difficult due to the computational complexity of EMD for 2D (EMD2) and absence of a suitable Hilbert Transform for 2D spectral analysis (HSA2). The real-time EMD2 and HSA2 comprise the real-time HHT2. Kizhner completed the real-time EMD2 and the HSA2 reference digital implementations respectively in 2013 & 2014. Still, the HHT2 outcome synthesis remains an active research area. This paper presents the initial concepts and preliminary results of HHT2-based synthesis and its application to processing of signals contaminated by Radio-Frequency Interference (RFI), as well as optical systems' fringe detection and mitigation at design stage. The Soil Moisture Active Passive (SMAP mission (SMAP) carries a radiometer instrument that measures Earth soil moisture at L1 frequency (1.4 GHz polarimetric - H, V, 3rd and 4th Stokes parameters). There is abundant RFI at L1 and because soil moisture is a strategic parameter, it is important to be able to recover the RFI-contaminated measurement samples (15% of telemetry). State-of-the-art only allows RFI detection and removes RFI-contaminated measurements. The HHT-based analysis and synthesis facilitates recovery of measurements contaminated by all kinds of RFI, including jamming [7-8]. The fringes are inherent in optical systems and multi-layer complex contour expensive coatings are employed to remove the unwanted fringes. HHT2-based analysis allows test image decomposition to analyze and detect fringes, and HHT2-based synthesis of useful image.

Spacecraft Electronics Equipment↗
NASA NTRS2005

Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations - A New Approach based on the Space-Time CE/SE Method

This paper reports on a significant advance in the area of nonreflecting boundary conditions (NRBCs) for unsteady flow computations. As a part of t he development of t he space-time conservation element and solution element (CE/SE) method, sets of NRBCs for 1D Euler problems are developed without using any characteristics- based techniques. These conditions are much simpler than those commonly reported in the literature, yet so robust that they are applicable to subsonic, transonic and supersonic flows even in the presence of discontinuities. In addition, the straightforward multidimensional extensions of the present 1D NRBCs have been shown numerically to be equally simple and robust. The paper details the theoretical underpinning of these NRBCs, and explains t heir unique robustness and accuracy in terms of t he conservation of space-time fluxes. Some numerical results for an extended Sod's shock-tube problem, illustrating the effectiveness of the present NRBCs are included, together with an associated simple Fortran computer program. As a preliminary to the present development, a review of the basic CE/SE schemes is also included.

Chang, S.-C.↗