Search results

This paper aims to solve inhomogeneous dielectric problems by matching boundary conditions at the interfaces among homogeneous subdomains. The capabilities of Hilbert transform computations are deeply investigated in the case of limited numbers of samples, and a refined model is presented by means of investigating accuracies in a case study with three subdomains.

The accuracies, refined by Richardson extrapolation to zero error, are compared to finite element (FEM) and finite difference methods. The boundary matching procedures can be easily applied to the results of a previous Schwarz–Christoffel (SC) conformal mapping stage in SC + BC procedures, to cope with field singularities or with open boundary problems.

The proposed field computations are of general interest both for electrostatic and magnetostatic field analysis and optimization. They can be useful as comparison tools for FEM results or when severe field singularities can impair the accuracies of other methods.

This static field methodology, of course, can be used to analyse transverse electro magnetic (TEM) or quasi-TEM propagation modes. It is possible that, in some case, these may make a contribution to the analysis of axis symmetrical problems.

The most relevant result is the possible introduction of SC + BC computations as a standard tool for solving inhomogeneous dielectric field problems.

This paper aims to analyze the pressure distribution of rectangular aerostatic thrust bearing with a single air supply inlet using the complex potential theory and conformal mapping.

The Möbius transform is used to map the interior of a rectangle onto the interior of a unit circle, from which the pressure distribution and load carrying capacity are obtained. The calculation results are verified by finite difference method.

The constructed Möbius formula is very effective for the performance characteristics researches for the rectangular thrust bearing with a single air supply inlet. In addition, it is also noted that to obtain the optimized load carrying capacity, the square thrust bearing can be adopted.

The Möbius transform is found suitable to describe the pressure distribution of the rectangular thrust bearing with a single air supply inlet.

This study aims to simulate the dendritic growth in Stokes flow by iteratively coupling a domain and boundary type meshless method.

A preconditioned phase-field model for dendritic solidification of a pure supercooled melt is solved by the strong-form space-time adaptive approach based on dynamic quadtree domain decomposition. The domain-type space discretisation relies on monomial augmented polyharmonic splines interpolation. The forward Euler scheme is used for time evolution. The boundary-type meshless method solves the Stokes flow around the dendrite based on the collocation of the moving and fixed flow boundaries with the regularised Stokes flow fundamental solution. Both approaches are iteratively coupled at the moving solid–liquid interface. The solution procedure ensures computationally efficient and accurate calculations. The novel approach is numerically implemented for a 2D case.

The solution procedure reflects the advantages of both meshless methods. Domain one is not sensitive to the dendrite orientation and boundary one reduces the dimensionality of the flow field solution. The procedure results agree well with the reference results obtained by the classical numerical methods. Directions for selecting the appropriate free parameters which yield the highest accuracy and computational efficiency are presented.

A combination of boundary- and domain-type meshless methods is used to simulate dendritic solidification with the influence of fluid flow efficiently.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

The rostrum of a paddlefish provides hydrodynamic stability during feeding process in addition to detect the food using receptors that are randomly distributed in the rostrum. The exterior tissue of the rostrum covers the cartilage that surrounds the bones forming interlocking star shaped bones.

The aim of this work is to assess the mechanical behavior of four finite element models varying the type of formulation as follows: linear-reduced integration, linear-full integration, quadratic-reduced integration and quadratic-full integration. The paper also presents the load transfer mechanisms of the bone structure of the rostrum. The base material used in the study was steel with elastic–plastic behavior as a homogeneous material before applying materials properties that represents the behavior of bones, cartilages and tissues.

Conclusions are based on comparison among the four models. There is no significant difference between integration orders for similar type of elements. Quadratic-reduced integration formulation resulted in lower structural stiffness compared with linear formulation as seen by higher displacements and stresses than using linearly formulated elements. It is concluded that second-order elements with reduced integration are the alternative to analyze biological structures as they can better adapt to the complex natural contours and can model accurately stress concentrations and distributions without over stiffening their general response.

The use of advanced computational mechanics techniques to analyze the complex geometry and components of the paddlefish rostrum provides a viable avenue to gain fundamental understanding of the proper finite element formulation needed to successfully obtain the system behavior and hot spot locations.

This paper aims to consider a multiscale electromagnetic wave problem for a housing with a ventilation grill. Using the standard finite element method to discretise the apertures leads to an unduly large number of unknowns. An efficient approach to simulate the multiple scales is introduced. The aim is to significantly reduce the computational costs.

A domain decomposition technique with upscaling is applied to cope with the different scales. The idea is to split the domain of computation into an exterior domain and multiple non-overlapping sub-domains. Each sub-domain represents a single aperture and uses the same finite element mesh. The identical mesh of the sub-domains is efficiently exploited by the hybrid discontinuous Galerkin method and a Schur complement which facilitates the transition from fine meshes in the sub-domains to a coarse mesh in the exterior domain. A coarse skeleton grid is used on the interface between the exterior domain and the individual sub-domains to avoid large dense blocks in the finite element discretisation matrix.

Applying a Schur complement to the identical discretisation of the sub-domains leads to a method that scales very well with respect to the number of apertures.

The error compared to the standard finite element method is negligible and the computational costs are significantly reduced.

We consider, as discretization in space, the nonconforming mesh developed in SUSHI (Scheme Using Stabilization and Hybrid Interfaces) developed in Eymard et al. (2010) for a semi-linear heat equation. The time discretization is performed using a uniform mesh. We are concerned with a nonlinear scheme that has been studied in Bradji (2016) in the context of the general framework GDM (Gradient Discretization Method) (Droniou et al., 2018) which includes SUSHI. We provide sufficient conditions on the size of the spatial mesh and the time step which allow to prove a W1,∞(L2)-error estimate. This error estimate can be viewed as an improvement for the W1,2(L2)-error estimate proved in Bradji (2016). The W1,∞(L2)-error estimate we want to prove in this note was stated without proof in Bradji (2016, Remark 7.2, Page 1302). Its proof is based on a comparison with an appropriately chosen auxiliary finite volume scheme along with the derivation of some new estimates on its solution.

This paper shows a new methodology for evaluating the value and sensitivity of autocall knock-in type equity-linked securities. While the existing evaluation methods, Monte Carlo simulation and finite difference method, have limitations in underestimating the knock-in effect, which is one of the important characteristics of this type, this paper presents a precise joint probability formula for multiple autocall chances and knock-in events. Based on this, the calculation results obtained by utilizing numerical and Monte Carlo integration are presented and compared with those of existing models. The results of the proposed model show notable improvements in terms of accuracy and calculation time.

The goal of this work is to create a computational finite element model to perform thermoelastic stress analysis (TSA) with the usage of a non-ideal load frequency, containing the effects of the material thermal properties.

Throughout this document, the methodology of the model is presented first, followed by the procedure and results. The last part is reserved to results, discussion and conclusions.

This work had the main goal to create a model to perform TSA with the usage of non-ideal loading frequencies, considering the materials’ thermal properties. Loading frequencies out of the ideal range were applied and the model showed capable of good results. The created model reproduced acceptably the TSA, with the desired conditions.

This work creates a model to perform TSA with the usage of non-ideal loading frequencies, considering the materials’ thermal properties.

The purpose of this study is to develop stable, convergent and accurate numerical method for solving singularly perturbed differential equations having both small and large delay.

This study introduces a fitted nonpolynomial spline method for singularly perturbed differential equations having both small and large delay. The numerical scheme is developed on uniform mesh using fitted operator in the given differential equation.

The stability of the developed numerical method is established and its uniform convergence is proved. To validate the applicability of the method, one model problem is considered for numerical experimentation for different values of the perturbation parameter and mesh points.

In this paper, the authors consider a new governing problem having both small delay on convection term and large delay. As far as the researchers' knowledge is considered numerical solution of singularly perturbed boundary value problem containing both small delay and large delay is first being considered.