Search results

Non‐linear reaction‐diffusion processes with cross‐diffusion in two‐dimensional, anisotropic media are analyzed by means of an implicit, iterative, time‐linearized approximate factorization technique as functions of the anisotropy of the heat and species diffusivity tensors, the Soret and Dufour cross‐diffusion effects, and five types of boundary conditions. It is shown that anisotropy and cross‐diffusion deform the reaction front and affect the front velocity, and the magnitude of these effects increases as the magnitude of the off‐diagonal components of the heat and species diffusivity tensors is increased. It is also shown that the five types of boundary conditions employed in this study produce similar results except when there is either strong anisotropy in the species or heat diffusivity tensors and there are no Soret and Dufour effects, or the species and heat diffusivity tensors are isotropic, but the anisotropy of the Soret and Dufour effects is important. If the species and heat diffusivity tensors are isotropic, the effects of either the Soret or the Dufour cross‐diffusion effects are small for the cases considered in this study. The time required to achieve steady state depends on the anisotropy of the heat and diffusivity tensors, the cross‐diffusion effects, and the boundary conditions.

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.

A method has been developed for the solution of inverse heat diffusion problems to find the initial condition, boundary condition, and the source and sink function in the heat diffusion equation. The method has been used in the development of a source‐and‐sink method to find the boundary conditions in inverse Stefan problems. Green's functions have been used in the solution, and the problems are solved by using two approaches: a series solution approach, and a time incremental approach. Both can be used to find the boundary conditions without reliance on the flux information to be supplied at both sides of the interface. The methods are efficient in that they require less equations to be solved for the conditions. The numerical results have shown to be accurate, convergent, and stable. Most of all, the results do not degrade with time as in other time marching schemes reported in the literature. Algorithms can also be easily developed for the solution of the conditions.

The numerical solution of problems involving two‐dimensional flow in an infinite or a semi‐infinite channel is considered. Beyond a certain finite region, where the flow and geometry may be general, a “tail” region is assumed where the flow is potential and the channel is uniform. This situation is typical in many cases of fluid‐structure interaction and flow around obstacles in a channel. The unbounded domain is truncated by means of an artificial boundary B, which separates between the finite computational domain and the “tail.” On B, special boundary conditions are devised. In the finite computational domain, the problem is solved using a finite element scheme. Both non‐local and local artificial boundary conditions are considered on B, and their performance is compared via numerical examples.

The purpose of this work is to apply the Fourier–Ritz method to study the vibration behavior of the moderately thick functionally graded (FG) parabolic and circular panels and shells of revolution with general boundary conditions.

The modified Fourier series is chosen as the basis function of the admissible functions of the structure to eliminate all the relevant discontinuities of the displacements and their derivatives at the edges, and the vibration behavior is solved by means of the Ritz method. The complete shells of revolution can be achieved by using the coupling spring technique to imitate the kinematic compatibility and physical compatibility conditions of FG parabolic and circular panels at the common meridian of θ = 0 and 2π. The convergence and accuracy of the present method are verified by other contributors.

Some new results of FG panels and shells with elastic restraints, as well as different geometric and material parameters, are presented and the effects of the elastic restraint parameters, power-law exponent, circumference angle and power-law distributions on the free vibration characteristic of the panels are also presented, which can be served as benchmark data for the designers and engineers to avoid the unpleasant, inefficient and structurally damaging resonant.

The paper could provide the reference for the research about the moderately thick FG parabolic and circular panels and shells of revolution with general boundary conditions. In addition, the change of the boundary conditions can be easily achieved by just varying the stiffness of the boundary restraining springs along all the edges of panels without making any changes in the solution procedure.

This paper aims to present a novel scheme for imposing periodic boundary conditions with downscaled macroscopic strain measures of gradient Cosserat continuum on the representative volume element (RVE) of discrete particle assembly in the frame of the second-order computational homogenization methods for granular materials.

The proposed scheme is based on the generalized Hill’s lemma of gradient Cosserat continuum and the incremental non-linear constitutive relation condensed to the peripheral particles of the RVE of discrete particle assembly. The generalized Hill’s lemma conducts to downscale the macroscopic strain or stress measures and to impose the periodic boundary conditions on the RVE boundary so that the Hill-Mandel energy equivalence condition is ensured. Because of the incremental non-linear constitutive relation condensed to the peripheral particles of the RVE, the periodic boundary displacement and traction constraints together with the downscaled macroscopic strains and strain gradients, micro-rotations and curvatures are imposed in the point-wise sense without the need of introducing the Lagrange multipliers for enforcing the periodic boundary displacement and traction constraints in a weak sense.

Numerical results demonstrate that the applicability and effectiveness of the proposed scheme in imposing the periodic boundary conditions on the RVE. The results of the RVE subjected to the periodic boundary conditions together with the displacement boundary conditions in the second-order computational homogenization for granular materials provide the desired estimations, which lie between the upper and the lower bounds provided by the displacement and the traction boundary conditions imposed on the RVE respectively.

Each grain in the particulate system under consideration is assumed to be rigid and circular.

The proposed scheme for imposing periodic boundary conditions on the RVE can be adopted solely for estimating the effective mechanical properties of granular materials and/or integrated into the frame of the second-order computational homogenization method with a nested finite element method-discrete element method solution procedure for granular materials. It will tend to provide, at least theoretically, more reasonable results for effective material properties and solutions of a macroscopic boundary value problem simulated by the computational homogenization method.

This paper presents a novel scheme for imposing periodic boundary conditions with downscaled macroscopic strain measures of gradient Cosserat continuum on the RVE of discrete particle assembly for granular materials without need of introducing Lagrange multipliers for enforcing periodic boundary conditions in a weak (integration) sense.

The purpose of this paper is to extend the penalty concept to treat partial slip, free surface, contact and related boundary conditions in viscous flow simulation.

The penalty partial‐slip formulation is analysed and related to the classical Navier slip condition. The same penalty scheme also allows partial penetration through a boundary, hence the implementation of porous wall boundaries. The finite element method is used for investigating and interpreting penalty approaches to boundary conditions.

The generalised penalty approach is verified by means of a novel variant of the circular‐Couette flow problem, having partial slip on one of the cylindrical boundaries, for which an analytic solution is derived. Further verificationis provided by consideration of viscous flow over a sphere with partial slip on the surface, and comparison of numerical and classical solutions. Numerical studies illustrate the versatility of the approach.

The penalty approach is applied to some different boundaries: partial slip and partial penetration with no/full slip/penetration as limiting cases; free surface; space‐ and time‐varying boundary conditions which allow progressive contact over time. Application is made to curved and inclined boundaries. Sensitivity of flow to penalty parameters is an avenue for continued research, as is application of the penalty approach for non‐Newtonian flows.

This is the first work to show the relation between penalty formulation of boundary conditions and physical boundary conditions. It provides a method that overcomes past difficulties in implementing partial slip on boundaries of general shape, and which handles progressive contact. It also provides useful benchmark problems for future studies.

The boundary integral method (BIM) provides unparalleled computational efficiency for solving problems wherever it is applicable. For Stokes flows, the BIM in its current form can only be applied to a limited class of problems that generally comprises boundaries with either specified velocity or stress. This study aims to radically extend the applicability by developing a general method within the BIM framework that can handle periodic, symmetry, zero normal-velocity gradient and the specified pressure boundary conditions. This study is limited in scope to steady-state flows.

The proposed method introduces a set of points near the boundary for the symmetry, zero normal-velocity gradient and specified pressure boundary conditions. The formulation for the first two boundary conditions use a spatial discretization procedure within the BIM framework to arrive at a set of equations for the unknowns. The specified pressure boundary condition warrants the decomposition of the unknown traction term into simpler components before the discretization procedure can be executed. Though the new methodology is illustrated in detail for two-dimensional rectangular domains, it can be generalized to more complex three-dimensional cases. This will be the subject for future investigations.

The current endeavor has successfully demonstrated the incorporation of the above boundary conditions through simple Stokes flow problems like plane channel flow, flow through ribbed duct and plane wall jet. The predicted results matched adequately with either analytical solutions or with available literature data.

To the best of the author’s knowledge, this is the first time that the exit boundary conditions like zero normal-velocity gradient and specified pressure have been formulated within the BIM for Stokes flows. These boundary conditions are extremely powerful and the current research initiative has the potential to dramatically increase the range of applicability of the BIM for Stokes flow simulations.

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

Numerical simulation of a fluid flow involves the specification of boundary conditions along all or part of the boundary. Designs a means of handling outflow boundary conditions for the incompressible Navier‐Stokes equations. Addresses through‐flow problems involving the specification of outflow conditions at the synthetic boundary. This outflow boundary condition is applicable to a developing flow problem. The underlying objectives behind designing the boundary condition at the truncated boundary are three‐fold, namely: matching with Navier‐Stokes equations inside the domain; taking both non‐linear and diffusive contributions into account; and ensuring the discrete divergence‐free condition. In order to meet these requirements, follows the concept of a free boundary condition by taking the outflow nodal values of u, v and p as unknowns, which are coupled with the interior unknowns through the surface integrals in the momentum equations. The computed solutions can be legitimately regarded as solutions to conservation equations under consideration when both components of the surface traction vector approach zero. With the convergent property accommodated in the present mixed finite element analysis, the task remains to simply improve the accuracy. Demonstrates the capability of the proposed non‐linear outflow boundary conditions through several benchmark tests.