Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Rayleigh Scattering Density Measurements, Cluster Theory, and Nucleation Calculations at Mach 10

In an exploratory investigation, quantitative unclustered laser Rayleigh scattering measurements of density were performed in the air in the NASA Langley Research Center's 31 in. Mach 10 wind tunnel. A review of 20 previous years of data in supersonic and Mach 6 hypersonic flows is presented where clustered signals typically overwhelmed molecular signals. A review of nucleation theory and accompanying nucleation calculations are also provided to interpret the current observed lack of clustering. Data were acquired at a fixed stagnation temperature near 990Kat five stagnation pressures spanning 2.41 to 10.0 MPa (350 to 1454 psi) using a pulsed argon fluoride excimer laser and double-intensified charge-coupled device camera. Data averaged over 371 images and 210 pixels along a 36.7mmline measured freestream densities that agree with computed isentropic-expansion densities to less than 2% and less than 6% at the highest and lowest densities, respectively. Cluster-free Mach 10 results are compared with previous clustered Mach 6 and condensation-free Mach 14 results. Evidence is presented indicating vibrationally excited oxygen and nitrogen molecules are absorbed as the clusters form, release their excess energy, and inhibit or possibly reverse the clustering process. Implications for delaying clustering and condensation onset in hypersonic and hypervelocity facilities are discussed.

Balla, R. Jeffrey↗
NASA NTRSDate not supplied

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗
NASA NTRSDate not supplied

Overview of the Subscale RAVEN Flight Controls and Modeling Testbed

The Research Aircraft for eVTOL Enabling TechNologies (RAVEN) Subscale Wind-Tunnel and Flight Test (SWFT) model is a subscale aircraft built for flight dynamics and controls research demonstrated in wind-tunnel and flight-test experiments. The intent of this paper is to provide a summary of past, current, and future efforts being pursued by the RAVEN-SWFT project. Initially, vehicle development guidelines were crafted by a multidisciplinary team to ensure that the RAVEN-SWFT vehicle was well suited for research in multiple areas, including aero-propulsive modeling, flight controls, and autonomy, among others. The vehicle has been used to obtain extensive wind-tunnel data, enabling aero-propulsive model development across the transition flight envelope and validation of computational tools. The vehicle will be used to conduct flight testing in order to evaluate modeling strategies and flight control logic. The RAVEN-SWFT model also serves as a risk reduction activity for a conceptual, full-scale vehicle in the 1000-lb class. The next steps in the project are to successfully demonstrate free flight in hover, transition, forward flight, and the reverse thereof, utilizing custom control laws integrated onto the RAVEN-SWFT avionics hardware. The project intends to publicize all of the geometry, data, and methods in future reports.

eVTOL↗
DOE OSTI2024

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗
DOE OSTI2024

Strain Fluctuations Unlock Ferroelectricity in Wurtzites

Ferroelectrics are of practical interest for non-volatile data storage due to their reorientable, crystallographically defined polarization. Yet efforts to integrate conventional ferroelectrics into ultrathin memories have been frustrated by film-thickness limitations, which impede polarization reversal under low applied voltage. Wurtzite materials, including magnesium-substituted zinc oxide (Zn,Mg)O, have been shown to exhibit scalable ferroelectricity as thin films. In this work, the origins of ferroelectricity in (Zn,Mg)O are explained, showing that large strain fluctuations emerge locally in (Zn,Mg)O and can reduce local barriers to ferroelectric switching by more than 40%. Concurrent experimental and computational evidence of these effects are provided by demonstrating polarization switching in ZnO/(Zn,Mg)O/ZnO heterostructures featuring built-in interfacial strain gradients. These results open up an avenue to develop scalable ferroelectrics by controlling strain fluctuations atomistically.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CO 2 Hydrogenation to Methanol over Inverse ZrO 2 /Cu(111) Catalysts: The Fate of Methoxy under Dry and Wet Conditions

Understanding the surface chemistry of CH 3 O species is essential for the production of methanol by CO 2 hydrogenation over Cu-based heterogeneous catalysts, as it facilitates the rational design of more efficient conversion processes. Recent research has identified inverse ZrO 2 /Cu catalysts as highly active and selective systems for the transformation of CO 2 to methanol with a performance that can be better than that of commercial Cu/ZnO catalysts. Here, we employed synchrotron-based ambient pressure X-ray photoelectron spectroscopy (AP-XPS) and calculations based on density functional theory (DFT) to understand the fate of CH 3 O groups under dry and wet environments. AP-XPS spectra revealed that under CO 2 hydrogenation conditions, formate and methoxy are two key intermediates to produce methanol. Furthermore, there are three different types of reactive sites on the surface: One is active for methoxy adsorption, which is stable and responsible for the methanol synthesis; Another one transforms CO 2 into CO; and a third one is active for CO 2 and methoxy dissociation, leading to C and methane formation. The theoretical calculations indicate that CH 3 OH readily dissociates to CH 3 O species following a highly exothermic (ΔE = -20.99 kcal/mol) and barrierless process. The water produced by the reverse water-gas shift reaction (CO 2 + H 2 → H 2 O + CO) can prevent the decomposition of CH 3 O species. We discovered that by introducing a tiny amount of water vapor (2 × 10 -6 Torr) into the reaction chamber, the energy barrier for the reaction CH 3 O(ads) + H(ads) → CH 3 OH(gas) is dramatically reduced. AP-XPS and computational modelling showed that water is quite capable of extracting adsorbed methoxy to form gaseous methanol. With this in mind, one could boost the methanol selectivity by adding appropriate amounts of water or steam, which is an inexpensive and feasible solution for industrial operations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of Organic Cations on Hydrogen Oxidation Reaction of Carbon Supported Platinum

Effect of organic cations on hydrogen oxidation reaction (HOR) of carbon supported platinum (Pt/C) is investigated using three 0.1 M alkaline electrolytes, tetramethylammonium hydroxide (TMAOH), tetrabutylammonium hydroxide (TBAOH) and tetrabutylphosphonium hydroxide (TBPOH). Rotating disk electrode experiments indicate that the HOR of Pt/C is adversely impacted by time-dependent and potential-driven chemisorption of organic cations. In-situ infrared reflection adsorption spectroscopy experiments indicated that the specific chemisorption of organic cations drives the hydroxide co-adsorption on Pt surface. The co-adsorption of TMA + and hydroxide at 0.1 V vs. reversible hydrogen electrode is the strongest; consequently, complete removal of the co-adsorbed layer from Pt surface is difficult even after exposure the Pt surface to 1.2 V. Conversely, the chemisorption of TBP+ is the weakest, yet notable decrease of HOR current density is still observed. The adsorption energies, ΔE, for TMA + , TBA + , and TBP + on Pt (111) surface from density functional theory are computed to be -2.79, -2.42 and -2.00 eV, respectively. The relatively low adsorption energy of TBP + is explained by the steric hindrance and electronic effect. This study emphasizes the importance of cationic group on HOR activity of alkaline anion exchange membrane fuel cells.

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

Global Characteristics of Observable Foreshocks for Large Earthquakes

Abstract Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.

Geochemistry & Geophysics↗
NASA NTRS1985

Detailed measurements of the flowfield in the vicinity of an airfoil with glaze ice

An experimental study has been conducted in the OSU subsonic tunnel to measure the characteristics of the separation bubble on an airfoil with glaze ice. A measured glaze ice accretion on a NACA 0012 airfoil was simulated in wood for this dry tunnel test. The 21 inch chord model was pressure belted and the ice shape internally tapped to obtain surface pressures, lift and moment coefficients. A wake survey probe was used to obtain airfoil drag. The separation bubble was explored by measuring the time averaged velocities using a split film probe. The probe was positioned using a computer controlled two-dimensional traversing system. In this paper, airfoil lift, drag, and moment coefficient data are compared for the airfoil with and without glaze ice. Velocity profiles in the separation bubble are presented for several chordwise stations at three angles of attack. The ice shape caused a severe lift and drag penalty. The velocity profiles show clearly the large bubble geometry, regions of reversed flow, and bubble reattachment.

Bragg, M. B.↗
NASA NTRS1942

A method for determining the camber and twist of a surface to support a given distribution of lift

A graphical method is described for finding the shape (camber and twist) of an airfoil having an arbitrary distribution of lift. The method consists in replacing the lifting surface and its wake with an equivalent arrangement of vortices and in finding the associated vertical velocities. By a division of the vortex pattern into circular strips concentric about the downwash point instead of into the usual rectangular strips, the lifting surface is reduced for each downwash point to an equivalent loaded line for which the induced velocity is readily computed. The ratio of the vertical velocity to the stream velocity is the slope of the surface in the free-stream direction. As an illustration, the shape of the wing consistent with the pressure distribution derived from the two-dimensional theories is found for two wings: a straight elliptical wing and one with 30 degree sweepback. Application of the method to solve the reverse problem - finding the lift distribution over a given surface - is briefly discussed.

Cohen, Doris↗
NASA NTRS2012

Application of the NASCAP Spacecraft Simulation Tool to Investigate Electrodynamic Tether Current Collection in LEO

Recent interest in using electrodynamic tethers (EDTs) for orbital maneuvering in Low Earth Orbit (LEO) has prompted the development of the Marshall ElectroDynamic Tether Orbit Propagator (MEDTOP) model. The model is comprised of several modules which address various aspects of EDT propulsion, including calculation of state vectors using a standard orbit propagator (e.g., J2), an atmospheric drag model, realistic ionospheric and magnetic field models, space weather effects, and tether librations. The natural electromotive force (EMF) attained during a radially-aligned conductive tether results in electrons flowing down the tether and accumulating on the lower-altitude spacecraft. The energy that drives this EMF is sourced from the orbital energy of the system; thus, EDTs are often proposed as de-orbiting systems. However, when the current is reversed using satellite charged particle sources, then propulsion is possible. One of the most difficult challenges of the modeling effort is to ascertain the equivalent circuit between the spacecraft and the ionospheric plasma. The present study investigates the use of the NASA Charging Analyzer Program (NASCAP) to calculate currents to and from the tethered satellites and the ionospheric plasma. NASCAP is a sophisticated set of computational tools to model the surface charging of three-dimensional (3D) spacecraft surfaces in a time-varying space environment. The model's surface is tessellated into a collection of facets, and NASCAP calculates currents and potentials for each one. Additionally, NASCAP provides for the construction of one or more nested grids to calculate space potential and time-varying electric fields. This provides for the capability to track individual particles orbits, to model charged particle wakes, and to incorporate external charged particle sources. With this study, we have developed a model of calculating currents incident onto an electrodynamic tethered satellite system, and first results are shown here.

Adams, Mitzi↗
NASA NTRS2019

Comparative Study of Active Flow Control Strategies for Lift Enhancement of a Simplified High-Lift Configuration

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of the conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-AFC configuration to recover the aerodynamic performance (lift) of the CRM-HL configuration. In the current paper, a Lattice Boltzmann method-based computational fluid dynamics (CFD) code, known as PowerFLOWQ® is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with several different types of AFC devices. The transonic version of the PowerFLOWQ® code that has been validated for high speed flows is used to accurately simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient AFC systems and actuator arrangements emerged based on the parametric studies performed prior to a Fall 2018 wind tunnel test. Preliminary comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations.

Vatsa, Veer N.↗
DOE OSTI2015

Peculiarities of plastic deformation nucleation in copper under nanoindentation

The computer simulation results on the atomic structure of the copper crystallite and its behavior in nanoindentation demonstrate the key role of local structural transformations in nucleation of plasticity. The generation of local structural transformations can be considered as an elementary event during the formation of higher scale defects, including partial dislocations and stacking faults. The cause for local structural transformations, both direct fcc-hcp and reverse hcp-fcc, is an abrupt local increase in atomic volume. A characteristic feature is that the values of local volume jumps in direct and reverse structural transformations are comparable with that in melting and lie in the range 5–7%.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Miniaturized Magnetoelastic Sensor System

This article describes the design, assembly, and implementation of a hand-held, magnetic-field-based sensor system that can be adapted for a variety of sensing applications. The miniaturized system is based on Chemical Identification by Magneto-Elastic Sensing (ChIMES) technology, which uses three concentric solenoid coils to wirelessly interrogate a sensor body comprised of a response material coupled to a magnetoelastic wire. The response material expands when it encounters a target, imposing mechanical stress on the wire and altering its magnetic permeability. The sensor bodies are passive, requiring no external power source, and they are small, measuring about 15 mm in length and 3.0 mm in diameter. Up to four sensor bodies can be configured as an evenly-spaced linear array. The sensor system operates by applying a low-frequency, current-stabilized, filtered triangle wave to a uniform-density excitation coil to switch the magnetic domains within the wire. Further, the responses from the sensors are picked up by a detection coil as stress-induced changes in the Faraday voltage, and the strong magnetic field induced by the excitation coil in the detection coil is nullified by a cancellation coil reverse-wound in series with the detection coil. The responses of the sensors in an array are separated in time by a linear gradient dc biasing coil. The sensors can be interrogated through metallic and nonmetallic barriers. The signals from the detection coil and the excitation coil are digitized by a pair of bipolar analog-to-digital converters (ADCs). A Raspberry Pi single-board computer (SBC) and associated software perform data acquisition and control all aspects of the sensor system hardware. The program allows the user to select the number of sensors in the array, the type of signal that is being collected, and the number of samples to take. The program also allows for signal processing of the sensor data, such as baseline correction. The program can differentiate sensor peaks from each other and calculate the magnitude of each sensor response with less than 1% error. The data are then displayed along with a graph of the signal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Modelling electron transfer in photosystem I: limits and perspectives

Uncovering the parameters underlying the electron transfer (ET) in photosynthetic reaction centres is of importance for understanding the molecular mechanisms underpinning their functionality. The reductive nature of most cofactors involved in photosynthetic ET makes the direct estimation of their properties difficult. Photosystem I (PSI) operates in a highly reducing regime, making the assessment of cofactor properties even more difficult. Kinetic modelling coupled to a non-adiabatic description of ET is a useful approach in overcoming this hindrance. Here we review the theory and modelling approaches that have been used in assessing parameters associated with ET reactions in PSI, with particular attention to ET reactions involving the phylloquinones and the iron–sulphur clusters. In most modelling studies, the goal is to estimate the driving force of ET, which is usually associated with the cofactor midpoint potentials. The driving force is sensitive to many factors, which define the ET rate, i.e. the reorganisation energy, the coupling with nuclear modes and the electronic matrix elements, which are explored and discussed here. The importance of an inclusive modelling of both forward and reverse ET processes is discussed and highlighted. It is shown that although estimates are indeed sensitive to the exact parameter sets employed in the modelling, a general consensus is still attained, pointing to a scenario where $ΔG^0_{A_{tA→Fx}}/ΔG^0_{A_{tB→Fx}}$ is weakly endergonic/exergonic, respectively. It is emphasised that to further refine those estimates, it will require a joint effort between computational modelling and more wide-ranging experimental studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Programmable Dynamic Self-Assembly of DNA Nanostructures (Final Technical Report)

This technical report summarizes our research toward the synthesis of active biomolecular materials that take advantage of the biological properties of DNA and RNA, but can operate in non-biological contexts. Our approach takes inspiration from the dynamic assembly and disassembly of cytoskeletal filaments inside cells, and focuses on the generation of artificial filaments with comparable adaptation and responsiveness. We report the construction of DNA nanotube systems that assemble and disassemble reversibly in response to a variety of molecular stimuli including other nucleic acids (DNA and RNA), enzymes, and pH, as well as our progress on developing a reaction network framework to regulate self-assembly. Our DNA polymers are active in that they rely on the programmable energetics of hybridization and on the kinetics of enzymatic reactions to perform assembly and disassembly. Results include experiments and computational models. The research described here was performed between June 15, 2016 and June 15, 2019 at the University of California at Riverside with the support of the collaborative award DE-SC0010595 to PI Elisa Franco at UC Riverside and PI Rebecca Schulman at Johns Hopkins. The report focuses on the research performed by the team of PI Franco. The report also details work performed during a no-cost extension (June 15, 2019 - December 15, 2019) awarded to PI Guillermo Aguilar and PI Franco.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Revolutionizing Turbine Cooling with Micro-Architectures Enabled by Direct Metal Laser Sintering

The objective of this research effort was to explore innovative cooling architectures enabled by additive manufacturing techniques for improved turbine cooling performance. The ability to create complex internal geometries was leveraged to better distribute coolant as well as to integrate inherently unstable flow devices to enhance internal and external heat transfer. This was accomplished with a multi-faceted approach including analytical, experimental, and computational components. This final report documents progress during the total 4.25 year effort that was extended (at no additional cost) from the original 3-year cooperative agreement. Prior to starting the project, significant effort was invested in investigating innovative cooling designs from the open literature. When the project started, we quickly down-selected to 4 configurations: impingement jet fluidic oscillators (for internal leading edge cooling), reverse film cooling (for pressure surface), sweeping film cooling jets (for suction surface), and trailing edge slot cooling with microchannels. During the 2nd year, the reverse film cooling design was shown to be sub-optimal and work on that topic was halted. The three remaining technologies were integrated into a large scale nozzle guide vane and installed into a low-speed linear cascade facility for further interrogation. Initial results showed the sweeping film cooling jets (on the suction surface) to be the most promising technology compared to a baseline diffusion shaped film hole (777 design). The benefit was particularly evident at high blowing ratios (>1.5) when the 777 coolant film separated from the downstream surface. Benefits were also evident at elevated freestream turbulence levels. The leading edge cooling with unsteady jets was less effective in terms of peak or average cooling – however it was superior to round hole impingement cooling in terms of spatial uniformity. Finally, the trailing edge cooling design with micro-channels was scrapped in favor of a pinned arrangement with centerbody. During the 3rd year, the same three technologies were integrated into a transonic linear cascade for an assessment of compressibility effects. Again, the sweeping film cooling jets proved superior at high blowing ratios while the unsteady leading edge impingement and trailing edge pin-fins with centerbody designs yielded mixed results. Finally, during the 4th and final year, these 3 technologies were integrated into a direct metal laser sintered (DMLS) nozzle guide vane for testing in a high temperature, transonic annular vane cascade. The facility matches the flow temperature, Mach number, and coolant pressure ratios of an actual gas turbine. Before doing this, an additional series of tests were completed to validate a redesigned trailing edge cooling architecture. This final design included a centerbody with triangular pins between the centerbody and the vane external skin. The majority of these triangular pins are fabricated with a 30% gap to the centerbody – reducing their pressure drop considerably while still providing excellent heat transfer augmentation. This additional test campaign required an additional 3-month extension request to complete testing in the high temperature NGV test facility (Turbine Reacting Flow Rig – TuRFR). Testing of the DMLS vane in TuRFR pitted the 3 innovative cooling technologies against more traditional technologies (shaped 777 film hole, round impingement jet, and full pin arrangement in the trailing edge. The sweeping film cooling jets (on the suction surface) showed approximately 15% improvement in overall effectiveness compared to the 777 film hole while the trailing edge showed up to 20% improvement. The leading edge impingement was disappointing with a substantially lower effectiveness than traditional direct impingement. In summary, the study was successful in demonstrating that DMLS-enabled cooling technologies can yield significant gains in cooling performance.

03 NATURAL GAS↗
NASA NTRS1984

A unified N-body and statistical treatment of stellar dynamics. I - The hybrid code. II - Applications to globular cluster cores

A unified N-body and statistical treatment of stellar dynamics is developed and applied to the late stages of core collapse and early stages of post collapse evolution in globular clusters. A 'hybrid' computer code is joined to a direct N-body code which is used to calculate exactly the behavior of particles in the inner spatial region, and the combination is used to follow particles statistically in the outer spatial region. A transition zone allows the exchange of particles and energy between the two regions. The main application results include: formation of a hard central binary system, reversal of core collapse and expansion due to the heat input from this binary, ejection of the binary from the core, and recollapse of the core; density profiles that form a one-parameter sequence during the core oscillations; and indications that these oscillations will eventually cease.

Mcmillan, S. L. W.↗