Byla provedena měření turbulentní Rayleighovy-Bénardovy konvekce v kryogenním heliu v rozsahu Rayleighových čísel 1e8

Metoda koneÄnch prvk implementovan v programu COMSOL Multiphysics je vyu­vna k analze tzv. free-standing kmitoÄtovÄ selektivn­ch povrch ve 3D. Tyto modely jsou nslednÄ doplnÄny o periodick© okrajov© podm­nky. Dle jsou free-standing povrchy doplnÄny o vrstvy dielektrika a je zkoumn jejich vliv na modul Äinitele odrazu. V analytick© Ästi jsou vyhodnoceny vlivy poÄtu element diskretizaÄn­ m­ky na pesnost vsledku a d©lku vpoÄt. Vsledky jsou srovnvny vzhledem k vsledkm uvedenm v literatue [5]. V zvÄreÄn© Ästi prce je vysvÄtlen postup pi generovn­ m-file pro obd©ln­kov element a pouit­ globln­ho optimalizaÄn­ho algoritmu PSO, kter automaticky upravuje rozmÄry vodiv©ho motivu tak, aby bylo dosaeno prbÄhu modulu Äinitele odrazu podle poadovan©ho prbÄhu.

La dynamique des flammes de prémélange est étudiée par deux approches numériques différentes. La première résout les équations compressibles de Navier-Stokes avec une chimie simplifiée (DNS). Afin de réduire les coûts de calcul, nous analysons et développons un schéma numérique à grille décalée. Le traitement des ondes acoustiques aux sorties est connu pour rendre les flammes cylindriques légèrement carrées. Ces déformations non-physiques sont expliquées en mettant en évidence la modélisation insuffisamment précise de l'accélération du fluide lorsque l'écoulement est oblique à la sortie. Une étude paramétrique et statistique de flammes turbulentes est menée en 2D et une simulation parallèle 3D est réalisée dans un domaine de (3cm)3. En considérant la flamme infiniment mince, l'approche EEM diminue considérablement les coûts de calcul. Les mêmes simulations sont réalisées et comparées aux résultats de DNS pour tester la capacité du modèle EEM à fournir des résultats quantitatifs. / We study premixed flame dynamics using DNS and EEM approaches. DNS solves compressible Navier-Stokes equations with simplified chemistry. To reduce computational costs and increase efifciency, we analyse and develop a modified staggered scheme. Treatment of acoustic waves at boundaries is known to slightly square cylindrical flames. We try to explain these unphysical distortions by highlighting the poor modeling of fluid acceleration when mainstrean is transverse to outow. A parametric and statistical study of expanding flames is carried out in 2D and also for an expanding (3cm)3 flame. The EEM approach models the flame as an infinitely thin interface. This perturbative strategy dramatically decreases cpu costs. Simulations are carried out and compared to DNS results to check the ability of EEM modeling to give quantitative results.

Non disponible / Non disponible

Title: Bifurcation of Ordinary Differential Equations from Points of Fučík Spectrum Author: Vendula Exnerová Department: Department of Mathematical Analysis Supervisor: doc. RNDr. Jana Stará, CSc., Department of Mathematical Analysis MFF UK, Prague Abstract: The main subject of the thesis is the Fučík spectrum of a system of two differential equations of the second order with mixed boundary conditions. In the first part of the thesis there are described Fučík spectra of problems of a differential equation with Dirichlet, mixed and Neumann boundary conditions. The other part deals with systems of two differential equations. It attends to basic properties of systems and their nontrivial solutions, to a possibility of a reduction of number of parameters and to a dependance of a problem with mixed boundary condition on one with Dirichlet boundary conditions. The thesis takes up the results of E. Massa and B. Ruff about the Dirichlet problem and improves some of their proofs. In the end the Fučík spectrum of a problem with mixed boundary conditions is described as the union of countably many continuously differentiable surfaces and there is proven that this spectrum is closed.

In this thesis is presented mathematical model of stratified 2D flow of viscous incopressible fluid and its program realization. Basic equations of fluid flow in Boussinesq approximation were solved by finite volume method on structured nonortogonal grid. Discretization was done by the principle of semi-discretisation. The space derivative was solved by AUSM me- thod with MUSCL velocity reconstruction. The viscid terms were solved through auxiliary grids. During time discretization artificial compressibility method was used in dual time. The resulting system of ODEs is integrated in time by a suitable Runge-Kutta multistage scheme. Numerical experiments were calculated for flow with Reynolds number equals 1000. Further 3 numerical experiments are presented with different boundary conditions. 1

This thesis outlines an iterative approach for determining the scattered wave for two dimensional multiple acoustic scattering problems using high order local absorbing boundary conditions and second order finite difference. We seek to approximate the total wave as it is scattered off of multiple arbitrarily shaped obstacles. This is done by decomposing the scattered wave into the superposition of single scattered waves. We then repeatedly solve the single scattering system for each obstacle, while updating the boundary conditions based off the incident wave and the scattered wave off the other obstacles. We solve each single scattering by enclosing the obstacle in a circular artificial boundary and generating a curvilinear coordinate system for the computational region between the obstacle and the artificial boundary. We impose an absorbing boundary condition, specifically Karp's Farfield Expansion ABC, on the artificial boundary. We use a finite difference method to discretize the governing equations and to discretize the absorbing boundary conditions. This will create a linear system whose solution will approximate the single scattered wave. The forcing vector of the linear system is determined from the total influence on the obstacle boundary from the incident wave and the scattered waves from the other obstacles. In each iteration, we solve the singular acoustic scattering problem for each obstacle by using the scattered wave approximations from the other obstacles obtained from the previous iteration. The iterations continue until the solutions converge. This iterative method scales well to multiple scattering configurations with many obstacles, and achieves errors on the order of 1E-5 in less than five minutes. This is due to using LU factorization to solve the linear systems, paired with parallelization. I will include numerical results which demonstrate the accuracy and advantages of this iterative technique.

No description available.

In this investigation, a single phase fluid is used to study the coupling between natural convection heat transfer within an enclosure and forced convection through computer covering case to cool the electronic chip. Two working fluids are used (water and air) within a rectangular enclosure and the air flow through the computer case is created by an exhaust fan installed at the back of the computer case. The optimum enclosure size configuration that keeps a maximum temperature of the heat source at a safe temperature level (85℃) is determined. The cooling system is tested for varying values of applied power in the range of 15−40𝑊. The study is based on both numerical models and experimental observations. The numerical work was developed using the commercial software (ANSYS-Icepak) to simulate the flow and temperature fields for the desktop computer and the cooling system. The numerical simulation has the same physical geometry as those used in the experimental investigations. The experimental work was aimed to gather the details for temperature field and use them in the validation of the numerical prediction. The results showed that, the cavity size variations influence both the heat transfer process and the maximum temperature. Furthermore, the experimental results ii compared favourably with those obtained numerically, where the maximum deviation in terms of the maximum system temperature, is within 3.5%. Moreover, it is seen that using water as the working fluid within the enclosure is capable of keeping the maximum temperature under 77℃ for a heat source of 40𝑊, which is below the recommended electronic chips temperature of not exceeding 85℃. As a result, the noise and vibration level is reduced. In addition, the proposed cooling system saved about 65% of the CPU fan power.

yes / There is an increasing need for accurate models describing the electrical behaviour of individual biological cells exposed to electromagnetic fields. In this area of solving linear problem, the most frequently used technique for computing the EM field is the Finite-Difference Time-Domain (FDTD) method. When modelling objects that are small compared with the wavelength, for example biological cells at radio frequencies, the standard Finite-Difference Time-Domain (FDTD) method requires extremely small time-step sizes, which may lead to excessive computation times. The problem can be overcome by implementing a quasi-static approximate version of FDTD, based on transferring the working frequency to a higher frequency and scaling back to the frequency of interest after the field has been computed. An approach to modeling and analysis of biological cells, incorporating the Hodgkin and Huxley membrane model, is presented here. Since the external medium of the biological cell is lossy material, a modified Berenger absorbing boundary condition is used to truncate the computation grid. Linear assemblages of cells are investigated and then Floquet periodic boundary conditions are imposed to imitate the effect of periodic replication of the assemblages. Thus, the analysis of a large structure of cells is made more computationally efficient than the modeling of the entire structure. The total fields of the simulated structures are shown to give reasonable and stable results at 900MHz, 1800MHz and 2450MHz. This method will facilitate deeper investigation of the phenomena in the interaction between EM fields and biological systems. Moreover, the nonlinear response of biological cell exposed to a 0.9GHz signal was discussed on observing the second harmonic at 1.8GHz. In this, an electrical circuit model has been proposed to calibrate the performance of nonlinear RF energy conversion inside a high quality factor resonant cavity with known nonlinear device. Meanwhile, the first and second harmonic responses of the cavity due to the loading of the cavity with the lossy material will also be demonstrated. The results from proposed mathematical model, give good indication of the input power required to detect the weakly effects of the second harmonic signal prior to perform the measurement. Hence, this proposed mathematical model will assist to determine how sensitivity of the second harmonic signal can be detected by placing the required specific input power.