Application of ultraspherical polynomials to forced oscillations of a third-order non-linear system

In this paper the method of ultraspherical polynomial approximation is applied to study the steady-state response in forced oscillations of a third-order non-linear system. The non-linear function is expanded in ultraspherical polynomials and the expansion is restricted to the linear term. The equation for the response curve is obtained by using the linearized equation and the results are presented graphically. The agreement between the approximate solution and the analog computer solution is satisfactory. The problem of stability is not dealt with in this paper.

Flow field and traverse times for fan forced injection of fumigant via circular or annular inlet into stored grain

Fan forced injection of phosphine gas fumigant into stored grain is a common method to treat infestation by insects. For low injection velocities the transport of fumigant can be modelled as Darcy flow in a porous medium where the gas pressure satisfies Laplace's equation. Using this approach, a closed form series solution is derived for the pressure, velocity and streamlines in a cylindrically stored grain bed with either a circular or annular inlet, from which traverse times are numerically computed. A leading order closed form expression for the traverse time is also obtained and found to be reasonable for inlet configurations close to the central axis of the grain storage. Results are interpreted for the case of a representative 6m high farm wheat store, where the time to advect the phosphine to almost the entire grain bed is found to be approximately one hour.

Automatic Path Planning and Control Design for Autonomous Landing of UAVs using Dynamic Inversion

In this paper a nonlinear control has been designed using the dynamic inversion approach for automatic landing of unmanned aerial vehicles (UAVs), along with associated path planning. This is a difficult problem because of light weight of UAVs and strong coupling between longitudinal and lateral modes. The landing maneuver of the UAV is divided into approach, glideslope and flare. In the approach UAV aligns with the centerline of the runway by heading angle correction. In glideslope and flare the UAV follows straight line and exponential curves respectively in the pitch plane with no lateral deviations. The glideslope and flare path are scheduled as a function of approach distance from runway. The trajectory parameters are calculated such that the sink rate at touchdown remains within specified bounds. It is also ensured that the transition from the glideslope to flare path is smooth by ensuring C-1 continuity at the transition. In the outer loop, the roll rate command is generated by assuring a coordinated turn in the alignment segment and by assuring zero bank angle in the glideslope and flare segments. The pitch rate command is generated from the error in altitude to control the deviations from the landing trajectory. The yaw rate command is generated from the required heading correction. In the inner loop, the aileron, elevator and rudder deflections are computed together to track the required body rate commands. Moreover, it is also ensured that the forward velocity of the UAV at the touch down remains close to a desired value by manipulating the thrust of the vehicle. A nonlinear six-DOF model, which has been developed from extensive wind-tunnel testing, is used both for control design as well as to validate it.

Coupled Solution for Forced Recharge in Confined Aquifers

Analytical solutions for forced well recharge currently in use were initially developed for pumping scenarios and applied for recharge cases assuming that radial flow in the recharge well replicates a mirror image of that in to a pumping well. Moreover these solutions were not extended to multiaquifer systems. Well bore numerical solutions were generally not considering the effect of well bore interaction, which has a significant effect in the case of a recharge well. In the present paper, improved analytical solutions are developed for a well fully penetrating either single or multiaquifers in respect.to of well storage, well loss, and interactions between the individual aquifers through well bore. The solution developed for constant and variable rates of injection and well loss is applied to the experimental data of the Hansol well injection project near the city of Ahmedabad in the Gujarat state in India. The paper also discusses the difference in well hydraulics of recharge and recovery wells.

Spirometric studies on the adult general population of Helsinki : bronchodilation responses, determinants, and intrasession repeatability of FEV1, FEV6, FVC, and forced expiratory time

Spirometry is the most widely used lung function test in the world. It is fundamental in diagnostic and functional evaluation of various pulmonary diseases. In the studies described in this thesis, the spirometric assessment of reversibility of bronchial obstruction, its determinants, and variation features are described in a general population sample from Helsinki, Finland. This study is a part of the FinEsS study, which is a collaborative study of clinical epidemiology of respiratory health between Finland (Fin), Estonia (Es), and Sweden (S). Asthma and chronic obstructive pulmonary disease (COPD) constitute the two major obstructive airways diseases. The prevalence of asthma has increased, with around 6% of the population in Helsinki reporting physician-diagnosed asthma. The main cause of COPD is smoking with changes in smoking habits in the population affecting its prevalence with a delay. Whereas airway obstruction in asthma is by definition reversible, COPD is characterized by fixed obstruction. Cough and sputum production, the first symptoms of COPD, are often misinterpreted for smokers cough and not recognized as first signs of a chronic illness. Therefore COPD is widely underdiagnosed. More extensive use of spirometry in primary care is advocated to focus smoking cessation interventions on populations at risk. The use of forced expiratory volume in six seconds (FEV6) instead of forced vital capacity (FVC) has been suggested to enable office spirometry to be used in earlier detection of airflow limitation. Despite being a widely accepted standard method of assessment of lung function, the methodology and interpretation of spirometry are constantly developing. In 2005, the ATS/ERS Task Force issued a joint statement which endorsed the 12% and 200 ml thresholds for significant change in forced expiratory volume in one second (FEV1) or FVC during bronchodilation testing, but included the notion that in cases where only FVC improves it should be verified that this is not caused by a longer exhalation time in post-bronchodilator spirometry. This elicited new interest in the assessment of forced expiratory time (FET), a spirometric variable not usually reported or used in assessment. In this population sample, we examined FET and found it to be on average 10.7 (SD 4.3) s and to increase with ageing and airflow limitation in spirometry. The intrasession repeatability of FET was the poorest of the spirometric variables assessed. Based on the intrasession repeatability, a limit for significant change of 3 s was suggested for FET during bronchodilation testing. FEV6 was found to perform equally well as FVC in the population and in a subgroup of subjects with airways obstruction. In the bronchodilation test, decreases were frequently observed in FEV1 and particularly in FVC. The limit of significant increase based on the 95th percentile of the population sample was 9% for FEV1 and 6% for FEV6 and FVC; these are slightly lower than the current limits for single bronchodilation tests (ATS/ERS guidelines). FEV6 was proven as a valid alternative to FVC also in the bronchodilation test and would remove the need to control duration of exhalation during the spirometric bronchodilation test.

Turbulence-induced melting of a nonequilibrium vortex crystal in a forced thin fluid film

To gain a better understanding of recent experiments on the turbulence-induced melting of a periodic array of vortices in a thin fluid film, we perform a direct numerical simulation of the two-dimensional Navier-Stokes equations forced such that, at low Reynolds numbers, the steady state of the film is a square lattice of vortices. We find that as we increase the Reynolds number, this lattice undergoes a series of nonequilibrium phase transitions, first to a crystal with a different reciprocal lattice and then to a sequence of crystals that oscillate in time. Initially, the temporal oscillations are periodic; this periodic behaviour becoming more and more complicated with increasing Reynolds number until the film enters a spatially disordered nonequilibrium statistical steady state that is turbulent. We study this sequence of transitions using fluid-dynamics measures, such as the Okubo-Weiss parameter that distinguishes between vortical and extensional regions in the flow, ideas from nonlinear dynamics, e.g. Poincare maps, and theoretical methods that have been developed to study the melting of an equilibrium crystal or the freezing of a liquid and that lead to a natural set of order parameters for the crystalline phases and spatial autocorrelation functions that characterize short- and long-range order in the turbulent and crystalline phases, respectively.

On the ignition and extinction problems in forced convection systems

In this paper the problem of ignition and extinction has been formulated for the flow of a compressible fluid with Prandtl and Schmidt numbers taken as unity. In particular, the problems of (i) a jet impinging on a wall of combustible material and (ii) the opposed jet diffusion flame have been studied. In the wall jet case, three approximations in the momentum equation namely, (i) potential flow, (ii) viscous flow, (ii) viscous incompressible with k = 1 and (iii) Lees' approximation (taking pressure gradient terms zero) are studied. It is shown that the predictions of the mass flow rates at extinction are not very sensitive to the approximations made in the momentum equation. The effects of varying the wall temperature in the case (i) and the jet temperature in the case (ii) on the extinction speeds have been studied. The effects of varying the activation energy and the free stream oxidant concentration in case (ii), have also been investigated.

Unsteady MHD forced flow due to a point sink

An analysis is performed to study the unsteady laminar incompressible boundary-layer flow of an electrically conducting fluid in a cone due to a point sink with an applied magnetic field. The unsteadiness in the flow is considered for two types of motion, viz. the motion arising due to the free stream velocity varying continuously with time and the transient motion occurring due to an impulsive change either in the strength of the point sink or in the wall temperature. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The magnetic field increases the skin friction but reduces heat transfer. The heat transfer and temperature field are strongly influenced by the viscous dissipation and Prandtl number. The velocity field is more affected at the early stage of the transient motion, caused by an impulsive change in the strength of the point sink, as compared to the temperature field. When the transient motion is caused by a sudden change in the wall temperature, both skin friction and heat transfer take more time to reach a new steady state. The transient nature of the flow and heat transfer is active for a short time in the case of suction and for a long time in the case of injection. The viscous dissipation prolongs the transient behavior of the flow.

Combined free and forced convection along a rotating vertical cylinder

The steady incompressible laminar mixed convection boundary layer flow along a rotating slender vertical cylinder with an isothermal wall has been studied. The transformed coupled nonlinear partial differential equations have been solved numerically using the Keller box method. In general, the rotation of the cylinder, the buoyancy forces and the curvature parameter are found to significantly affect the skin friction, heat transfer, velocity and temperature profiles as well as the pressure distribution. The buoyancy forces cause an overshoot in the axial velocity profile but the rotation and curvature parameters reduce it.

Wind loading over the top portions of tall stacks with and without external landing platforms

The variation of the drag force near the top portions of tall stacks with and without external landing platforms, and with the exit open and closed, has been examined by model studies in a wind tunnel at Reynolds numbers of about 10(5). Pressure measurements on three models of different height to diameter ratios have been supplemented by flow visualisation studies. Observations confirm that when there is no platform, significant load enhancement over the top three to four diameters occurs, due to the high suction caused by the sharp separation of the flow over the top from the rim, in the aft regions of the stack. The enhanced loading is found to be greater if the exit is closed. A platform at the top, of less than twice the exit diameter, further increases the drag force near the top, but a still larger platform at the top, of about three times the exit diameter, decreases the drag force to values less than those much further below, effectively nullifying the enhanced drag force. It was found that such a reduction of the enhanced drag force in the top regions can also be achieved by a smaller platform of 1.1 to 1.3 times the local diameter, located at about three to five diameters below the top.

Simulation of ENSO-Related Surface Winds in the Tropical Pacific by an Atmospheric General Circulation Model Forced by Observed Sea Surface Temperatures

The authors present the simulation of the tropical Pacific surface wind variability by a low-resolution (R15 horizontal resolution and 18 vertical levels) version of the Center for Ocean-Land-Atmosphere Interactions, Maryland, general circulation model (GCM) when forced by observed global sea surface temperature. The authors have examined the monthly mean surface winds acid precipitation simulated by the model that was integrated from January 1979 to March 1992. Analyses of the climatological annual cycle and interannual variability over the Pacific are presented. The annual means of the simulated zonal and meridional winds agree well with observations. The only appreciable difference is in the region of strong trade winds where the simulated zonal winds are about 15%-20% weaker than observed, The amplitude of the annual harmonics are weaker than observed over the intertropical convergence zone and the South Pacific convergence zone regions. The amplitudes of the interannual variation of the simulated zonal and meridional winds are close to those of the observed variation. The first few dominant empirical orthogonal functions (EOF) of the simulated, as well as the observed, monthly mean winds are found to contain a targe amount of high-frequency intraseasonal variations, While the statistical properties of the high-frequency modes, such as their amplitude and geographical locations, agree with observations, their detailed time evolution does not. When the data are subjected to a 5-month running-mean filter, the first two dominant EOFs of the simulated winds representing the low-frequency EI Nino-Southern Oscillation fluctuations compare quite well with observations. However, the location of the center of the westerly anomalies associated with the warm episodes is simulated about 15 degrees west of the observed locations. The model simulates well the progress of the westerly anomalies toward the eastern Pacific during the evolution of a warm event. The simulated equatorial wind anomalies are comparable in magnitude to the observed anomalies. An intercomparison of the simulation of the interannual variability by a few other GCMs with comparable resolution is also presented. The success in simulation of the large-scale low-frequency part of the tropical surface winds by the atmospheric GCM seems to be related to the model's ability to simulate the large-scale low-frequency part of the precipitation. Good correspondence between the simulated precipitation and the highly reflective cloud anomalies is seen in the first two EOFs of the 5-month running means. Moreover, the strong correlation found between the simulated precipitation and the simulated winds in the first two principal components indicates the primary role of model precipitation in driving the surface winds. The surface winds simulated by a linear model forced by the GCM-simulated precipitation show good resemblance to the GCM-simulated winds in the equatorial region. This result supports the recent findings that the large-scale part of the tropical surface winds is primarily linear.

Alfven resonances and forced reconnection

The correspondence between the forced magnetic reconnection induced by perturbing the boundary of the simple Taylor model and the surface-wave-induced magnetic reconnection given by Alfven resonance theory is pointed out explicitly by showing that the theory of forced magnetic reconnection is actually embedded in the Alfven resonance theory. The advantages of viewing the forced reconnection as surface-wave-induced reconnection are briefly discussed in the context of the formation of small-scale structures at the magnetospheric boundary and solar coronal heating.