975 resultados para open-boundary-conditions
Resumo:
The primary goal of this work is related to the extension of an analytic electro-optical model. It will be used to describe single-junction crystalline silicon solar cells and a silicon/perovskite tandem solar cell in the presence of light-trapping in order to calculate efficiency limits for such a device. In particular, our tandem system is composed by crystalline silicon and a perovskite structure material: metilammoniumleadtriiodide (MALI). Perovskite are among the most convenient materials for photovoltaics thanks to their reduced cost and increasing efficiencies. Solar cell efficiencies of devices using these materials increased from 3.8% in 2009 to a certified 20.1% in 2014 making this the fastest-advancing solar technology to date. Moreover, texturization increases the amount of light which can be absorbed through an active layer. Using Green’s formalism it is possible to calculate the photogeneration rate of a single-layer structure with Lambertian light trapping analytically. In this work we go further: we study the optical coupling between the two cells in our tandem system in order to calculate the photogeneration rate of the whole structure. We also model the electronic part of such a device by considering the perovskite top cell as an ideal diode and solving the drift-diffusion equation with appropriate boundary conditions for the silicon bottom cell. We have a four terminal structure, so our tandem system is totally unconstrained. Then we calculate the efficiency limits of our tandem including several recombination mechanisms such as Auger, SRH and surface recombination. We focus also on the dependence of the results on the band gap of the perovskite and we calculare an optimal band gap to optimize the tandem efficiency. The whole work has been continuously supported by a numerical validation of out analytic model against Silvaco ATLAS which solves drift-diffusion equations using a finite elements method. Our goal is to develop a simpler and cheaper, but accurate model to study such devices.
Resumo:
The present work studies a km-scale data assimilation scheme based on a LETKF developed for the COSMO model. The aim is to evaluate the impact of the assimilation of two different types of data: temperature, humidity, pressure and wind data from conventional networks (SYNOP, TEMP, AIREP reports) and 3d reflectivity from radar volume. A 3-hourly continuous assimilation cycle has been implemented over an Italian domain, based on a 20 member ensemble, with boundary conditions provided from ECMWF ENS. Three different experiments have been run for evaluating the performance of the assimilation on one week in October 2014 during which Genova flood and Parma flood took place: a control run of the data assimilation cycle with assimilation of data from conventional networks only, a second run in which the SPPT scheme is activated into the COSMO model, a third run in which also reflectivity volumes from meteorological radar are assimilated. Objective evaluation of the experiments has been carried out both on case studies and on the entire week: check of the analysis increments, computing the Desroziers statistics for SYNOP, TEMP, AIREP and RADAR, over the Italian domain, verification of the analyses against data not assimilated (temperature at the lowest model level objectively verified against SYNOP data), and objective verification of the deterministic forecasts initialised with the KENDA analyses for each of the three experiments.
Resumo:
Sinotubular junction dilation is one of the most frequent pathologies associated with aortic root incompetence. Hence, we create a finite element model considering the whole root geometry; then, starting from healthy valve models and referring to measures of pathological valves reported in the literature, we reproduce the pathology of the aortic root by imposing appropriate boundary conditions. After evaluating the virtual pathological process, we are able to correlate dimensions of non-functional valves with dimensions of competent valves. Such a relation could be helpful in recreating a competent aortic root and, in particular, it could provide useful information in advance in aortic valve sparing surgery.
Resumo:
Argillaceous formations generally act as aquitards because of their low hydraulic conductivities. This property, together with the large retention capacity of clays for cationic contaminants, has brought argillaceous formations into focus as potential host rocks for the geological disposal of radioactive and other waste. In several countries, programmes are under way to characterise the detailed transport properties of such formations at depth. In this context, the interpretation of profiles of natural tracers in pore waters across the formations can give valuable information about the large-scale and long-term transport behaviour of these formations. Here, tracer-profile data, obtained by various methods of pore-water extraction for nine sites in central Europe, are compiled. Data at each site comprise some or all of the conservative tracers: anions (Cl(-), Br(-)), water isotopes (delta(18)O, delta(2)H) and noble gases (mainly He). Based on a careful evaluation of the palaeo-hydrogeological evolution at each site, model scenarios are derived for initial and boundary pore-water compositions and an attempt is made to numerically reproduce the observed tracer distributions in a consistent way for all tracers and sites, using transport parameters derived from laboratory or in situ tests. The comprehensive results from this project have been reported in Mazurek et al. (2009). Here the results for three sites are presented in detail, but the conclusions are based on model interpretations of the entire data set. In essentially all cases, the shapes of the profiles can be explained by diffusion acting as the dominant transport process over periods of several thousands to several millions of years and at the length scales of the profiles. Transport by advection has a negligible influence on the observed profiles at most sites, as can be shown by estimating the maximum advection velocities that still give acceptable fits of the model with the data. The advantages and disadvantages of different conservative tracers are also assessed. The anion Cl(-) is well suited as a natural tracer in aquitards, because its concentration varies considerably in environmental waters. It can easily be measured, although the uncertainty regarding the fraction of the pore space that is accessible to anions in clays remains an issue. The stable water isotopes are also well suited, but they are more difficult to measure and their values generally exhibit a smaller relative range of variation. Chlorine isotopes (delta(37)Cl) and He are more difficult to interpret because initial and boundary conditions cannot easily be constrained by independent evidence. It is also shown that the existence of perturbing events such as the activation of aquifers due to uplift and erosion, leading to relatively sharp changes of boundary conditions, can be considered as a pre-requisite to obtain well-interpretable tracer signatures. On the other hand, gradual changes of boundary conditions are more difficult to parameterise and so may preclude a clear interpretation.
Resumo:
A main field in biomedical optics research is diffuse optical tomography, where intensity variations of the transmitted light traversing through tissue are detected. Mathematical models and reconstruction algorithms based on finite element methods and Monte Carlo simulations describe the light transport inside the tissue and determine differences in absorption and scattering coefficients. Precise knowledge of the sample's surface shape and orientation is required to provide boundary conditions for these techniques. We propose an integrated method based on structured light three-dimensional (3-D) scanning that provides detailed surface information of the object, which is usable for volume mesh creation and allows the normalization of the intensity dispersion between surface and camera. The experimental setup is complemented by polarization difference imaging to avoid overlaying byproducts caused by inter-reflections and multiple scattering in semitransparent tissue.
Resumo:
We present examples of isospectral operators that do not have the same heat content. Several of these examples are planar polygons that are isospectral for the Laplace operator with Dirichlet boundary conditions. These include examples with infinitely many components. Other planar examples have mixed Dirichlet and Neumann boundary conditions. We also consider Schrodinger operators acting in L-2[0,1] with Dirichlet boundary conditions, and show that an abundance of isospectral deformations do not preserve the heat content.
Resumo:
In a recent study of the self-adjoint extensions of the Hamiltonian of a particle confined to a finite region of space, in which we generalized the Heisenberg uncertainty relation to a finite volume, we encountered bound states localized at the wall of the cavity. In this paper, we study this situation in detail both for a free particle and for a hydrogen atom centered in a spherical cavity. For appropriate values of the self-adjoint extension parameter, the bound states localized at the wall resonate with the standard hydrogen bound states. We also examine the accidental symmetry generated by the Runge–Lenz vector, which is explicitly broken in a spherical cavity with general Robin boundary conditions. However, for specific radii of the confining sphere, a remnant of the accidental symmetry persists. The same is true for an electron moving on the surface of a finite circular cone, bound to its tip by a 1/r1/r potential.
Resumo:
The problem of re-sampling spatially distributed data organized into regular or irregular grids to finer or coarser resolution is a common task in data processing. This procedure is known as 'gridding' or 're-binning'. Depending on the quantity the data represents, the gridding-algorithm has to meet different requirements. For example, histogrammed physical quantities such as mass or energy have to be re-binned in order to conserve the overall integral. Moreover, if the quantity is positive definite, negative sampling values should be avoided. The gridding process requires a re-distribution of the original data set to a user-requested grid according to a distribution function. The distribution function can be determined on the basis of the given data by interpolation methods. In general, accurate interpolation with respect to multiple boundary conditions of heavily fluctuating data requires polynomial interpolation functions of second or even higher order. However, this may result in unrealistic deviations (overshoots or undershoots) of the interpolation function from the data. Accordingly, the re-sampled data may overestimate or underestimate the given data by a significant amount. The gridding-algorithm presented in this work was developed in order to overcome these problems. Instead of a straightforward interpolation of the given data using high-order polynomials, a parametrized Hermitian interpolation curve was used to approximate the integrated data set. A single parameter is determined by which the user can control the behavior of the interpolation function, i.e. the amount of overshoot and undershoot. Furthermore, it is shown how the algorithm can be extended to multidimensional grids. The algorithm was compared to commonly used gridding-algorithms using linear and cubic interpolation functions. It is shown that such interpolation functions may overestimate or underestimate the source data by about 10-20%, while the new algorithm can be tuned to significantly reduce these interpolation errors. The accuracy of the new algorithm was tested on a series of x-ray CT-images (head and neck, lung, pelvis). The new algorithm significantly improves the accuracy of the sampled images in terms of the mean square error and a quality index introduced by Wang and Bovik (2002 IEEE Signal Process. Lett. 9 81-4).
Resumo:
Purpose: Development of an interpolation algorithm for re‐sampling spatially distributed CT‐data with the following features: global and local integral conservation, avoidance of negative interpolation values for positively defined datasets and the ability to control re‐sampling artifacts. Method and Materials: The interpolation can be separated into two steps: first, the discrete CT‐data has to be continuously distributed by an analytic function considering the boundary conditions. Generally, this function is determined by piecewise interpolation. Instead of using linear or high order polynomialinterpolations, which do not fulfill all the above mentioned features, a special form of Hermitian curve interpolation is used to solve the interpolation problem with respect to the required boundary conditions. A single parameter is determined, by which the behavior of the interpolation function is controlled. Second, the interpolated data have to be re‐distributed with respect to the requested grid. Results: The new algorithm was compared with commonly used interpolation functions based on linear and second order polynomial. It is demonstrated that these interpolation functions may over‐ or underestimate the source data by about 10%–20% while the parameter of the new algorithm can be adjusted in order to significantly reduce these interpolation errors. Finally, the performance and accuracy of the algorithm was tested by re‐gridding a series of X‐ray CT‐images. Conclusion: Inaccurate sampling values may occur due to the lack of integral conservation. Re‐sampling algorithms using high order polynomialinterpolation functions may result in significant artifacts of the re‐sampled data. Such artifacts can be avoided by using the new algorithm based on Hermitian curve interpolation
Resumo:
This dissertation concerns convergence analysis for nonparametric problems in the calculus of variations and sufficient conditions for weak local minimizer of a functional for both nonparametric and parametric problems. Newton's method in infinite-dimensional space is proved to be well-defined and converges quadratically to a weak local minimizer of a functional subject to certain boundary conditions. Sufficient conditions for global converges are proposed and a well-defined algorithm based on those conditions is presented and proved to converge. Finite element discretization is employed to achieve an implementable line-search-based quasi-Newton algorithm and a proof of convergence of the discretization of the algorithm is included. This work also proposes sufficient conditions for weak local minimizer without using the language of conjugate points. The form of new conditions is consistent with the ones in finite-dimensional case. It is believed that the new form of sufficient conditions will lead to simpler approaches to verify an extremal as local minimizer for well-known problems in calculus of variations.
Resumo:
Scaphoid is one of the 8 carpal bones found adjacent to the thumb supported proximally by Radius bone. During the free fall, on outstretched hand, the impact load gets transferred to the scaphoid at its free anterior end. Unique arrangement of other carpal bones in the palm is also one of the reasons for the load to get transferred to scaphoid. About half of the total load acting upon carpal bone gets transferred to scaphoid at its distal pole. There are about 10 to 12 clinically observed fracture pattern in the scaphoid due to free fall. The aim of the study is to determine the orientation of the load, magnitude of the load and the corresponding fracture pattern. This study includes both static and dynamic finite element models validated by experiments. The scaphoid model has been prepared from CT scans of a 27 year old person. The 2D slices of the CT scans have been converted to 3D model by using MIMICS software. There are four cases of loading studied which are considered to occur clinically more frequently. In case (i) the load is applied at the posterior end at distal pole whereas in case (ii), (iii) and (iv), the load is applied at anterior end at different directions. The model is given a fixed boundary condition at the region which is supported by Radius bone during the impact. Same loading and boundary conditions have been used in both static and dynamic explicit finite element analysis. The site of fracture initiation and path of fracture propagation have been identified by using max principal stress / gradient and max principal strain / gradient criterion respectively in static and dynamic explicit finite element analysis. Static and dynamic impact experiments were performed on the polyurethane foam specimens to validate the finite element results. Experimental results such as load at fracture, site of fracture initiation and path of fracture propagation have been compared with the results of finite element analysis. Four different types of fracture patterns observed in clinical studies have been identified in this study.
Resumo:
This technical report discusses the application of Lattice Boltzmann Method (LBM) in the fluid flow simulation through porous filter-wall of disordered media. The diesel particulate filter (DPF) is an example of disordered media. DPF is developed as a cutting edge technology to reduce harmful particulate matter in the engine exhaust. Porous filter-wall of DPF traps these soot particles in the after-treatment of the exhaust gas. To examine the phenomena inside the DPF, researchers are looking forward to use the Lattice Boltzmann Method as a promising alternative simulation tool. The lattice Boltzmann method is comparatively a newer numerical scheme and can be used to simulate fluid flow for single-component single-phase, single-component multi-phase. It is also an excellent method for modelling flow through disordered media. The current work focuses on a single-phase fluid flow simulation inside the porous micro-structure using LBM. Firstly, the theory concerning the development of LBM is discussed. LBM evolution is always related to Lattice gas Cellular Automata (LGCA), but it is also shown that this method is a special discretized form of the continuous Boltzmann equation. Since all the simulations are conducted in two-dimensions, the equations developed are in reference with D2Q9 (two-dimensional 9-velocity) model. The artificially created porous micro-structure is used in this study. The flow simulations are conducted by considering air and CO2 gas as fluids. The numerical model used in this study is explained with a flowchart and the coding steps. The numerical code is constructed in MATLAB. Different types of boundary conditions and their importance is discussed separately. Also the equations specific to boundary conditions are derived. The pressure and velocity contours over the porous domain are studied and recorded. The results are compared with the published work. The permeability values obtained in this study can be fitted to the relation proposed by Nabovati [8], and the results are in excellent agreement within porosity range of 0.4 to 0.8.
Resumo:
Die Topologieoptimierung hat sich in den letzten Jahren zu einer sehr praktischen und vielseitig ein-setzbaren Design- und Entwicklungsmethode entwickelt. Diese Methode soll nun an einer Seilscheibe, die aus der Industrie nicht mehr wegzudenken ist, angewendet werden. Im Vordergrund steht vor allem die Reduzierung der Masse sowie die Anpassung der Speichenform an die unterschiedlichen Randbedingungen.
Resumo:
In July and August 2010 floods of unprecedented impact afflicted Pakistan. The floods resulted from a series of intense multi-day precipitation events in July and early August. At the same time a series of blocking anticyclones dominated the upper-level flow over western Russia and breaking waves i.e. equatorward extrusions of stratospheric high potential vorticity (PV) air formed along the downstream flank of the blocks. Previous studies suggested that these extratropical upper-level breaking waves were crucial for instigating the precipitation events in Pakistan. Here a detailed analysis is provided of the extratropical forcing of the precipitation. Piecewise PV inversion is used to quantify the extratropical upper-level forcing associated with the wave breaking and trajectories are calculated to study the pathways and source regions of the moisture that precipitated over Pakistan. Limited-area model simulations are carried out to complement the Lagrangian analysis. The precipitation events over Pakistan resulted from a combination of favourable boundary conditions with strong extratropical and monsoonal forcing factors. Above-normal sea-surface temperatures in the Indian Ocean led to an elevated lower-tropospheric moisture content. Surface monsoonal depressions ensured the transport of moist air from the ocean towards northeastern Pakistan. Along this pathway the air parcel humidity increased substantially (60–90% of precipitated moisture) via evapotranspiration from the land surface. Extratropical breaking waves influenced the surface wind field substantially by enhancing the wind component directed towards the mountains which reinforced the precipitation.
Resumo:
The goal of this paper is to establish exponential convergence of $hp$-version interior penalty (IP) discontinuous Galerkin (dG) finite element methods for the numerical approximation of linear second-order elliptic boundary-value problems with homogeneous Dirichlet boundary conditions and piecewise analytic data in three-dimensional polyhedral domains. More precisely, we shall analyze the convergence of the $hp$-IP dG methods considered in [D. Schötzau, C. Schwab, T. P. Wihler, SIAM J. Numer. Anal., 51 (2013), pp. 1610--1633] based on axiparallel $\sigma$-geometric anisotropic meshes and $\bm{s}$-linear anisotropic polynomial degree distributions.