Search results

Some results concerning the well‐posedness of the hydrodynamic model of semiconductor devices in two dimensions are given. We show the non‐ellipticity of the stationary model; give representations which exhibit its elliptic and hyperbolic components, and obtain some appropriate boundary conditions from an examination of the time‐dependent problem.

A mathematical analysis of the time‐dependent multi‐dimensional Hydrodynamic model is performed to determine the well‐posed boundary conditions for semiconductor device simulation. The number of independent boundary conditions that need to be specified at electrical contacts of a semi‐conductor device are derived. Using the classical energy method, a mathematical relation among the physical parameters is established to define the well‐posed boundary conditions for the problem. Several possible sets of boundary conditions are given to illustrate the proper boundary conditions. Natural boundary conditions that can be specified are obtained from the boundary integrals of the weak‐form finite element formulations. An example is included to illustrate the importance of well‐posedness of the boundary conditions for device simulation.

The initial value problem for a semi-linear high-order heat equation is investigated. In the focusing case, global well-posedness and exponential decay are obtained. In the focusing sign, global and non global existence of solutions are discussed via the potential well method.

In this paper, we solve by a finite difference upwinded method an extended hydrodynamic model for semiconductors, with viscous terms in the momentum equation. In particular, we consider the simulation of a one‐dimensional n⁺‐n ‐n⁺ diode, whose solution exhibits at low temperatures strong discontinuities, and investigate the effect of the momentum viscosity on the shock waves. Numerical experiments, performed also on a two‐dimensional test case, demonstrate that the numerical scheme, working on non‐uniform grids, is suitable to describe solutions with strong variations in time and space. Well‐posedness for the boundary conditions is discussed, and a linear stability estimate is established for the one‐dimensional n⁺‐n ‐n⁺ diode benchmark problem.

To provide sufficient conditions for existence, uniqueness and finite element approximability of the solution of time‐harmonic electromagnetic boundary value problems involving metamaterials.

The objectives are achieved by analysing the most simple conditions under which radiation, scattering and cavity problems are well posed and can be reliably solved by the finite element method. The above “most simple conditions” refer to the hypotheses allowing the exploitation of the simplest mathematical tools dealing with the well posedness of variationally formulated problems, i.e. Lax‐Milgram and first Strang lemmas.

The results of interest are found to hold true whenever the effective dielectric permittivity is uniformly positive definite on the regions where no losses are modelled in it and, moreover, the effective magnetic permeability is uniformly negative definite on the regions where no losses are modelled in it. The same good features hold true if “positive” is replaced by “negative” and vice versa in the previous sentence.

It is a priori known that more sophisticated mathematical tools, like Fredholm alternative and compactness results, can provide more general results. However this would require a more complicated analysis and could be considered in a future research.

The design of practical devices involving metamaterials requires the use of reliable electromagnetic simulators. The finite element method is shown to be reliable even when metamaterials are involved, provided some simple conditions are satisfied.

For the first time to the best of authors' knowledge a numerical method is shown to be reliable in problems involving metamaterials.

This paper aims to considered the problem of identification of temperature-dependent thermal conductivity in the nonstationary, nonlinear heat equation. To describe the heat transfer in the furnace charge occupied by a homogeneous porous material, the heat equation is formulated. The inverse problem consists in finding the heat conductivity parameter, which depends on the temperature, from the measurements of the temperature in fixed points of the material.

A numerical method based on the finite-difference scheme and the least squares approach for numerical solution of the direct and inverse problems has been recently developed.

The influence of different numerical scheme parameters on the accuracy of the identified conductivity coefficient is studied. The results of the experiment carried out on real measurements are presented. Their results confirm the ones obtained earlier by using other methods.

Novelty is in a new, easy way to identify thermal conductivity by known temperature measurements. This method is based on special finite-difference scheme, which gives a resolvable system of algebraic equations. The results sensitivity on changes in the method parameters was studies. The algorithms of identification in the case of a purely mathematical experiment and in the case of real measurements, their differences and the practical details are presented.

A generalized variational principle of 2D unsteady compressible flow around oscillating airfoils is established directly from the governing equations and boundary/initial conditions via the semi‐inverse method proposed by He. In this method, an energy integral with an unknown F is used as a trial‐functional. The identification of the unknown F is similar to the identification of the Lagrange multiplier. Based on the variational theory with variable domain, a variational principle for the inverse problem (given as the time‐averaged pressure over the airfoil contour, while the corresponding airfoil shape is unknown) is constructed, and all the boundary/initial conditions are converted into natural ones, leading to well‐posedness and the unique solution of the inverse problems.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from excessive mesh dependence when strain‐softening models are used in numerical analyses and cannot reproduce the size effect commonly observed in quasi‐brittle failure. In this contribution three different approaches will be scrutinized which may be used to remedy these two intimately related deficiencies of the classical theory, namely (i) the addition of higher‐order deformation gradients, (ii) the use of micropolar continuum models, and (iii) the addition of rate dependence. By means of a number of numerical simulations it will be investigated under which conditions these enriched continuum theories permit localization of deformation without losing ellipticity for static problems and hyperbolicity for dynamic problems. For the latter class of problems the crucial role of dispersion in wave propagation in strain‐softening media will also be highlighted.

This work is concerned with the development and the numerical investigation of a hybridizable discontinuous Galerkin (HDG) method for the simulation of two‐dimensional time‐harmonic electromagnetic wave propagation problems.

The proposed HDG method for the discretization of the two‐dimensional transverse magnetic Maxwell equations relies on an arbitrary high order nodal interpolation of the electromagnetic field components and is formulated on triangular meshes. In the HDG method, an additional hybrid variable is introduced on the faces of the elements, with which the element‐wise (local) solutions can be defined. A so‐called conservativity condition is imposed on the numerical flux, which can be defined in terms of the hybrid variable, at the interface between neighbouring elements. The linear system of equations for the unknowns associated with the hybrid variable is solved here using a multifrontal sparse LU method. The formulation is given, and the relationship between the considered HDG method and a standard upwind flux‐based DG method is also examined.

The approximate solutions for both electric and magnetic fields converge with the optimal order of p+1 in L₂ norm, when the interpolation order on every element and every interface is p and the sought solution is sufficiently regular. The presented numerical results show the effectiveness of the proposed HDG method, especially when compared with a classical upwind flux‐based DG method.

The work described here is a demonstration of the viability of a HDG formulation for solving the time‐harmonic Maxwell equations through a detailed numerical assessment of accuracy properties and computational performances.

The purpose of this paper is to study the coupled fixed point problem and the coupled best proximity problem for single-valued and multi-valued contraction type operators defined on cyclic representations of the space. The approach is based on fixed point results for appropriate operators generated by the initial problems.