Search results

The purpose of this paper is to improve the accuracy and stability of the existing solutions to 1D Stefan problem with time-dependent Dirichlet boundary conditions. The accuracy improvement should come with respect to both temperature distribution and moving boundary location.

The variable space grid method based on mixed finite element/finite difference approach is applied on 1D Stefan problem with time-dependent Dirichlet boundary conditions describing melting process. The authors obtain the position of the moving boundary between two phases using finite differences, whereas finite element method is used to determine temperature distribution. In each time step, the positions of finite element nodes are updated according to the moving boundary, whereas the authors map the nodal temperatures with respect to the new mesh using interpolation techniques.

The authors found that computational results obtained by proposed approach exhibit good agreement with the exact solution. Moreover, the results for temperature distribution, moving boundary location and moving boundary speed are more accurate than those obtained by variable space grid method based on pure finite differences.

The authors’ approach clearly differs from the previous solutions in terms of methodology. While pure finite difference variable space grid method produces stable solution, the mixed finite element/finite difference variable space grid scheme is significantly more accurate, especially in case of high alpha. Slightly modified scheme has a potential to be applied to 2D and 3D Stefan problems.

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.

Introduces papers from this area of expertise from the ISEF 1999 Proceedings. States the goal herein is one of identifying devices or systems able to provide prescribed performance. Notes that 18 papers from the Symposium are grouped in the area of automated optimal design. Describes the main challenges that condition computational electromagnetism’s future development. Concludes by itemizing the range of applications from small activators to optimization of induction heating systems in this third chapter.

Detailed results of numerical calculations of transient, 2D incompressible flow around and in the wake of a square prism at Re = 100, 200 and 500 are presented. An implicit finite‐difference operator‐splitting method, a version of the known SIMPLEC‐like method on a staggered grid, is described. Appropriate theoretical results are presented. The method has second‐order accuracy in space, conserving mass, momentum and kinetic energy. A new modification of the multigrid method is employed to solve the elliptic pressure problem. Calculations are performed on a sequence of spatial grids with up to 401 × 321 grid points, at sequentially halved time steps to ensure grid‐independent results. Three types of flow are shown to exist at Re = 500: a steady‐state unstable flow and two which are transient, fully periodic and asymmetric about the centre line but mirror symmetric to each other. Discrete frequency spectra of drag and lift coefficients are presented.

The purpose of this study is to develop a new method of lines for one-dimensional (1D) advection-reaction-diffusion (ADR) equations that is conservative and provides piecewise analytical solutions in space, compare it with other finite-difference discretizations and assess the effects of advection and reaction on both 1D and two-dimensional (2D) problems.

A conservative method of lines based on the piecewise analytical integration of the two-point boundary value problems that result from the local solution of the advection-diffusion operator subject to the continuity of the dependent variables and their fluxes at the control volume boundaries is presented. The method results in nonlinear first-order, ordinary differential equations in time for the nodal values of the dependent variables at three adjacent grid points and triangular mass and source matrices, reduces to the well-known exponentially fitted techniques for constant coefficients and equally spaced grids and provides continuous solutions in space.

The conservative method of lines presented here results in three-point finite difference equations for the nodal values, implicitly treats the advection and diffusion terms and is unconditionally stable if the reaction terms are implicitly treated. The method is shown to be more accurate than other three-point, exponentially fitted methods for nonlinear problems with interior and/or boundary layers and/or source/reaction terms. The effects of linear advection in 1D reacting flow problems indicates that the wave front steepens as it approaches the downstream boundary, whereas its back corresponds to a translation of the initial conditions; for nonlinear advection, the wave front exhibits steepening but the wave back shows a linear dependence on space. For a system of two nonlinearly coupled, 2D ADR equations, it is shown that a counter-clockwise rotating vortical field stretches the spiral whose tip drifts about the center of the domain, whereas a clock-wise rotating one compresses the wave and thickens its arms.

A new, conservative method of lines that implicitly treats the advection and diffusion terms and provides piecewise-exponential solutions in space is presented and applied to some 1D and 2D advection reactions.

The purpose of this paper is to develop a new finite difference method for solving the seismic wave propagation in fluid-solid media, which can be described by the acoustic and viscoelastic wave equations for the fluid and solid parts, respectively.

In this paper, the authors introduced a coordinate transformation method for seismic wave simulation method. In the new method, the irregular fluid–solid interface is transformed into a horizontal interface. Then, a multi-block coordinate transformation method is proposed to mesh every layer to curved grids and transforms every interface to horizontal interface. Meanwhile, a variable grid size is used in different regions according to the shape and the velocity within each region. Finally, a Lebedev-standard staggered coupled grid scheme for curved grids is applied in the multi-block coordinate transformation method to reduce the computational cost.

The instability in the auxiliary coordinate system caused by the standard staggered grid scheme is resolved using a curved grid viscoelastic wave field separation strategy. Several numerical examples are solved using this new method. It has been shown that the new method is stable, efficient and highly accurate in solving the seismic wave equation defined on domain with irregular fluid–solid interface.

First, the irregular fluid–solid interface is transformed into a horizontal interface by using the coordinate transformation method. The conversion between pressures and stresses is easy to implement and adaptive to different irregular fluid–solid interface models, because the normal stress and shear stress vanish when the normal angle is 90° in the interface. Moreover, in the new method, the strong false artificial boundary reflection and instability caused by ladder-shaped grid discretion are resolved as well.

We describe a novel mathematical approach to deriving and solving covolume models of the incompressible 2‐D Navier‐Stokes flow equations. The approach integrates three technical components into a single modelling algorithm: 1. Automatic Grid Generation. An algorithm is described and used to automatically discretize the flow domain into a Delaunay triangulation and a dual Voronoi polygonal tessellation. 2. Covolume Finite Difference Equation Generation. Three covolume discretizations of the Navier‐Stokes equations are presented. The first scheme conserves mass over triangular control volumes, the second scheme over polygonal control volumes and the third scheme conserves mass over both. Simple consistent finite difference equations are derived in terms of the primitive variables of velocity and pressure. 3. Dual Variable Reduction. A network theoretic technique is used to transform each of the finite difference systems into equivalent systems which are considerably smaller than the original primitive finite difference system.

The purpose of this paper is to propose an accurate and efficient technique for computing flow sensitivities by finite differences of perturbed flow fields. It relies on computing the perturbed flows on coarser grid levels only: to achieve the same fine-grid accuracy, the approximate value of the relative local truncation error between coarser and finest grids unperturbed flow fields, provided by a standard multigrid method, is added to the coarse grid equations. The gradient computation is introduced in a hybrid genetic algorithm (HGA) that takes advantage of the presented method to accelerate the gradient-based search. An application to a classical transonic airfoil design is reported.

Genetic optimization algorithm hybridized with classical gradient-based search techniques; usage of fast and accurate gradient computation technique.

The new variant of the prolongation operator with weighting terms based on the volume of grid cells improves the accuracy of the MAFD method for turbulent viscous flows. The hybrid GA is capable to efficiently handle and compensate for the error that, although very limited, is present in the multigrid-aided finite-difference (MAFD) gradient evaluation method.

The proposed new variants of HGA, while outperforming the simple genetic algorithm, still require tuning and validation to further improve performance.

Significant speedup of CFD-based optimization loops.

Introduction of new multigrid prolongation operator that improves the accuracy of MAFD method for turbulent viscous flows. First application of MAFD evaluation of flow sensitivities within a hybrid optimization framework.

The purpose of this paper is to both determine the effects of the nonlinearity on the wave dynamics and assess the temporal and spatial accuracy of five finite difference methods for the solution of the inviscid generalized regularized long-wave (GRLW) equation subject to initial Gaussian conditions.

Two implicit second- and fourth-order accurate finite difference methods and three Runge-Kutta procedures are introduced. The methods employ a new dependent variable which contains the wave amplitude and its second-order spatial derivative. Numerical experiments are reported for several temporal and spatial step sizes in order to assess their accuracy and the preservation of the first two invariants of the inviscid GRLW equation as functions of the spatial and temporal orders of accuracy, and thus determine the conditions under which grid-independent results are obtained.

It has been found that the steepening of the wave increase as the nonlinearity exponent is increased and that the accuracy of the fourth-order Runge-Kutta method is comparable to that of a second-order implicit procedure for time steps smaller than 100th, and that only the fourth-order compact method is almost grid-independent if the time step is on the order of 1,000th and more than 5,000 grid points are used, because of the initial steepening of the initial profile, wave breakup and solitary wave propagation.

This is the first study where an accuracy assessment of wave breakup of the inviscid GRLW equation subject to initial Gaussian conditions is reported.

The overall pressure drop in an axisymmetric contraction is minimised using two different grid sizes. The transition region was parameterised with only two design variables to make it possible to create surface plots of the objective function in the design space, which were based on 121 CFD calculations for each grid. The coarse grid showed to have significant numerical noise in the objective function while the finer grid had less numerical noise. The optimisation was performed with two methods, a Response Surface Model (RSM) and a gradient‐based method (the Method of Feasible Directions) to study the influence from numerical noise. Both optimisation methods were able to find the global optimum with the two different grid sizes (the search path for the gradient‐based method on the coarse grid was able to avoid the region in the design space containing local minima). However, the RSM needed fewer iterations in reaching the optimum. From a grid convergence study at two points in the design space the level of noise appeared to be sufficiently low, when the relative step size is 10^–4 for the finite difference calculations, to not influence the convergence if the errors are below 5 per cent for this contraction geometry.