Original Boundary Conditions Research Articles

A Shifted Boundary Method for the compressible Euler equations

The Shifted Boundary Method (SBM) is applied to compressible Euler flows, with and without shock discontinuities. The SBM belongs to the class of unfitted (or immersed, or embedded) finite element methods and avoids integration over cut cells (and the associated implementation/stability issues) by reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. We specifically discuss the advantages the proposed method offers in avoiding spurious numerical artifacts in two scenarios: (a) when curved boundaries are represented by body-fitted polygonal approximations and (b) when the Kutta condition needs to be imposed in immersed simulations of airfoils. An extensive suite of numerical tests is included.

Investigation of the influence of different boundary conditions in SPH on ship dynamics

This paper explores ship dynamics through numerical investigation using Smoothed Particle Hydrodynamics (SPH). To enhance the performance of SPH simulations, we implement the Dummy Particle Condition (DPC) floating-fluid interaction scheme in our numerical model, substituting the original Dynamic Boundary Condition (DBC). The accuracy of the updated model is validated by conducting a series of two-dimensional (2D) benchmark tests. Subsequently, based on both DBC and DPC, we simulate three-dimensional (3D) ship motion in a series of regular waves, including low and high waves, and compare the simulation results with model test data. It is demonstrated that, compared with the original DBC, the DPC floating-fluid interaction scheme effectively enhances SPH simulation performance in terms of ship dynamics for all wave conditions. Then, the correlation between simulation accuracy, wave height, and particle size in SPH is quantitatively discussed, and a clear relationship is observed for DPC results.

A time viscosity-splitting method for incompressible flows with temperature-dependent viscosity and thermal conductivity

A fractional-step method is proposed and analyzed for solving the incompressible thermal Navier–Stokes equations coupled to the convection–conduction equation for heat transfer with a generalized source term for which the viscosity and thermal conductivity are temperature-dependent under the Boussinesq assumption. The proposed method consists of four steps all based on a viscosity-splitting algorithm where the convection and diffusion terms of both velocity and temperature solutions are separated while a viscosity term is kept in the correction step at all times. This procedure preserves the original boundary conditions on the corrected velocity and it removes any pressure inconsistencies. As a main feature, our method allows the temperature to be transported by a non-divergence-free velocity, in which case we show how to handle the subtle temperature convection term in the error analysis and establish full first-order error estimates for the velocity and the temperature solutions and 1/2-order estimates for the pressure solution in their appropriate norms. The theoretical results are examined by an accuracy test example with known analytical solution and using a benchmark problem of Rayleigh–Bénard convection with temperature-dependent viscosity and thermal conductivity. We also apply the method for solving a problem of unsteady flow over a heated airfoil. The obtained results demonstrate the convergence, accuracy and applicability of the proposed time viscosity-splitting method.

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

Numerical-analytical approach to solving problems of non-stationary thermal conductivity of a non-thin annular plate

This paper considers the first stage of calculating the initial boundary value problem of non-stationary thermal conductivity of cylindrical bodies using a modified method of lines, namely dimension reduction of the original differential equations, initial and boundary conditions. The original equations of thermal conductivity are defined in a cylindrical coordinate system in a spatial setting. An object is a cylindrical body with commensurate dimensions. This area of research is relevant, because when calculating the load bearing elements of structures to thermal effects, the first step is to determine the distribution of temperature fields. Boundary conditions are considered as conditions of convective heat transfer, which by means of boundary transition are transformed into boundary conditions of the first and second types. Dimension reduction with respect to spatial coordinate is performed by the Bubnov-Galorkin-Petrov projection method using local basis functions. These functions are called "cover" functions, which are related to the lines drawn on the domain of the task. Normalized trigonometric series are used to reduce the dimensionality of the equations with respect to circular coordinate. All transformations are performed in index form. In addition to differential equations, the projection method reduces the dimensionality of the initial and boundary conditions. In this paper, the most optimal form of writing reduced equations is determined, which provides the ease of reducing the dimensionality of the original differential equations. Initial and boundary conditions take into account the impact of the environment. All this makes it possible to set a reduced initial boundary value problem for further calculation by numerical finite difference methods.

Математичні моделі локального нагрівання елементів електронних пристроїв

Linear and non-linear mathematical models for the determination of the temperature field, and subsequently for the analysis of temperature regimes in isotropic spatial heat-active media sub-jected to internal and external local heat load, have been developed. In the case of nonlinear boundary-value problems, the Kirchhoff transformation was applied, using which the original nonlinear heat conduction equations and nonlinear boundary conditions were linearized, and as a result, linearized second-order differential equations with partial derivatives and a discontinu-ous right-hand side and partially linearized boundary conditions were obtained. For the final linearization of the partially linearized differential equation and boundary conditions, the ap-proximation of the temperature according to one of the spatial coordinates on the boundary sur-faces of the inclusion was performed by piecewise constant functions. To solve linear bounda-ry-value problems, as well as obtained linearized boundary-value problems with respect to the Kirchhoff transformation, the Henkel integral transformation method was used, as a result of which analytical solutions of these problems were obtained. For a heat-sensitive environment, as an example, a linear dependence of the coefficient of thermal conductivity of the structural material of the structure on temperature, which is often used in many practical problems, was chosen. As a result, analytical relations for determining the temperature distribution in this envi-ronment were obtained. Numerical analysis of temperature behavior as a function of spatial co-ordinates for given values of geometric and thermophysical parameters was performed. The in-fluence of the power of internal heat sources and environmental materials on the temperature distribution was studied. To determine the numerical values of the temperature in the given structure, as well as to analyze the heat exchange processes in the middle of these structures, caused by the internal and external heat load, software tools were developed, using which a ge-ometric image of the temperature distribution depending on the spatial coordinates was made.

Mathematical models for analysis of temperature regimes in vehicle braking systems

Linear and non-linear mathematical models for the determination of the temperature field, and subsequently for the analysis of temperature regimes in the braking systems of vehicles, which are geometrically depicted as isotropic spatial heat-active media that are subject to internal local thermal heating, have been developed. With the use of classical methods, it is not possible to obtain analytical solutions of linear and nonlinear boundary value problems of mathematical physics in a closed form. This is especially the case when the right-hand sides of differential equations with partial derivatives and boundary conditions are discontinuous functions. The given approach is based on the application of the apparatus of generalized functions to describe the local concentration of thermal influence. This made it possible to apply the integral transformation and, on this basis, to obtain analytical solutions of both linear and nonlinear boundary value problems. In the case of a nonlinear boundary value problem, the Kirchhoff transformation was applied, using which the original nonlinear heat conduction equation and nonlinear boundary conditions were linearized, and as a result, a linearized second-order differential equation with partial derivatives and boundary conditions with a discontinuous right-hand side were obtained. To solve the linear boundary value problem, as well as the obtained linearized boundary value problem with respect to the Kirchhoff transformation, the Henkel integral transformation method was used, as a result of which analytical solutions of these problems were obtained. For a heat-sensitive environment, as an example, a linear dependence of the coefficient of thermal conductivity of the structural material of the structure on temperature, which is often used in many practical problems, was chosen. As a result, an analytical relationship was obtained for determining the temperature distribution in this medium. On the basis of the developed mathematical models, a computational algorithm was created and on this basis, software tools were created, using which the heat exchange processes in the middle of the brake structures for the selected materials of the brake pads were analyzed in terms of their effectiveness, as well as the determination of the optimal temperature values for the effective operation of the braking system of vehicles. The developed linear and nonlinear mathematical models for determining the temperature field in spatial heat-active media with internal heating make it possible to analyze their thermal stability. As a result, it becomes possible to increase it and protect it from overheating, which can cause the destruction of not only individual nodes and individual elements, but also the entire structure.

Error estimates for a viscosity-splitting scheme in time applied to non-Newtonian fluid flows

A time fractional-step method is presented for numerical solutions of the incompressible non-Newtonian fluids for which the viscosity is non-linear depending on the shear-rate magnitude according to a generic model. The method belongs to a class of viscosity-splitting procedures and it consists of separating the convection term and incompressibility constraint into two time steps. Unlike the conventional projection methods, the viscosity is not dropped in the last step allowing to enforce the full original boundary conditions on the end-of-step velocity which eliminates any concerns about the numerical boundary layers. We carry out a rigorous error analysis and provide a full first-order error estimate for both the velocity and pressure solutions in the relevant norms. Numerical results are presented for two test examples of non-Newtonian fluid flows to demonstrate the theoretical analysis and confirm the reliability of this viscosity-splitting scheme.

Investigating hygrothermal bending behavior of FG nanobeams via local/nonlocal stress gradient theory of elasticity with general boundary conditions

The bending response of functionally graded (FG) nanobeams under hygrothermal loading was investigated to emphasize the different scenarios that arise when using simplified and original boundary conditions. The governing equations were derived by using the principle of virtual work on the basis of the local/nonlocal stress gradient theory of elasticity. A Wolfram language code in Mathematica was written by the authors to develop a numerical investigation for different values of the material gradient index, the gradient length parameter, the nonlocal parameter, and considering two distinct types of thermal loading, that is, uniform temperature rise and heat conduction across the thickness of FG nanobeam cross-section.

Математичні моделі температурного поля у елементах електронних пристроїв із чужорідним включенням

Linear and non-linear mathematical models for the determination of the temperature field, and subsequently for the analysis of temperature regimes in isotropic media with semi-transparent foreign inclusions that are also subject to external heat load, have been developed. To solve the nonlinear boundary-value problem, a linearizing function was introduced, using which the original nonlinear heat conduction equation and nonlinear boundary conditions were linearized, and as a result, a partially linearized second-order differential equation with partial derivatives and discontinuous and singular coefficients relative to the linearizing function and partially linearized boundary conditions was obtained. For the final linearization of the partially linearized differential equation and boundary conditions, the approximation of the temperature according to one of the spatial coordinates on the boundary surfaces of the inclusion was performed by piecewise constant functions. To solve the obtained linear boundary value problem, the Hankel integral transformation method was used, as a result of which an analytical solution was obtained, which determines the introduced linearizing function. For the numerical analysis of temperature behavior and heat exchange processes caused by external heat load, software tools have been developed, which are used to create a geometric image of temperature distribution depending on spatial coordinates. The obtained results indicate the correspondence of the developed mathematical models of the analysis of heat exchange processes in heterogeneous media with external heating to a real physical process. The obtained results indicate the correspondence of the developed mathematical models of the analysis of heat exchange processes in heterogeneous media with external heating to a real physical process. The developed models make it possible to analyze environments with foreign elements under external thermal loads regarding their thermal resistance. As a result, it becomes possible to increase it and protect it from overheating, which can cause the destruction of not only individual elements, but also the entire structure.

New Difference Schemes for Impedance Boundary Conditions of Wide-Angle Shift-Map PE

New impedance boundary conditions (IBC) difference schemes are proposed for wide-angle shift-map PE. A simple experiment is set up to check the accuracy of the PE. We explain the drawbacks of original boundary conditions (BCs) difference schemes. Some problems are caused by the PEC boundary. We use one-sided difference schemes for the BCs as substitutes. The proposed schemes not only meet the physical property but also perform well in numerical calculation. The numerical experiments are done on PEC and media boundaries to study the characteristic of the new BCs. The reflection calculated by PE are compared with Fresnel reflection. The numerical analysis reveals that the new BCs scheme errors are smaller than those of the original BCs. The new 1st-order PE is accuracy for the slope up to 14° with the reflection angle offset being less than 0.5°, where a horizontally spreading, 0.5° beamwidth levy source illuminates the slope. Under the same condition, the 2nd-order Claerbout PE is accurate when the slope is less than 16° and the 2nd-order Greene PE is accurate when the slope is less than 19°. In addition, the ripples caused by the 2nd-order PE with new BCs is smaller than the original one.

An improved partial differential equation filter scheme for topology optimization of additively manufactured coated structure

In this paper, an improved partial differential equation (PDE) filter is presented for density-based topology optimization of coated structures with both regular and irregular design domains. The proposed filter efficiently addresses the boundary effect problems that causes identification failure of the material interface for modeling coating layer during topology optimization by replacing the original Dirichlet boundary conditions with non-homogeneous Neumann boundary conditions. One attractive feature of the proposed approach is that it can serve as an ideal alternative to the common domain extension approach, particularly for designing geometrically complex engineering structures with irregular design domain meshed by unstructured mesh. Based on the improved PDE filter, a new topology optimization design framework for coated structures with complete coating layer and controllable coating thickness is achieved. Examples of designing regular and irregular two-dimensional and three-dimensional coated structures are provided to demonstrate the effectiveness of the proposed approach.

Mathematical Modeling of Single and Phase Autowaves in a Ferrocolloid

This paper describes a mathematical model of an autowave process in a cell with a ferrocolloid. The model is a system of differential coupled equations of the second order and differs from the previously presented model in terms of its original boundary conditions. The mathematical modeling of autowaves presented in this work constitutes an innovative approach, since the characteristics of the wave process are not initially included in the model but the model demonstrates a wave motion. A 2D solution of the model, which shows the correctness of the described mechanism of the autowave process, i.e., the recharging of magnetic particles in dense near-electrode layers formed near the electrodes under the influence of an electric field, is obtained. The propagation of single and phase autowaves is demonstrated in a computer experiment.

МАТЕМАТИЧНІ МОДЕЛІ ВИЗНАЧЕННЯ ТЕМПЕРАТУРНИХ ПОЛІВ У НЕОДНОРІДНИХ ЕЛЕМЕНТАХ ЦИФРОВИХ ПРИСТРОЇВ ІЗ УРАХУВАННЯМ ТЕРМОЧУТЛИВОСТІ Primary tabs View(active tab) Edit Translate

Linear and nonlinear mathematical models for determining the temperature field and subsequently analyzing temperature regimes in isotropic spatial media with semi-through foreign inclusions subjected to internal and external thermal loads are developed. For this purpose, the heat transfer coefficient for such structures is described as a single unit using asymmetric unit functions, which makes it possible to consider boundary value problems of heat transfer with one linear and nonlinear differential equations of heat transfer with discontinuous and singular coefficients and linear and nonlinear boundary conditions on the boundary surfaces of the media. In the case of a nonlinear boundary value problem, the introduced linearizing function is used to linearize the original nonlinear heat conduction equation and nonlinear boundary conditions, and as a result, a partially linearized second-order differential equation with partial derivatives and discontinuous and singular coefficients is obtained relative to the linearizing function with partially linearized boundary conditions. For the final linearization of the partially linearized differential equation and boundary conditions, the temperature is approximated by one of the spatial coordinates on the boundary surfaces of the inclusion by piecewise linear functions, as a result of which both the differential equation and boundary conditions become fully linearized. To solve the resulting linear boundary value problem, the Hankel integral transformation method is used, which results in an analytical solution that determines the introduced linearizing function. As an example, the linear dependence of the thermal conductivity coefficient of structural materials of a structure on temperature, which is often used in many practical problems, is chosen. As a result, analytical relations in the form of quadratic equations were obtained to determine the temperature distribution in a thermally sensitive layer with a foreign semi-through inclusion under external heating in the form of a heat flux. A numerical analysis of the temperature behavior as a function of spatial coordinates for given values of geometric and thermophysical parameters is performed. The influence of a foreign inclusion on the temperature distribution is investigated if the VK94-I ceramic is chosen as the material of the medium and the inclusion is silver. To determine the numerical values of temperature in the above structures, as well as to analyze heat transfer processes inside these structures caused by internal and external thermal loads, software tools have been developed that have been used to perform a geometric image of the temperature distribution depending on spatial coordinates. The obtained numerical temperature values indicate that the developed mathematical models for analyzing heat transfer processes in spatially heterogeneous environments with internal and external heating correspond to a real physical process. The software also makes it possible to analyze such environments subjected to internal and external thermal loads in terms of their thermal resistance. As a result, it becomes possible to increase it and protect it from overheating, which can cause the destruction of not only individual elements but also the entire structure.

Математичні моделі визначення температурних полів у теплоактивних елементах цифрових пристроїв з локальним внутрішнім нагріванням та із урахуванням термочутливості

Linear and non-linear mathematical models for the determination of the temperature field, and subsequently for the analysis of temperature regimes in isotropic spatial heat-active media subjected to internal local heat load, have been developed. In the case of a nonlinear boundary-value problem, the Kirchhoff transformation is applied, using which the original nonlinear heat conduction equation and nonlinear boundary conditions are linearized, and as a result, a linearized second-order differential equation with partial derivatives and a discontinuous right-hand side and partially linearized boundary conditions is obtained. For the final linearization of the partially linearized boundary conditions, the approximation of the temperature by the radial spatial coordinate on the boundary surface of the thermosensitive medium was performed by a piecewise constant function, as a result of which the boundary value problem was obtained completely linearized. To solve the linear boundary value problem, as well as the obtained linearized boundary value problem with respect to the Kirchhoff transformation, the Henkel integral transformation method was used, as a result of which analytical solutions of these problems were obtained. For a heat-sensitive environment, as an example, a linear dependence of the coefficient of thermal conductivity of the structural material of the structure on temperature, which is often used in many practical problems, was chosen. As a result, an analytical relationship was obtained for determining the temperature distribution in this medium. Numerical analysis of temperature behavior as a function of spatial coordinates for given values of geometric and thermophysical parameters was performed. The influence of the power of internal heat sources and environmental materials on the temperature distribution was studied. To determine the numerical values of the temperature in the given structure, as well as to analyze the heat exchange processes in the middle of these structures, caused by the internal heat load, software tools were developed, using which a geometric image of the temperature distribution depending on the spatial coordinates was made. The developed linear and nonlinear mathematical models for determining the temperature field in spatial heat-active environments with internal heating testify to their adequacy to a real physical process. They make it possible to analyze such environments for their thermal stability. As a result, it becomes possible to increase it and protect it from overheating, which can cause the destruction of not only individual nodes and their elements, but also the entire structure.

A Diffuse-Domain Phase-Field Lattice Boltzmann Method for Two-Phase Flows in Complex Geometries

This study constructs a new diffuse-domain (DD) lattice Boltzmann (LB) method for two-phase flows in complex geometries. We first coupled the DD method with the consistent and conservative phase-field (CCPF) method which can be described by the DD-CC Navier–Stokes–Cahn–Hilliard (NSCH) equations with the consistency of reduction, the consistency of mass and momentum transport, and the consistency of mass conservation. Then we conducted a matched asymptotic analysis and found that the DD-CCNSCH equations would converge to the CCNSCH equations as the thickness of the DD interface approaches to zero. Since the boundary conditions imposed on the complex geometries are incorporated into the DD-CCNSCH equations as the source terms, the system can be implemented more readily in a large and regular domain, and there is no need to treat the complex boundary conditions specially. We further developed a DD-CCPFLB method, and through the Chapman–Enskog analysis, the DD-CCPFLB method can correctly recover the DD-CCNSCH equations for two-phase flows in complex geometries. To quantitatively validate the DD-CCPFLB method, three benchmark problems are considered, and the numerical results are in agreement with those obtained through directly solving the CCNSCH equations combined with the original boundary conditions, the analytical solution, or the reported data. Finally, the present method is used to study the problems of two-phase flows in complex geometries, including the droplet dynamics in the Y-shaped channel and a droplet passing through cylindrical obstacles. The effects of the wettability and viscosity ratio are mainly investigated, and it is found that these two factors have significant impacts on the droplet dynamics.

Extension of Babinet Principle to Complementary Metallic Elements on the Interface of Different Substrates

In this letter, the existing Babinet principle is extended to two-dimensional complementary metallic elements on the interface of different substrates. Specifically, the existing Babinet principle requests homogenized substrate to cover complementary metallic elements. In order to break this restriction, the concept of effective relative permittivity is introduced to keep the original boundary condition unchanged. Further, the extended Babinet principle is expressed as a mathematical equation including vector fields of complementary metallic elements on the substrate interface. Especially, the proposed theory is related to the accuracy of the effective permittivity. In order to verify the proposed theory, the extended Babinet principle is applied to complementary metasurfaces with the approximate effective permittivity. Great consistence exists between theoretical and simulated results. Thus, the extended Babinet principle not only provides a theoretical approach to analyze complementary metallic elements on the substrate interface but also greatly enriches the physical connotation of existing Babinet principle.

Penalization method for the Navier–Stokes–Fourier system

We apply the method of penalization to the Dirichlet problem for the Navier–Stokes–Fourier system governing the motion of a general viscous compressible fluid confined to a bounded Lipschitz domain. The physical domain is embedded into a large cube on which the periodic boundary conditions are imposed. The original boundary conditions are enforced through a singular friction term in the momentum equation and a heat source/sink term in the internal energy balance. The solutions of the penalized problem are shown to converge to the solution of the limit problem. In particular, we extend the available existence theory to domains with rough (Lipschitz) boundary. Numerical experiments are performed to illustrate the efficiency of the method.

Thermal stresses in an elastic parallelepiped

A semi-analytical algorithm is presented for solving a three-dimensional thermoelasticity problem for a parallelepiped with free edges. By implementing the direct integration method, the problem is reduced to the determination of three Vihak key functions introduced as the integrals of the equilibrium equations. The governing equations for the Vihak functions are derived based on the compatibility equations in terms of stresses. Based on original boundary conditions, the integral conditions are derived for the Vihak functions verbalizing their resultant “force” and “momentum.” An approach for separating variables in the governing equations and integral conditions is suggested by implementing special sets of the associated- and eigenfunctions. Numerical case studies are discussed with special attention given to the effect of the Poisson ratio.

Schwarz methods by domain truncation

Schwarz methods use a decomposition of the computational domain into subdomains and need to impose boundary conditions on the subdomain boundaries. In domain truncation one restricts the unbounded domain to a bounded computational domain and must also put boundary conditions on the computational domain boundaries. In both fields there are vast bodies of literature and research is very active and ongoing. It turns out to be fruitful to think of the domain decomposition in Schwarz methods as a truncation of the domain onto subdomains. Seminal precursors of this fundamental idea are papers by Hagstrom, Tewarson and Jazcilevich (1988), Després (1990) and Lions (1990). The first truly optimal Schwarz method that converges in a finite number of steps was proposed by Nataf (1993), and used precisely transparent boundary conditions as transmission conditions between subdomains. Approximating these transparent boundary conditions for fast convergence of Schwarz methods led to the development of optimized Schwarz methods – a name that has become common for Schwarz methods based on domain truncation. Compared to classical Schwarz methods, which use simple Dirichlet transmission conditions and have been successfully used in a wide range of applications, optimized Schwarz methods are much less well understood, mainly due to their more sophisticated transmission conditions.A key application of Schwarz methods with such sophisticated transmission conditions turned out to be time-harmonic wave propagation problems, because classical Schwarz methods simply do not work in this case. The past decade has given us many new Schwarz methods based on domain truncation. One review from an algorithmic perspective (Gander and Zhang 2019) showed the equivalence of many of these new methods to optimized Schwarz methods. The analysis of optimized Schwarz methods, however, is lagging behind their algorithmic development. The general abstract Schwarz framework cannot be used for the analysis of these methods, and thus there are many open theoretical questions about their convergence. Just as for practical multigrid methods, Fourier analysis has been instrumental for understanding the convergence of optimized Schwarz methods and for tuning their transmission conditions. Similar to local Fourier mode analysis in multigrid, the unbounded two-subdomain case is used as a model for Fourier analysis of optimized Schwarz methods due to its simplicity. Many aspects of the actual situation,e.g.boundary conditions of the original problem and the number of subdomains, were thus neglected in the unbounded two-subdomain analysis. While this gave important insight, new phenomena beyond the unbounded two-subdomain models were discovered.This present situation is the motivation for our survey: to give a comprehensive review and precise exploration of convergence behaviours of optimized Schwarz methods based on Fourier analysis, taking into account the original boundary conditions, many-subdomain decompositions and layered media. We consider as our model problem the operator$-\Delta + \eta $in the diffusive case$\eta>0$(screened Laplace equation) or the oscillatory case$\eta <0$(Helmholtz equation), in order to show the fundamental difference in behaviour of Schwarz solvers for these problems. The transmission conditions we study include the lowest-order absorbing conditions (Robin), and also more advanced perfectly matched layers (PMLs), both developed first for domain truncation. Our intensive work over the last two years on this review has led to several new results presented here for the first time: in the bounded two-subdomain analysis for the Helmholtz equation, we see strong influence of the original boundary conditions imposed on the global problem on the convergence factor of the Schwarz methods, and the asymptotic convergence factors with small overlap can differ from the unbounded two-subdomain analysis. In the many-subdomain analysis, we find the scaling with the number of subdomains,e.g.when the subdomain size is fixed, robust convergence of the double-sweep Schwarz method for the free-space wave problem, either with fixed overlap and zeroth-order Taylor conditions or with a logarithmically growing PML, and we find that Schwarz methods with PMLs work like smoothers that converge faster for higher Fourier frequencies; in particular, for the free-space wave problem, plane waves (in the error) passing through interfaces at a right angle converge more slowly. In addition to our main focus on analysis in Sections 2 and 3, we start in Section 1 with an expository historical introduction to Schwarz methods, and in Section 4 we give a brief interpretation of the recently proposed optimal Schwarz methods for decompositions with cross-points from the viewpoint of transmission conditions. We conclude in Section 5 with a summary of open research problems. In Appendix A we provide a Matlab program for a block LU form of an optimal Schwarz method with cross-points, and in Appendix B we give the Maple program for the two-subdomain Fourier analysis.