Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local mixing”

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 271 records · Page 15

On the thermal structure of Triton's thermosphere

The analysis of the Voyager 2 Ultraviolet Spectrometer solar occultation data obtained at Triton is consistent with a spherically symmetric, isothermal thermosphere above 400 km at T(infinity) = 96 K. A detailed calculation of energy loss processes in a pure N2 atmosphere indicates that solar heating, with calculated T(infinity) = 70 K, is insufficient to account for the inferred T(infinity) = 96 K. The magnetosphere must deposit twice as much power as the sun to heat the thermosphere to 96 K and generate the observed N2 tangential column densities above 450 km. The thermal escape of H and N atoms and the downward diffusion of N atoms to recombine below 130 km results in local ionospheric heating efficiency of 24 percent. An upper limit on the tropopause CO mixing ratio of 2 x 10 exp -4 is inferred in the absence of aerosol heating to balance its efficient cooling by LTE rotational line emission.

Stevens, Michael H.↗

Optimization of Orifice Geometry for Cross-Flow Mixing in a Cylindrical Duct

Mixing of gaseous jets in a cross-flow has significant applications in engineering, one example of which is the dilution zone of a gas turbine combustor. Despite years of study, the design of jet injection in combustors is largely based on practical experience. A series of experiments was undertaken to delineate the optimal mixer orifice geometry. A cross-flow to core-flow momentum-flux ratio of 40 and a mass flow ratio of 2.5 were selected as representative of an advanced design. An experimental test matrix was designed around three variables: the number of orifices, the orifice aspect ratio (long-to-short dimension), and the orifice angle. A regression analysis was performed on the data to arrive at an interpolating equation that predicted the mixing performance of orifice geometry combinations within the range of the test matrix parameters. Results indicate that mixture uniformity is a non-linear function of the number of orifices, the orifice aspect ratio, and the orifice angle. Optimum mixing occurs when the asymptotic mean jet trajectories are in the range of 0.35 less than r/R less than 0.5 (where r = 0 is at the mixer wall) at z/R = 1.0. At the optimum number of orifices, the difference between shallow-angled slots with large aspect ratios and round holes is minimal and either approach will lead to good mixing performance. At the optimum number of orifices, it appears possible to have two local optimums where one corresponds to an aspect ratio of 1.0 and the other to a high aspect ratio.

Sowa, W. A.↗

Experimental Study of Vane Heat Transfer and Film Cooling at Elevated Levels of Turbulence

This report documents the results of an experimental study on the influence of high level turbulence on vane film cooling and the influence of film cooling on vane heat transfer. Three different cooling configurations were investigated which included one row of film cooling on both pressure and suction surfaces, two staggered rows of film cooling on both suction and pressure surfaces, and a shower-head cooling array. The turbulence had a strong influence on film cooling effectiveness, particularly on the pressure surface where local turbulence levels were the highest. For the single row of holes, the spanwise mixing quickly reduced centerline effectiveness levels while mixing in the normal direction was more gradual. The film cooling had a strong influence on the heat transfer in the laminar regions of the vane. The effect of film cooling on heat transfer was noticeable in the turbulent regions but augmentation ratios were significantly lower. In addition to heat transfer and film cooling, velocity profiles were taken downstream of the film cooling rows at three spanwise locations. These profile comparisons documented the strong spanwise mixing due to the high turbulence. Total pressure exit measurements were also documented for the three configurations.

Ames, Forrest E.↗

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗

Modeling the Gas Dynamics Environment in a Subscale Solid Rocket Test Motor

Subscale test motors are often used for the evaluation of solid rocket motor component materials such as internal insulation. These motors are useful for characterizing insulation performance behavior, screening insulation material candidates and obtaining material thermal and ablative property design data. One of the primary challenges associated with using subscale motors however, is the uncertainty involved when extrapolating the results to full-scale motor conditions. These uncertainties are related to differences in such phenomena as turbulent flow behavior and boundary layer development, propellant particle interactions with the wall, insulation off-gas mixing and thermochemical reactions with the bulk flow, radiation levels, material response to the local environment, and other anomalous flow conditions. In addition to the need for better understanding of physical mechanisms, there is also a need to better understand how to best simulate these phenomena using numerical modeling approaches such as computational fluid dynamics (CFD). To better understand and model interactions between major phenomena in a subscale test motor, a numerical study of the internal flow environment of a representative motor was performed. Simulation of the environment included not only gas dynamics, but two-phase flow modeling of entrained alumina particles like those found in an aluminized propellant, and offgassing from wall surfaces similar to an ablating insulation material. This work represents a starting point for establishing the internal environment of a subscale test motor using comprehensive modeling techniques, and lays the groundwork for improving the understanding of the applicability of subscale test data to full-scale motors. It was found that grid resolution, and inclusion of phenomena in addition to gas dynamics, such as two-phase and multi-component gas composition are all important factors that can effect the overall flow field predictions.

Eaton, Andrew M.↗

Properties, Propagation, and Excitation of EMIC Waves Properties, Propagation, and Excitation of EMIC Waves

Electromagnetic ion cyclotron (EMIC) waves (0.1-5 Hz) play an important role in particle dynamics in the Earth's magnetosphere. EMIC waves are preferentially excited in regions where hot anisotropic ions and cold dense plasma populations spatially overlap. While the generation region of EMIC waves is usually on or near the magnetic equatorial plane in the inner magnetosphere, EMIC waves have both equatorial and off-equator source regions on the dayside in the compressed outer magnetosphere. Using field and plasma measurements from the Magnetospheric Multiscale (MMS) mission, we perform a case study of EMIC waves and associated local plasma conditions observed on 19 October 2015. From 0315 to 0810 UT, before crossing the magnetopause into the magnetosheath, all four MMS spacecraft detected long-lasting He(exp +)-band EMIC wave emissions around local noon (MLT = 12.7 - 14.0) at high L-shells (L = 8.8 - 15.2) and low magnetic latitudes (MLAT = -21.8deg - -30.3deg). Energetic (greater than 1 keV) and anisotropic ions were present throughout this event that was in the recovery phase of a weak geomagnetic storm (min. Dst = -48 nT at 1000 UT on 18 October 2015). The testing of linear theory suggests that the EMIC waves were excited locally. Although the wave event is dominated by small normal angles, its polarization is mixed with right- and left-handedness and its propagation is bi-directional with regard to the background magnetic field. The short inter-spacecraft distances (as low as ~15 km) of the MMS mission make it possible to accurately determine the k vector of the waves using the phase difference technique. Preliminary analysis finds that the k vector magnitude, phase speed, and wavelength of the 0.3-Hz wave packet at 0453:55 UT are 0.005 km(exp -1), 372.9 km/s, and 1242.9 km, respectively.

Zhang, Jichun↗

Partially Premixed Flame (PPF) Research for Fire Safety

Incipient fires typically occur after the partial premixing of fuel and oxidizer. The mixing of product species into the fuel/oxidizer mixture influences flame stabilization and fire spread. Therefore, it is important to characterize the impact of different levels of fuel/oxidizer/product mixing on flame stabilization, liftoff and extinguishment under different gravity conditions. With regard to fire protection, the agent concentration required to achieve flame suppression is an important consideration. The initial stage of an unwanted fire in a microgravity environment will depend on the level of partial premixing and the local conditions such as air currents generated by the fire itself and any forced ventilation (that influence agent and product mixing into the fire). The motivation of our investigation is to characterize these impacts in a systematic and fundamental manner.

Puri, Ishwar K.↗

Tidal Energy Available for Deep Ocean Mixing: Bounds from Altimetry Data

Maintenance of the large-scale thermohaline circulation has long presented a problem to oceanographers. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon (T/P) satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Egbert, Gary D.↗

Global seasonal variations of the near-surface relative humidity levels on present-day Mars

We investigate the global seasonal variations of near-surface relative humidity and relevant attributes, like temperature and water vapor volume mixing ratio on Mars using calculations from modelled and measurement data. We focus on 2 AM local time snapshots to eliminate daily effects related to differences in insolation, and to be able to compare calculations based on modelling data from the Laboratoire de Météorologie Dynamique Mars General Circulation Model with the observations of Mars Global Surveyor Thermal Emission Spectrometer. We study the seasonal effects by examining four specific dates in the Martian year, the northern spring equinox, summer solstice, autumn equinox, and winter solstice. We identify three specific zones, where the near-surface relative humidity levels are systematically higher than in their vicinity regardless of season. We find that these areas coincide with low thermal inertia features, which control surface temperatures on the planet, and are most likely covered with unconsolidated fine dust with grain sizes smaller than ∼ 40 μm. By comparing the data of relative humidity, temperature and water vapor volume mixing ratio at three different heights (near-surface, ∼ 4m and ∼ 23m above the surface), we demonstrate that the thermal inertia could play an important role in determining near-surface humidity levels. We also notice that during the night the water vapor levels drop at ∼ 4m above the surface. This, together with the temperature and thermal inertia values, shows that water vapor likely condenses in the near-surface atmosphere and on the ground during the night at the three aforementioned regions. This condensation may be in the form of brines, wettening of the fine grains by adsorption or deliquescence. This study specifies areas of interest on the surface of present day Mars for the proposed condensation, which may be examined by in-situ measurements in the future.

Bernadett Pál↗

Localization of Ad-Hoc Lunar Constellations in Communication Failure Modes for Distributed Spacecraft Autonomy

As lunar missions increase in complexity inspired by NASA’s Artemis Program, they will require reliable and sufficient capability of the Position, Navigation, and Timing (PNT) system to support their scientific objectives. In addition, NASA's Commercial Lunar Payload Services (CLPS) program initiates the proliferation of public and private exploration partnerships using small satellites from commercial and private organizations, expanding traditionally confined low Earth orbit to be used for missions beyond geosynchronous orbit (Zucherman et al., 2022). Therefore, the Lunar PNT system is also required to provide navigation services compatible with the smaller platforms being sent by the public and private sectors, like CubeSats. However, traditional approaches to deep space missions’ navigation based on ground radio facilities have difficulties in providing sufficient support for the increasing number of users and communication at a distance from the Earth (Kaplev et al., 2022). In particular, the existing Lunar navigation technologies such as weak signal global positioning system (GPS) and deep space network (DSN) are not able to ensure operations of the upcoming small-scale Lunar missions due to their limitations in localization performance as well as capacity aspects. Another way to provide Lunar PNT service is to create a dedicated Lunar global navigation satellite system (GNSS) constellation, like GNSS systems on Earth. Space agencies like NASA, ESA, and JAXA are now developing the lunar communications relay and navigation systems (LCRNS) and Lunar navigation satellite systems (LNSS). In their systems, satellites will be deployed in moon orbits to provide the communication, positioning, navigation, and timing (CPNT) service at the lunar south pole region where the Artemis base camp will be expected (Murata et al., 2022). Meanwhile, common challenges considered in lunar PNT research arise from poor geometry of the terrestrial GNSS satellites when seen from the lunar user, highly perturbed lunar orbits, and limitations in power, size, and cost of the equipment on lunar satellites (Iiyama et al., 2023). It is also not clear if there will be enough Lunar users to support the cost and resources this would require as the Low-cost surface missions may not be able to support the large power, mass, and weight requirements that these navigation solutions entail (Niemoeller et al., 2022). As an alternative, existing Lunar science and exploration assets could be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions (Hagenau et al., 2021). Introducing the non-dedicated and ad-hoc Lunar navigation constellation gives a way to provide PNT services on-demand. The non-dedicated swarm assets of Lunar constellations are designed to localize themselves with minimal interaction with Earth by adding cooperative autonomous localization to lunar missions, freeing up valuable bandwidth and ground segment resources. An autonomous localization of Lunar constellations is based on the concept of the decentralized PNT system with a distributed extended Kalman filter (DEKF) approach to state estimation for minimal onboard operating costs. In the distributed data processing algorithm, computation is broken down and assigned to each satellite, resulting in a considerably decreased computational amount while maintaining the accuracy of the orbit ephemeris and clock offsets as the result of centralized data processing (Wen et al., 2019). The DEKF requires spacecraft to perform two-way ranging operations with each other to communicate simultaneously, leveraging neighbor two-way intersatellite link (ISL) measurements such as pseudoranges to, and relative velocities between, visible satellites as sensor values (Frank et al., 2021). The Lunar autonomous PNT simulation (LAPS) demonstrated the feasibility of orbital asset localization among ad-hoc Lunar small-sat constellations based on the DEKF in Hagenau et al. (2021) and evaluated the matching algorithm proposed by Frank et al. (2021) in scheduling position estimation updates. In previous papers, all assets and measurements are assumed to be always available without consideration of the impact of intermittent and permanent communication failure. This study presents localization performance with increasing levels of network degradation for swarm assets and users to demonstrate the robustness of the decentralized Lunar PNT service in more realistic scenarios. Main issues arising from communication failure include spacecraft permanent or transient loss, antenna failures, message delays, etc. We tested four possible reasons for network degradation for 7 days in 21 satellites frozen with an altitude of 5500 km, evenly spaced around 3 circular, 40 inclination orbital planes where each spacecraft has two directional antennas. As anchor nodes with an independent estimate of their position are required in the DEKF approach, two ground nodes in each pole and one node in the gateway were implemented in the simulation. First, the most probable failure scenario involves the loss of a single spacecraft due to solar interference and technical malfunctions of the assets. Losing the availability of a single spacecraft means losing the two-way ISL measurement of the asset in the DEKF update. In order to provide the best possible quality of PNT service with limited time and resources, the distributed Lunar constellations must schedule the communication activities. The scheduler leverages mixed-integer linear programming (MILP) for the coordination and scheduling of the desired “as-needed” localization service (Niemoeller et al., 2022). We assume the scheduler has completely excluded the spacecraft information before the DEKF update in the failure scenario. When a random spacecraft has been turned off at a specific time, the robustness of the autonomous Lunar PNT system is evaluated. The simulation results give an 11.5% degradation in median position accuracy compared to the idealized performance excluding the asset loss. Second, a large number of assets may vanish due to major hardware problems or meteor strikes around the moon. A multiple spacecraft loss can degrade the localization performance very fast by losing the communication ability to do cross-plane measurements and in-plane measurements in a 3-plane constellation. When the matching-based scheduler is aware of ISL availability, we investigate a large number of in-plane and cross-plane asset vanishments both in close proximity and equally spaced throughout the orbital plane. According to the simulations, the loss of in-plane measurements gives 40.2% degradation while cross-plane measurements degrade 50.5% of asset localization performance among available assets. Therefore, it is concluded that cross-plane measurements are more important in improving the position estimation accuracy. Third, spacecraft failure information can be lost due to the internal message delay, resulting in the DEKF update scheduler to solve the matching problem with unavailable assets. The DEKF update cycle is comprised of network setup, communication, and computations where a global broadcast network and a 2-way ISL network setup take 6 minutes in total (Frank et al., 2021). Once the broadcast network successfully transmits and receives information, a random spacecraft may lose its availability right before solving the matching problem. This means the matching solution is no longer optimal, resulting in degradation in the localization performance. A numerical assessment shows the matching-based scheduler with knowing failure holds 11.5% of position accuracy degradation, whereas the scheduler without knowing failure gives 34% degraded localization performance without asset loss. Fourth, a transient loss of a single or multiple spacecraft may occur due to their antenna outages. After losing the two-way ISL availability for a few DEKF update cycles, the availability of spacecraft can easily be recovered as their states have been independently updated using measurements from anchor nodes. It is likely that the longer failure will result in worse localization performance. We have tested the transient failure of a random single asset for 30 min in the simulation, which is losing 3 update cycles in the DEKF system. From the simulation results, the position accuracy has been degraded to 4.84% which is better than the degraded localization performance of 11.5% from the permanent loss scenario among available assets. In conclusion, the autonomous Lunar PNT system based on the DEKF approach shows the ability to maintain resilience and robustness in the possible communication failure scenarios, ensuring that localization accuracy is preserved across various network degradation and outages. Future studies on investigating user localization performance near the South Pole and the broadcast network system will be continued in the following months.

Yeji Kim↗

The Generation of Mixing Across the Stratospheric Polarvortex Edge By Breaking Gravity Waves

Modern high resolution global atmospheric analyses are capable of resolving larger gravity waves, including their generation, propagation, and breaking. Orographically generated gravity waves that break in the middle stratosphere often appear as disturbances in maps of potential vorticity, a key stratospheric dynamical quantity. Using analyses from the NASA Global Modeling and Assimilation Office, we show how gravity waves not only disturb the local flow above them but can generate flow instabilities leading to eddies that propagate around the stratospheric polar vortex, creating significant mixing across the vortex edge and into the polar region. These disturbances and mixing events are especially noticeable during NH winter months with low planetary scale wave activity such as in January 2022. During January 2022, the initially circular flow around the polar vortex was disrupted by breaking gravity waves, leading to unstable flow that created a chain of eddies around the vortex edge. The mixing associated with this process led to increased ozone in the upper stratosphere polar region.

Polar Vortex↗

The Western North American Cretaceous-Tertiary (K-T) boundary interval and its content of shock-metamorphosed minerals: Implications concerning the K-T boundary impact-extinction theory

At 20 sites in the Raton Basin of Colorado and New Mexico, and at several other sites in Wyoming, Montana, and Canada, a pair of claystone units, an Ir abundance anomaly, and a concentration of shock-metamorphosed minerals mark the palynological K-T boundary. The K-T boundary claystone, which is composed of kaolinite and small amounts of illite/smectite mixed-layer clay, is similar in most respects to kaolinite tonstein layers in coal beds. At some, but not all, K-T boundary localities, the boundary claystone contains solid kaolinite and hollow and solid goyazite spherules, 0.05 to 1.2 mm in diameter. The upper unit, the K-T boundary impact layer, consists chiefly of kaolinite and various amounts of illite/smectite mixed-layer clay. The impact layer and boundary claystone are similar chemically, except that the former has slightly more Fe, K, Ba, Cr, Cu, Li, V, and Zn than the latter. The facts that the boundary claystone and impact layer contain anomalous amounts of Ir, comprise a stratigraphic couplet at Western North American sites, and form thin, discrete layers, similar to air-fall units (volcanic or impact), suggest that the claystone units are of impact origin. Significantly, the impact layer contains as much as 2 percent clastic mineral grains, about 30 percent of which contain multiple sets of shock lamellae. Only one such concentration of shocked minerals has been found near the K-T boundary. The type of K-T boundary shock-metamorphosed materials (quartzite and metaquartzite) in the impact layer and the lack of shock lamellae in quartz and feldspar of pumice lapilli and granitic xenoliths in air-fall pumice units of silicic tuffs, such as the Bishop Tuff, eliminate the possibility that the shock-metamorphosed minerals in the K-T impact layer are of volcanic origin. The global size distribution and abundance of shock-metamorphosed mineral grains suggest that the K-T impact occurred in North America.

Izett, G. A.↗

Dynamo magnetic field modes in thin astrophysical disks - An adiabatic computational approximation

An adiabatic approximation is applied to the calculation of turbulent MHD dynamo magnetic fields in thin disks. The adiabatic method is employed to investigate conditions under which magnetic fields generated by disk dynamos permeate the entire disk or are localized to restricted regions of a disk. Two specific cases of Keplerian disks are considered. In the first, magnetic field diffusion is assumed to be dominated by turbulent mixing leading to a dynamo number independent of distance from the center of the disk. In the second, the dynamo number is allowed to vary with distance from the disk's center. Localization of dynamo magnetic field structures is found to be a general feature of disk dynamos, except in the special case of stationary modes in dynamos with constant dynamo number. The implications for the dynamical behavior of dynamo magnetized accretion disks are discussed and the results of these exploratory calculations are examined in the context of the protosolar nebula and accretion disks around compact objects.

Stepinski, T. F.↗

A simple apparatus for the experimental study of non-steady flow thrust-augmenter ejector configurations

A flexible, and easily modified, test rig is described which allows a one dimensional nonsteady flow stream to be generated, economically from a steady flow source of compressed air. This nonsteady flow is used as the primary stream in a nonsteady flow ejector constituting part of the test equipment. Standard piezo-electric pressure transducers etc. allow local pressures to be studied, as functions of time, in both the primary and secondary (mixed) flow portions of the apparatus. Provision is also made for measuring the primary and secondary mass flows and the thrust generated. Sample results obtained with the equipment are presented.

Khare, J. M.↗

Comparison of soil dielectric mixing models for Soil Moisture Retrieval using SMAP Brightness Temperature over croplands in India

The accurate estimation of soil moisture (SM) using microwave remote sensing depends mostly on careful selection of retrieval parameters among which the soil dielectric mixing model is the important one. These models are often categorized into empirical, semi-empirical or volumetric based on their methodologies and input data requirements. To study in detail, the comparative performance of four dielectric mixing models -- Wang & Schmugge model, Hallikainen model, Dobson model and Mironov model were used with Soil Moisture Active Passive (SMAP) L-band brightness temperature and Single Channel Algorithm for SM retrieval over agricultural landscapes in India. The highest performance statistics combination in terms of Root Mean Square Error (RMSE), correlation coefficient (R^(2)) and percentage bias (PBIAS) against the concurrent in-situ SM measurements were calculated at the selected validation sites. The overall results indicate that the best performance was given by the Mironov model (RMSE = 0.07 cu. m/cu. m), followed by Wang & Schmugge model (RMSE = 0.08 cu. m3/cu. m), Hallikainen model (RMSE = 0.09 cu. m/cu. m), Dobson model (RMSE = 0.10 cu. m/cu. m) and original SMAP radiometer SM (RMSE = 0.12 cu. m/cu. m). Findings of this study provides important insights into application and performance of dielectric mixing models in mapping surface SM variations. This study also underlines the pivotal role of local conditions for SM retrieval which should be carefully included in the algorithms.

Swati Suman↗

Adaptable Constrained Genetic Programming: Extensions and Applications

An evolutionary algorithm applies evolution-based principles to problem solving. To solve a problem, the user defines the space of potential solutions, the representation space. Sample solutions are encoded in a chromosome-like structure. The algorithm maintains a population of such samples, which undergo simulated evolution by means of mutation, crossover, and survival of the fittest principles. Genetic Programming (GP) uses tree-like chromosomes, providing very rich representation suitable for many problems of interest. GP has been successfully applied to a number of practical problems such as learning Boolean functions and designing hardware circuits. To apply GP to a problem, the user needs to define the actual representation space, by defining the atomic functions and terminals labeling the actual trees. The sufficiency principle requires that the label set be sufficient to build the desired solution trees. The closure principle allows the labels to mix in any arity-consistent manner. To satisfy both principles, the user is often forced to provide a large label set, with ad hoc interpretations or penalties to deal with undesired local contexts. This unfortunately enlarges the actual representation space, and thus usually slows down the search. In the past few years, three different methodologies have been proposed to allow the user to alleviate the closure principle by providing means to define, and to process, constraints on mixing the labels in the trees. Last summer we proposed a new methodology to further alleviate the problem by discovering local heuristics for building quality solution trees. A pilot system was implemented last summer and tested throughout the year. This summer we have implemented a new revision, and produced a User's Manual so that the pilot system can be made available to other practitioners and researchers. We have also designed, and partly implemented, a larger system capable of dealing with much more powerful heuristics.

Janikow, Cezary Z.↗

Small hydrogen/oxygen rocket flowfield behavior from heat flux measurements

The mixing and heat transfer phenomena in small rocket flow fields with fuel film cooling is not well understood. An instrumented, water-cooled chamber with a gaseous hydrogen/gaseous oxygen injector was used to gather steady-state inner and outer wall temperature profiles. The chamber was tested at 414 kPa (60 psia) chamber pressure, from mixture ratios of 3.41 to 8.36. Sixty percent of the fuel was used for film cooling. These temperature profiles were used as boundary conditions in a finite element analysis program, MSC/NASTRAN, to calculate the local radial and axial heat fluxes in the chamber wall. The normal heat fluxes were then calculated and used as a diagnostic of the rocket's flow field behavior. The normal heat fluxes determined were on the order of 1.0 to 3.0 MW/meters squared (0.6 to 1.8 Btu/sec-inches squared). In the cases where mixture ratio was 5 or above, there was a sharp local heat flux maximum in the barrel section of the chamber. This local maximum seems to indicate a reduction or breakdown of the fuel film cooling layer, possibly due to increased mixing in the shear layer between the film and core flows. However, the flow was thought to be completely laminar, as the throat Reynolds numbers were below 50,000 for all the cases. The increased mixing in the shear layer in the higher mixture ratio cases appeared not to be due to the transition of the flow from laminar to turbulent, but rather due to increased reactions between the hydrogen film and oxidizer-rich core flows.

Reed, Brian D.↗

A Shell/3D Modeling Technique for Delaminations in Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗