Reply to comment by A. G. J. Hilberts and P. A. Troch on "Influence of capillarity on a simple harmonic oscillating water table: Sand column experiments and modeling"

We thank Hilberts and Troch [2006] for their comment on our paper [Cartwright et al, 2005]. Before proceeding with our specific replies to the comments we would first like to clarify the definitions and meanings of equations (1)-(3) as presented by Hilberts and Troch [2006]. First, equation (1) is the fundamental definition of the (complex) effective porosity as derived by Nielsen and Perrochet [2000]. Equations (2) and (3), however, represent the linear frequency response function of the water table in the sand column responding to simple harmonic forcing. This function, which was validated by Nielsen and Perrochet [2000], provides an alternative method for estimating the complex effective porosity from the experimental sand column data in the absence of direct measurements of h_(tot) (which are required if equation (1) is to be used).

MODELING OF TRANSFER PATH FOR DETERMINATION OF COMBUSTION AND NOISE METRICS ON DIESEL ENGINES

Determination of combustion metrics for a diesel engine has the potential of providing feedback for closed-loop combustion phasing control to meet current and upcoming emission and fuel consumption regulations. This thesis focused on the estimation of combustion metrics including start of combustion (SOC), crank angle location of 50% cumulative heat release (CA50), peak pressure crank angle location (PPCL), and peak pressure amplitude (PPA), peak apparent heat release rate crank angle location (PACL), mean absolute pressure error (MAPE), and peak apparent heat release rate amplitude (PAA). In-cylinder pressure has been used in the laboratory as the primary mechanism for characterization of combustion rates and more recently in-cylinder pressure has been used in series production vehicles for feedback control. However, the intrusive measurement with the in-cylinder pressure sensor is expensive and requires special mounting process and engine structure modification. As an alternative method, this work investigated block mounted accelerometers to estimate combustion metrics in a 9L I6 diesel engine. So the transfer path between the accelerometer signal and the in-cylinder pressure signal needs to be modeled. Depending on the transfer path, the in-cylinder pressure signal and the combustion metrics can be accurately estimated - recovered from accelerometer signals. The method and applicability for determining the transfer path is critical in utilizing an accelerometer(s) for feedback. Single-input single-output (SISO) frequency response function (FRF) is the most common transfer path model; however, it is shown here to have low robustness for varying engine operating conditions. This thesis examines mechanisms to improve the robustness of FRF for combustion metrics estimation. First, an adaptation process based on the particle swarm optimization algorithm was developed and added to the single-input single-output model. Second, a multiple-input single-output (MISO) FRF model coupled with principal component analysis and an offset compensation process was investigated and applied. Improvement of the FRF robustness was achieved based on these two approaches. Furthermore a neural network as a nonlinear model of the transfer path between the accelerometer signal and the apparent heat release rate was also investigated. Transfer path between the acoustical emissions and the in-cylinder pressure signal was also investigated in this dissertation on a high pressure common rail (HPCR) 1.9L TDI diesel engine. The acoustical emissions are an important factor in the powertrain development process. In this part of the research a transfer path was developed between the two and then used to predict the engine noise level with the measured in-cylinder pressure as the input. Three methods for transfer path modeling were applied and the method based on the cepstral smoothing technique led to the most accurate results with averaged estimation errors of 2 dBA and a root mean square error of 1.5dBA. Finally, a linear model for engine noise level estimation was proposed with the in-cylinder pressure signal and the engine speed as components.

Analisi vibroacustica di un sistema di scarico per monoposto Formula SAE

Lo scopo del presente lavoro è quello di determinare una metodologia di analisi vibroacustica di validità generale applicabile a tutti i casi nei quali la forzante sia di tipo vibrazionale. Nello specifico si è analizzato il comportamento strutturale e acustico di un impianto di scarico di una monoposto Formula SAE. Ricorrendo all’analisi FEM (Finite Element Method) è possibile determinare e quantificare gli effetti dannosi causati dalle vibrazioni già nella fase di prototipazione permettendo una sostanziale riduzione dei tempi e costi. La determinazione del comportamento strutturale del modello alle vibrazioni è iniziata dall’analisi modale, grazie alla quale sono state determinate le frequenze naturali e i modi propri dell’impianto di scarico. Successivamente, l’analisi FRF (Frequency Response Function) ha permesso di conoscere la risposta del nostro sistema ad una forzante imposta mettendo in luce le diverse criticità strutturali. Con il presupposto di ottenere delle condizioni di carico che fossero il più vicine possibili alle normali condizioni operative si è impostata un’analisi PSD (Power Spectral Density). Per concludere la prima parte dell’analisi si è reso necessario indagare anche il comportamento a fatica vibrazionale, valutando in questo modo le zone soggette a vita finita e quindi le prime a cedere in fase di esercizio. La parte finale è stata dedicata all’analisi del rumore. Dall’analisi FRF si è determinata la SPL (Sound Pressure Level) ottenendo come output un valore di pressione sonora prodotto dall’effetto della propagazione delle onde di pressione generate dalla vibrazione strutturale dell’impianto di scarico. Infine, l’analisi di Transmission Loss ha permesso di valutare l’efficacia della geometria del silenziatore sulla riduzione del livello acustico generato dal transito dei gas di scarico alle diverse frequenze.

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Frequency responses of hearing aids and their effects on the speech perception of eight hard-of-hearing subjects

This dissertation examines the relationship between frequency response and word-discrimination performance of hearing-impaired persons. Three questions are addressed: does the restoration of the normal field-to-eardrum transfer function improve word discrimination; is the restoration of the normal shape of the audibility curve (uniform hearing level at all frequencies) beneficial to hearing-impaired listeners; and can speech discrimination be improved by an extension of the present narrow-band response in hearing aids.

Piezoelectric and electroacoustic properties of Ti-doped AlN thin films as a function of Ti content

In this work we present the assessment of the structural and piezoelectric properties of Al(0.5-x)TixN0.5 compounds (titanium content menor que6% atomic), which are expected to possess improved properties than conventional AlN films, such as larger piezoelectric activity, thermal stability of frequency and temperature resistance. Al:Ti:N films were deposited from a twin concentric target of Al and Ti by reactive AC sputtering, which provided films with a radial gradient of the Ti concentration. The properties of the films were investigated as a function of their composition, which was measured by electron dispersive energy dispersive X-ray spectroscopy and Rutherford backscattering spectrometry. The microstructure and morphology of the films were assessed by X-ray diffraction and infrared reflectance. Their electroacoustic properties and dielectric constant were derived from the frequency response of BAW test resonators. Al:Ti:N films properties appear to be strongly dependent on the Ti content, which modifies the AlN wurtzite crystal structure leading to greater dielectric constant, lower sound velocities, lower electromechanical factor and moderately improved temperature coefficient of the resonant frequency.

Monte Carlo Simulation Of The Response Functions Of Cdte Detectors To Be Applied In X-ray Spectroscopy.

In this work, the energy response functions of a CdTe detector were obtained by Monte Carlo (MC) simulation in the energy range from 5 to 160keV, using the PENELOPE code. In the response calculations the carrier transport features and the detector resolution were included. The computed energy response function was validated through comparison with experimental results obtained with (241)Am and (152)Eu sources. In order to investigate the influence of the correction by the detector response at diagnostic energy range, x-ray spectra were measured using a CdTe detector (model XR-100T, Amptek), and then corrected by the energy response of the detector using the stripping procedure. Results showed that the CdTe exhibits good energy response at low energies (below 40keV), showing only small distortions on the measured spectra. For energies below about 80keV, the contribution of the escape of Cd- and Te-K x-rays produce significant distortions on the measured x-ray spectra. For higher energies, the most important correction is the detector efficiency and the carrier trapping effects. The results showed that, after correction by the energy response, the measured spectra are in good agreement with those provided by a theoretical model of the literature. Finally, our results showed that the detailed knowledge of the response function and a proper correction procedure are fundamental for achieving more accurate spectra from which quality parameters (i.e., half-value layer and homogeneity coefficient) can be determined.

A hybrid Particle Swarm Optimization - Simplex algorithm (PSOS) for structural damage identification

This study proposes a new PSOS-model based damage identification procedure using frequency domain data. The formulation of the objective function for the minimization problem is based on the Frequency Response Functions (FRFs) of the system. A novel strategy for the control of the Particle Swarm Optimization (PSO) parameters based on the Nelder-Mead algorithm (Simplex method) is presented; consequently, the convergence of the PSOS becomes independent of the heuristic constants and its stability and confidence are enhanced. The formulated hybrid method performs better in different benchmark functions than the Simulated Annealing (SA) and the basic PSO (PSO(b)). Two damage identification problems, taking into consideration the effects of noisy and incomplete data, were studied: first, a 10-bar truss and second, a cracked free-free beam, both modeled with finite elements. In these cases, the damage location and extent were successfully determined. Finally, a non-linear oscillator (Duffing oscillator) was identified by PSOS providing good results. (C) 2009 Elsevier Ltd. All rights reserved

Modeling and Analysis of Piezoelectric Energy Harvesting From Aeroelastic Vibrations Using the Doublet-Lattice Method

Multifunctional structures are pointed out as an important technology for the design of aircraft with volume, mass, and energy source limitations such as unmanned air vehicles (UAVs) and micro air vehicles (MAVs). In addition to its primary function of bearing aerodynamic loads, the wing/spar structure of an UAV or a MAV with embedded piezoceramics can provide an extra electrical energy source based on the concept of vibration energy harvesting to power small and wireless electronic components. Aeroelastic vibrations of a lifting surface can be converted into electricity using piezoelectric transduction. In this paper, frequency-domain piezoaeroelastic modeling and analysis of a canti-levered platelike wing with embedded piezoceramics is presented for energy harvesting. The electromechanical finite-element plate model is based on the thin-plate (Kirchhoff) assumptions while the unsteady aerodynamic model uses the doublet-lattice method. The electromechanical and aerodynamic models are combined to obtain the piezoaeroelastic equations, which are solved using a p-k scheme that accounts for the electromechanical coupling. The evolution of the aerodynamic damping and the frequency of each mode are obtained with changing airflow speed for a given electrical circuit. Expressions for piezoaeroelastically coupled frequency response functions (voltage, current, and electrical power as well the vibratory motion) are also defined by combining flow excitation with harmonic base excitation. Hence, piezoaeroelastic evolution can be investigated in frequency domain for different airflow speeds and electrical boundary conditions. [DOI:10.1115/1.4002785]

Bandwidth considerations for multilevel converters

Multilevel converters can achieve an overall effective switch frequency multiplication and consequent ripple reduction through the cancellation of the lowest order switch frequency terms. This paper investigates the harmonic content and the frequency response of these multimodulator converters. It is shown that the transfer function of uniformly sampled modulators is a bessel function associated with the inherent sampling process. Naturally sampled modulators have a flat transfer function, but multiple switchings per switch cycle will occur unless the input is slew-rate limited. Lower sideband harmonics of the effective carrier frequency and, in uniform converters, harmonics of the input signal also limit the useful bandwidth. Observations about the effect of the number of converters, their type (naturally or uniformly sampled), and the ratio of modulating frequency and switch frequency are made.

The preparation and readiness for voluntary movement: a high-field event-related fMRI study of the Bereitschafts-BOLD response

Activity within motor areas of the cortex begins to increase 1 to 2 s prior to voluntary self-initiated movement (termed the Bereitschaftspotential or readiness potential). There has been much speculation and debate over the precise source of this early premovement activity as it is important for understanding the roles of higher order motor areas in the preparation and readiness for voluntary movement. In this study, we use high-field (3-T) event-related fMRI with high temporal sampling (partial brain volumes every 250 ms) to specifically examine hemodynamic response time courses during the preparation, readiness, and execution of purely self-initiated voluntary movement. Five right-handed healthy volunteers performed a rapid sequential finger-to-thumb movement performed at self-determined times (12-15 trials). Functional images for each trial were temporally aligned and the averaged time series for each subject was iteratively correlated with a canonical hemodynamic response function progressively shifted in time. This analysis method identified areas of activation without constraining hemodynamic response timing. All subjects showed activation within frontal mesial areas, including supplementary motor area (SMA) and cingulate motor areas, as well as activation in left primary sensorimotor areas. The time courses of hemodynamic responses showed a great deal of variability in shape and timing between subjects; however, four subjects clearly showed earlier relative hemodynamic responses within SMA/cingulate motor areas compared with left primary motor areas. These results provide further evidence that the SMA and cingulate motor areas are major contributors to early stage premovement activity and play an important role in the preparation and readiness for voluntary movement. (C) 2003 Elsevier Inc. All rights reserved.

From EEG to BOLD: Brain mapping and estimating transfer functions in simultaneous EEG-fMRI acquisitions

Simultaneous acquisition of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) aims to disentangle the description of brain processes by exploiting the advantages of each technique. Most studies in this field focus on exploring the relationships between fMRI signals and the power spectrum at some specific frequency bands (alpha, beta, etc.). On the other hand, brain mapping of EEG signals (e.g., interictal spikes in epileptic patients) usually assumes an haemodynamic response function for a parametric analysis applying the GLM, as a rough approximation. The integration of the information provided by the high spatial resolution of MR images and the high temporal resolution of EEG may be improved by referencing them by transfer functions, which allows the identification of neural driven areas without strong assumptions about haemodynamic response shapes or brain haemodynamic`s homogeneity. The difference on sampling rate is the first obstacle for a full integration of EEG and fMRI information. Moreover, a parametric specification of a function representing the commonalities of both signals is not established. In this study, we introduce a new data-driven method for estimating the transfer function from EEG signal to fMRI signal at EEG sampling rate. This approach avoids EEG subsampling to fMRI time resolution and naturally provides a test for EEG predictive power over BOLD signal fluctuations, in a well-established statistical framework. We illustrate this concept in resting state (eyes closed) and visual simultaneous fMRI-EEG experiments. The results point out that it is possible to predict the BOLD fluctuations in occipital cortex by using EEG measurements. (C) 2010 Elsevier Inc. All rights reserved.

A hybrid MLS technique for room impulse response estimation

The measurement of room impulse response (RIR) when there are high background noise levels frequently means one must deal with very low signal-to-noise ratios (SNR). if such is the case, the measurement might yield unreliable results, even when synchronous averaging techniques are used. Furthermore, if there are non-linearities in the apparatus or system time variances, the final SNR can be severely degraded. The test signals used in RIR measurement are often disturbed by non-stationary ambient noise components. A novel approach based on the energy analysis of ambient noise - both in the time and in frequency - was considered. A modified maximum length sequence (MLS) measurement technique. referred to herein as the hybrid MLS technique, was developed for use in room acoustics. The technique consists of reducing the noise energy of the captured sequences before applying the averaging technique in order to improve the overall SNRs and frequency response accuracy. Experiments were conducted under real conditions with different types of underlying ambient noises. Results are shown and discussed. Advantages and disadvantages of the hybrid MLS technique over standard MLS technique are evaluated and discussed. Our findings show that the new technique leads to a significant increase in the overall SNR. (C) 2008 Elsevier Ltd. All rights reserved.

Fractional-order impulse response of the respiratory system

This paper presents the measurement, frequency-response modeling and identification, and the corresponding impulse time response of the human respiratory impedance and admittance. The investigated adult patient groups were healthy, diagnosed with chronic obstructive pulmonary disease and kyphoscoliosis, respectively. The investigated children patient groups were healthy, diagnosed with asthma and cystic fibrosis, respectively. Fractional order (FO) models are identified on the measured impedance to quantify the respiratory mechanical properties. Two methods are presented for obtaining and simulating the time-domain impulse response from FO models of the respiratory admittance: (i) the classical pole-zero interpolation proposed by Oustaloup in the early 90s, and (ii) the inverse discrete Fourier Transform (DFT). The results of the identified FO models for the respiratory admittance are presented by means of their average values for each group of patients. Consequently, the impulse time response calculated from the frequency response of the averaged FO models is given by means of the two methods mentioned above. Our results indicate that both methods provide similar impulse response data. However, we suggest that the inverse DFT is a more suitable alternative to the high order transfer functions obtained using the classical Oustaloup filter. Additionally, a power law model is fitted on the impulse response data, emphasizing the intrinsic fractal dynamics of the respiratory system.