Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 523 records · Page 29

COMPASS Final Report: Lunar Communications Terminal (LCT)

The Lunar Communications Terminal (LCT) COllaborative Modeling and Parametric Assessment of Space Systems (COMPASS) session designed a terminal to provide communications between lunar South Pole assets, communications relay to/from these assets through an orbiting Lunar Relay Satellite (LRS) and navigation support. The design included a complete master equipment list, power requirement list, configuration design, and brief risk assessment and cost analysis. The Terminal consists of a pallet containing the communications and avionics equipment, surrounded by the thermal control system (radiator), an attached, deployable 10-m tower, upon which were mounted locally broadcasting and receiving modems and a deployable 1 m diameter Ka/S band dish which provides relay communications with the lunar relay satellites and, as a backup, Earth when it is in view. All power was assumed to come from the lunar outpost Habitat. Three LCT design options were explored: a stand-alone LCT servicing the manned outpost, an integrated LCT (into the Habitat or Lunar Lander), and a mini-LCT which provides a reduced level of communication for primarily robotic areas dealing as in situ resource utilization (ISRU) and remote science. Where possible all the designs assumed single fault tolerance. Significant mass savings were found when integrating the LCT into the Habitat or Lander but increases in costs occurred depending upon the level of man rating required for such designs.

Oleson, Steven R.↗

Process-Based Cost Estimation for Ramjet/Scramjet Engines

Process-based cost estimation plays a key role in effecting cultural change that integrates distributed science, technology and engineering teams to rapidly create innovative and affordable products. Working together, NASA Glenn Research Center and Boeing Canoga Park have developed a methodology of process-based cost estimation bridging the methodologies of high-level parametric models and detailed bottoms-up estimation. The NASA GRC/Boeing CP process-based cost model provides a probabilistic structure of layered cost drivers. High-level inputs characterize mission requirements, system performance, and relevant economic factors. Design alternatives are extracted from a standard, product-specific work breakdown structure to pre-load lower-level cost driver inputs and generate the cost-risk analysis. As product design progresses and matures the lower level more detailed cost drivers can be re-accessed and the projected variation of input values narrowed, thereby generating a progressively more accurate estimate of cost-risk. Incorporated into the process-based cost model are techniques for decision analysis, specifically, the analytic hierarchy process (AHP) and functional utility analysis. Design alternatives may then be evaluated not just on cost-risk, but also user defined performance and schedule criteria. This implementation of full-trade study support contributes significantly to the realization of the integrated development environment. The process-based cost estimation model generates development and manufacturing cost estimates. The development team plans to expand the manufacturing process base from approximately 80 manufacturing processes to over 250 processes. Operation and support cost modeling is also envisioned. Process-based estimation considers the materials, resources, and processes in establishing cost-risk and rather depending on weight as an input, actually estimates weight along with cost and schedule.

Singh, Brijendra↗

Perceived Noise Analysis for Offset Jets Applied to Commercial Supersonic Aircraft

A systems analysis was performed with experimental jet noise data, engine/aircraft performance codes and aircraft noise prediction codes to assess takeoff noise levels and mission range for conceptual supersonic commercial aircraft. A parametric study was done to identify viable engine cycles that meet NASA's N+2 goals for noise and performance. Model scale data from offset jets were used as input to the aircraft noise prediction code to determine the expected sound levels for the lateral certification point where jet noise dominates over all other noise sources. The noise predictions were used to determine the optimal orientation of the offset nozzles to minimize the noise at the lateral microphone location. An alternative takeoff procedure called "programmed lapse rate" was evaluated for noise reduction benefits. Results show there are two types of engines that provide acceptable mission range performance; one is a conventional mixed-flow turbofan and the other is a three-stream variable-cycle engine. Separate flow offset nozzles reduce the noise directed toward the thicker side of the outer flow stream, but have less benefit as the core nozzle pressure ratio is reduced. At the systems level for a three-engine N+2 aircraft with full throttle takeoff, there is a 1.4 EPNdB margin to Chapter 3 noise regulations predicted for the lateral certification point (assuming jet noise dominates). With a 10% reduction in thrust just after clearing the runway, the margin increases to 5.5 EPNdB. Margins to Chapter 4 and Chapter 14 levels will depend on the cumulative split between the three certification points, but it appears that low specific thrust engines with a 10% reduction in thrust (programmed lapse rate) can come close to meeting Chapter 14 noise levels. Further noise reduction is possible with engine oversizing and derated takeoff, but more detailed mission studies are needed to investigate the range impacts as well as the practical limits for safety and takeoff regulations.

Turbofans↗

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗

Large Solar-Rejection Filter

NASA utilized Image Intensified Video Cameras for ATV data acquisition from a jet flying at 12.8 km. Afterwards the video was digitized and then analyzed with a modified commercial software package, Image Systems Trackeye. Astrometric results were limited by saturation, plate scale, and imposed linear plate solution based on field reference stars. Time-dependent fragment angular trajectories, velocities, accelerations, and luminosities were derived in each video segment. It was evident that individual fragments behave differently. Photometric accuracy was insufficient to confidently assess correlations between luminosity and fragment spatial behavior (velocity, deceleration). Use of high resolution digital video cameras in future should remedy this shortcoming. A scenically accurate description of matter interpreted as a substance made up of corpuscular constituents was established during the course of the 19th century. In this description, atoms--the building blocks of the matter--form molecules. The properties of the molecules were described by chemistry or thermodynamics depending on what characteristics of the matter were investigated. In both theories, the molecules can dissociate to atoms when the kinetic energies of the atoms exceed the strength of the chemical bonds. The number of atoms is always preserved in a closed system. This is not true, however, when the matter takes up much higher energies at relativistic scales. New particles can be produced at the expense of the kinetic energy. The number of particles is no longer preserved. There are other conserved quantities, however, these quantities, the charge, baryon number, lepton number, are associated with particles that are considered elementary today. The properties and behavior of these elementary particles is the subject of Particle Physics or High Energy Physics. Practice Page: A crewmember touching a positively charged surface was thought to be galvanically isolated from the vehicle ground analogous to a bird on a high voltage power wire. Recent analysis confirms that positive floating potentials, ionospheric currents to the EVA suit, can be hazardous. The analysis is wrong in that the ionospheric plasma itself can close the circuit. Parametric analysis of very low voltage exposures (2 to 15 volts) could cause pain and/or involuntary muscle tetani or spinal cord shock. NASA worked with the Naval Health Research Center Detachment Directed Energy Bioeffects Laboratory to examine the affects electrical hazards could have on extravehicular activity using two models. The results of the two computational models were combined to predict areas of the body in which neurons of different diameters would be excited. They predicted that physiologically active current could be conducted across the crew member causing catastrophic hazards. Future work to analyze additional current paths was proposed. The FUSE spectrum of BB Dor, observed in a high state, is modeled with an accretion disk with a very low inclination (possibly lower than 10 degrees). Assuming an average WD mass of 0.8 solar mass leads to a distance of the order of approximately 650pc, consistent with the extremely low galactic reddening in its direction, and a mass accretion rate of 10 (exp -9) solar mass a year. The spectrum presents some broad and deep silicon and sulfur absorption lines, indicating that these elements are over-abundant: silicon is 3 times solar, and sulfur is 20 times solar. The FUSE spectrum of BB Dor, observed in a high state, is modeled with an accretion disk with a very low inclination (possibly lower than 10 degrees). Assuming an average WD mass of 0.8 solar mass leads to a distance of the order of approximately 650pc, consistent with the extremely low galactic reddening in its direction, and a mass accretion rate of 10 (exp -9) solar mass a year. The spectrum presents some broad and deep silicon and sulfur absorption lines, indicating that these elements are over-abundant: silicon is times solar, and sulfur is 20 times solar. From its previously observed optical emission lines, P831-57 (WD 0334 6400 or Ret 1 in A Catalog and Atlas of Cataclysmic Variables: Living Edition) has been suspected to contain an accretion disk associated with a companion star in orbit around a subdwarf star with a temperature T is greater than 21,000K. P831-57 has therefore been classified as a nova-like. However, our present observations show it to be a DA + dMe binary. The analysis of its Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum (continuum and lines) reveal an average mass white dwarf (Log(g) approximately equals 7.8 plus or minus 0.1)with a temperature T approximately equals 37,000 plus or minus 500K, an extremely low projected rotational velocity, and a distance of about 115 plus or minus 5pc. The photosphere contains C, N, Si, and S (at about 1% of solar abundances). The dMe star is seen as a flux excess in near-infrared photometry and appears to show occasional flaring of about one magnitude as seen in the Harvard plates. There is no evidence of periodic variability in the spectroscopic or photometric data. We find no evidence of a an accretion disk, instead we find evidence of wind accretion as the stellar carbon abundance (N(C)/N(H) = 2.5 x 10 (exp -6) is about ten times larger than predicted by radiative levitation for such a gravity and temperature. The power needs and solutions for the space exploration and lunar mobility program are discussed. Long term missions in space and on the lunar surface require high energy batteries. Rechargeable batteries for mobility systems and portable utility pallet are needed for successful exploration missions. Nanomaterial usage increases the energy density of the cells apart from increasing the power density. The symptoms and threats from acute mountain sickness (AMS) are discussed. The underlying assumptions concerning spacecraft atmosphere mean there is a potential risk to astronauts. The baseline worst case scenario is about 25% based on direct ascent to 8.0 psia with 32% O2. The practical mitigation approach is a staged depressurization scheme. Due to the uncertainty about potential AMS risk the flight surgeon should prepare, focused research should proceed, and current analytical efforts should continue. This slide presentation reviews the threat astronauts face from acute mountain sickness (AMS). The incidence of AMS is highly variable and the risk to astronauts is unclear. However, the potential risk is higher than current EAAs suggest and precautions should be implemented. Due to the uncertain nature of the risk flight surgeons should prepare, and focused research should proceed. Quasiparticle tunneling can be understood at low temperatures using a non-equilibrium kinetic theory. Low-temperature tunneling can be suppressed using existing techniques: island-lead gap engineering, lead cooling with SIN junctions, and quasiparticle traps. The slide presentation reviews the science goals and mission planning for the proposed Jupiter System Observer (JSO). The spacecraft would accommodate 9 instruments making the instrument payload more capable than any previously flown to Jupiter. The would-be mission includes study programs for Jupiter's satellites, interior, magnetosphere and atmosphere. The power needs and solutions for the space exploration and lunar mobility program are discussed. Long term missions in space and on the lunar surface require high energy batteries. Rechargeable batteries for mobility systems and portable utility pallet are needed for successful exploration missions. Nanomaterial usage increases the energy density of the cells apart from increasing the power density. This slide presentation reviews the power requirements for the space exploration and the lunar surface mobility programs. It includes information about the specifications for high energy batteries and the power requirements for lunar rovers,

Roberts, William↗

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗

Study of effect of PWR cold leg temperature gradient on reactor core condition

Effects of temperature and flow gradients in Westinghouse designed three-loop Pressurized Water Reactor (PWR) cold legs, the piping between the main coolant pump and the reactor vessel, were evaluated using Computational Fluid Dynamics (CFD) code STAR CCM+ and coupled neutronic and thermal-hydraulic (T/H) code system VERA. In the parametric study, several symmetric and asymmetric temperature gradients that were significantly larger than those observed from plant measurements were applied to the cold leg inlets for comparison with the base case without any temperature gradient. A CFD model using the STAR-CCM+ code was developed for a portion of the RCS region between the Reactor Coolant Pump (RCP) and the core inlet based on previously validated modeling approach. The CFD simulation results were processed for the temperature and flow rate distributions at the core inlet as input to the VERA calculations. The VERA code system consists of COBRA-TF (CTF) for thermal-hydraulics, MPACT for reactor physics and neutron transport, and ORIGEN for isotopic depletion. The VERA model was for depletion calculations of a high-burnup loading pattern with the reactor core in pin-by-pin and subchannel resolution. The results of the study indicate that the postulated temperature gradients within the PWR cold legs do not result in any significant changes in the core inlet temperature distributions and the core power distributions during the reactor operation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of Radiative Heating Uncertainty for Hyperbolic Earth Entry

This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.

Johnston, Christopher O.↗

Perceived Noise Analysis for Offset Jets Applied to Commercial Supersonic Aircraft

A systems analysis was performed with experimental jet noise data, engine/aircraft performance codes and aircraft noise prediction codes to assess takeoff noise levels and mission range for conceptual supersonic commercial aircraft. A parametric study was done to identify viable engine cycles that meet NASAs N+2 goals for noise and performance. Model scale data from offset jets was used as input to the aircraft noise prediction code to determine the expected sound levels for the lateral certification point where jet noise dominates over all other noise sources. The noise predictions were used to determine the optimal orientation of the offset nozzles to minimize the noise at the lateral microphone location. An alternative takeoff procedure called programmed lapse rate was evaluated for noise reduction benefits. Results show there are two types of engines that provide acceptable range performance; one is a standard mixed-flow turbofan with a single-stage fan, and the other is a three-stream variable-cycle engine with a multi-stage fan. The engine with a single-stage fan has a lower specific thrust and is 8 to 10 EPNdB quieter for takeoff. Offset nozzles reduce the noise directed toward the thicker side of the outer flow stream, but have less benefit as the core nozzle pressure ratio is reduced and the bypass-to-core area ratio increases. At the systems level for a three-engine N+2 aircraft with full throttle takeoff, there is a 1.4 EPNdB margin to Chapter 3 noise regulations predicted for the lateral certification point (assuming jet noise dominates). With a 10 reduction in thrust just after takeoff rotation, the margin increases to 5.5 EPNdB. Margins to Chapter 4 and Chapter 14 levels will depend on the cumulative split between the three certification points, but it appears that low specific thrust engines with a 10 reduction in thrust (programmed lapse rate) can come close to meeting Chapter 14 noise levels. Further noise reduction is possible with additional reduction in takeoff thrust using programmed lapse rate, but studies are needed to investigate the practical limits for safety and takeoff regulations.

Turbofans↗

Perceived Noise Analysis for Offset Jets Applied to Commercial Supersonic Aircraft

A systems analysis was performed with experimental jet noise data, engineaircraft performance codes and aircraft noise prediction codes to assess takeoff noise levels and mission range for conceptual supersonic commercial aircraft. A parametric study was done to identify viable engine cycles that meet NASAs N+2 goals for noise and performance. Model scale data from offset jets was used as input to the aircraft noise prediction code to determine the expected sound levels for the lateral certification point where jet noise dominates over all other noise sources. The noise predictions were used to determine the optimal orientation of the offset nozzles to minimize the noise at the lateral microphone location. An alternative takeoff procedure called programmed lapse rate was evaluated for noise reduction benefits. Results show there are two types of engines that provide acceptable range performance; one is a standard mixed-flow turbofan with a single-stage fan, and the other is a three-stream variable-cycle engine with a multi-stage fan. The engine with a single-stage fan has a lower specific thrust and is 8 to 10 EPNdB quieter for takeoff. Offset nozzles reduce the noise directed toward the thicker side of the outer flow stream, but have less benefit as the core nozzle pressure ratio is reduced and the bypass-to-core area ratio increases. At the systems level for a three-engine N+2 aircraft with full throttle takeoff, there is a 1.4 EPNdB margin to Chapter 3 noise regulations predicted for the lateral certification point (assuming jet noise dominates). With a 10 reduction in thrust just after takeoff rotation, the margin increases to 5.5 EPNdB. Margins to Chapter 4 and Chapter 14 levels will depend on the cumulative split between the three certification points, but it appears that low specific thrust engines with a 10 reduction in thrust (programmed lapse rate) can come close to meeting Chapter 14 noise levels. Further noise reduction is possible with additional reduction in takeoff thrust using programmed lapse rate, but studies are needed to investigate the practical limits for safety and takeoff regulations.

Jet Noise↗

Rapid Airplane Parametric Input Design (RAPID)

RAPID is a methodology and software system to define a class of airplane configurations and directly evaluate surface grids, volume grids, and grid sensitivity on and about the configurations. A distinguishing characteristic which separates RAPID from other airplane surface modellers is that the output grids and grid sensitivity are directly applicable in CFD analysis. A small set of design parameters and grid control parameters govern the process which is incorporated into interactive software for 'real time' visual analysis and into batch software for the application of optimization technology. The computed surface grids and volume grids are suitable for a wide range of Computational Fluid Dynamics (CFD) simulation. The general airplane configuration has wing, fuselage, horizontal tail, and vertical tail components. The double-delta wing and tail components are manifested by solving a fourth order partial differential equation (PDE) subject to Dirichlet and Neumann boundary conditions. The design parameters are incorporated into the boundary conditions and therefore govern the shapes of the surfaces. The PDE solution yields a smooth transition between boundaries. Surface grids suitable for CFD calculation are created by establishing an H-type topology about the configuration and incorporating grid spacing functions in the PDE equation for the lifting components and the fuselage definition equations. User specified grid parameters govern the location and degree of grid concentration. A two-block volume grid about a configuration is calculated using the Control Point Form (CPF) technique. The interactive software, which runs on Silicon Graphics IRIS workstations, allows design parameters to be continuously varied and the resulting surface grid to be observed in real time. The batch software computes both the surface and volume grids and also computes the sensitivity of the output grid with respect to the input design parameters by applying the precompiler tool ADIFOR to the grid generation program. The output of ADIFOR is a new source code containing the old code plus expressions for derivatives of specified dependent variables (grid coordinates) with respect to specified independent variables (design parameters). The RAPID methodology and software provide a means of rapidly defining numerical prototypes, grids, and grid sensitivity of a class of airplane configurations. This technology and software is highly useful for CFD research for preliminary design and optimization processes.

Smith, Robert E.↗

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Analysis and Modeling of a Two-Phase Jet Pump of a Flow Boiling Test Facility for Aerospace Applications

Jet pumps are devices capable of pumping fluids to a higher pressure employing a nozzle/diffuser/mixing chamber combination. A primary fluid is usually allowed to pass through a converging-diverging nozzle where it can accelerate to supersonic speeds at the nozzle exit. The relatively high kinetic energy that the primary fluid possesses at the nozzle exit is accompanied by a low pressure region in order to satisfy Bernoulli's equation. The low pressure region downstream of the nozzle exit permits a secondary fluid to be entrained into and mixed with the primary fluid in a mixing chamber located downstream of the nozzle. Several combinations may exist in terms of the nature of the primary and secondary fluids in so far as whether they are single or two-phase fluids. Depending on this, the jet pump may be classified as gas/gas, gas/liquid, liquid/liquid, two-phase/liquid, or similar combinations. The mixing chamber serves to create a homogeneous single-phase or two-phase mixture which enters a diffuser where the high kinetic energy of the fluid is converted into pressure energy. If the fluid mixture entering the diffuser is in the supersonic flow regime, a normal shock wave usually develops inside the diffuser. If the fluid mixture is one that can easily change phase, a condensation shock would normally develop. Because of the overall rise in pressure in the diffuser as well as the additional rise in pressure across the shock layer, condensation becomes more likely. Associated with the pressure rise across the shock is a velocity reduction from the supersonic to the subsonic range. If the two-phase flow entering the diffuser is predominantly gaseous with liquid droplets suspended in it, it will transform into a predominantly liquid flow containing gaseous bubbles (bubbly flow) somewhere in the diffuser. While past researchers have been able to model the two-phase flow jet pump using the one-dimensional assumption with no shock waves and no phase change, there is no research known to the authors apart from that of Anand (1992) which accounted for condensation shocks. One of the objectives of this research effort is to develop a comprehensive model in which the effects of phase slip and inter-phase heat transfer as well as the wall friction and shock waves are accounted for. While this modeling effort is predominantly analytical in nature and is primarily intended to provide a parametric understanding of the jet pump performance under different operating scenarios, another parallel effort employing a commercial CFD code is also implemented. The latter effort is primarily intended to model an axisymmetric counterpart of the problem in question. The viability of using the CFD code to model a two-phase flow jet pump will be assessed by attempting to recreate some of the existing performance data of similar jet pumps. The code will eventually be used to generate the jet pump performance characteristics of several scenarios involving jet pump geometries as well as flow regimes in order to be able to determine an optimum design which would be suitable for a two-phase flow boiling test facility at NASA-Marshall. Because of the extensive nature of the analytical model developed, the following section will only provide very brief highlights of it, while leaving the details to a more complete report submitted to the NASA colleague. This report will also contain some of the simulation results obtained using the CFD code.

Sherif, S. A.↗

Application of Parametrized Post-Newtonian Methods to the Gravitational IS of Satellite Energy Exchange Data

Project Satellite Energy Exchange (SEE) is a free-flying, high altitude satellite that utilizes space to construct a passive, low-temperature, nano-g environment in order to accurately measure the poorly known gravitational constant G plus other gravitational parameters that are difficult to measure in an earth-based laboratory. Eventually data received from SEE must be analyzed using a model of the gravitational interaction including parameters that describe deviations from general relativity and experiment. One model that can be used to fit tile data is the Parametrized post- Newtonian (PPN) approximation of general relativity (GR) which introduces ten parameters which have specified values in (GR). It is the lowest-order, consistent approximation that contains non linear terms. General relativity predicts that the Robertson parameters, gamma (light deflection), and beta (advance of the perihelion), are both 1 in GR. Another eight parameters, alpha(sub k), k=1,2,3 and zeta(sub k), k=1,2,3,4 and Xi are all zero in GR. Non zero values for alpha(sub k) parameters predict preferred frame effects; for zeta(sub k) violations of globally conserved quantities such as mass, momentum and angular momentum; and for Xi a contribution from the Whitehead theory of gravitation, once thought to be equivalent to GR. In addition, there is the possibility that there may be a preferred frame for the universe. If such a frame exists, then all observers must measure the velocity omega of their motion with respect to this universal rest frame. Such a frame is somewhat reminiscent of the concept of the ether which was supposedly the frame in which the velocity of light took the value c predicted by special relativity. The SEE mission can also look for deviations from the r(exp -2) law of Newtonian gravity, adding parameters alpha and lamda for non Newtonian behavior that describe the magnitude and range of the r(exp -2) deviations respectively. The foundations of the GR supposedly agree with Newtonian gravity to first order so that the parameters alpha and lamda are zero in GR. More important, however, GR subsequently depends on this Newtonian approximation to build up the non linear higher-order terms which forms the basis of the PPN frame work.

Smalley, Larry L.↗

NASA Microgravity Combustion Science Program

Combustion has been a subject of increasingly vigorous scientific research for over a century, not surprising considering that combustion accounts for approximately 85% of the world's energy production and is a key element of many critical technologies used by contemporary society. Although combustion technology is vital to our standard of living, it also poses great challenges to maintaining a habitable environment. A major goal of combustion research is production of fundamental (foundational) knowledge that can be used in developing accurate simulations of complex combustion processes, replacing current "cut-and-try" approaches and allowing developers to improve the efficiency of combustion devices, to reduce the production of harmful emissions, and to reduce the incidence of accidental uncontrolled combustion. With full understanding of the physics and chemistry involved in a given combustion process, including details of the unit processes and their interactions, physically accurate models which can then be used for parametric exploration of new combustion domains via computer simulation can be developed, with possible resultant definition of radically different approaches to accomplishment of various combustion goals. Effects of gravitational forces on earth impede combustion studies more than they impede most other areas of science. The effects of buoyancy are so ubiquitous that we often do not appreciate the enormous negative impact that they have had on the rational development of combustion science. Microgravity offers potential for major gains in combustion science understanding in that it offers unique capability to establish the flow environment rather than having it dominated by uncontrollable (under normal gravity) buoyancy effects and, through this control, to extend the range of test conditions that can be studied. It cannot be emphasized too strongly that our program is dedicated to taking advantage of microgravity to untangle complications caused by gravity, allowing major strides in our understanding of combustion processes and in subsequent development of improved combustion devices leading to improved quality of life on Earth. Fire and/or explosion events aboard spacecraft could be devastating to international efforts to expand the human presence in space. Testing to date has shown that ignition and flame spread on fuel surfaces (e.g., paper, wire insulation) behave quite differently under partial gravity and microgravity conditions. In addition, fire signatures-i.e., heat release, smoke production, flame visibility, and radiation-are now known to be quite different in reduced gravity environments; this research has provided data to improve the effectiveness of fire prevention practices, smoke and fire detectors, and fire extinguishment systems. The more we can apply our scientific and technological understanding to potential fire behavior in microgravity and partial gravity, the more assurance can be given to those people whose lives depend on the environment aboard spacecraft or eventually on habitats on the Moon or Mars.

King, Merrill K.↗

Turbulent Jet Flames Into a Vitiated Coflow

Examined is the vitiated coflow flame, an experimental condition that decouples the combustion processes of flows found in practical combustors from the associated recirculating fluid mechanics. The configuration consists of a 4.57 mm diameter fuel jet into a coaxial flow of hot combustion products from a lean premixed flame. The 210 mm diameter coflow isolates the jet flame from the cool ambient, providing a hot environment similar to the operating conditions of advanced combustors; this important high temperature element is lacking in the traditional laboratory experiments of jet flames into cool (room) air. A family of flows of increasing complexity is presented: 1) nonreacting flow, 2) all hydrogen flame (fuel jet and premixed coflow), and 3) set of methane flames. This sequence of experiments provides a convenient ordering of validation data for combustion models. Laser Raman-Rayleigh-LIF diagnostics at the Turbulent Diffusion Flame laboratory of Sandia National Laboratories produced instantaneous multiscalar point measurements. These results attest to the attractive features of the vitiated coflow burner and the well-defined boundary conditions provided by the coflow. The coflow is uniform and steady, isolating the jet flame from the laboratory air for a downstream distance ranging from z/d = 50-70. The statistical results show that differential diffusion effects in this highly turbulent flow are negligible. Complementing the comprehensive set of multiscalar measurements is a parametric study of lifted methane flames that was conducted to analyze flame sensitivity to jet and coflow velocity, as well as coflow temperature. The linear relationship found between the lift-off height and the jet velocity is consistent with previous experiments. New linear sensitivities were found correlating the lift-off height to coflow velocity and temperature. A blow-off study revealed that the methane flame blows off at a common coflow temperature (1260 K), regardless of coflow or jet velocity. An explanation for this phenomenon is that entrainment of ambient air at the high lift-off heights prevents autoignition. Analysis of the results suggests that flame stabilization occurs through a combination of flame propagation, autoignition, and localized extinction processes. Proposed is an expanded view of distributed reaction combustion based on analysis of the distributions of probe volume conditions at the stabilization region of the lifted hydrogen and methane flames. Turbulent eddies the size of the flame thickness mix fuel and hot coflow across the flame front, thereby enhancing the reaction zone with autoignition of reactants at elevated temperatures; this is the reverse effect of turbulent flames in ambient air, where intense turbulence in cool mixtures result in localized extinction. Each of the three processes (i.e., flame propagation, autoignition and localized extinction) contributes to flame stabilization in varying degrees, depending on flow conditions.

Holdeman, James D.↗

Initial Circulation and Peak Vorticity Behavior of Vortices Shed from Airfoil Vortex Generators

An extensive parametric study of vortices shed from airfoil vortex generators has been conducted to determine the dependence of initial vortex circulation and peak vorticity on elements of the airfoil geometry and impinging flow conditions. These elements include the airfoil angle of attack, chord length, span, aspect ratio, local boundary layer thickness, and free stream Mach number. In addition, the influence of airfoil-to-airfoil spacing on the circulation and peak vorticity has been examined for pairs of co-rotating and counter-rotating vortices. The vortex generators were symmetric airfoils having a NACA-0012 cross-sectional profile. These airfoils were mounted either in isolation, or in pairs, on the surface of a straight pipe. The turbulent boundary layer thickness to pipe radius ratio was about 17 percent. The circulation and peak vorticity data were derived from cross-plane velocity measurements acquired with a seven-hole probe at one chord-length downstream of the airfoil trailing edge location. The circulation is observed to be proportional to the free-stream Mach number, the angle-of-attack, and the span-to-boundary layer thickness ratio. With these parameters held constant, the circulation is observed to fall off in monotonic fashion with increasing airfoil aspect ratio. The peak vorticity is also observed to be proportional to the free-stream Mach number, the airfoil angle-of-attack, and the span-to-boundary layer thickness ratio. Unlike circulation, however, the peak vorticity is observed to increase with increasing aspect ratio, reaching a peak value at an aspect ratio of about 2.0 before falling off again at higher values of aspect ratio. Co-rotating vortices shed from closely spaced pairs of airfoils have values of circulation and peak vorticity under those values found for vortices shed from isolated airfoils of the same geometry. Conversely, counter-rotating vortices show enhanced values of circulation and peak vorticity when compared to values obtained in isolation. The circulation may be accurately modeled with an expression based on Prandtl's relationship between finite airfoil circulation and airfoil geometry. A correlation for the peak vorticity has been derived from a conservation relationship equating the moment at the airfoil tip to the rate of angular momentum production of the shed vortex, modeled as a Lamb (ideal viscous) vortex. This technique provides excellent qualitative agreement to the observed behavior of peak vorticity for low aspect ratio airfoils typically used as vortex generators.

Wendt, Bruce J.↗