Computational method to predict thermodynamic, transport, and flow properties for the modified Langley 8-foot high-temperature tunnel

The Langley 8 foot high temperature tunnel (8 ft HTT) is used to test components of hypersonic vehicles for aerothermal loads definition and structural component verification. The test medium of the 8 ft HTT is obtained by burning a mixture of methane and air under high pressure; the combustion products are expanded through an axisymmetric conical contoured nozzle to simulate atmospheric flight at Mach 7. This facility was modified to raise the oxygen content of the test medium to match that of air and to include Mach 4 and Mach 5 capabilities. These modifications will facilitate the testing of hypersonic air breathing propulsion systems for a wide range of flight conditions. A computational method to predict the thermodynamic, transport, and flow properties of the equilibrium chemically reacting oxygen enriched methane-air combustion products was implemented in a computer code. This code calculates the fuel, air, and oxygen mass flow rates and test section flow properties for Mach 7, 5, and 4 nozzle configurations for given combustor and mixer conditions. Salient features of the 8 ft HTT are described, and some of the predicted tunnel operational characteristics are presented in the carpet plots to assist users in preparing test plans

An IoT-Based Beehive Monitoring System for Real-Time Monitoring of Apis cerana indica Colonies

A study was conducted to monitor the bee activity in the colonies of diferente strengths in real time using an IoT-based device. The in-hive temperature and relative humidity were measured in the colonies of Apis cerana indica Fabricius of different strengths using the sensor-laden IoT device that was correlated with the movement of foragers into and out of the hive. A significantly higher movement of foragers was recorded at an in-hive temperature and relative humidity of 27.84 ºC and 61.47% at 5-6 p.m. with an observed activity of 9,638 bees/hive/hour in the strong colonies. In the weak colonies, the mean forager activity was 1,436.3 bees/hive/hour, which was recorded at an in-hive temperature of 26.52 ºC and 61.42% relative humidity. The mean honey area in the strong and weak colonies were 1,300.80±177.61 cm2 and 508.80±156.84 cm2, respectively. Pollen area in the strong and weak colonies were 447.60±112.08 cm2 and 116.20±66.43 cm2, respectively. In the strong and weak colonies, the area under egg brood was 470±53.06 cm2 and 88.20±36.85 cm2, larvae brood was 583.40±11.04 cm2 and 80.00±24.67 cm2 and sealed brood was 684.20±57.98 cm2 and 102.80±16.59 cm2, respectively. The real-time data on the movement of foragers in the colonies of different strengths enabled us to undertake timely intervention in the maintenance of the bee colonies

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors

Validation of two-equation turbulence models for propulsion flowfields

The objective of the study is to assess the capability of two-equation turbulence models for simulating propulsion-related flowfields. The standard kappa-epsilon model with Chien's low Reynolds number formulation for near-wall effects is used as the baseline turbulence model. Several experimental test cases, representative of rocket combustor internal flowfields, are used to catalog the performance of the baseline model. Specific flowfields considered here include recirculating flow behind a backstep, mixing between coaxial jets and planar shear layers. Since turbulence solutions are notoriously dependent on grid and numerical methodology, the effects of grid refinement and artificial dissipation on numerical accuracy are studied. In the latter instance, computational results obtained with several central-differenced and upwind-based formulations are compared. Based on these results, improved turbulence modes such as enhanced kappa-epsilon models as well as other two-equation formulations (e.g., kappa-omega) are being studied. In addition, validation of swirling and reacting flowfields are also currently underway

Development of a computational testbed for numerical simulation of combustion instability

A synergistic hierarchy of analytical and computational fluid dynamic techniques is used to analyze three-dimensional combustion instabilities in liquid rocket engines. A mixed finite difference/spectral procedure is employed to study the effects of a distributed vaporization zone on standing and spinning instability modes within the chamber. Droplet atomization and vaporization are treated by a variety of classical models found in the literature. A multi-zone, linearized analytical solution is used to validate the accuracy of the numerical simulations at small amplitudes for a distributed vaporization region. This comparison indicates excellent amplitude and phase agreement under both stable and unstable operating conditions when amplitudes are small and proper grid resolution is used. As amplitudes get larger, expected nonlinearities are observed. The effect of liquid droplet temperature fluctuations was found to be of critical importance in driving the instabilities of the combustion chamber

CFD modeling of microwave electrothermal thrusters

Microwave-heated plasmas in convergent nozzles are analyzed using a coupled Maxwell and Navier-Stokes solver to examine relevant issues associated with microwave thermal propulsion. Parametric studies are conducted to understand the effect of power, pressure, and plasma location with respect to the nozzle throat. For nozzles in the 0.5 to 3 N range with helium flow, results show that specific impulses up to 550-650 seconds are possible, with further increases being limited by severe wall-heating. Coupling efficiencies of over 90 percent are consistently obtained, with overall efficiencies ranging from 40 percent to 80 percent. Size scale-up studies-done by scaling the frequency from 2.45 GHz to 0.91 GHz-indicate that plasma migration toward the walls occurs more frequently for the lower frequency. Increasing the cavity aspect ratio and detuning the cavity are found to be effective ways of keeping the plasma on axis

Convergence acceleration of implicit schemes in the presence of high aspect ratio grid cells

The performance of Navier-Stokes codes are influenced by several phenomena. For example, the robustness of the code may be compromised by the lack of grid resolution, by a need for more precise initial conditions or because all or part of the flowfield lies outside the flow regime in which the algorithm converges efficiently. A primary example of the latter effect is the presence of extended low Mach number and/or low Reynolds number regions which cause convergence deterioration of time marching algorithms. Recent research into this problem by several workers including the present authors has largely negated this difficulty through the introduction of time-derivative preconditioning. In the present paper, we employ the preconditioned algorithm to address convergence difficulties arising from sensitivity to grid stretching and high aspect ratio grid cells. Strong grid stretching is particularly characteristic of turbulent flow calculations where the grid must be refined very tightly in the dimension normal to the wall, without a similar refinement in the tangential direction. High aspect ratio grid cells also arise in problems that involve high aspect ratio domains such as combustor coolant channels. In both situations, the high aspect ratio cells can lead to extreme deterioration in convergence. It is the purpose of the present paper to address the reasons for this adverse response to grid stretching and to suggest methods for enhancing convergence under such circumstances. Numerical algorithms typically possess a maximum allowable or optimum value for the time step size, expressed in non-dimensional terms as a CFL number or vonNeumann number (VNN). In the presence of high aspect ratio cells, the smallest dimension of the grid cell controls the time step size causing it to be extremely small, which in turn results in the deterioration of convergence behavior. For explicit schemes, this time step limitation cannot be exceeded without violating stability restrictions of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right. Primary among these are the proper definition of the time step size, proper selection of viscous preconditioner and the precise treatment of boundary conditions. These algorithmic improvements are then applied to a variety of test cases to demonstrate uniform convergence at all aspect ratios

Efficiency and reliability enhancements in propulsion flowfield modeling

The implementation of traditional CFD algorithms in practical propulsion related flowfields often leads to dramatic reductions in efficiency and/or robustness. The present research is directed at understanding the reasons for this deterioration and finding methods to circumvent it. Work to date has focussed on low Mach number regions, viscous dominated regions, and high grid aspect ratios. Time derivative preconditioning, improved definition of the local time stepping, and appropriate application of boundary conditions are employed to decrease the required time to obtain a solution, while maintaining accuracy. A number of cases having features typical of rocket engine flowfields are computed to demonstrate the improvement over conventional methods. These cases include laminar and turbulent high Reynolds number flat plate boundary layers, flow over a backward-facing step, a diffusion flame, and wall heat-flux calculations in a turbulent converging-diverging nozzle. Results from these cases show convergence that is virtually independent of the local Mach number and the grid aspect ratio, which translates to a convergence speed-up of up to several orders of magnitude over conventional algorithms. Current emphasis is in extending these results to three-dimensional flows with highly stretched grids

Microbial Community Structures of Novel Icelandic Hot Spring Systems Revealed by PhyloChip G3 Analysis

Microbial community profiles of recently formed hot spring systems ranging in temperatures from 57°C to 100°C and pH values from 2 to 4 in Hveragerði (Iceland) were analyzed with PhyloChip G3 technology. In total, 1173 bacterial operational taxonomic units (OTUs) spanning 576 subfamilies and 38 archaeal OTUs covering 32 subfamilies were observed. As expected, the hyperthermophilic (100°C) spring system exhibited both low microbial biomass and diversity when compared to thermophilic (60°C) springs. Ordination analysis revealed distinct bacterial and archaeal diversity in geographically distinct hot springs. Slight variations in temperature (from 57°C to 64°C) within the interconnected pools led to a marked fluctuation in microbial abundance and diversity. Correlation and PERMANOVA tests provided evidence that temperature was the key environmental factor responsible for microbial community dynamics, while pH, H_(2)S, and SO_2 influenced the abundance of specific microbial groups. When archaeal community composition was analyzed, the majority of detected OTUs correlated negatively with temperature, and few correlated positively with pH. Key Words: Microbial diversity—PhyloChip G3—Acidophilic—Thermophilic—Hot springs—Iceland. Astrobiology 14, xxx–xxx

Low frequency Raman studies of multi-wall carbon nanotubes: experiments and theory

In this paper, we investigate the low frequency Raman spectra of multi-wall carbon nanotubes (MWNT) prepared by the electric arc method. Low frequency Raman modes are unambiguously identified on purified samples thanks to the small internal diameter of the MWNT. We propose a model to describe these modes. They originate from the radial breathing vibrations of the individual walls coupled through the Van der Waals interaction between adjacent concentric walls. The intensity of the modes is described in the framework of bond polarization theory. Using this model and the structural characteristics of the nanotubes obtained from transmission electron microscopy allows to simulate the experimental low frequency Raman spectra with an excellent agreement. It suggests that Raman spectroscopy can be as useful regarding the characterization of MWNT as it is in the case of single-wall nanotubes.Comment: 4 pages, 2 eps fig., 2 jpeg fig., RevTex, submitted to Phys. Rev.