用激波管驱动水喷雾

为了了解用激波管产生大流量水喷雾的流体力学原理，进行了实验研究。用可视化方法观察了管内两相流及所产生的喷雾的发展过程。 用PVDF压力传感器测量了液波管内的压力波形。结果表明，水柱内的压力波明显地不同于空气中的压力波；水柱被经过加速后，其雾化形态发生了新的变化。

酚醛树脂高温热解的激波管实验研究

采用激波管方法研究了酚醉树脂在温度1100K～1800K范围之间的热解过程.在激波管高温和短实验时间条件下,分析了实验样品颗粒在高温气相中的传热过程,讨论了样品颗粒达到热平衡的条件.通过色谱和质谱方法检测热解产物,获得了酚醛树脂高温热解产物分布和热解速率常数.酚醛树脂高温热解最主要的产物为水、一氧化碳、氢、乙炔和苯.温度1400 K将酚醛树脂热解分为高温和低温区,分别表现出不同的热解速率常数与温度的关系.

Boundary-layer transition on a slender cone in hypervelocity flow with real gas effects

The laminar to turbulent transition process in boundary layer flows in thermochemical nonequilibrium at high enthalpy is measured and characterized. Experiments are performed in the T5 Hypervelocity Reflected Shock Tunnel at Caltech, using a 1 m length 5-degree half angle axisymmetric cone instrumented with 80 fast-response annular thermocouples, complemented by boundary layer stability computations using the STABL software suite. A new mixing tank is added to the shock tube fill apparatus for premixed freestream gas experiments, and a new cleaning procedure results in more consistent transition measurements. Transition location is nondimensionalized using a scaling with the boundary layer thickness, which is correlated with the acoustic properties of the boundary layer, and compared with parabolized stability equation (PSE) analysis. In these nondimensionalized terms, transition delay with increasing CO₂ concentration is observed: tests in 100% and 50% CO₂, by mass, transition up to 25% and 15% later, respectively, than air experiments. These results are consistent with previous work indicating that CO₂ molecules at elevated temperatures absorb acoustic instabilities in the MHz range, which is the expected frequency of the Mack second-mode instability at these conditions, and also consistent with predictions from PSE analysis. A strong unit Reynolds number effect is observed, which is believed to arise from tunnel noise. N_Tr for air from 5.4 to 13.2 is computed, substantially higher than previously reported for noisy facilities. Time- and spatially-resolved heat transfer traces are used to track the propagation of turbulent spots, and convection rates at 90%, 76%, and 63% of the boundary layer edge velocity, respectively, are observed for the leading edge, centroid, and trailing edge of the spots. A model constructed with these spot propagation parameters is used to infer spot generation rates from measured transition onset to completion distance. Finally, a novel method to control transition location with boundary layer gas injection is investigated. An appropriate porous-metal injector section for the cone is designed and fabricated, and the efficacy of injected CO₂ for delaying transition is gauged at various mass flow rates, and compared with both no injection and chemically inert argon injection cases. While CO₂ injection seems to delay transition, and argon injection seems to promote it, the experimental results are inconclusive and matching computations do not predict a reduction in N factor from any CO₂ injection condition computed.

Resolução numérica de escoamentos compressíveis empregando um método de partículas livre de malhas e o processamento em paralelo (CUDA)

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem).

Autoignition of alcohols and ethers in a rapid compression machine

The autoignition characteristics of methanol, ethanol and MTBE (methyl tert-butyl ether) have been investigated in a rapid compression machine at pressures in the range 20-40 atm and temperatures within 750-1000 K. All three oxygenated fuels tested show higher autoignition temperatures than paraffins, a trend consistent with the high octane number of these fuels. The autoignition delay time for methanol was slightly lower than predicted values using reported reaction mechanisms. However, the experimental and measured values for the activation energy are in very good agreement around 44 kcal/mol. The measured activation energy for ethanol autoignition is in good agreement with previous shock tube results (31 kcal/mol), although ignition times predicted by the shock tube correlation are a factor of three lower than the measured values. The measured activation energy for MTBE, 41.4 kcal/mol, was significantly higher than the value previously observed in shock tubes (28.1 kcal/mol). The onset of preignition, characterized by a slow energy release prior to early ignition was observed in some instances. Possible reasons for these ocurrences are discussed. © Copyright 1993 Society of Automotive engineers, Inc.

Cabaret finite-difference schemes for the one-dimensional Euler equations

In the present paper we consider second order compact upwind schemes with a space split time derivative (CABARET) applied to one-dimensional compressible gas flows. As opposed to the conventional approach associated with incorporating adjacent space cells we use information from adjacent time layer to improve the solution accuracy. Taking the first order Roe scheme as the basis we develop a few higher (i.e. second within regions of smooth solutions) order accurate difference schemes. One of them (CABARET3) is formulated in a two-time-layer form, which makes it most simple and robust. Supersonic and subsonic shock-tube tests are used to compare the new schemes with several well-known second-order TVD schemes. In particular, it is shown that CABARET3 is notably more accurate than the standard second-order Roe scheme with MUSCL flux splitting.

激波管壁AlO自由基B~2∑~+-X~2∑~+和C~2∏_r-X~2∑~+带系辐射光谱的研究

爆轰驱动过程中产生的高温高压气流对铝质膜片、激波管壁产生烧蚀和冲刷作用,以致激波管壁、端盖上附有氧化铝等杂质,而高温下AlO自由基在气体分子的高速碰撞下被激发并产生强烈的辐射,从而干扰了高温气体辐射光谱的分析。用爆轰驱动加热技术将空气加热到4 000~7 000 K,利用多通道光学分析仪对AlO自由基辐射光谱进行分析,实验发现在460~530 nm波长范围内有多支辐射非常强烈的AlO自由基B2Σ+-X2Σ+(T00=20 689 cm-1)带系辐射谱带,且每支谱带都由多个带头组成,带头间隔约为2 nm,带头处于高频位置并向低频方向伸延。通过实验与理论计算相结合,重点分析了AlO自由基B2Σ+-X2Σ+带系辐射光谱的结构特征。AlO自由基C2Πr-X2Σ+(T00=33 047 cm-1)带系辐射光谱处于270~335 nm波长范围内,其辐射强度相对于B2Σ+-X2Σ+带系较弱,并且与OH基A2Σ+-X2Π(T00=32 682 cm-1)带系辐射光谱互相干扰而难以分辨,对该波段高温空气的辐射光谱分析产生不利的影响。

An efficient numerical scheme for the Euler equations

A finite difference scheme based on flux difference splitting is presented for the solution of the Euler equations for the compressible flow of an ideal gas. A linearised Riemann problem is defined, and a scheme based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency, leading to arithmetic averaging. This is in contrast to the usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. The scheme is applied to a shock tube problem and a blast wave problem. Each approximate solution compares well with those given by other schemes, and for the shock tube problem is in agreement with the exact solution.

MATERIAL MODELING AND ANALYSIS FOR THE DEVELOPMENT OF A REALISTIC BLAST HEADFORM

Blast traumatic brain injury (BTBI) has become an important topic of study because of the increase of such incidents, especially due to the recent growth of improvised explosive devices (IEDs). This thesis discusses a project in which laboratory testing of BTBI was made possible by performing blast loading on experimental models simulating the human head. Three versions of experimental models were prepared – one having a simple geometry and the other two having geometry similar to a human head. For developing the head models, three important parts of the head were considered for material modeling and analysis – the skin, skull and brain. The materials simulating skin, skull and brain went through many testing procedures including dynamic mechanical analysis (DMA). For finding a suitable brain simulant, several materials were tested under low and high frequencies. Step response analysis, rheometry and DMA tests were performed on materials such as water based gels, oil based mixtures and silicone gels cured at different temperatures. The gelatins and silicone gels showed promising results toward their use as brain surrogate materials. Temperature degradation tests were performed on gelatins, indicating the fast degradation of gelatins at room temperature. Silicone gels were much more stable compared to the water based gels. Silicone gels were further processed using a thinner-type additive gel to bring the dynamic modulus values closer to those of human brain matter. The obtained values from DMA were compared to the values for human brain as found in literature. Then a silicone rubber brain mold was prepared to give the brain model accurate geometry. All the components were put together to make the entire head model. A steel mount was prepared to attach the head for testing at the end of the shock tube. Instrumentation was implemented in the head model to obtain effective results for understanding more about the possible mechanisms of BTBI. The final head model was named the Realistic Explosive Dummy Head or the “RED Head.” The RED Head offered potential for realistic experimental testing in blast loading conditions by virtue of its material properties and geometrical accuracy.

Implementation of monotonic higher order upwind scheme in KIVA 4

KIVA is a FORTRAN code developed by Los Alamos national lab to simulate complete engine cycle. KIVA is a flow solver code which is used to perform calculation of properties in a fluid flow field. It involves using various numerical schemes and methods to solve the Navier-Stokes equation. This project involves improving the accuracy of one such scheme by upgrading it to a higher order scheme. The numerical scheme to be modified is used in the critical final stage calculation called as rezoning phase. The primitive objective of this project is to implement a higher order numerical scheme, to validate and verify that the new scheme is better than the existing scheme. The latest version of the KIVA family (KIVA 4) is used for implementing the higher order scheme to support handling the unstructured mesh. The code is validated using the traditional shock tube problem and the results are verified to be more accurate than the existing schemes in reference with the analytical result. The convection test is performed to compare the computational accuracy on convective transfer; it is found that the new scheme has less numerical diffusion compared to the existing schemes. A four valve pentroof engine, an example case of KIVA package is used as application to ensure the stability of the scheme in practical application. The results are compared for the temperature profile. In spite of all the positive results, the numerical scheme implemented has a downside of consuming more CPU time for the computational analysis. The detailed comparison is provided. However, in an overview, the implementation of the higher order scheme in the latest code KIVA 4 is verified to be successful and it gives better results than the existing scheme which satisfies the objective of this project.

Enhancing SPH using moving least-squares and radial basis functions

In this paper we consider two sources of enhancement for the meshfree Lagrangian particle method smoothed particle hydrodynamics (SPH) by improving the accuracy of the particle approximation. Namely, we will consider shape functions constructed using: moving least-squares approximation (MLS); radial basis functions (RBF). Using MLS approximation is appealing because polynomial consistency of the particle approximation can be enforced. RBFs further appeal as they allow one to dispense with the smoothing-length - the parameter in the SPH method which governs the number of particles within the support of the shape function. Currently, only ad hoc methods for choosing the smoothing-length exist. We ensure that any enhancement retains the conservative and meshfree nature of SPH. In doing so, we derive a new set of variationally-consistent hydrodynamic equations. Finally, we demonstrate the performance of the new equations on the Sod shock tube problem.

Pathogenesis of Fallopian tube damage caused by Chlamydia trachomatis infections

Chlamydia trachomatis is the leading cause of bacterial sexually transmitted disease worldwide resulting in 4–5 million new cases of Chlamydia annually and an estimated 100 million cases per annum. Infections of the lower female genital tract (FGT) frequently are asymptomatic so they often remain undiagnosed or untreated. If infections are either not resolved, or are left untreated, chlamydia can ascend to the upper FGT and infect the fallopian tubes (FTs) causing salpingitis that may lead to functional damage of the FTs and tubal factor infertility (TFI). Clinical observations and experimental data have indicated a role for antibodies against C. trachomatis proteins such as the 60 kDa heat-shock protein 60 (cHSP60) in the immunopathogenesis of TFI. When released from infected cells cHSP60 can induce pro-inflammatory immune responses that may functionally impair the FTs leading to fibrosis and luminal occlusion. Chlamydial pathogenesis of irreversible and permanent tubal damage is a consequence of innate and adaptive host immune responses to ongoing or repeated infections. The extracellular matrix (ECM) that is regulated by metalloproteinases (MMPs) may also be modified by chlamydial infections of the FGT. This review will highlight protective and pathogenic immune responses to ongoing and repeated chlamydial infections of the FGT. It will also present two recent hypotheses to explain mechanisms that may contribute to FT damage during a C. trachomatis infection. If Chlamydia immunopathology can be controlled it might yield a method of inducing fibrosis and thus provide a means of non-surgical permanent contraception for women.