Discrete Vortex Method-Based Model for Ground-Effect Studies

A discrete vortex method-based model has been proposed for two-dimensional/three-dimensional ground-effect prediction. The model merely requires two-dimensional sectional aerodynamics in free flight. This free-flight data can be obtained either from experiments or a high-fidelity computational fluid dynamics solver. The first step of this two-step model involves a constrained optimization procedure that modifies the vortex distribution on the camber line as obtained from a discrete vortex method to match the free-flight data from experiments/computational fluid dynamics. In the second step, the vortex distribution thus obtained is further modified to account for the presence of the ground plane within a discrete vortex method-based framework. Whereas the predictability of the lift appears as a natural extension, the drag predictability within a potential flow framework is achieved through the introduction of what are referred to as drag panels. The need for the use of the generalized Kutta-Joukowski theorem is emphasized. The extension of the model to three dimensions is by the way of using the numerical lifting-line theory that allows for wing sweep. The model is extensively validated for both two-dimensional and three-dimensional ground-effect studies. The work also demonstrates the ability of the model to predict lift and drag coefficients of a high-lift wing in ground effect to about 2 and 8% accuracy, respectively, as compared to the results obtained using a Reynolds-averaged Navier-Stokes solver involving grids with several million volumes. The model shows a lot of promise in design, particularly during the early phase.

Investigation of the separated region ahead of three-dimensional protuberances on plates and cones in hypersonic flows with laminar boundary layers

Heat transfer rate and pressure measurements were made upstream of surface pro-tuberances on a flat plate and a sharp cone subjected to hypersonic flow in a conventional shock tunnel. Heat flux was measured using platinum thin-film sensors deposited on macor substrate and the pressure measurements were made using fast acting piezoelectric sensors. A distinctive hot spot with highest heat flux was obtained near the foot of the protuberance due to heavy vortex activity in the recirculating region. Schlieren flow visualization was used to capture the shock structures and the separation distance ahead of the protrusions was quantitatively measured for varying protuberance heights. A computational analysis was conducted on the flat plate model using commercial computational fluid dynamics software and the obtained trends of heat flux and pressure were compared with the experimental observation. Experiments were also conducted by physically disturbing the laminar boundary layer to check its effect on the magnitude of the hot spot heat flux. In addition to air, argon was also used as test gas so that the Reynolds number can be varied. (C) 2014 AIP Publishing LLC.

Stability of the flow in a soft tube deformed due to an applied pressure gradient

A linear stability analysis is carried out for the flow through a tube with a soft wall in order to resolve the discrepancy of a factor of 10 for the transition Reynolds number between theoretical predictions in a cylindrical tube and the experiments of Verma and Kumaran J. Fluid Mech. 705, 322 (2012)]. Here the effect of tube deformation (due to the applied pressure difference) on the mean velocity profile and pressure gradient is incorporated in the stability analysis. The tube geometry and dimensions are reconstructed from experimental images, where it is found that there is an expansion and then a contraction of the tube in the streamwise direction. The mean velocity profiles at different downstream locations and the pressure gradient, determined using computational fluid dynamics, are found to be substantially modified by the tube deformation. The velocity profiles are then used in a linear stability analysis, where the growth rates of perturbations are calculated for the flow through a tube with the wall modeled as a neo-Hookean elastic solid. The linear stability analysis is carried out for the mean velocity profiles at different downstream locations using the parallel flow approximation. The analysis indicates that the flow first becomes unstable in the downstream converging section of the tube where the flow profile is more pluglike when compared to the parabolic flow in a cylindrical tube. The flow is stable in the upstream diverging section where the deformation is maximum. The prediction for the transition Reynolds number is in good agreement with experiments, indicating that the downstream tube convergence and the consequent modification in the mean velocity profile and pressure gradient could reduce the transition Reynolds number by an order of magnitude.

CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

An Optimizing Code Generator for a Class of Lattice-Boltzmann Computations

The Lattice-Boltzmann method (LBM), a promising new particle-based simulation technique for complex and multiscale fluid flows, has seen tremendous adoption in recent years in computational fluid dynamics. Even with a state-of-the-art LBM solver such as Palabos, a user has to still manually write the program using library-supplied primitives. We propose an automated code generator for a class of LBM computations with the objective to achieve high performance on modern architectures. Few studies have looked at time tiling for LBM codes. We exploit a key similarity between stencils and LBM to enable polyhedral optimizations and in turn time tiling for LBM. We also characterize the performance of LBM with the Roofline performance model. Experimental results for standard LBM simulations like Lid Driven Cavity, Flow Past Cylinder, and Poiseuille Flow show that our scheme consistently outperforms Palabos-on average by up to 3x while running on 16 cores of an Intel Xeon (Sandybridge). We also obtain an improvement of 2.47x on the SPEC LBM benchmark.

Implicit scheme for meshless compressible Euler solver

In this paper, an implicit scheme is presented for a meshless compressible Euler solver based on the Least Square Kinetic Upwind Method (LSKUM). The Jameson and Yoon's split flux Jacobians formulation is very popular in finite volume methodology, which leads to a scalar diagonal dominant matrix for an efficient implicit procedure (Jameson & Yoon, 1987). However, this approach leads to a block diagonal matrix when applied to the LSKUM meshless method. The above split flux Jacobian formulation, along with a matrix-free approach, has been adopted to obtain a diagonally dominant, robust and cheap implicit time integration scheme. The efficacy of the scheme is demonstrated by computing 2D flow past a NACA 0012 airfoil under subsonic, transonic and supersonic flow conditions. The results obtained are compared with available experiments and other reliable computational fluid dynamics (CFD) results. The present implicit formulation shows good convergence acceleration over the RK4 explicit procedure. Further, the accuracy and robustness of the scheme in 3D is demonstrated by computing the flow past an ONERA M6 wing and a clipped delta wing with aileron deflection. The computed results show good agreement with wind tunnel experiments and other CFD computations.

On the control of rapidly rotating convection by an axially varying magnetic field

The magnetic field in rapidly rotating dynamos is spatially inhomogeneous. The axial variation of the magnetic field is of particular importance because tall columnar vortices aligned with the rotation axis form at the onset of convection. The classical picture of magnetoconvection with constant or axially varying magnetic fields is that the Rayleigh number and wavenumber at onset decrease appreciably from their non-magnetic values. Nonlinear dynamo simulations show that the axial lengthscale of the self-generated azimuthal magnetic field becomes progressively smaller as we move towards a rapidly rotating regime. With a small-scale field, however, the magnetic control of convection is different from that in previous studies with a uniform or large-scale field. This study looks at the competing viscous and magnetic mode instabilities when the Ekman number E (ratio of viscous to Coriolis forces) is small. As the applied magnetic field strength (measured by the Elsasser number Lambda) increases, the critical Rayleigh number for onset of convection initially increases in a viscous branch, reaches an apex where both viscous and magnetic instabilities co-exist, and then falls in the magnetic branch. The magnetic mode of onset is notable for its dramatic suppression of convection in the bulk of the fluid layer where the field is weak. The viscous-magnetic mode transition occurs at Lambda similar to 1, which implies that small-scale convection can exist at field strengths higher than previously thought. In spherical shell dynamos with basal heating, convection near the tangent cylinder is likely to be in the magnetic mode. The wavenumber of convection is only slightly reduced by the self-generated magnetic field at Lambda similar to 1, in agreement with previous planetary dynamo models. The back reaction of the magnetic field on the flow is, however, visible in the difference in kinetic helicity between cyclonic and anticyclonic vortices.

Thermal Modeling of Organic Light-Emitting Diode Display Panels

Power densities required to operate active-matrix organic-light-emitting diode (AMOLED) based displays for high luminance applications, lead to temperature rise due to self heating. Temperature rise leads to significant degradation and consequent reduction in life time. In this work numerical techniques based computational fluid dynamics (CFD) is used to determine the temperature rise and its distribution for an AMOLED based display for a given power density and size. Passive cooling option in form of protruded rectangular fins is implemented to reduce the display temperature.

Numerical Investigations for Leakage and Windage Heating in Straight-Through Labyrinth Seals

The ability to quantify leakage flow and windage heating for labyrinth seals with honeycomb lands is critical in understanding gas turbine engine system performance and predicting its component life. Variety of labyrinth seal configurations (number of teeth, stepped or straight, honeycomb cell size) are in use in gas turbines, and for each configuration, there are many geometric factors that can impact a seal's leakage and windage characteristics. This paper describes the development of a numerical methodology aimed at studying the effect of honeycomb lands on leakage and windage heating. Specifically, a three-dimensional computational fluid dynamics (CFD) model is developed utilizing commercial finite volume-based software incorporating the renormalization group (RNG) k-epsilon turbulence model with modified Schmidt number. The modified turbulence model is benchmarked and fine-tuned based on several experiments. Using this model, a broad parametric study is conducted by varying honeycomb cell size, pressure ratio (PR), and radial clearance for a four-tooth straight-through labyrinth seal. The results show good agreement with available experimental data. They further indicate that larger honeycomb cells predict higher seal leakage and windage heating at tighter clearances compared to smaller honeycomb cells and smooth lands. However, at open seal clearances larger honeycomb cells have lower leakage compared to smaller honeycomb cells.

DESIGN REQUIREMENTS FOR DIRECT SUPERCRITICAL CARBON DIOXIDE RECEIVER DEVELOPMENT AND TESTING

This paper establishes the design requirements for the development and testing of direct supercritical carbon dioxide (sCO2) solar receivers. Current design considerations are based on the ASME Boiler and Pressure Vessel Code (BPVC). Section I (BPVC) considers typical boilers/superheaters (i.e. fired pressure vessels) which work under a constant low heat flux. Section VIII (BPVC) considers pressure vessels with operating pressures above 15 psig 2 bar] (i.e. unfired pressure vessels). Section III, Division I - Subsection NH (BPVC) considers a more detailed stress calculation, compared to Section I and Section VIII, and requires a creep-fatigue analysis. The main drawback from using the BPVC exclusively is the large safety requirements developed for nuclear power applications. As a result, a new set of requirements is needed to perform detailed thermal-structural analyses of solar thermal receivers subjected to a spatially-varying, high-intensity heat flux. The last design requirements document of this kind was an interim Sandia report developed in 1979 (SAND79-8183), but it only addresses some of the technical challenges in early-stage steam and molten-salt solar receivers but not the use of sCO2 receivers. This paper presents a combination of the ASME BPVC and ASME B31.1 Code modified appropriately to achieve the reliability requirements in sCO(2) solar power systems. There are five main categories in this requirements document: Operation and Safety, Materials and Manufacturing, Instrumentation, Maintenance and Environmental, and General requirements. This paper also includes the modeling guidelines and input parameters required in computational fluid dynamics and structural analyses utilizing ANSYS Fluent, ANSYS Mechanical, and nCode Design Life. The main purpose of this document is to serve as a reference and guideline for design and testing requirements, as well as to address the technical challenges and provide initial parameters for the computational models that will be employed for the development of sCO(2) receivers.

流体动力学的研究:面向国家需求

流体动力学这门科学描述液体和气体的运动及其和固体的相互作用.它几乎涉及日常生活的方方面面,并在大部分科学和工程领域中处于核心地位,所以,这是一门宽广的多学科的学科领域.流体力学影响着国防、国土安全、交通、制造业、医学、生物学、能源和环境科学.预测人体内血液的流动,微流动器件装置的行为,飞机、汽车、船舶的流体动力学性能,电子元件的冷却,或者灾害性天气和气候的预报等,都需要对流体动力学有细致的了解和深入的研究.

Shock wave diffraction by a square cavity filled with dusty gas

A simple two-dimensional square cavity model is used to study shock attenuating effects of dust suspension in air. The GRP scheme for compressible flows was extended to simulate the fluid dynamics of dilute dust suspensions, employing the conventional two-phase approximation. A planar shock of constant intensity propagated in pure air over Aat ground and diffracted into a square cavity filled with a dusty quiescent suspension. Shock intensities were M-s = 1.30 and M-s = 2.032, dust loading ratios were alpha = 1 and alpha = 5, and particle diameters were d = 1, 10 and 50 mum. It was found that the diffraction patterns in the cavity were decisively attenuated by the dust suspension, particularly for the higher loading ratio. The particle size has a pronounced effect on the flow and wave pattern developed inside the cavity. Wall pressure historics were recorded for each of the three cavity walls, showing a clear attenuating effect of the dust suspension.

Pulverized-coal combustion in Wall-Protecting-Jets combustor (WPJC)

A new type of pulverized-coal combustor, called "Wall-Protecting-Jets Combustor" (hereafter, WPJC has been proposed, designed and studied with both CFD (Computational Fluid Dynamics) and experimental methods. The WPJC is based on a novel concept in which all inlet jets are along the combustor wall. Pilot combustion experiments were conducted to investigate the combustion performance of WPJC. Two-phase flows and pulverized-coal combustion were simulated to study the mechanism of),WPJC using the commercial software FLUENT. The results show that the WPJC has many remarkable advantages: wall-protection by the cold jets without the use of refractory materials; low-temperature and three-stage combustion with low NOx emission; negligible ash/slag-deposition; multiple functions with convenient switching between them; effective adjustment of the combustion intensity and the ignition position.

Modeling of Ammonothermal Growth of Nitrides

Bulk single crystals of GaN and AlN can be grown from supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. In this process, a mineralizer such as amide, imide or nitride is used to attack a bulk nitride feedstock at temperatures from 200°C to 500°C and pressures from 1 to 4 kbar. Ammonothermal systems have been modeled here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the fluid flow is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of particle size on flow pattern and temperature distribution in an autoclave are analyzed.

Characterisation of cell adhesion to substrate materials and the resistance to enzymatic and mechanical cell-removal

Cell-implant adhesive strength is important for prostheses. In this paper, an investigation is described into the adhesion of bovine chondrocytes to Ti6Al4V-based substrates with different surface roughnesses and compositions. Cells were cultured for 2 or 5 days, to promote adhesion. The ease of cell removal was characterised, using both biochemical (trypsin) and mechanical (accelerated buoyancy and liquid flow) methods. Computational fluid dynamics (CFD) modelling has been used to estimate the shear forces applied to the cells by the liquid flow. A comparison is presented between the ease of cell detachment indicated using these methods, for the three surfaces investigated. © 2008 Materials Research Society.