Search results

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 polymer injection processes has become increasingly common in mould design. In industry, such a task is accomplished mainly by using commercial packages. Owing to the complexities inherent of this class of problems, most commercial codes attempt to combine realistic rheological descriptions with simplified numerical models. In spite of the apparent success, such approaches are not able to capture important aspects of the flow topology. The present work aims to describe a more elaborate mathematical model based on finite volumes which is able to provide both accurate solutions and further insights on the physics of the polymer flow.

The mathematical model comprises the momentum and energy equations and a Poisson equation for pressure to impose the incompressibility constraint. The governing equations are discretized using the finite volume method based on central, second‐order accurate formulas for both convection and diffusion terms. Artificial dissipation terms are added externally in order to control the odd‐even decoupling problem.

The numerical model was conceived within the framework of a generalized Newtonian formulation. The capability of the numerical scheme is illustrated by simulations using three distinct constitutive relations to approach the non‐Newtonian behaviour of the polymer melt: isothermal power‐law, modified Arrhenius power‐law and cross models.

This paper extends the computational strategies previously developed to Newtonian fluids to account for more complex constitutive relations. The velocity and temperature coupled solution for polymer melts using only second‐order accurate formulas constitute also a relevant contribution.

A software package is developed for the modelling and animation of viscous incompressible fluids. The full time‐dependent Navier‐Stokes equations are used to simulate 2D and 3D incompressible fluid phenomena which include shallow and deep fluid flow, transient dynamic flow, vorticity and splashing in simulated physical environments. The package also allows the inclusion of variously shaped and spaced static or moving obstacles that are fully submerged or penetrate the fluid surface. Stable numerical analysis techniques based on finite‐differences are used for the solution of the Navier‐Stokes equations. To model free‐surface fluids, a technique based on the Marker‐and‐Cell method is presented. Based on the fluid’s pressure and velocities obtained from the solution of the Navier‐Stokes equations this technique allows modelling of the fluid’s free surface either by solving a surface equation of by tracking the motion of marker particles. The latter technique is suitable for visualization of splashing and vorticity. Furthermore, an editing tool is developed for easy definition of a physical‐world which includes obstacles, boundaries and fluid properties such as viscosity, initial velocity and pressure. Using the editor, complex fluid simulations can be performed without prior knowledge of the underlying fluid dynamics equations. Finally, depending on the application fluid rendering techniques are developed using standard Silicon Graphics workstation hardware routines.

Describes further development of a 3D finite difference code written to model turbulent flows in an open channel with a moving free surface. The code has been developed so that the computational domain can have side‐walls and/or periodic directions and that the flow may also be buoyancy driven. Either a full simulation or large eddy simulation (LES) of the turbulence can be performed. Results are presented of a simulation of periodic streamwise flow in an open channel with parallel side‐walls and also of a thermal jet into an open tank. Both simulations were carried out on a UNIX workstation using resolutions that enable the results to be viewed within an “engineering context”. The LES application demands numerical approximations which conserve mass, momentum and total energy with high precision, and which permit wave motion with very little numerical dispersion or dissipation. The free surface is tracked using a split‐merge technique which combines the volume of fluid (VOF) and height function methods in a way that is conservative.

The purpose of this paper is to develop a counter-extrapolation approach for computational heat and mass transfer with the interfacial discontinuity considered at conjugate interfaces.

By applying finite-difference approximations for the interfacial gradients along the local normal direction, the conjugate system can be simplified to the Dirichlet boundary problems for individual domains. A suitable method for the Dirichlet boundary value condition can then be used. The lattice Boltzmann method has been used to demonstrate the method. The model has been carefully validated by comparing the simulation results and theoretical solutions for steady and unsteady systems with flat or circular interfaces. Furthermore, the cooling process of a hot cylinder in a cold flow, which involves unsteady flow and heat transfer across a curved interface, has been simulated as an example to illustrate the practical usefulness of this model.

Good agreement has been observed in comparisons of simulations and theoretical solutions. The convergence and stability of the method have also been examined and satisfactory results have been obtained. Results of the cylinder cooling process show that a surface insulation layer can effectively reduce the heat transfer process and slow down the cooling process.

This method possesses several technical advantages, including the simple and straightforward algorithm, and accurate representation of the interface geometry. The basic idea and algorithm of the counter-extrapolation procedure presented here can be readily extended to other lattice Boltzmann models and even other computational technologies for heat and mass transfer systems with interface discontinuity.

This paper aims to improve the computational efficiency of higher-order accurate Noye–Hayman [NH (9,9)] implicit finite difference scheme for the solution of electromagnetic scattering problems in tunnel environments.

The proposed method consists of two major steps: First, the higher-order NH (9,9) scheme is numerically discretized using the finite-difference method. The second step is to use an algorithm based on hierarchical interpolative factorization (HIF) to accelerate the solution of this scheme.

It is observed that the simulation results obtained from the numerical tests illustrate very high accuracy of the NH (9,9) method in typical tunnel environments. HIF algorithm makes the NH (9,9) method computationally efficient for two-dimensional (2D) or three-dimensional (3D) problems. The proposed method could help in reducing the computational cost of the NH (9,9) method very close to O(n) usual O(n3) for a full matrix.

For simplicity, in this study, perfect electric conductor boundary conditions are considered. Future research may also include the utilization of meteorological techniques, including the effects of backward traveling waves, and make comparisons with the experimental data.

This study is directly applicable to typical problems in the field of tunnel propagation modeling for both national commercial and military applications.

This paper aims to present the development of a highly parallel finite-difference computational fluid dynamics code in generalized curvilinear coordinates system. The objectives are to handle internal and external flows in fairly complex geometries including shock waves, compressible turbulence and heat transfer.

The code is equipped with high-order discretization schemes to improve the computational accuracy of the solution algorithm. Besides, a new method to deal with the geometrical singularities, so-called domain decomposition method (DDM), is implemented. The DDM consists of using two different meshes communicating with each other, where the base mesh is Cartesian and the overlapped one a hollow cylinder.

The robustness of the present implemented code is appraised through several numerical test cases including a vortex advection, supersonic compressible flow over a cylinder, Poiseuille flow, turbulent channel and pipe flows. The results obtained here are in an excellent agreement when compared to the experimental data and the previous direct numerical simulation (DNS). As for the DDM strategy, it was successful as simulation time is clearly decreased and the connection between the two subdomains does not create spurious oscillations.

In sum, the developed solver was capable of solving, accurately and with high-precision, two- and three-dimensional compressible flows including fairly complex geometries. It is noted that the data provided by the DNS of supersonic pipe flows are not abundant in the literature and therefore will be available online for the community.

The paper aims to deal with the exact computation of the Jacobian of a time criteria from a numerical simulation of power electronics structures, for the sizing by gradient-based optimization algorithm.

Runge Kutta 44 is used to solve the state equations. The generic approach combines numerical and symbolic approaches. The modelling of the static converter is based on ideal switches.

The paper extends the state equations to derivate any state variable according a sizing parameter. The integral expressions used for some sizing performances (e.g. average or RMS values) mix symbolic and numerical approaches. Choices are made for the derivatives of the extrema of which the search is not a continuous process. The use of an object-oriented implementation allows to have generic formulation of some design performances.

The paper aims to propose and to test formulations of sizing criteria and their gradients; so, the modelling of the study case is carried out manually. Due to generic modelling approach used for the power electronics, the model is not completely continuous. So, the derivatives according some parameters (e.g. switch controls) must be carried out by finite differences. However, as the global behaviour is continuous, it is not critical.

The proposed formulations can be easily applied on simple static converter applications. For applications with large state equations, it should be possible to use the basic model of switches used in simulation tools of power electronics. The solving process and the sizing criteria formulation (with their derivatives) are generic and can be instantiate for any study.

The approach proposes formulations giving a numerical sizing dynamic model with a Jacobian computed, if possible, by an exact derivation useful for optimization studies. The approach gives fast simulation and fast computation of the derivatives by combining numerical and analytical approaches.

To present a numerical method for the solution of the unsteady incompressible Navier‐Stokes equations in a generic setting.

The equations are discretized in space by the finite element method, and in time by a semi‐implicit finite difference scheme, using a fractional‐step method to enforce incompressibility.

The presented results demonstrate the satisfactory accuracy of the method in the simulation of vortical flows in laminar regime and the stability of the solution in presence of a strong boundary layer.

The successful integration of the CFD into the industrial design depends on its capability to produce accurate and reliable simulations of real life applications. These considerations drive the development of the proposed method: it can be used in conjunction with finite elements of any order of accuracy, providing accurate and numerically stable results for complex flows. Moreover, the computational requirements are low when compared with other similar strategies.

This paper aims to present a novel hybrid time integration approach for efficient numerical simulations of multiscale problems involving interactions of electromagnetic fields with fine structures.

The entire computational domain is discretized with a coarse grid and a locally refined subgrid containing the tiny objects. On the coarse grid, the time integration of Maxwell’s equations is realized by the conventional finite-difference technique, while on the subgrid, the unconditionally stable Krylov-subspace-exponential method is adopted to breakthrough the Courant–Friedrichs–Lewy stability condition.

It is shown that in contrast with the conventional finite-difference time-domain method, the proposed approach significantly reduces the memory costs and computation time while providing comparative results.

An efficient hybrid time integration approach for numerical simulations of multiscale electromagnetic problems is presented.