Search

Analysis of the slipstream development around the high-speed trains in tunnels would provide references for assessing the transient gust loads on trackside workers and trackside furniture in tunnels. This paper focuses on the computational analysis of the slipstream caused by high-speed trains passing through double-track tunnels with a cross-sectional area of 100 m2. Three-dimensional unsteady compressible Reynolds-averaged Navier-Stokes equations and a realizable k-ε turbulence model were used to describe the airflow characteristics around a high-speed train in the tunnel. The moving boundary problem was treated using the sliding mesh technology. Three cases were simulated in this paper, including two tunnel lengths and two different configurations of the train. The train speed in these three cases was 250 km/h. The accuracy of the numerical method was validated by the experimental data from full-scale tests, and reasonable consistency was obtained. The results show that the flow field around the high-speed trains can be divided into three distinct regions: the region in front of the train nose, the annular region and the wake region. The slipstream development along the two sides of train is not in balance and offsets to the narrow side in the double-track tunnels. Due to the piston effect, the slipstream has a larger peak value in the tunnel than in open air. The tunnel length, train length and length ratio affect the slipstream velocities; in particular, the velocities increase with longer trains. Moreover, the propagation of pressure waves also induces the slipstream fluctuations: substantial velocity fluctuations mainly occur in front of the train, and weaken with the decrease in amplitude of the pressure wave. [ABSTRACT FROM AUTHOR]

An ordinary differential equation (ODE) can be split into simpler sub equations and each of the sub equations is solved subsequently by a numerical method. Such a procedure involves splitting error and numerical error caused by the time stepping methods applied to sub equations. The aim of the paper is to present an integral formula for the global error expansion of a splitting procedure combined with any numerical ODE solver. [ABSTRACT FROM AUTHOR]

In this paper, the authors study the forward and inverse problems for a fractional boundary value problem with Dirichlet boundary conditions. The existence and uniqueness of solutions for the forward problem is first proved. Then an inverse source problem is considered. [ABSTRACT FROM AUTHOR]

This paper is devoted to the analytical study of initial–boundary value problems for a system of hyperbolic balance laws derived from a repulsive chemotaxis model with logarithmic sensitivity. In the first part of the paper we show that, subject to the Dirichlet boundary conditions, classical solutions exist globally in time for large initial data. Asymptotically in time, the solutions are shown to converge to their boundary data at an exponential rate as time goes to infinity. Numerical simulations are supplied to corroborate the analytical results. The analytic approach developed herein can be utilized to handle a family of initial and initial–boundary value problems of the model and related models with similar mathematical structure. As a demonstration of the effectiveness of our approach, in the second part of the paper we show that, subject to the Neumann–Dirichlet boundary conditions, classical solutions exist and converge to constant equilibrium states for large initial data and for arbitrary values of the chemical diffusion coefficient. This improves a previous result obtained in [30] where the smallness of the chemical diffusion coefficient was required. [ABSTRACT FROM AUTHOR]

In this paper, the appearing probability of truck is introduced and an extended car-following model is presented to analyze the traffic flow based on the consideration of driver's characteristics, under honk environment. The stability condition of this proposed model is obtained through linear stability analysis. In order to study the evolution properties of traffic wave near the critical point, the mKdV equation is derived by the reductive perturbation method. The results show that the traffic flow will become more disorder for the larger appearing probability of truck. Besides, the appearance of leading truck affects not only the stability of traffic flow, but also the effect of other aspects on traffic flow, such as: driver's reaction and honk effect. The effects of them on traffic flow are closely correlated with the appearing probability of truck. Finally, the numerical simulations under the periodic boundary condition are carried out to verify the proposed model. And they are consistent with the theoretical findings. [ABSTRACT FROM AUTHOR]

The paper presents the Finite Element formulation of the equations proposed by Tsiampousi et al. (2016) for coupled consolidation in unsaturated soils. Their coupling is discussed in relation to a conceptual model which divides soil behaviour into zones ranging from fully saturated to dry states. The numerical simulation of a laboratory experiment involving drainage of water from a vertical column of sand is used to validate the equations. Finally, the example of rainfall infiltration into a cut slope highlights how aspects of the conceptual model are reflected in the numerical analysis of boundary value problems involving unsaturated soils. [ABSTRACT FROM AUTHOR]

This paper is a review of the use of particulate discrete element modeling (DEM) in geomechanics. The overall objective of the paper is to serve as an introduction to researchers and practitioners in geomechanics who are considering adopting DEM in their work or using the results of DEM simulations to guide other studies, for example, the development of constitutive models for continuum-based numerical analysis. It is hoped that prior converts to the use of DEM will also benefit from a relatively objective overview of current DEM use in geomechanics. The introductory sections present the background to the method and give an overview of the evolution of the use of particulate DEM in recent geotechnical research. The general principals of the algorithm are then presented, considering the types of particles typically used, the calculation of contact forces, and formulation of simulation boundary conditions. Some techniques available to interpret and postprocess of DEM results and provide the information to link the particle scale and overall response are outlined. Approaches used to validate and calibrate DEM models to verify that DEM simulation results representative of physical reality are discussed. An overview of the application of DEM modeling to field-scale problems is then presented. Finally the conclusions consider future developments in the area and emphasize the need to maintain quality in DEM simulations. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, numerical results of a drop weight impact test examining the dynamic response of a clamped steel beams struck transversely at the centre by a mass with a rectangular indenter are presented. The numerical simulation showed that the results are very sensitive to the way in which the supports are modelled and thus this paper studies the modelling of the structural supports used to clamp the specimen beam in the experiment. The first approach simplifies the structural support as boundary condition on the nodes; the second model represents the entire support as the real approximation; the last approach uses simplified support plate model acting with the beam to simulate the experimental set up. The complete force–displacement curve is used to compare the plastic deformation of the impacted beams, thus the best approximation is selected for further calculations. The finite element model was performed using the LS-DYNA non-linear, dynamic finite element software. The numerical calculations can predict accurately the response of deflections, forces and absorbed energies, when proper boundary conditions are applied, in this case with shell elements representing the support plates. [Copyright &y& Elsevier]

Abstract: Taking correctly into account long ranged interactions in molecular simulations of molecular models based on classical atomistic representations are essential to obtain reliable results on model systems and in simulations of biological systems. A lot of numerical methods have been developed to this end; the most important of them are reviewed in this paper. Particular attention is paid to the analytical relations between the methods, this allows comparisons on efficiency and accuracy between them and also to proceed to precise implementations of these techniques. While most of the methods have been developed for Coulomb interactions, we give also some analytical details to apply these methods to screened Coulomb (Yukawa interactions), inverse power law and dipolar interactions. Particular types of systems, the quasi-2D systems, are also considered in this paper. Quasi-2D systems represent a large class of physical systems where the spatial extension in one direction of the space is very small by comparison to the other two. These representations are very useful to describe the properties of interfaces, surfaces, fluids confined in slab geometry, etc. In computer simulations, these systems are studied with partial periodic boundary conditions: periodic boundary conditions are taken in directions where spatial extensions are large and some other boundary conditions are taken in directions with smaller extensions. In this review, we describe also the numerical methods developed to handle long ranged interactions in numerical simulations of quasi-2D systems. The properties of quasi-2D systems depend strongly on interactions between components; more specifically electrostatic and magnetic interactions and interactions with external fields are of particular interest in these systems. [Copyright &y& Elsevier]

In this paper, we study a one-dimensional morphogenesis model considered by C. Stinner et al. (Math. Meth. Appl. Sci. 2012; 35: 445-465). Under homogeneous boundary conditions, we prove the existence of nonconstant positive steady states through local bifurcation theories. Then we rigorously study the stability of these nonconstant solutions when the sensitivity functions are chosen to be linear and logarithmic, respectively. Finally, we present numerical solutions to illustrate the formation of stable inhomogeneous spatial patterns. Our numerical simulations show that this model can develop very complicated and interesting structures even over one-dimensional finite domains. [ABSTRACT FROM AUTHOR]

This paper is concerned with the observer design for one-dimensional linear parabolic partial differential equations whose output is a weighted spatial average of the state over the entire spatial domain. We focus on the backstepping approach, which provides a systematic procedure to design an observer gain for systems with boundary measurement. If the output is not a boundary value of the state, the backstepping approach is not directly applicable to obtaining an observer gain that stabilizes the error dynamics. Therefore, we attempt to convert the error system into another system to which backstepping is applicable. The conversion is successfully achieved for a class of weighting functions, and the resultant observer realizes exponential convergence of the estimation error with an arbitrary decay rate in terms of the L 2 norm. In addition, an explicit expression of the observer gain is available in a special case. The effectiveness of the proposed observer is also confirmed by numerical simulations. [ABSTRACT FROM AUTHOR]

This paper presents a two-dimensional numerical investigation of deformation and motion of erythrocytes in stenotic microvessels using the immersed boundary-fictitious domain method. The erythrocytes were modeled as biconcave-shaped closed membranes filled with cytoplasm. We studied the biophysical characteristics of human erythrocytes traversing constricted microchannels with the narrowest cross-sectional diameter as small as 3 μm. The effects of essential parameters, namely, stenosis severity, shape of the erythrocytes, and erythrocyte membrane stiffness, were simulated and analyzed in this study. Moreover, simulations were performed to discuss conditions associated with the shape transitions of the cells along with the relative effects of radial position and initial orientation of erythrocytes, membrane stiffness, and plasma environments. The simulation results were compared with existing experiment findings whenever possible, and the physical insights obtained were discussed. The proposed model successfully simulated rheological behaviors of erythrocytes in microscale flow and thus is applicable to a large class of problems involving fluid flow with complex geometry and fluid-cell interactions. Our study would be helpful for further understanding of pathology of malaria and some other blood disorders. [ABSTRACT FROM AUTHOR]

Helical self-organized reversed-field pinch (RFP) regimes emerge both numerically--in 3D visco-resistive magnetohydrodynamic (MHD) simulations--and experimentally, as in the RFX-mod device at high current (IP above 1 MA). These states, called quasi-single helicity (QSH) states, are characterized by the action of a MHD mode that impresses a quasi-helical symmetry to the system, thus allowing a high degree of magnetic chaos healing. This is in contrast with the multiple helicity (MH) states, where magnetic fluctuations create a chaotic magnetic field degrading the confinement properties of the RFP. This paper reports an extensive numerical study performed in the frame of 3D visco-resistive MHD which considers the effect of helical magnetic boundary conditions, i.e. of a finite value of the radial magnetic field at the edge (magnetic perturbation, MP). We show that the system can be driven to a selected QSH state starting from both spontaneous QSH and MH regimes. In particular, a high enough MP can force a QSH helical self-organization with a helicity different from the spontaneous one. Moreover, MH states can be turned into QSH states with a selected helicity. A threshold in the amplitude of MP is observed above which is able to influence the system. Analysis of the magnetic topology of these simulations indicates that the dominant helical mode is able to temporarily sustain conserved magnetic structures in the core of the plasma. The region occupied by conserved magnetic surfaces increases reducing secondary modes' amplitude to experimental-like values. [ABSTRACT FROM AUTHOR]

This paper presents a hybrid model for opinion formation in a large group of agents exposed to the persuasive action of a small number of strong opinion leaders. The model is defined by coupling a finite difference equation for the dynamics of leaders opinion with a continuous integro-differential equation for the dynamics of the others. Such a definition stems from the idea that the leaders are few and tend to retain original opinions, so that their dynamics occur on a longer time scale with respect to the one of the other agents. A general well-posedness result is established for the initial value problem linked to the model. The asymptotic behavior in time of the related solution is characterized for some general parameter settings, which mimic distinct social scenarios, where different emerging behaviors can be observed. Analytical results are illustrated and extended through numerical simulations. [ABSTRACT FROM AUTHOR]

Abstract: The lattice Boltzmann method (LBM) for fluid flow has already proven itself a viable alternative to classical CFD (computational fluid dynamics) methods based on discrete scheme application to Navier Stokes equations governing the fluid flow. In this paper, we aim to verify the applicability of the LBM in a special case of a Couette–Taylor photobioreactor—a device comprised of two coaxial cylinders with a rotating inner wall. An appropriate numerical approach to the curved moving wall boundary condition is presented and verified by velocity flow field comparison with analytical solution. The parallelism potential of the LBM is exploited on the parallel platform of CUDA (Compute Unified Device Architecture). Microalgae growth in the flow field simulated by the LBM is based on the photosynthetic factory (PSF) model treated in a parallel stochastic manner and implemented also on the parallel platform of CUDA. The parallel stochastic PSF solver has been validated by comparison with an analytical solution of the PSF model at constant irradiance. [Copyright &y& Elsevier]

Abstract: Many novel reactor concepts based on the integration of membranes in fluidized beds have been proposed. For the detailed numerical investigation of such systems, a method is required describing the motion of the particles while accounting for (i) the gas motion and (ii) the encounters between particles and confining walls/membranes. In this paper a novel hybrid computational method is presented, combining the Discrete Particle Model (DPM) and the Immersed Boundary Method (IBM). The first part focuses on the implementation and validation of these models. The IBM component of the model has been verified by comparing the computed drag coefficients for Simple Cubic (SC) arrays of cylinders with data of Sangani and Acrivos (1982) at creeping flow conditions. For higher Reynolds numbers the code has been compared to Ferziger and Peric (2002) for a single off-center cylinder, after which simulations with addition and extraction of gas through the membrane surface have been conducted. In the second part of this work we will focus on actual fluidized bed membrane systems, and will report results of simulations for several membrane geometries as well as permeation ratios. [Copyright &y& Elsevier]

Abstract: The advent of isogeometric analysis (IGA) using the same basis functions for design and analysis constitutes a milestone in the unification of geometric modeling and numerical simulation. However, an important class of geometric models based on the CSG (Constructive Solid Geometry) concept such as trimmed NURBS surfaces do not fully support the isogeometric paradigm, since basis functions do not explicitly represent the boundary. The finite cell method (FCM) is a high-order fictitious domain method, which offers simple meshing of potentially complex domains into a structured grid of cuboid cells, while still achieving exponential rates of convergence for smooth problems. In the present paper, we first discuss the possibility to directly couple the finite cell method to CSG, without any necessity for meshing the three-dimensional domain, and then explore a combination of the best of the two approaches IGA and FCM, closely following ideas of the recently introduced shell FCM. The resulting finite cell extension to isogeometric analysis achieves a truly straightforward transfer of a trimmed NURBS surface into an analysis suitable NURBS basis, while benefiting from the favorable properties of the high-order and high-continuity basis functions. Accuracy and efficiency of the new approach are demonstrated by a numerical benchmark, and its versatility is outlined by the analysis of different trimmed design variants of a brake disk. [Copyright &y& Elsevier]

In this paper, we apply a regularized sinc method to compute the eigenvalues of a discontinuous regular Dirac system with transmission conditions at the point of discontinuity. The regularized technique allows us to insert some parameters to the well-known sinc method, strengthening the existing technique, and to avoid the aliasing error. The error analysis is established considering both the truncation and amplitude errors associated with the sampling theorem. Numerical examples together with tables and illustrative figures are given. [ABSTRACT FROM PUBLISHER]

Abstract: A numerical procedure based on the method of lines for time-dependent electrodiffusion transport has been developed. Two types of boundary conditions (Neumann and Dirichlet) are considered. Finite difference space discretization with suitably selected weights based on a non-uniform grid is applied. Consistency of this method and the method put forward by Brumleve and Buck are analysed and compared. The resulting stiff system of ordinary differential equations is effectively solved using the RADAU5, RODAS and SEULEX integrators. The applications to selected electrochemical systems: liquid junction, bi-ionic case, ion selective electrodes and electrochemical impedance spectroscopy have been demonstrated. In the paper we promote the use of the full form of the Nernst–Planck and Poisson (NPP) equations, that is including explicitly the electric field as an unknown variable with no simplifications like electroneutrality or constant field assumptions. An effective method of the numerical solution of the NPP problem for arbitrary number of ionic species and valence numbers either for a steady state or a transient state is shown. The presented formulae – numerical solutions to the NPP problem – are ready to be implemented by anyone. Moreover, we make the resulting software freely available to anybody interested in using it. [Copyright &y& Elsevier]

Abstract: This paper presents a numerical simulation and simplified methods for estimating the residual strength of a damaged double bottom. A typical double-bottom structure from a shuttle tanker was modelled. The structure was idealised to a small degree. Imperfections were introduced to the whole structure according to design rules. The most severe situation of the fully loaded condition was investigated. Hull girder bending was considered. Damage was caused by a variety of indenters that were specially designed to obtain a desired damage profile. A total of 16 numerical-calculation cases were investigated using the explicit commercial code LS-DYNA 971. A simplified method was derived. This method includes elastic plus rigid-plastic analyses of a 3-span, single-stiffener model. Unlike Smith''s method, the present approach takes damaged stiffened panels into account in the residual strength assessment. The elastic and rigid-plastic analyses were combined to provide the load-end shortening curve for the damaged stiffeners. A simplified damage mechanism for the tripping of stiffeners is presented, as is a 3-span single stiffener model with proper boundary conditions. Both the analytical and single stiffener models were used to estimate the residual strength of a damaged double bottom. Both methods were validated with a non-linear, finite-element analysis (NLFEA) simulation carried out using LS-DYNA 971. [Copyright &y& Elsevier]

Abstract: In this paper, natural convective flow within a rectangular enclosure has been investigated numerically. All the walls of the enclosure are adiabatic except the bottom wall, which is partially heated and cooled by sinusoidal temperature profile. Both situations: medium is hydro-dynamically isotropic and anisotropic are considered. The governing equations are written under assumption of Brinkman-extended non-Darcy model, including material derivative, and then solved by numerically using spectral element method (SEM). Main emphasize is given on effect of periodicity parameter (N) on local heat transfer rate (Nu x ) as well as flow mechanism in the enclosure. The result shows that, as the periodicity is decreased on increasing N, the absolute value of Nu x at the bottom left corner point increases in both situations. For odd values of N, the local heat transfer profile is symmetric about the line X =0.5. The entire flow is governed by two type convective cells: (i) rotating clockwise (ii) rotating anticlockwise. Furthermore, increasing of N increases the multiple cellular structures in the form of N +1. The relative impact of media permeability on Nu x at (0,0) shows that enhancement of media permeability along Y-direction is much more effective than the same along X-direction. In contrast to isotropic porous media, the effect of the permeability and permeability orientation angle on the flow is complex and often non intuitive. In particular the present analysis shows that, different periodicity of temperature boundary condition has significant effect on the flow mechanism and consequently on the heat transfer rate for both situations. [Copyright &y& Elsevier]

Two-dimensional numerical simulations are carried out to examine the problem of transient electroconvection stability of dielectric liquids subjected to unipolar injection. The entire set of electrohydrodynamics equations associated with the electroconvective phenomena that occur in a layer of a dielectric liquid between two parallel electrodes subjected to a potential difference are solved numerically. We first validate the numerical simulation by comparing our linear stability electroconvection criteria with those obtained by other authors with a stability approach. In this paper, we restrict the study to the strong injection case, which corresponds to values of the non-dimensional injection parameter C greater than or equal to 10. The numerical solution of the electroconvective problem is then presented for rigid lateral boundary conditions. A detailed analysis of the scenario that occurs for different characteristic values of the stability parameter T is provided. The flow structure and its behaviour highlight the existence of different regimes, from laminar to chaotic. The development of charged plumes has been observed in particular. We compute the electrical Nusselt number for different values of the stability parameter and ion mobility. The electrical Nusselt number saturates with increasing T, a fact that it is in agreement with available experimental data. Finally, a spectral analysis is conducted for different aspect ratios of the computational domain. The spectral analysis gives an insight into the physical origin of the velocity and current oscillations. [ABSTRACT FROM AUTHOR]

Open boundaries (OBs) are usually unavoidable in numerical coastal circulation simulations. At OBs, an appropriate open boundary condition (OBC) is required so that outgoing waves freely pass to the exterior without creating reflections back into the interior of the computational domain. In this paper, the authors derive, based on the shallow-water equations including bottom friction and neglecting Coriolis effect and by means of nonlinear characteristic analysis, an OBC formulation with two predictive parameters, phase speed cr, and decay time Tf. Simple idealized tests are performed to demonstrate the proposed OBC's excellent skills in elimination of unwanted reflections at OBs when the motion is periodic, as assumed in its theoretical derivation. It turns out that the formulas for the two OBC parameters become independent of period in the limit of small friction and/or short period. This feature is used to derive an OBC applicable when information about the typical period of the motion to be simulated is unavailable. Simple, idealized tests of this period independent OBC demonstrate its ability to afford excellent results, even when the limitations inherent in its derivation are exceeded. Finally, the OBC is applied in more realistic simulations, including Coriolis effects of 2D tidal flows, and is shown to yield excellent results, especially for residual flows. [ABSTRACT FROM AUTHOR]

This paper deals with a diffusive logistic equation on one dimensional isotropically growing domain. The model equation on growing domains is first presented, and the comparison principle is then proved. The asymptotic behavior of temporal solutions to the reaction-diffusion problem is given by constructing upper and lower solutions. Our result shows that when the domain grows slowly, the species successfully spreads to the whole habitat and stabilizes at a positive steady state, while it dies out in the long run if the domain grows fast. Numerical simulations are also presented to illustrate the analytical result. [ABSTRACT FROM AUTHOR]

Abstract: During installation of a displacement pile, the soil around the pile is heavily distorted. The resulting changes in soil density and stress state around the pile determine the ultimate pile capacity. In most finite element models, the installation phase is not explicitly modelled. In this paper the full installation phase is modelled in two different ways, in a numerical code capable of large deformations. In the first approach the pile is fixed and the soil flows around the pile. This modelling approach requires somewhat unrealistic boundary conditions and requires that the results for the pile installation are for a non-stationary phase of the calculation, whereas formally only the values at the stationary full penetration phase are reliable. A second modelling approach has been introduced to overcome these limitations. In the new approach the initial conditions are set as a stress level, increasing linearly with depth, resulting from a gravity loading stage, and a stepwise penetration of the pile into the soil is prescribed. This approach updates the geometry of the problem domain and keeps updating the convective terms in a fixed mesh. Results are compared with experimental results from centrifuge tests. Although both models show the development of large effective vertical stress below the pile base and porosity change near the pile shaft, also differences with the experimental results are found. Especially the stiffness response during pile installation is difficult to model correctly. [Copyright &y& Elsevier]

Numerical modelling plays an important role in helping us understand the characteristics of seismic wave propagation. The presence of spurious reflections from the boundaries of the truncated computational domain is a prominent problem in finite difference computations. The nearly perfectly matched layer has been proven to be a very effective boundary condition to absorb outgoing waves in both electromagnetic and acoustic media. In this paper, the nearly perfectly matched layer technique is applied to elastic isotropic media to further test the method's absorbing ability. The staggered-grid finite-difference method (fourth-order accuracy in space and second-order accuracy in time) is used in the numerical simulation of seismic wave propagation in 2D Cartesian coordinates. In the numerical tests, numerical comparisons between the nearly perfectly matched layer and the convolutional perfectly matched layer, which is considered the best absorbing layer boundary condition, is also provided. Three numerical experiments demonstrate that the nearly perfectly matched layer has a similar performance to the convolutional perfectly matched layer and can be a valuable alternative to other absorbing layer boundary conditions. [ABSTRACT FROM AUTHOR]

Our motivation is a mathematical model describing the spatial propagation of an epidemic disease through a population. In this model, the pathogen diversity is structured into two clusters and then the population is divided into eight classes which permits to distinguish between the infected/uninfected population with respect to clusters. In this paper, we prove the weak and the global existence results of the solutions for the considered reaction-diffusion system with Neumann boundary. Next, mathematical Turing formulation and numerical simulations are introduced to show the pattern formation for such systems. [ABSTRACT FROM AUTHOR]

This paper is concerned with a two-species predator-prey reaction-diffusion system with Beddington-DeAngelis functional response and subject to homogeneous Neumann boundary conditions. By linearizing the system at the positive constant steady-state solution and analyzing the associated characteristic equation in detail, the asymptotic stability of the positive constant steady-state solution and the existence of local Hopf bifurcations are investigated. Also, it is shown that the appearance of the diffusion and homogeneous Neumann boundary conditions can lead to the appearance of codimension two Bagdanov-Takens bifurcation. Moreover, by applying the normal form theory and the center manifold reduction for partial differential equations (PDEs), the explicit algorithm determining the direction of Hopf bifurcations and the stability of bifurcating periodic solutions is given. Finally, numerical simulations supporting the theoretical analysis are also included. [ABSTRACT FROM AUTHOR]

Abstract: Thermoplastic composite materials offer new perspective in the mechanical industry, especially in aeronautic industry. The assembly of huge structures by welding is made possible by the ability of the matrix to melt. This paper focuses on ultrasonic welding process where heating is confined at the welding interface. This is achieved thanks to a local mechanical dissipation in triangles specially located at the interface, called energy director. In order to better understand the phenomena that occur at the energy director scale, we propose to model and simulate the polymer flow at the interface. Based on a previous work, the flow under vibration is modeled using three coupled boundary value problems. A specific simulation tool is then developed for solving those three problems. It entails specific numerical methods: a level set method allows to handle the large geometry change, and an iterative solver manages the multiphysical aspects. The novel simulation obtained is validated with a qualitative comparison to experiments. Then, an analysis of the numerical results allows to understand the phenomena that enable welding. A thermomechanical localization heats the tip of the energy director. This initiates a fold of polymer that progressively fills the gap between the two plates to weld, and ensures conditions for adhesion. [Copyright &y& Elsevier]

Abstract: This paper explored numerical techniques to simulate the movement of an applied tracer plume in a laboratory beach subjected to tide. The MARUN model () was used. The movement of the tracer plume in the beach and the development of a transient seepage on the beach surface were modeled. By calibrating the results of the numerical model against the observed data, the laboratory beach was found to have two different zones: a top layer (thickness 9–17cm) with a saturated hydraulic conductivity, K 0, of 0.2cm/s and a lower layer with a K 0 of 0.12cm/s. The simulations revealed that the plume moved seaward during falling tides and downward during rising tides, which is consistent with previous studies. It was also observed that the plume developed a tail that extended in the landward direction, and a freshwater pool became entrapped between the main plume and its extended tail. A sensitivity analysis was conducted to evaluate the effects of boundary conditions, seepage face development, capillarity, and dispersivity values. It was found that a Dirichlet boundary condition landward of the beach causes the plume of applied tracer to sink deeper into the beach in comparison with a Neumann boundary condition. The suppression of the seepage face by forcing the water table to match the tide level did not have great effects on beach hydrodynamics or on the movement of the tracer, which was probably due to the small slope of the beach (10%). However, seepage flow was around half of the seaward flow leaving the beach through the submerged face. Investigation on the capillarity through the van Genuchten parameter revealed that a decrease in the pore size of the beach greatly affects solute hydrodynamics; we noted that the effects on solute transport increase with depth. Thus, we conclude that using a saturated flow model for tidally influenced beaches should be discouraged, especially for fine textured media (i.e., small α or strong capillarity). The simulations also showed that dispersivity does not affect the travel time of the peak of solute concentration at a specific location but affects its magnitude. [Copyright &y& Elsevier]

This paper presents the determination of virtual micro-forming limit diagrams from two types of numerical simulations based on the finite element method: one with modelling of the full tool for micro-deep drawing and a thin blank with geometric imperfections, based on a defined roughness; and a second called "reduced simulation" where different deformation paths were simulated with appropriate boundary conditions by introducing the same type of geometric imperfections on aluminium 1050A (99.5%). A new test for detecting the onset of necking, called the "change of slope" criterion has been defined. Several methods based on histograms to represent the distributions of major and minor strains, have been used to determine the strain at the onset of necking. The numerical micro-forming limit diagrams (MFLD) were then compared to one obtained experimentally. [ABSTRACT FROM AUTHOR]

In this paper, we construct a group of Saul'yev type asymmetric difference formulas for the dispersive equation. Based on these formulas we derive a new alternating 6-point group algorithm to solve dispersive equations with periodic boundary conditions. The algorithm has a high-order accuracy in space and an unconditional stability. The theoretical results are conformed to the numerical simulation. A comparison of this algorithm with the previous Alternating Group Explicit method is presented. [ABSTRACT FROM AUTHOR]

A time-derivative preconditioned system of equations suitable for the numerical simulation of inviscid multicomponent and multiphase flows at all speeds is described. The system is shown to be hyperbolic in time and remains well conditioned in the incompressible limit, allowing time marching numerical methods to remain an efficient solution strategy. It is well known that the application of conservative numerical methods to multicomponent flows containing sharp fluid interfaces will generate nonphysical pressure and velocity oscillations across the component interface. These oscillations may lead to stability problems when the interface separates fluids with large density ratio, such as water and air. The effect of which may lead to the requirement of small physical time steps and slow subiteration convergence for implicit time marching numerical methods. At low speeds the use of nonconservative methods may be considered In this paper a characteristic-based preconditioned nonconservative method is described. This method preserves pressure and velocity equilibrium across fluid interfaces, obtains density ratio independent stability and convergence, and remains well conditioned in the incompressible limit of the equations. To extend the method to transonic and supersonic flows containing shocks, a hybrid formulation is described, which combines a conservative preconditioned Roe method with the nonconservative preconditioned characteristic-based method. The hybrid method retains the pressure and velocity equilibrium at component interfaces and converges to the physically correct weak solution. To demonstrate the effectiveness of the nonconservative and hybrid approaches, a series of one-dimensional multicomponent Riemann problems is solved with each of the methods. The solutions are compared with the exact solution to the Riemann problem, and stability of the numerical methods are discussed. [ABSTRACT FROM AUTHOR]

This paper introduces a new spontaneous potential log model for the case in which formation resistivity is not piecewise constant. The spontaneous potential satisfies an elliptic boundary value problem with jump conditions on the interfaces. It has been shown that the elliptic interface problem has a unique weak solution. Furthermore, a jump condition capturing finite difference scheme is proposed and applied to solve such elliptic problems. Numerical results show validity and effectiveness of the proposed method. [ABSTRACT FROM AUTHOR]

Abstract: A Direct Numerical Simulation (DNS) method has been developed to solve the heat transfer equations for the computation of thermal convection in particulate flows. This numerical method makes use of a finite difference method in combination with the Immersed Boundary (IB) method for treating the particulate phase. A regular Eulerian grid is used to solve the modified momentum and energy equations for the entire flow region simultaneously. In the region that is occupied by the solid particles, a second particle-based Lagrangian grid is used, which tracks particles, and a force density function or an energy density function is introduced to represent the momentum interaction or thermal interaction between particle and fluid. The numerical methods developed in this paper have been validated extensively by comparing the present simulation results with those obtained by others. [Copyright &y& Elsevier]

This paper presents computer simulation of the forging process using the finite volume method (FVM). The process of forging is highly non-linear, where both large deformations and continuously changing boundary conditions occur. In most practical cases, the initial billet shape is relatively simple, but the final shape of the end product is often geometrically complex, to the extent that it is commonly obtained using multiple forming stages. Examples of the numerical simulation of the forged pieces provided were created using Msc/SuperForge computer code. The main results of the analysis are deformed shape, temperature, pressure, effective plastic strain, effective stress and forces acting on the die. [ABSTRACT FROM AUTHOR]

In this paper, we will extend the lattice Boltzmann model (LBM) to deal with the complex coupled nonlinear Schrödinger equations. A D1Q5 nonstandard complex lattice Boltzmann model (CLBM) with amending term for the coupled nonlinear Schrödinger equations (CNLSEs) is proposed. Complex equilibrium distribution functions and amending terms are selected, and the CNLSEs can be recovered by using the Chapman–Enskog multi-scale analysis and Taylor expansion technique. Some initial boundary value problems of the coupled system are simulated to verify the effectiveness and stability of our model. [ABSTRACT FROM AUTHOR]

This paper presents a review of computational uncertainties in scientific computing, as well as quantification of these uncertainties in the context of numerical simulations for thermo-fluid problems. The need for defining a measure of the numerical error that takes into account errors arising from different numerical building blocks of the simulation methods is discussed. In the above context, the effects of grid resolution, initial and boundary conditions, numerical discretization, and physical modeling constraints are presented. [ABSTRACT FROM AUTHOR]