A novel method to produce small droplets from large nozzles.

This work presents a new method to generate droplets with diameters significantly smaller than the nozzle from which they emerge. The electrical waveform used to produce the jetting consists of a single square negative pulse. The negative edge of the pressure wave pulls the meniscus in, overturning the surface in such a way that a cavity is created. This cavity is then forced to collapse under the action of the positive edge of the pressure wave. This violent collapse produces a thin jet that eventually breaks up and produces droplets. Four droplet generator prototypes that demonstrate the capabilities of this novel mechanism are described. It is also shown that the proposed mechanism extends the existing limits of the commonly accepted inkjet operating regime.

An analysis of the changes to the mechanical and acoustical properties of paper caused by the damage mechanisms present in the fibre structure

This paper will analyse two of the likely damage mechanisms present in a paper fibre matrix when placed under controlled stress conditions: fibre/fibre bond failure and fibre failure. The failure process associated with each damage mechanism will be presented in detail focusing on the change in mechanical and acoustic properties of the surrounding fibre structure before and after failure. To present this complex process mathematically, geometrically simple fibre arrangements will be chosen based on certain assumptions regarding the structure and strength of paper, to model the damage mechanisms. The fibre structures are then formulated in terms of a hybrid vibro-acoustic model based on a coupled mass/spring system and the pressure wave equation. The model will be presented in detail in the paper. The simulation of the simple fibre structures serves two purposes; it highlights the physical and acoustic differences of each damage mechanism before and after failure, and also shows the differences in the two damage mechanisms when compared with one another. The results of the simulations are given in the form of pressure wave contours, time-frequency graphs and the Continuous Wavelet Transform (CWT) diagrams. The analysis of the results leads to criteria by which the two damage mechanisms can be identified. Using these criteria it was possible to verify the results of the simulations against experimental acoustic data. The models developed in this study are of specific practical interest in the paper-making industry, where acoustic sensors may be used to monitor continuous paper production. The same techniques may be adopted more generally to correlate acoustic signals to damage mechanisms in other fibre-based structures.

Simulation of the Unsteady Gas Flow through a Three-way Automotive Catalyst

This paper describes a model of a 1.8-litre four-cylinder four-stroke gasoline engine fitted with a close-coupled three-way catalyst (TWC). Designed to meet EURO 3 emissions standards, the engine includes some advanced emission control features in addition to the TWC, namely: variable valve timing (VVT), swirl control plates, and exhaust gas recirculation (EGR). Gas flow is treated as one-dimensional (1D) and unsteady in the engine ducting and in the catalyst. Reflection and transmission of pressure waves at the boundaries of the catalyst monolith are modelled. In-cylinder combustion is represented by a two-zone burn model with dissociation and reaction kinetics. A single Wiebe analysis of measured in-cylinder pressure data is used to determine the mass fraction burned as a function of crank angle (CA) at each engine speed. Measured data from steady-state dynamometer tests are presented for operation at wide open throttle (WOT) over a range of engine speeds. These results include CA-resolved traces of pressure at various locations throughout the engine together with cycle-averaged traces of gas composition entering the catalyst as indicated by a fast-response emissions analyser. Simulated engine performance and pressure wave action throughout the engine are well validated by the measured data.

Modelling of Unsteady Gas-Dynamic Flow in a Pipe and at its Exit Using CFD

Traditionally the simulation of the thermodynamic aspects of the internal combustion engine has been undertaken using one-dimensional gas-dynamic models to represent the intake and exhaust systems. CFD analysis of engines has been restricted to modelling of in-cylinder flow structures. With the increasing accessibility of CFD software it is now worth considering its use for complete gas-dynamic engine simulation. This paper appraises the accuracy of various CFD models in comparison to a 1D gas-dynamic simulation. All of the models are compared to experimental data acquired on an apparatus that generates a single gas-dynamic pressure wave. The progress of the wave along a constant area pipe and its subsequent reflection from the open pipe end are recorded with a number of high speed pressure transducers. It was found that there was little to choose between the accuracy of the 1D model and the best CFD model. The CFD model did not require experimentally derived loss coefficients to accurately represent the open pipe end; however, it took several hundred times longer to complete its analysis. The best congruency between the CFD models and the experimental data was achieved using the RNG k-e turbulence model. The open end of the pipe was most effectively represented by surrounding it with a relatively small volume of cells connected to the rest of the environment using a pressure boundary.

Modeling Frequency-Dependent Boundaries as Digital Impedance Filters in FDTD and K-DWM Room Acoustics Simulations

A new method for modeling-frequency-dependent boundaries in finite-difference time-domain (FDTD) and Kirchhoff variable digital waveguide mesh (K-DWM) room acoustics simulations is presented. The proposed approach allows the direct incorporation of a digital impedance filter (DIF) in the Multidimensional (2D or 3D) FDTD boundary model of a locally reacting surface. An explicit boundary update equation is obtained by carefully constructing a Suitable recursive formulation. The method is analyzed in terms of pressure wave reflectance for different wall impedance filters and angles of incidence. Results obtained from numerical experiments confirm the high accuracy of the proposed digital impedance filter boundary model, the reflectance of which matches locally reacting surface (LRS) theory closely. Furthermore a numerical boundary analysis (NBA) formula is provided as a technique for an analytic evaluation of the numerical reflectance of the proposed digital impedance filter boundary formulation.

Formulation of Locally Reacting Surfaces in FDTD/K-DWM Modelling of Acoustic Spaces

In this paper, we present new methods for constructing and analysing formulations of locally reacting surfaces that can be used in finite difference time domain (FDTD) simulations of acoustic spaces. Novel FDTD formulations of frequency-independent and simple frequency-dependent impedance boundaries are proposed for 2D and 3D acoustic systems, including a full treatment of corners and boundary edges. The proposed boundary formulations are designed for virtual acoustics applications using the standard leapfrog scheme based on a rectilinear grid, and apply to FDTD as well as Kirchhoff variable digital waveguide mesh (K-DWM) methods. In addition, new analytic evaluation methods that accurately predict the reflectance of numerical boundary formulations are proposed. numerical experiments and numerical boundary analysis (NBA) are analysed in time and frequency domains in terms of the pressure wave reflectance for different angles of incidence and various impedances. The results show that the proposed boundary formulations structurally adhere well to the theoretical reflectance. In particular, both reflectance magnitude and phase are closely approximated even at high angles of incidence and low impedances. Furthermore, excellent agreement was found between the numerical boundary analysis and the experimental results, validating both as tools for researching FDTD boundary formulations.

Pre- to Post-discharge Behavior in Saline Solution

The development of a plasma discharge at low voltage (200-600 V) in saline solution is characterized using fast and standard CCD camera imaging. Vapor formation, plasma formation, and vapor collapse and subsequent pressure wave propagation are observed. If, with increasing voltage, the total energy deposited is kept approximately constant, the sequence and nature of events are similar but develop faster and more reproducibly at the higher voltages. This is attributed to the slower temporal evolution of the vapor layer at lower voltages which means a greater sensitivity to hydrodynamic instabilities at the vapor-liquid interface.

Oscillatory Behavior of Supersonic Impinging Jet Flows

Oscillatory flows of a choked underexpanded supersonic impinging jet issuing from a convergent nozzle have been computed using the axisymmetric unsteady Navier–Stokes system. This paper focuses on the oscillatory flow features associated with the variation of the nozzle-to-plate distance and nozzle pressure ratio. Frequencies of the surface pressure oscillation and flow structural changes from computational results have been analyzed. Staging behaviour of the oscillation frequency has been observed for both cases of nozzle-to-plate distance variation and pressure ratio variation. However, the staging behavior for each case exhibits different features. These two distinct staging behaviors of the oscillation frequency are found to correlate well if the frequency and the distance are normalized by the length of the shock cell. It is further found that the staging behaviour is strongly correlated with the change of the pressure wave pattern in the jet shear layer, but not with the shock cell structure.

Étude par dynamique moléculaire de l'ablation par impulsions laser ultrabrèves de cibles nanocristallines

L’ablation de cibles d’Al nanocristallines (taille moyenne des cristallites d = 3,1 et 6,2 nm) par impulsions laser ultrabrèves (200 fs) a été étudiée par l’entremise de si- mulations combinant la dynamique moléculaire et le modèle à deux températures (two- temperature model, TTM) pour des fluences absorbées allant de 100 à 1300 J/m2. Nos simulations emploient un potentiel d’interaction de type EAM et les propriétés électro- niques des cibles en lien avec le TTM sont représentées par un modèle réaliste possédant une forme distincte dans le solide monocristallin, le solide nanocristallin et le liquide. Nous avons considéré l’effet de la taille moyenne des cristallites de même que celui de la porosité et nous avons procédé à une comparaison directe avec des cibles mono- cristallines. Nous avons pu montrer que le seuil d’ablation des métaux nanocristallins est significativement plus bas, se situant à 400 J/m2 plutôt qu’à 600 J/m2 dans le cas des cibles monocristallines, l’écart étant principalement dû à l’onde mécanique plus im- portante présente lors de l’ablation. Leur seuil de spallation de la face arrière est aussi significativement plus bas de par la résistance à la tension plus faible (5,40 GPa contre 7,24 GPa) des cibles nanocristallines. Il est aussi apparu que les contraintes résiduelles accompagnant généralement l’ablation laser sont absentes lors de l’ablation de cibles d’aluminium nanocristallines puisque la croissance cristalline leur permet d’abaisser leur volume spécifique. Nos résultats indiquent aussi que le seuil de fusion des cibles nano- cristallines est réduit de façon marquée dans ces cibles ce qui s’explique par la plus faible énergie de cohésion inhérente à ces matériaux. Nos simulations permettent de montrer que les propriétés structurelles et électroniques propres aux métaux nanocristallins ont toutes deux un impact important sur l’ablation.

Ablation laser femtoseconde de verres métalliques de Cu_x Zr_(1−x) : une étude par dynamique moléculaire

L’ablation laser de verres métalliques de CuxZr1−x (x = 0.33, 0.50 et 0.67) et d’un alliage métallique cristallin de CuZr2 dans la structure C11b a été étudiée par dynamique moléculaire (DM) combinée à un modèle à deux températures (TTM). Le seuil d’ablation (Fth) a été déterminé pour chacun des 4 échantillons et s'est avéré plus bas pour les échantillons plus riches en Cu étant donné que la cohésion du Cu est plus faible que celle du Zr dans tous les échantillons. Pour x=0.33, Fth est plus bas pour le cristal que pour l’amorphe car le couplage électron-phonon est plus faible dans ce dernier, ce qui implique que l’énergie est transférée plus lentement du système électronique vers le système ionique pour le a-CuZr2 que le c-CuZr2. La vitesse de l’onde de pression créée par l’impact du laser croît avec la fluence dans l’échantillon cristallin, contrairement aux échantillons amorphes dans lesquels sa vitesse moyenne est relativement constante avec la fluence. Ceci est expliqué par le fait que le module de cisaillement croît avec la pression pour le cristal, ce qui n’est pas le cas pour les verres métalliques étudiés. Finalement, la zone affectée par la chaleur (HAZ) a été étudiée via la profondeur de fusion et les déformations plastiques. La plus faible température de fusion des échantillons amorphes implique que la profondeur de fusion est plus importante dans ceux-ci que dans l’échantillon cristallin. Dans les verres métalliques, les déformations plastiques ont été identifiées sous forme de zones de transformation par cisaillement (STZ) qui diffusent et fusionnent à plus haute fluence. Aucune déformation plastique importante n’a été identifiée dans le c-CuZr2 mis à part de légères déformations près du front de fusion causées par les contraintes résiduelles. Ce travail a ainsi permis d’améliorer notre compréhension de l’ablation laser sur les verres métalliques et de l’étendue des dommages qu’elle peut entraîner.

Some numerical experiments on the effect of the variation of static stability in two-layer quasi-geostrophic models

A method to solve a quasi-geostrophic two-layer model including the variation of static stability is presented. The divergent part of the wind is incorporated by means of an iterative procedure. The procedure is rather fast and the time of computation is only 60–70% longer than for the usual two-layer model. The method of solution is justified by the conservation of the difference between the gross static stability and the kinetic energy. To eliminate the side-boundary conditions the experiments have been performed on a zonal channel model. The investigation falls mainly into three parts: The first part (section 5) contains a discussion of the significance of some physically inconsistent approximations. It is shown that physical inconsistencies are rather serious and for these inconsistent models which were studied the total kinetic energy increased faster than the gross static stability. In the next part (section 6) we are studying the effect of a Jacobian difference operator which conserves the total kinetic energy. The use of this operator in two-layer models will give a slight improvement but probably does not have any practical use in short periodic forecasts. It is also shown that the energy-conservative operator will change the wave-speed in an erroneous way if the wave-number or the grid-length is large in the meridional direction. In the final part (section 7) we investigate the behaviour of baroclinic waves for some different initial states and for two energy-consistent models, one with constant and one with variable static stability. According to the linear theory the waves adjust rather rapidly in such a way that the temperature wave will lag behind the pressure wave independent of the initial configuration. Thus, both models give rise to a baroclinic development even if the initial state is quasi-barotropic. The effect of the variation of static stability is very small, qualitative differences in the development are only observed during the first 12 hours. For an amplifying wave we will get a stabilization over the troughs and an instabilization over the ridges.

Photoacoustic spectroscopy of aromatic amino acids in proteins

This paper concerns the use of photoacoustic spectroscopy (PAS) to study the presence of aromatic amino acid in proteins. We examined the aromatic amino acids in six proteins with well-known structures using absorption spectra of near ultraviolet PAS over the wavelength range 240-320 nm. The fundamental understanding of the physical and chemical properties that govern the absorption of light and a subsequent release of heat to generate a transient pressure wave was used to test the concept of monitoring aromatic amino acids with this method. Second derivative spectroscopy in the ultraviolet region of proteins was also used to study the regions surrounding the aromatics and the percentage area in each band was related in order to determine the contribution in function of the respective molar extinction coefficients for each residue. Further investigation was conducted into the interaction between sodium dodecyl sulphate (SDS) and bothropstoxin-I (BthTx-I), with the purpose of identifying the aromatics that participate in the interaction. The clear changes in the second derivative and curve-fitting procedures suggest that initial SDS binding to the tryptophan located in the dimer interface and above 10 SDS an increased intensity between 260 and 320 nm, demonstrating that the more widespread tyrosine and phenylalanine residues contribute to the SDS/BthTx-I interactions. These results demonstrate the potential of near UV-PAS for the investigation of membrane proteins/detergent complexes in which light scattering is significant.

Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

Studio del metodo "systolic volume balance" per la stima non invasiva della portata media cardiaca

La portata media cardiaca, (cardiac output “CO”) è un parametro essenziale per una buona gestione dei pazienti o per il monitoraggio degli stessi durante la loro permanenza nell’unità di terapia intensiva. La stesura di questo elaborato prende spunto sull’articolo di Theodore G. Papaioannou, Orestis Vardoulis, and Nikos Stergiopulos dal titolo “ The “systolic volume balance” method for the non invasive estimation of cardiac output based on pressure wave analysis” pubblicato sulla rivista American Journal of Physiology-Heart and Circulatory Physiology nel Marzo 2012. Nel sopracitato articolo si propone un metodo per il monitoraggio potenzialmente non invasivo della portata media cardiaca, basato su principi fisici ed emodinamici, che usa l’analisi della forma d’onda di pressione e un metodo non invasivo di calibrazione e trova la sua espressione ultima nell’equazione Qsvb=(C*PPao)/(T-(Psm,aorta*ts)/Pm). Questa formula è stata validata dagli autori, con buoni risultati, solo su un modello distribuito della circolazione sistemica e non è ancora stato validato in vivo. Questo elaborato si pone come obiettivo quello di un’analisi critica di questa formula per la stima della portata media cardiaca Qsvb. La formula proposta nell'articolo verrà verificata nel caso in cui la circolazione sistemica sia approssimata con modelli di tipo windkessel. Dallo studio svolto emerge il fatto che la formula porta risultati con errori trascurabili solo se si approssima la circolazione sistemica con il modello windkessel classico a due elementi (WK2) e la portata aortica con un’onda rettangolare. Approssimando la circolazione sistemica con il modello windkessel a tre elementi (WK3), o descrivendo la portata aortica con un’onda triangolare si ottengono risultati con errori non più trascurabili che variano dal 7%-9% nel caso del WK2 con portata aortica approssimata con onda triangolare ad errori più ampi del 20% nei i casi del WK3 per entrambe le approssimazioni della portata aortica.

Ventricular-arterial coupling in a rat model of reduced arterial compliance provoked by hypervitaminosis D and nicotine

The vitamin D(3) and nicotine (VDN) model is one of isolated systolic hypertension (ISH) in which arterial calcification raises arterial stiffness and vascular impedance. The effects of VDN treatment on arterial and cardiac hemodynamics have been investigated; however, a complete analysis of ventricular-arterial interaction is lacking. Wistar rats were treated with VDN (VDN group, n = 9), and a control group (n = 10) was included without the VDN. At week 8, invasive indexes of cardiac function were obtained using a conductance catheter. Simultaneously, aortic pressure and flow were measured to derive vascular impedance and characterize ventricular-vascular interaction. VDN caused significant increases in systolic (138 +/- 6 vs. 116 +/- 13 mmHg, P < 0.01) and pulse (42 +/- 10 vs. 26 +/- 4 mmHg, P < 0.01) pressures with respect to control. Total arterial compliance decreased (0.12 +/- 0.08 vs. 0.21 +/- 0.04 ml/mmHg in control, P < 0.05), and pulse wave velocity increased significantly (8.8 +/- 2.5 vs. 5.1 +/- 2.0 m/s in control, P < 0.05). The arterial elastance and end-systolic elastance rose significantly in the VDN group (P < 0.05). Wave reflection was augmented in the VDN group, as reflected by the increase in the wave reflection coefficient (0.63 +/- 0.06 vs. 0.52 +/- 0.05 in control, P < 0.05) and the amplitude of the reflected pressure wave (13.3 +/- 3.1 vs. 8.4 +/- 1.0 mmHg in control, P < 0.05). We studied ventricular-arterial coupling in a VDN-induced rat model of reduced arterial compliance. The VDN treatment led to development of ISH and provoked alterations in cardiac function, arterial impedance, arterial function, and ventricular-arterial interaction, which in many aspects are similar to effects of an aged and stiffened arterial tree.