On the forced matching numbers of bipartite graphs

Let G be a graph that admits a perfect matching. A forcing set for a perfect matching M of G is a subset S of M, such that S is contained in no other perfect matching of G. This notion has arisen in the study of finding resonance structures of a given molecule in chemistry. Similar concepts have been studied for block designs and graph colorings under the name defining set, and for Latin squares under the name critical set. There is some study of forcing sets of hexagonal systems in the context of chemistry, but only a few other classes of graphs have been considered. For the hypercubes Q(n), it turns out to be a very interesting notion which includes many challenging problems. In this paper we study the computational complexity of finding the forcing number of graphs, and we give some results on the possible values of forcing number for different matchings of the hypercube Q(n). Also we show an application to critical sets in back circulant Latin rectangles. (C) 2003 Elsevier B.V. All rights reserved.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (2nd report, direct numerical simulation using townsend's eddy forcing)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method using 'simple model eddies' (Townsend 1976) for DNS of stationary homogeneous isotropic turbulence is proposed. A force field is obtained in real space by sprinkling many space-filling 'simple model eddies' whose centers are randomly but uniformly distributed in space and whose axes of rotation are random. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects are investigated. The results show that stationary homogeneous isotropic turbulence is generated in real space using the present method. By using different model eddies with different sizes and rotation speeds, we could change the turbulence properties such as the integral and micro scales, the turbulent Reynolds number and the isotropy of turbulence. Turbulence intensities near the wall showed good agreements with the previous measurement and the linear analysis based on a rapid distortion theory (RDT). The splat effect (i.e., turbulence intensities of the components parallel to the boundary are amplified) occurs near the boundary and the viscous effect prohibits the splat effect at the quasi steady state at low Reynolds number.

Natural convection and heat transfer in attics subject to periodic thermal forcing

The heat transfer through the attics of buildings under realistic thermal forcing has been considered in this study. A periodic temperature boundary condition is applied on the sloping walls of the attic to show the basic flow features in the attic space over diurnal cycles. The numerical results reveal that, during the daytime heating stage, the flow in the attic space is stratified; whereas at the night-time cooling stage, the flow becomes unstable. A symmetrical solution is seen for relatively low Rayleigh numbers. However, as the Ra gradually increases, a transition occurs at a critical value of Ra. Above this critical value, an asymmetrical solution exhibiting a pitchfork bifurcation arises at the night-time. It is also found that the calculated heat transfer rate at the night-time cooling stage is much higher than that during the daytime heating stage.

Natural convection in a triangular enclosure subject to periodic thermal forcing

The effects of periodic thermal forcing on the flow field and heat transfer through an attic space are examined numerically in this paper. We consider the case with a fixed aspect ratio of 0.5 and a fixed Grashof number of 1.33×106. The numerical results reveal that, during the daytime, the flow is stratified; whereas at the night-time, the flow becomes unstable. A number of regular plumes and vortices are observed in the contours of isotherms and stream functions respectively. Moreover, the flow appears to be symmetric during the daytime, and becomes asymmetric at the night-time. It is also found that the flow is weaker during the daytime than that at the night-time in the present case, and the calculated heat transfer rate at the night-time is approximately three times greater than the heat transfer rate during the daytime.

Natural convection adjacent to an inclined flat plate and in an attic space under various thermal forcing conditions

The effect of viscous dissipation on natural convection from a vertical plate placed in a thermally stratified environment has been investigated numerically. The reduced equations are integrated by employing the implicit finite difference scheme or Ke1ler-box method and obtained the effect of heat due to viscous dissipation on the local skin-friction and loca1 Nusselt number at various stratification levels, for fluids having Prandtl number equals 10, 50, and 100. Solutions are also obtained using the perturbation technique for small values of viscous dissipation parameters and compared with the Finite Difference solutions. Effect of the heat transfer due to viscous dissipation and the temperature stratification are also shown on the velocity and temperature distributions in the boundary layer region. A numerical study of laminar doubly diffusive free convection flows adjacent to a vertical surface in a stable thermally stratified medium is also considered for this study. Solutions are obtained using the implicit Finite Difference method and compared with the local non-similarity method. The velocity and temperature distributions for different values of stratification parameter are shown graphically. The results show many interesting aspects of complex interaction of the two buoyant mechanisms.

A closed-form optimal tuning of mass dampers for one degree-of-freedom systems under rotating unbalance forcing

This paper is focused on the study of a vibrating system forced by a rotating unbalance and coupled to a tuned mass damper (TMD). The analysis of the dynamic response of the entire system is used to define the parameters of such device in order to achieve optimal damping properties. The inertial forcing due to the rotating unbalance depends quadratically on the forcing frequency and it leads to optimal tuning parameters that differ from classical values obtained for pure harmonic forcing. Analytical results demonstrate that frequency and damping ratios, as a function of the mass parameter, should be higher than classical optimal parameters. The analytical study is carried out for the undamped primary system, and numerically investigated for the damped primary system. We show that, for practical applications, proper TMD tuning allows to achieve a reduction in the steady-state response of about 20% with respect to the response achieved with a classically tuned damper. Copyright © 2015 by ASME.

Aerosol size distribution and its connection to cloud droplet number concentration

In order to predict the current state and future development of Earth s climate, detailed information on atmospheric aerosols and aerosol-cloud-interactions is required. Furthermore, these interactions need to be expressed in such a way that they can be represented in large-scale climate models. The largest uncertainties in the estimate of radiative forcing on the present day climate are related to the direct and indirect effects of aerosol. In this work aerosol properties were studied at Pallas and Utö in Finland, and at Mount Waliguan in Western China. Approximately two years of data from each site were analyzed. In addition to this, data from two intensive measurement campaigns at Pallas were used. The measurements at Mount Waliguan were the first long term aerosol particle number concentration and size distribution measurements conducted in this region. They revealed that the number concentration of aerosol particles at Mount Waliguan were much higher than those measured at similar altitudes in other parts of the world. The particles were concentrated in the Aitken size range indicating that they were produced within a couple of days prior to reaching the site, rather than being transported over thousands of kilometers. Aerosol partitioning between cloud droplets and cloud interstitial particles was studied at Pallas during the two measurement campaigns, First Pallas Cloud Experiment (First PaCE) and Second Pallas Cloud Experiment (Second PaCE). The method of using two differential mobility particle sizers (DMPS) to calculate the number concentration of activated particles was found to agree well with direct measurements of cloud droplet. Several parameters important in cloud droplet activation were found to depend strongly on the air mass history. The effects of these parameters partially cancelled out each other. Aerosol number-to-volume concentration ratio was studied at all three sites using data sets with long time-series. The ratio was found to vary more than in earlier studies, but less than either aerosol particle number concentration or volume concentration alone. Both air mass dependency and seasonal pattern were found at Pallas and Utö, but only seasonal pattern at Mount Waliguan. The number-to-volume concentration ratio was found to follow the seasonal temperature pattern well at all three sites. A new parameterization for partitioning between cloud droplets and cloud interstitial particles was developed. The parameterization uses aerosol particle number-to-volume concentration ratio and aerosol particle volume concentration as the only information on the aerosol number and size distribution. The new parameterization is computationally more efficient than the more detailed parameterizations currently in use, but the accuracy of the new parameterization was slightly lower. The new parameterization was also compared to directly observed cloud droplet number concentration data, and a good agreement was found.

Dust aerosols over India and adjacent continents retrieved using METEOSAT infrared radiance Part II: quantification of wind dependence and estimation of radiative forcing

Long-range transport of continental dust makes these particles a significant constituent even at locations far from their sources. It is important to study the temporal variations in dust loading over desert regions and the role of meteorology, in order to assess its radiative impact. In this paper, infrared radiance (10.5-12.5 mu m), acquired by the METEOSAT-5 satellite (similar to 5-km resolution) during 1999 and 2003 was used to quantify wind dependence of dust aerosols and to estimate the radiative forcing. Our analysis shows that the frequency of occurrence of dust events was higher during 2003 compared to 1999. Since the dust production function depends mainly on the surface wind speed over regions which are dry and without vegetation, the role of surface wind on IDDI was examined in detail. It was found that an increase of IDDI with wind speed was nearly linear and the rate of increase in IDDI with surface wind was higher during 2003 compared to 1999. It was also observed that over the Indian desert, when wind speed was the highest during monsoon months (June to August), the dust production rate was lower because of higher soil moisture (due to monsoon rainfall). Over the Arabian deserts, when the wind speed is the highest during June to August, the dust production rate is also highest, as soil moisture is lowest during this season. Even though nothing can be said precisely on the reason why 2003 had a greater number of dust events, examination of monthly mean soil moisture at source regions indicates that the occurrence of high winds simultaneous with high soil moisture could be the reason for the decreased dust production efficiency in 1999. It appears that the deserts of Northwest India are more efficient dust sources compared to the deserts of Saudi Arabia and Northeast Africa (excluding Sahara). The radiative impact of dust over various source regions is estimated, and the regionally and annually averaged top of the atmosphere dust radiative forcing (short wave, clear-sky and over land) over the entire study region (0-35 degrees N; 30 degrees-100 degrees E) was in the range of -0.9 to +4.5 W m(-2). The corresponding values at the surface were in the range of -10 to -25 W m(-2). Our studies demonstrate that neglecting the diurnal variation of dust can cause errors in the estimation of long wave dust forcing by as much as 50 to 100%, and nighttime retrieval of dust can significantly reduce the uncertainties. A method to retrieve dust aerosols during nighttime is proposed. The regionally and annually averaged long wave dust radiative forcing was +3.4 +/- 1.6 W m(-2).

Forcing of self-excited round jet diffusion flames

In this experimental and numerical study, two types of round jet are examined under acoustic forcing. The first is a non-reacting low density jet (density ratio 0.14). The second is a buoyant jet diffusion flame at a Reynolds number of 1100 (density ratio of unburnt fluids 0.5). Both jets have regions of strong absolute instability at their base and this causes them to exhibit strong self-excited bulging oscillations at welldefined natural frequencies. This study particularly focuses on the heat release of the jet diffusion flame, which oscillates at the same natural frequency as the bulging mode, due to the absolutely unstable shear layer just outside the flame. The jets are forced at several amplitudes around their natural frequencies. In the non-reacting jet, the frequency of the bulging oscillation locks into the forcing frequency relatively easily. In the jet diffusion flame, however, very large forcing amplitudes are required to make the heat release lock into the forcing frequency. Even at these high forcing amplitudes, the natural mode takes over again from the forced mode in the downstream region of the flow, where the perturbation is beginning to saturate non-linearly and where the heat release is high. This raises the possibility that, in a flame with large regions of absolute instability, the strong natural mode could saturate before the forced mode, weakening the coupling between heat release and incident pressure perturbations, hence weakening the feedback loop that causes combustion instability. © 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Advancements in jet turbulence and noise modeling: accurate one-way solutions and empirical evaluation of the nonlinear forcing of wavepackets

Jet noise reduction is an important goal within both commercial and military aviation. Although large-scale numerical simulations are now able to simultaneously compute turbulent jets and their radiated sound, lost-cost, physically-motivated models are needed to guide noise-reduction efforts. A particularly promising modeling approach centers around certain large-scale coherent structures, called wavepackets, that are observed in jets and their radiated sound. The typical approach to modeling wavepackets is to approximate them as linear modal solutions of the Euler or Navier-Stokes equations linearized about the long-time mean of the turbulent flow field. The near-field wavepackets obtained from these models show compelling agreement with those educed from experimental and simulation data for both subsonic and supersonic jets, but the acoustic radiation is severely under-predicted in the subsonic case. This thesis contributes to two aspects of these models. First, two new solution methods are developed that can be used to efficiently compute wavepackets and their acoustic radiation, reducing the computational cost of the model by more than an order of magnitude. The new techniques are spatial integration methods and constitute a well-posed, convergent alternative to the frequently used parabolized stability equations. Using concepts related to well-posed boundary conditions, the methods are formulated for general hyperbolic equations and thus have potential applications in many fields of physics and engineering. Second, the nonlinear and stochastic forcing of wavepackets is investigated with the goal of identifying and characterizing the missing dynamics responsible for the under-prediction of acoustic radiation by linear wavepacket models for subsonic jets. Specifically, we use ensembles of large-eddy-simulation flow and force data along with two data decomposition techniques to educe the actual nonlinear forcing experienced by wavepackets in a Mach 0.9 turbulent jet. Modes with high energy are extracted using proper orthogonal decomposition, while high gain modes are identified using a novel technique called empirical resolvent-mode decomposition. In contrast to the flow and acoustic fields, the forcing field is characterized by a lack of energetic coherent structures. Furthermore, the structures that do exist are largely uncorrelated with the acoustic field. Instead, the forces that most efficiently excite an acoustic response appear to take the form of random turbulent fluctuations, implying that direct feedback from nonlinear interactions amongst wavepackets is not an essential noise source mechanism. This suggests that the essential ingredients of sound generation in high Reynolds number jets are contained within the linearized Navier-Stokes operator rather than in the nonlinear forcing terms, a conclusion that has important implications for jet noise modeling.

Climate forcing of Pacific herring recruitment and growth in southern British Columbia

Over the last 50 years, much of the variability in ocean climate and herring recruitment has occurred at two dominant periods centered around 5 and 16 years. Herring growth has also exhibited a dominant 5- and 18-year periodicity. A recent analysis of a number of relevant time series suggests that interannual variations in oceanic conditions off the west coast of Vancouver Island affect survival of herring and their principal predator, Pacific hake, which also exhibits a marked 16-year oscillation in abundance. Thus the dynamics of the herring stock are modulated by a combination of climate and predator forcing. Much of the interannual variation in herring growth is centered around the 5-year (moderate ENSO period) and 16-year (strong ENSO period) ocean climate oscillations and the 16-year recruitment oscillation.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (3rd report, direct numerical simulation on wall blocking effects for anisotropic turbulence)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method proposed in the second report using Townsend's "simple model eddies" for DNS was extended to generate axisymmetric anisotropic turbulence. A force field is obtained in real space by sprinkling many space-filling "simple model eddies" whose centers are randomly but uniformly distributed in space. The axes of rotation are controlled in this study to generate axisymmetric anisotropic turbulence. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects for anisotropic turbulence are investigated. The results show that stationary axisymmetric anisotropic turbulence is generated using the present method. Turbulence intensities near the wall showed good agreements with the rapid distortion theory (RDT) for small t (t ≪ TL), where TL. is the eddy turnover time. The splat effect (i. e. turbulence intensities of the components parallel to the surface are amplified) occurs near the boundary and the viscous effect attenuates the splat effect at the quasi steady state at low Reynolds number as for Isotropic turbulence. Prandtl's secondary flow of the second kind does not occur for low Reynolds number flows, which qualitatively agrees with previous observetion in a mixing-box.

Updated correlation between aircraft smoke number and black carbon concentration

Aircraft emissions of black carbon (BC) contribute to anthropogenic climate forcing and degrade air quality. The smoke number (SN) is the current regulatory measure of aircraft particulate matter emissions and quantifies exhaust plume visibility. Several correlations between SN and the exhaust mass concentration of BC (CBC) have been developed, based on measurements relevant to older aircraft engines. These form the basis of the current standard method used to estimate aircraft BC emissions (First Order Approximation version 3 [FOA3]) for the purposes of environmental impact analyses. In this study, BC with a geometric mean diameter (GMD) of 20, 30, and 60 nm and filter diameters of 19 and 35 mm are used to investigate the effect of particle size and sampling variability on SN measurements. For BC with 20 and 30 nm GMD, corresponding to BC emitted by modern aircraft engines, a smaller SN results from a given CBC than is the case for BC with 60 nm GMD, which is more typical of older engines. An updated correlation between CBC and SNthat accounts for typical size of BC emitted by modern aircraft is proposed. An uncertainty of ±25% accounts for variation in GMD in the range 20-30 nm and for the range of filter diameters. The SN-CBC correlation currently used in FOA3 underestimates by a factor of 2.5-3 for SN <15, implying that current estimates of aircraft BC emissions derived from SN are underestimated by the same factor. Copyright © American Association for Aerosol Research.

Climate change and decadal to centennial-scale periodicities recorded in a late Holocene NE Pacific marine record: Examining the role of solar forcing

We present a decadal-scale late Holocene climate record based on diatoms, biogenic silica, and grain size from a 12-m sediment core (VEC02A04) obtained from Frederick Sound in the Seymour-Belize Inlet Complex of British Columbia, Canada. Sediments are characterized by graded, massive, and laminated intervals. Laminated intervals are most common between c. 2948–2708 cal. yr BP and c. 1992–1727 cal. yr BP. Increased preservation of laminated sediments and diatom assemblage changes at this time suggest that cli- mate became moderately drier and cooler relative to the preceding and succeeding intervals. Spectral and wavelet analyses are used to test for statistically significant periodicities in time series of proxies of primary production (total diatom abundance, biogenic silica) and hydrology (grain size) preserved in the Frederick Sound record. Periodicities of c. 42–53, 60–70, 82–89, 241–243, and 380 yrs are present. Results are com- pared to reconstructed sunspot number data of Solanki et al. (2004) using cross wavelet transform to evalu- ate the role of solar forcing on NE Pacific climate. Significant common power of periodicities between c. 42– 60, 70–89, 241–243, and of 380 yrs occur, suggesting that celestial forcing impacted late Holocene climate at Frederick Sound. Replication of the c. 241–243 yr periodicity in sunspot time series is most pronounced be- tween c. 2900 cal. yr BP and c. 2000 cal. yr BP, broadly correlative to the timing of maximum preservation of laminated sedimentary successions and diatom assemblage changes. High solar activity at the Suess/de Vries band may have been manifested as a prolonged westward shift and/or weakening of the Aleutian Low in the mid-late Holocene, which would have diverted fewer North Pacific storms and resulted in the relatively dry conditions reconstructed for the Seymour-Belize Inlet Complex.