Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current conservation”

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 397 records · Page 22

New York Ecological Forecasting: Utilizing NASA Earth Observations to Map Ash Distribution and Inform Emerald Ash Borer Control

Since their first sightings in the U.S. in 2002, emerald ash borer beetles (Agrilus planipennis; EAB) have killed millions of native ash (Fraxinus spp.) trees across 35 states. Infected ash stands frequently exhibit complete mortality, with the predicted result being the functional extinction of native ash in U.S. forests. In August of 2020, EAB was discovered in the 6.1-million-acre Adirondack Park. The team’s partners at the Adirondack Park Invasive Plant Program (APIPP) desired ash tree distribution and EAB susceptibility information to help improve EAB bio-control efficiency and apply the methodology to future invasive programs. To assist, the team mapped ash tree distribution using NASA Earth observations from Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Shuttle Radar Topography Mission (SRTM), along with hyperspectral imagery from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). Field data from the Monitoring and Managing Ash (MaMA) project, iMapInvasives and iNaturalist databases, and the New York State Department of Environmental Conservation (NYSDEC) provided ground truthing for mapping and modeling. Results indicate that for ash detection, the team’s Spectral Angle Mapping (SAM) hyperspectral classification is slightly more sensitive but less accurate than multispectral Random Forest (RF) classification, though neither method was above a ~20% detection rate. End products include maps of ash extent derived from both imagery types, a model forecasting future spread scenarios based on current EAB presence, and outreach materials. These products inform APIPP’s management decisions and facilitate public awareness of EAB’s threat to communities within the region.

Liam Megraw↗

Limb brightening on Uranus - An interpretation of the lambda 7300 A methane band

Measurements of limb brightening on the Uranus disk within the lambda 7300 A CH4 band are interpreted using an elementary inhomogeneous radiative transfer model to describe the atmosphere. A two layer model which consists of a finite, optically thin, region of conservatively scattering particles overlying a semi-infinite clear H2-CH4 atmosphere satisfactorily explains the observations. The maximum optical thickness of the upper layer appears to lie in the range 0.1 to 0.2. The CH4/H2 mixing ratio in the lower layer is larger than the corresponding solar value by a factor on the order of three or greater. The results are discussed briefly in terms of current models of the Uranus atmosphere.

Price, M. J.↗

Electromagnetic fluctuations for anisotropic media and the generalized Kirchhoff's law

In this paper the polarimetric emission parameters for anisotropic media are derived using the generalized Kirchhoff's law for media with a uniform temperature and the fluctuation-dissipation theory for media with a temperature profile. Both finite-size objects and half-space media are considered. When the object has a uniform temperature across its body, the Kirchhoff's law, based on the condition of energy conservation in thermal equilibrium is generalized to obtain the emission parameters of an anisotropic medium, which can be interpreted as the absorptivity or the absorption cross section of the complementary object with a permittivity that is the transpose of the original object. When the medium has a nonuniform temperature distribution, the fluctuation-dissipation theory is applied for deriving the covariances between vector components of the thermal currents and, consequently, the covariances of the polarizations of electric fields radiated by the thermal currents. To verify the formulas derived from the fluctuation-dissipation theory, we let the temperature of the object be a constant and show that the results reduce to those obtained from the generalized Kirchhoff's law.

Yueh, Simon H.↗

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗

Detonation to Deflagration-Mode Transitions in Pulsed Plasma Accelerators

PULSED plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas. The current in such discharges will ionize the gas and produce a strong magnetic field, which interacts with the flowing current to accelerate the plasma through the Lorentz body force. For the present work, two plasma accelerator types employing this general scheme are of interest: the gas-fed pulsed plasma thruster (PPT) 1 and the quasi-steady magnetoplasmadynamic (MPD) accelerator. 2 The gas-fed pulsed plasma accelerator is generally understood as a completely transient device discharging in ∼1-10 μs. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j × B body force, entraining and accelerating propellant it encounters. This process is sometimes referred to in literature as 'snowplowing' the propellant or accelerating the gas in a detonation-mode because the current sheet representation approximates that of a strong detonation shockwave propagating through the gas. For these devices, acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new 'crowbar' discharge at the breech and depriving the initial sheet of additional acceleration. In general, PPTs typically claim thrust efficiencies (ratio of jet kinetic energy to input electrical energy) registering in the teens or lower. 3 In the quasi-steady MPD accelerator, the pulse is lengthened to ∼1 ms or longer and maintained at an approximately constant level during discharge through the use of a pulse-forming network (PFN) of capacitors. After an initial transient discharge, which is typically short relative to the overall discharge period, the plasma assumes a relatively steady-state configuration, known as 'quasi-steady' MPD operation. 2 In this state, ionized gas flows through a stationary current channel in a manner that is sometimes referred to as deflagration-mode operation owing to the similarities to deflagration waves in gases. The plasma experiences electromagnetic acceleration as the plasma flows through the current channel towards the exit of the device. Quasi-steady MPD thrusters claim efficiencies up to 50% for certain propellants.4 There has been significant and sustained research over several decades on both gas-fed PPTs and quasi-steady MPD thrusters, however there have been pulsed thrusters that do not appear to exactly fit either classification, instead possessing a mixture of operational qualities characteristic of both thruster variants. The Coaxial High ENerGy (CHENG) thruster by Cheng, et al.5 operated on the short 10 μs timescales characteristic of PPTs, but claimed the high thrust densities, high efficiencies, and low electrode erosion rates that are more consistent with the MPD/deflagration mode of plasma acceleration. Gas-fed PPT research by Ziemer, et al. 3, 6 identified two separate regimes of performance in those thrusters. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency as a function of mass bit, while the second regime at very low mass bits (termed Mode II) exhibited an increase in efficiency with decreasing mass bit. Work by Poehlmann et al.7 and by Sitaraman and Raja8 sought to understand the performance of the CHENG thruster and the Mode I/Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes of plasma acceleration representing two distinct sets of solutions to the relevant conservation laws. In these works, it was proposed that the values of the various controllable parameters determined whether the accelerator would operate in detonation or deflagration mode. Our hypothesized view of the acceleration process in the CHENG thruster and in PPTs experiencing a transition from Mode I to Mode II is inspired by observations of the transition from the PPT mode of operation to the quasi-steady MPD mode. Specifically, the quasi-steady MPD was discovered by driving a PPT to extended pulse lengths. Above a certain pulse length threshold the transient plasma current sheet transitions into a stable plasma acceleration mode that 'replicates in every observable detail steady flow self-field magnetoplasmadynamic acceleration.' 9 In the present work, instead of treating the accelerator as if it were only operating in a single mode during a pulse, we consider the initial stage of the discharge in all cases as a current sheet forming at the breach of the accelerator and moving towards the exit as a detonation wave. If the current sheet reaches the exit of the accelerator before the discharge is completed, the view of the acceleration mode transitions to the deflagration mode-type found in quasi-steady MPD thrusters. In previous work10 we presented a modeling framework that first captured the time-evolution of the current sheet (detonation) mode of the thruster and then transitioned into the quasi-steady MPD (deflagration) mode of plasma acceleration. In the present work, variations of the controllable parameters - specifically the pulsed circuit properties, the amount of mass injected into the thruster, and the relative timing between the initial gas injection and the initiation of the plasma current sheet - will be used to explore the thruster performance. A range of parameters are explored to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗

The delayed rectifier, IKI, is the major conductance in type I vestibular hair cells across vestibular end organs

Hair cells were dissociated from the semicircular canal, utricle, lagena and saccule of white king pigeons. Type I hair cells were identified morphologically based on the ratios of neck width to cuticular plate width (NPR < 0.72) as well as neck width to cell body width (NBR < 0.64). The perforated patch variant of the whole-cell recording technique was used to measure electrical properties from type I hair cells. In voltage-clamp, the membrane properties of all identified type I cells were dominated by a predominantly outward potassium current, previously characterized in semicircular canal as IKI. Zero-current potential, activation, deactivation, slope conductance, pharmacologic and steady-state properties of the complex currents were not statistically different between type I hair cells of different vestibular end organs. The voltage dependence causes a significant proportion of this conductance to be active about the cell's zero-current potential. The first report of the whole-cell activation kinetics of the conductance is presented, showing a voltage dependence that could be best fit by an equation for a single exponential. Results presented here are the first data from pigeon dissociated type I hair cells from utricle, saccule and lagena suggesting that the basolateral conductances of a morphologically identified population of type I hair cells are conserved between functionally different vestibular end organs; the major conductance being a delayed rectifier characterized previously in semicircular canal hair cells as IKI.

NASA Discipline Number 40-99↗

The Evolution of Open Magnetic Flux Driven by Photospheric Dynamics

The coronal magnetic field is of paramount importance in solar and heliospheric physics. Two profoundly different views of the coronal magnetic field have emerged. In quasi-steady models, the predominant source of open magnetic field is in coronal holes. In contrast, in the interchange model, the open magnetic flux is conserved, and the coronal magnetic field can only respond to the photospheric evolution via interchange reconnection. In this view the open magnetic flux diffuses through the closed, streamer belt fields, and substantial open flux is present in the streamer belt during solar minimum. However, Antiochos and co-workers, in the form of a conjecture, argued that truly isolated open flux cannot exist in a configuration with one heliospheric current sheet (HCS) - it will connect via narrow corridors to the polar coronal hole of the same polarity. This contradicts the requirements of the interchange model. We have performed an MHD simulation of the solar corona up to 20R solar to test both the interchange model and the Antiochos conjecture. We use a synoptic map for Carrington Rotation 1913 as the boundary condition for the model, with two small bipoles introduced into the region where a positive polarity extended coronal hole forms. We introduce flows at the photospheric boundary surface to see if open flux associated with the bipoles can be moved into the closed-field region. Interchange reconnection does occur in response to these motions. However, we find that the open magnetic flux cannot be simply injected into closed-field regions - the flux eventually closes down and disconnected flux is created. Flux either opens or closes, as required, to maintain topologically distinct open and closed field regions, with no indiscriminate mixing of the two. The early evolution conforms to the Antiochos conjecture in that a narrow corridor of open flux connects the portion of the coronal hole that is nearly detached by one of the bipoles. In the later evolution, a detached coronal hole forms, in apparent violation of the Antiochos conjecture. Further investigation reveals that this detached coronal hole is actually linked to the extended coronal hole by a separatrix footprint on the photosphere of zero width. Therefore, the essential idea of the conjecture is preserved, if we modify it to state that coronal holes in the same polarity region are always linked, either by finite width corridors or separatrix footprints. The implications of these results for interchange reconnection and the sources of the slow solar wind are briefly discussed.

Linker, Jon A.↗

An Effort to Map and Monitor Baldcypress Forest Areas in Coastal Louisiana, Using Landsat, MODIS, and ASTER Satellite Data

This presentation discusses a collaborative project to develop, test, and demonstrate baldcypress forest mapping and monitoring products for aiding forest conservation and restoration in coastal Louisiana. Low lying coastal forests in the region are being negatively impacted by multiple factors, including subsidence, salt water intrusion, sea level rise, persistent flooding, hydrologic modification, annual insect-induced forest defoliation, timber harvesting, and conversion to urban land uses. Coastal baldcypress forests provide invaluable ecological services in terms of wildlife habitat, forest products, storm buffers, and water quality benefits. Before this project, current maps of baldcypress forest concentrations and change did not exist or were out of date. In response, this project was initiated to produce: 1) current maps showing the extent and location of baldcypress dominated forests; and 2) wetland forest change maps showing temporary and persistent disturbance and loss since the early 1970s. Project products are being developed collaboratively with multiple state and federal agencies. Products are being validated using available reference data from aerial, satellite, and field survey data. Results include Landsat TM- based classifications of baldcypress in terms of cover type and percent canopy cover. Landsat MSS data was employed to compute a circa 1972 classification of swamp and bottomland hardwood forest types. Landsat data for 1972-2010 was used to compute wetland forest change products. MODIS-based change products were applied to view and assess insect-induced swamp forest defoliation. MODIS, Landsat, and ASTER satellite data products were used to help assess hurricane and flood impacts to coastal wetland forests in the region.

Spruce, Joseph P.↗

Drones for aerodynamic and structural testing /DAST/ - A status report

A program for providing research data on aerodynamic loads and active control systems on wings with supercritical airfoils in the transonic speed range is described. Analytical development, wind tunnel tests, and flight tests are included. A Firebee II target drone vehicle has been modified for use as a flight test facility. The program currently includes flight experiments on two aeroelastic research wings. The primary purpose of the first flight experiment is to demonstrate an active control system for flutter suppression on a transport-type wing. Design and fabrication of the wing are complete and after installing research instrumentation and the flutter suppression system, flight testing is expected to begin in early 1979. The experiment on the second research wing - a fuel-conservative transport type - is to demonstrate multiple active control systems including flutter suppression, maneuver load alleviation, gust load alleviation, and reduce static stability. Of special importance for this second experiment is the development and validation of integrated design methods which include the benefits of active controls in the structural design.

Murrow, H. N.↗

Population Viability Analysis (PVA) as a Platform for Predicting Outcomes of Management Options for the Florida Scrub-Jay in Brevard County

The Florida scrub-jay (FSJ) is a species in decline because of extinction debt caused by habitat fragmentation and degradation. Understanding and managing the species and its habitats is challenging due to the complex interactions among the social system, a dynamic mosaic of scrub habitat, and active management. To provide insights into the possible fates of the FSJ populations of mainland, cape, and island sections of Brevard County, we modeled the population dynamics using the Vortex population viability analysis (PVA) software. Vortex is an individual-based simulation that allowed us to include such factors as demographic rates dependent on habitat state, impact of helpers on breeding success, and helper to breeder transition probabilities responding to availability of nearby vacant optimal habitat and vacancies due to the death of breeders. Detailed modeling of the FSJ was possible only because a lot of data about the species and its habitats have been gathered over decades of intensive research. We followed a phased approach to constructing population models that incorporated the best available science and data to address a variety of conservation actions. The first phase focused on constructing a model that incorporates sociobiology and source-sink habitat dynamics. This model allowed us to address some of the most important questions about population size and habitat quality. Once the model framework was in place, we considered the real landscapes and actual local populations, rather than just generic representations of typical FSJ dynamics. After examining the viability of the metapopulations under current conditions, we explored the likely consequences of various management actions that might slow or reverse population declines.

Population Viability Analysis↗

Cape Hatteras Ecological Conservation: Delineating Shoreline and Mapping Change Along the Cape Hatteras National Seashore for Coastline Management and Transportation Corridor Adaptation Strategies

The National Park Service at Cape Hatteras National Seashore works to protect North Carolina’s Outer Banks where frequent storms bring heavy winds and flooding, leading to overwash events on the main highway, North Carolina Highway 12. Shorelines are susceptible to erosion directly affecting transportation and housing infrastructures. Storm events can disrupt transportation on NC-12 and ferry service from Ocracoke Island to Hatteras Island, leaving inhabitants stranded on Ocracoke Island for indefinite amounts of time. The National Park Service’s current decision-making practices involve mitigating infrastructure damage by enlisting the help of North Carolina’s Department of Transportation and finding ways to relocate this infrastructure, as well as dredging and sediment placement that support beach nourishment efforts. Our NASA DEVELOP team partnered with the National Park Service to explore the use of optical data to delineate shoreline and map coastline change from 2014 to 2024. We used Earth observations collected by Landsat 8 Operational Land Imager, Landsat 9 Operational Land Imager-2, and Sentinel-2 MultiSpectral instrument to support decision-making related to prioritizing strategic planning for transportation corridor adaptations including potential relocation of infrastructure, and beach nourishment efforts. We derived coastline maps by consolidating images from each year to create composite images for winter and summer coastline seasonal oscillation patterns. The results highlight shoreline loss over the ~10-year study period, and the difference in shoreline inundation in the winter and summer seasons. These observations create a better understanding of successful mitigation efforts and allow the National Park Service to continue to focus on plans for infrastructure updates.

Ella Haugen↗

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

Bellan, Josette↗

Studies on Somatic Embryogenesis in Sweetpotato

The purpose of this study was to improve the somatic embryo (SE) system for plant production of sweetpotato Ipomoea batatas L.(Lam)l. Explants isolated from SE-derived sweet potato plants were compared with control (non SE-derived) plants for their competency for SE production. Leaf explants were cultured on Murashige-Skoog (MS) medium with 2,4-dichlorophenoxy acetic acid (0.2 mg/L) and 6-benzylaminopurine (2.5 mg/L) for 2 weeks in darkness and transferred to MS medium with abscisic acid (2.5 Explants isolated from those plants developed through somatic embryo-genesis produced new somatic embryos rapidly and in higher frequency than those isolated from control plants. They also appeared to grow faster in tissue culture than the control plants. Current studies in the laboratory are examining whether plants derived from a cyclical embryogenesis system (five cycles) would have any further positive impact on the rapidity and frequency of somatic embryo development. More detailed studies using electron microscopy are expected to show the point of origin of the embryos and to allow determination of their quality throughout the cyclical process. This study may facilitate improved plant micropropagation, gene transfer and germplasm conservation in sweet potato.

Bennett, J. Rasheed↗

Studies for Somatic Embryogenesis in Sweet Potato

The purpose of this study was to improve the somatic embryo (SE) system for plant production of sweet potato (Ipomoea batatas L(Lam)). Explants isolated from SE-derived sweet potato plants were compared with control (non SE-derived) plants for their competency for SE production. Leaf explants were cultured on Murashige-Skoog (MS) medium with 2,4-dichlorophenoxy acetic acid (0.2 mg/L) and 6-benzylaminopurine (2.5 mg/L) for 2 weeks in darkness and transferred to MS medium with abscisic acid (2.5 mg/L). Explants isolated from those plants developed through somatic embryogenesis produced new somatic embryos rapidly and in higher frequency than those isolated from control plants They also appeared to grow faster in tissue culture than the control plants. Current studies in the laboratory are examining whether plants derived from a cyclical embryogenesis system (five cycles) would have any further positive impact on the rapidity and frequency of somatic embryo development. More detailed studies using electron microscopy are expected to show the point of origin of the embryos and to allow determination of their quality throughout the cyclical process. This study may facilitate improved plant micropropagation, gene transfer and germplasm conservation in sweet potato.

Bennett, J. Rasheed↗

Control-oriented modeling of combustion and flow processes in liquid propellant rocket engines

This paper presents a control-oriented model of the flow, reaction, and transport processes in liquid propellant rocket combustion chambers, based on the multicomponent conservation laws of gas dynamics. This model provides a framework for the inclusion of detailed chemical kinetic relations, viscous and other dissipative effects, a variety of actuators and sensors, as well as process and measurement disturbances. In addition to its potential usefulness to the designer in understanding the dynamical complexity of the system and the sources of model uncertainty, the model provides a rigorous basis for control system design. An appraisal of current and feasible actuators and sensors, and their mathematical representation are included.

Bentsman, Joseph↗

New developments in the method of space-time conservation element and solution element: Applications to the Euler and Navier-Stokes equations

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods--i.e., finite difference, finite volume, finite element, and spectral methods. It is conceptually simple and designed to avoid several key limitations to the above traditional methods. An explicit model scheme for solving a simple 1-D unsteady convection-diffusion equation is constructed and used to illuminate major differences between the current method and those mentioned above. Unexpectedly, its amplification factors for the pure convection and pure diffusion cases are identical to those of the Leapfrog and the DuFort-Frankel schemes, respectively. Also, this explicit scheme and its Navier-Stokes extension have the unusual property that their stabilities are limited only by the CFL condition. Moreover, despite the fact that it does not use any flux-limiter or slope-limiter, the Navier-Stokes solver is capable of generating highly accurate shock tube solutions with shock discontinuities being resolved within one mesh interval. An accurate Euler solver also is constructed through another extension. It has many unusual properties, e.g., numerical diffusion at all mesh points can be controlled by a set of local parameters.

Chang, Sin-Chung↗

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗

Aerothermal Analysis of the Project Fire II Afterbody Flow

Computational fluid dynamics (CFD) is used to simulate the wake flow and afterbody heating of the Project Fire II ballistic reentry to Earth at 11.4 km/sec. Laminar results are obtained over a portion of the trajectory between the initial heat pulse and peak afterbody heating. Although non-catalytic forebody convective heating results are in excellent agreement with previous computations, initial predictions of afterbody heating were about a factor of two below the experimental values. Further analysis suggests that significant catalysis may be occurring on the afterbody heat shield. Computations including finite-rate catalysis on the afterbody surface are in good agreement with the data over the early portion of the trajectory, but are conservative near the peak afterbody heating point, especially on the rear portion of the conical frustum. Further analysis of the flight data from Fire II shows that peak afterbody heating occurs before peak forebody heating, a result that contradicts computations and flight data from other entry vehicles. This result suggests that another mechanism, possibly pyrolysis, may be occurring during the later portion of the trajectory, resulting in less total heat transfer than the current predictions.

Wright, Michael J.↗