Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Direct Numerical Simulation of Shock/Turbulent Boundary Layer Interaction in a Supersonic Compression Ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

Effects of Magnetohydrodynamic Interaction-Zone Position on Shock-Wave/Boundary-Layer Interaction

A third-order weighted essentially nonoscillatory and non-free-parameter difference scheme magnetohydrodynamic solver has been established to investigate the mechanisms of magnetohydrodynamics controlling separation induced by an oblique shock wave impinging on a flat plate. The effects of magnetohydrodynamic interaction-zone location on the separation point, reattachment point, separation-bubble size, and boundary-layer velocity profiles are analyzed. The results show that there exists a best location for the magnetohydrodynamic zone to be applied, where the separation point is delayed the farthest, and the separation bubble is decreased up to about 50% in size compared to the case without magnetohydrodynamic control, which demonstrated the promising of magnetohydrodynamics suppressing the separation induced by shock-wave/boundary-layer interactions.

Boundary Characterization Experiment Series Overview

Ocean acoustic propagation and reverberation in continental shelf regions is often controlled by the seabed and sea surface boundaries. A series of three multi-national and multi-disciplinary experiments was conducted between 2000-2002 to identify and measure key ocean boundary characteristics. The frequency range of interest was nominally 500-5000 Hz with the main focus on the seabed, which is generally considered as the boundary of greatest importance and least understood. Two of the experiments were conducted in the Mediterranean in the Strait of Sicily and one experiment in the North Atlantic with sites on the outer New Jersey Shelf (STRATAFORM area) and on the Scotian Shelf. Measurements included seabed reflection, seabed, surface, and biologic scattering, propagation, reverberation, and ambient noise along with supporting oceanographic, geologic, and geophysical data. This paper is primarily intended to provide an overview of the experiments and the strategies that linked the various measurements together, with detailed experiment results contained in various papers in this volume and other sources

extending the boundary of spreadsheet programming: lessons learned from chinese governmental projects

ICSE

Effects of Pr dopant on grain boundary and electrical properties of Ce5.2Sm0.8MoO15-delta

Material formulated as Ce5.2Sm0.8-xPrxMo15-(delta) (x=0.08) was prepared by adding small amounts of Pr dopant in oxide Ce5.2SM0.8-xPrxMoO15-delta. Structural and electrical properties were investigated by means of X-ray diffraction (XRD), Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), field emission scanning electron microscopy (FE-SEM) and AC impedance spectroscopy. The effect of small amounts of Pr on microstructure and electrical conductivity was discussed. It was showed that the material doped with Pr has a lot of dents and small openings, which provide channels for oxygen ions, resulting in lower grain boundary and total conductivity activation energy. Thus the corresponding grain boundary conductivity and total conductivity of the material were improved notably. The grain boundary conductivity of the material doped with Pr is 6.79 X 10(-3) S center dot cm(-1) at 500 degrees C, which is twice as large as that without Pr (5.61 X 10(-5) S center dot cm(-1)).

Grain boundary conductivity of high purity neodymium-doped ceria nanosystem with and without the doping of molybdenum oxide

Solid solutions of Ce1-xNdxO2-x/2 (0.05 <= x <= 0.2) and (Ce1-xNdx)(0.95)MO0.05O2-delta (0.05 <= x <= 0.2) have been synthesized by a modified sol-gel method. Both materials have very low content of SiO2 (similar to 27 ppm). Their structures and ionic conductivities were characterized by X-ray diffraction (XRD), field-emission scanning electron microscopy (FESEM) and electrochemical impedance spectroscopy (M). The XRD patterns indicate that these materials are single phases with a cubic fluorite structure. The powders calcined at 300 degrees C with a crystal size of 5.7 nm have good sinterability, and the relative density could reach above 96% after being sintered at 1450 degrees C. With the addition Of MoO3, the sintering temperature could be decreased to 1250 degrees C. Impedance spectroscopy measurement in the temperature range of 250-800 degrees C indicates that a sharp increase of conductivity is observed when a small amount of Nd2O3 is added into ceria, of which Ce0.85Nd0.15O1.925 (15NDC) shows the highest conductivity. With the addition of a small amount Of MoO3, the grain boundary conductivity of 15NDC at 600 degrees C increases from 2.56 S m(-1) to 5.62 S m(-1).

Boundary effect of relief structure on crystallization of diblock copolymer in thin films

In this Letter, crystal growth of a symmetric crystalline-amorphous diblock copolymer, poly(styrene-b-epsilon-caprolactone) (PS-b-PCL), in thin films was investigated by atomic force microscopy (AFM), Relief structures of holes and islands were formed during annealing the film at the molten state, and the in situ observation of subsequent crystal growth at room temperature indicated that the crystals were preferred to occur at the edge of holes or islands and grew into the interior area. It was concluded that the stretched PCL blocks at the edge of relief structures, caused by material transportation or deformation of the interface, could act as nucleation agents during polymer crystallization. The crystal growth rate of individual lamellae varied both from lamellae to lamellae and in time, but the area occupied by crystals increased constantly with time. At 22 degreesC, the growth rate was 1.2 x 10(-2) mum(2)/min with the scan size 2 x 2 mum(2).

Growth rates of different polymorphs from interspherulitic boundary profiles

Quantitative data on the crystallization kinetics of polymorphic polymers can be derived from the investigation of gross spherulitic morphology formed in isothermal conditions. Depending on distance between centers, and the time lag between their formation and relative growth rates, various types of boundary lines can be generated by the impinging of two spherical bodies whose radii increase linearly with time, In polymorphic polymers, different types of spherulites often develop simultaneously at different rates from sporadic or predetermined nuclei. In same cases, the so-called growth transformation, in which a nucleus of the fast growing specie is formed at the tip of an advancing lamella of the slower crystal form, provides an alternative mode of nucleation, It is shown that if only one event of growth transformation takes place at the front of a slow growing body, the fast growing spherulite swallows the parent one and the resultant shape of interspherulitic boundary is described by two symmetrical logarithmic spirals whose parameters can be extracted from micrographs taken at the end of crystallization. These concepts are applied to determine the radial growth rate of gamma form spherulites of polypivalolactone in a wide range of temperatures through analysis of the alpha/gamma interspherulitic profiles formed in isothermal conditions and direct measurement of the growth rate of the alpha counterparts at the same temperature.

EPITAXIAL BEHAVIOR OF HDPE ON THE BOUNDARY OF HIGHLY ORIENTED IPP SUBSTRATES

The epitaxial crystallization behavior of high-density polyethylene on the boundary of highly oriented isotactic polypropylene (iPP) substrates has been investigated by means of atomic force microscopy (AFM) and transmission electron microscopy (TEM). The results obtained from AFM and TEM indicate that the epitaxial nucleation of HDPE on the highly oriented iPP substrates occurs earlier than that in the pure HDPE phase, i.e., homogeneous nucleation. Therefore the epitaxially grown HDPE lamellae can grow across the boundary of the iPP substrate into the HDPE spherulitic phase with the epitaxial orientation relationship remaining.

On the western boundary currents in the Philippine Sea

In this study we describe the velocity structure and transport of the North Equatorial Current (NEC), the Kuroshio, and the Mindanao Current (MC) using repeated hydrographic sections near the Philippine coast. A most striking feature of the current system in the region is the undercurrent structure below the surface flow. Both the Luzon Undercurrent and the Mindanao Undercurrent appear to be permanent phenomena. The present data set also provides an estimate of the mean circulation diagram (relative to 1500 dbar) that involves a NEC transport of 41 Sverdrups (Sv), a Kuroshio transport of 14 Sv, and a MC transport of 27 Sv, inducing a mass balance better than 1 Sv within the region enclosed by stations. The circulation diagram is insensitive to vertical displacements of the reference level within the depth range between 1500 and 2500 dbar. Transport fluctuations are, in general, consistent with earlier observations; that is, the NEC and the Kuroshio vary in the same phase with a seasonal signal superimposed with interannual variations, and the transport of the MC is dominated by a quasi-biennial oscillation. Dynamic height distributions are also examined to explore the dynamics of the current system.

Roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean

An ocean general circulation model (OGCM) is used to study the roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean. The western boundary reflection is defined as the total Kelvin waves leaving the western boundary, which include the reflection of the equatorial Rossby waves as well as the effects of alongshore winds, off-equatorial Rossby waves, and nonlinear processes near the western boundary. The evaluation of the reflection is based on a wave decomposition of the OGCM results and experiments with linear models. It is found that the alongshore winds along the east coast of Africa and the Rossby waves in the off-equatorial areas contribute significantly to the annual harmonics of the equatorial Kelvin waves at the western boundary. The semiannual harmonics of the Kelvin waves, on the other hand, originate primarily from a linear reflection of the equatorial Rossby waves. The dynamics of a dominant annual oscillation of sea level coexisting with the dominant semiannual oscillations of surface zonal currents in the central equatorial Indian Ocean are investigated. These sea level and zonal current patterns are found to be closely related to the linear reflections of the semiannual harmonics at the meridional boundaries. Because of the reflections, the second baroclinic mode resonates with the semiannual wind forcing; that is, the semiannual zonal currents carried by the reflected waves enhance the wind-forced currents at the central basin. Because of the different behavior of the zonal current and sea level during the reflections, the semiannual sea levels of the directly forced and reflected waves cancel each other significantly at the central basin. In the meantime, the annual harmonic of the sea level remains large, producing a dominant annual oscillation of sea level in the central equatorial Indian Ocean. The linear reflection causes the semiannual harmonics of the incoming and reflected sea levels to enhance each other at the meridional boundaries. In addition, the weak annual harmonics of sea level in the western basin, resulting from a combined effect of the western boundary reflection and the equatorial zonal wind forcing, facilitate the dominance by the semiannual harmonics near the western boundary despite the strong local wind forcing at the annual period. The Rossby waves are found to have a much larger contribution to the observed equatorial semiannual oscillations of surface zonal currents than the Kelvin waves. The westward progressive reversal of seasonal surface zonal currents along the equator in the observations is primarily due to the Rossby wave propagation.