共查询到20条相似文献，搜索用时 15 毫秒
1.
The Laplace transform method (LTM) is introduced to solve Burgers' equation. Because of the nonlinear term in Burgers' equation, one cannot directly apply the LTM. Increment linearization technique is introduced to deal with the situation. This is a key idea in this paper. The increment linearization technique is the following: In time level t, we divide the solution u(x, t) into two parts: u(x, tk) and w(x, t), tk⩽t⩽tk+1, and obtain a time‐dependent linear partial differential equation (PDE) for w(x, t). For this PDE, the LTM is applied to eliminate time dependency. The subsequent boundary value problem is solved by rational collocation method on transformed Chebyshev points. To face the well‐known computational challenge represented by the numerical inversion of the Laplace transform, Talbot's method is applied, consisting of numerically integrating the Bromwich integral on a special contour by means of trapezoidal or midpoint rules. Numerical experiments illustrate that the present method is effective and competitive. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
2.
3.
Jian Guo Zhou 《国际流体数值方法杂志》2009,61(8):848-863
A lattice Boltzmann method is developed for solute transport. Proper expressions for the local equilibrium distribution functions enable the method to be formulated on rectangular lattice with the same simple procedure as that on a square lattice. This provides an additional advantage over a lattice Boltzmann method on a square lattice for problems characterized by dominant phenomenon in one direction and relatively weak in another such as solute transport in shear flow over a narrow channel, where the problems can efficiently be approached with fine and coarse meshes, respectively, resulting in more efficient algorithm. The stability conditions are also described. The proposed method on a square lattice is naturally recovered when a square lattice is used. It is verified by solving four tests and compared with the analytical/exact solutions. They are in good agreement, demonstrating that the method is simple, accurate and robust for solute transport. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
4.
We introduce a new submesh strategy for the two‐level finite element method. The numerical results show that the new submesh is able to better capture the boundary layer which is caused by the choice of bubble functions. The effect of an improved approximation of the residual free bubbles is studied for the advective–diffusive equation. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
5.
Operator splitting algorithms are frequently used for solving the advection–diffusion equation, especially to deal with advection dominated transport problems. In this paper an operator splitting algorithm for the three-dimensional advection–diffusion equation is presented. The algorithm represents a second-order-accurate adaptation of the Holly and Preissmann scheme for three-dimensional problems. The governing equation is split into an advection equation and a diffusion equation, and they are solved by a backward method of characteristics and a finite element method, respectively. The Hermite interpolation function is used for interpolation of concentration in the advection step. The spatial gradients of concentration in the Hermite interpolation are obtained by solving equations for concentration gradients in the advection step. To make the composite algorithm efficient, only three equations for first-order concentration derivatives are solved in the diffusion step of computation. The higher-order spatial concentration gradients, necessary to advance the solution in a computational cycle, are obtained by numerical differentiations based on the available information. The simulation characteristics and accuracy of the proposed algorithm are demonstrated by several advection dominated transport problems. © 1998 John Wiley & Sons, Ltd. 相似文献
6.
When transport is advection-dominated, classical numerical methods introduce excessive artificial diffusion and spurious oscillations. Special methods are required to overcome these phenomena. To solve the advection‒diffusion equation, a numerical method is developed using a discontinuous finite element method for the discretization of the advective terms. At the discontinuities of the approximate solution, numerical advective fluxes are calculated using one-dimensional approximate Riemann solvers. The method is stabilized with a multidimensional slope limiter which introduces small amounts of numerical diffusion when sharp concentration fronts occur. In addition, the diffusive term is discretized using a mixed hybrid finite element method. With this approach, numerical oscillations are completely avoided for a full range of cell Peclet numbers. The combination of discontinuous and mixed finite elements can be easily applied to 2D and 3D models using various types of elements in regular and irregular meshes. Numerical tests show good agreement with 1D and 2D analytical solutions. This approach is compared at the same time with two different numerical methods, a standard mixed finite method and a finite volume approach with high-resolution upwind terms. Regular and irregular meshes are used for the numerical tests to study the mesh effects on the numerical results. Our data show that in all cases this approach performs well. © 1997 by John Wiley & Sons, Ltd. 相似文献
7.
This paper is concerned with the solution of heterogeneous problems by the interface control domain decomposition (ICDD) method, a strategy introduced for the solution of partial differential equations in computational domains partitioned into subdomains that overlap. After reformulating the original boundary value problem by introducing new additional control variables, the unknown traces of the solution at internal subdomain interfaces; the latter are determined by requiring that the (a priori) independent solutions in each subdomain undergo the minimization of a suitable cost functional. We provide an abstract formulation for coupled heterogeneous problems and a general theorem of well‐posedness for the associated ICDD problem. Then, we illustrate and validate an efficient algorithm based on the solution of the Schur‐complement system restricted solely to the interface control variables by considering two kinds of heterogeneous boundary value problems: the coupling between pure advection and advection–diffusion equations and the coupling between Stokes and Darcy equations. In the latter case, we also compare the ICDD method with a classical approach based on the Beavers–Joseph–Saffman conditions. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
8.
Multidimensional residual distribution schemes for the convection–diffusion equation are described. Compact upwind cell vertex schemes are used for the discretization of the convective term. For the diffusive term, two approaches are compared: the classical finite element Galerkin formulation, which preserves the compactness of the stencil used for the convective part, and various residual-based approaches in which the diffusive term, evaluated after a reconstruction step, is upwinded along with the convective term. 相似文献
9.
Recently, we developed an explicit a posteriori error estimator especially suited for fluid dynamics problems solved with a stabilized method. The technology is based upon the theory that inspired stabilized methods, namely, the variational multiscale theory. The salient features of the formulation are that it can be readily implemented in existing codes, it is a very economical procedure, and it yields very accurate local error estimates uniformly from the diffusive to the advective regime. In this work, the variational multiscale error estimator is applied to develop adaptive strategies for the advection–diffusion‐reaction equation. The performance of L1 and L2 local error norms combined with three strategies to adapt the mesh is investigated. Emphasis is placed on flows with sharp boundary and interior layers but also attention is given to diffusion‐dominated flows. Computational results show that the method generates meshes with a smooth transition of the element size, which capture all the flow features. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
10.
A. Pestiaux S.A. Melchior J.F. Remacle T. Kärnä T. Fichefet J. Lambrechts 《国际流体数值方法杂志》2014,75(5):365-384
The discretization of a diffusion equation with a strong anisotropy by a discontinuous Galerkin finite element method is investigated. This diffusion term is implemented in the tracer equation of an ocean model, thanks to a symmetric tensor that is composed of diapycnal and isopycnal diffusions. The strong anisotropy comes from the difference of magnitude order between both diffusions. As the ocean model uses interior penalty terms to ensure numerical stability, a new penalty factor is required in order to correctly deal with the anisotropy of this diffusion. Two penalty factors from the literature are improved and established from the coercivity property. One of them takes into account the diffusion in the direction normal to the interface between the elements. After comparison, the latter is better because the spurious numerical diffusion is weaker than with the penalty factor proposed in the literature. It is computed with a transformed coordinate system in which the diffusivity tensor is diagonal, using its eigenvalue decomposition. Furthermore, this numerical scheme is validated with the method of manufactured solutions. It is finally applied to simulate the evolution of temperature and salinity due to turbulent processes in an idealized Arctic Ocean. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
11.
A. P. S. Selvadurai 《Transport in Porous Media》2004,56(1):51-60
This paper presents a proof of the uniqueness theorem for the initial boundary value problem governing advective–diffusive transport of a chemical in a fluid-saturated non-deformable isotropic, homogeneous porous medium. The advective Darcy flow in the porous medium results from the gradient of a hydraulic potential, which is derived from a well-posed problem in potential theory. The paper discusses the relevant set of consistent boundary conditions applicable to the potential inducing the advective flow and to the concentration field, which ensures uniqueness of the solution. 相似文献
12.
M. J. Martinez 《国际流体数值方法杂志》2006,50(3):347-376
The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass‐lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov–Galerkin streamline upwind schemes in the context of advection–diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection‐dominated problems. Published in 2005 by John Wiley & Sons, Ltd. 相似文献
13.
We derive a smoothed particle hydrodynamics (SPH) approximation for anisotropic dispersion that only depends upon the first derivative of the kernel function and study its numerical properties. In addition, we compare the performance of the newly derived SPH approximation versus an implementation of the particle strength exchange (PSE) method and a standard finite volume method for simulating multiple scenarios defined by different combinations of physical and numerical parameters. We show that, for regularly spaced particles, given an adequate selection of numerical parameters such as kernel function and smoothing length, the new SPH approximation is comparable with the PSE method in terms of convergence and accuracy and similar to the finite volume method. On other hand, the performance of both particle methods (SPH and PSE) decreases as the degree of disorder of the particle increases. However, we demonstrate that in these situations the accuracy and convergence properties of both particle methods can be improved by an adequate choice of some numerical parameters such as kernel core size and kernel function. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
14.
AbstractEnlightened by the Caputo fractional derivative, the present study deals with a novel mathematical model of generalized thermoelasticity to investigate the transient phenomena due to the influence of magnetic field and moving heat source in a rod in the context of three-phase lag (TPL) theory of thermoelasticity. Both ends of the rod are fixed and heat insulated. Employing Laplace transform as a tool, the problem has been transformed into the space-domain and solved analytically. Finally, solutions in the real-time domain are obtained by applying the inverse Laplace transform. Numerical calculation for stress, displacement, and temperature within the rod is carried out and displayed graphically. The effect of moving heat source speed on temperature, stress, and temperature is studied. It is found from the distributions that the temperature, thermally induced displacement and stress of the rod are found to decrease at large source speed. For the better understanding of the effect of moving heat source on all the distributions, three animations are added. 相似文献
15.
插值型重构核粒子法的形函数具有离散点插值特性和不低于核函数的高阶光滑性,因而不仅可以直接施加本质边界条件,同时也保证了较高的计算精度．本文将弹性动力学方程作拉氏变换后,在变换域内用插值型重构核粒子法求解,最后再借助Durbin数值反演方法求得时间域的解．针对典型的弹性动力学问题,给出了插值型重构核粒子法的数值算例,并验证了本文方法的有效性． 相似文献
16.
F. Navarrina I. Colominas M. Casteleiro L. Cueto‐Felgueroso H. Gómez J. Fe A. Soage 《国际流体数值方法杂志》2008,56(5):507-523
In this paper, a numerical model for the simulation of the hydrodynamics and of the evolution of the salinity in shallow water estuaries is presented. This tool is intended to predict the possible effects of Civil Engineering public works and other human actions (dredging, building of docks, spillages, etc.) on the marine habitat, and to evaluate their environmental impact in areas with high productivity of fish and of seafood. The prediction of these effects is essential in the decision making about the different options that could be implemented. The mathematical model consists of two coupled systems of differential equations: the shallow water hydrodynamic equations (that describe the evolution of the depth and of the velocity field) and the shallow water advective–diffusive transport equation (that describes the evolution of the salinity level). Some important issues that must be taken into account are the effects of the tides (including that the seabed could be exposed), the volume of fresh water provided by the rivers and the effects of the winds. Thus, different types of boundary conditions are considered. The numerical model proposed for solving this problem is a second‐order Taylor–Galerkin finite element formulation. The proposed approach is applied to a real case: the analysis of the possible effects of dredging Los Lombos del Ulla, a formation of sandbanks in the Arousa Estuary (Galicia, Spain). A number of simulations have been carried out to compare the actual salinity level with the predicted situation if the different dredging options were executed. Some of the obtained results are presented and discussed. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
17.
Federico Vismara;Tommaso Benacchio; 《国际流体数值方法杂志》2024,96(2):161-188
We introduce an extended discontinuous Galerkin discretization of hyperbolic–parabolic problems on multidimensional semi-infinite domains. Building on previous work on the one-dimensional case, we split the strip-shaped computational domain into a bounded region, discretized by means of discontinuous finite elements using Legendre basis functions, and an unbounded subdomain, where scaled Laguerre functions are used as a basis. Numerical fluxes at the interface allow for a seamless coupling of the two regions. The resulting coupling strategy is shown to produce accurate numerical solutions in tests on both linear and nonlinear scalar and vectorial model problems. In addition, an efficient absorbing layer can be simulated in the semi-infinite part of the domain in order to damp outgoing signals with negligible spurious reflections at the interface. By tuning the scaling parameter of the Laguerre basis functions, the extended DG scheme simulates transient dynamics over large spatial scales with a substantial reduction in computational cost at a given accuracy level compared to standard single-domain discontinuous finite element techniques. 相似文献
18.
Based on linearized 2-D Navier-Stokes equation, a Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis is proposed. Under assumption of incompressibility for the fluid, the corresponding equivalent boundary integral equation in terms of the potential function and stream function is first established by Lamb's transform in the Laplace transform domain. It enables us to solve impact water problems in frequency domain by the boundary element method, in which the effect of viscous flow on the dynamic response can be taken into account. Then a complete solution of the problem under consideration in time domain is obtained by means of Durbin's formulas for the numerical inversion of the Laplace transform. Finally, a practical example is given to confirm the validity of the present method. Project supported by the National Defence Foundation of Science & Technology of China (No. J14. 8. 1. JW0515). 相似文献
19.
We present a new stabilized method for advection–diffusion equations, which combines a control volume FEM formulation of the governing equations with a novel multiscale approximation of the total flux. The latter incorporates information about the exact solution that cannot be represented on the mesh. To define this flux, we solve the governing equations along suitable mesh segments under the assumption that the flux varies linearly along these segments. This procedure yields second‐order accurate fluxes on the edges of the mesh. Then, we use curl‐conforming elements of the same order to lift these edge fluxes into the mesh elements. In so doing, we obtain a stabilized control volume FEM formulation that is second‐order accurate and does not require mesh‐dependent stabilization parameters. Numerical convergence studies on uniform and nonuniform grids along with several standard advection tests illustrate the computational properties of the new method. Published 2015. This article is a U.S. Government work and is in the public domain in the USA. 相似文献
20.
In the last decade, the characterization of transport in porous media has benefited largely from numerical advances in applied mathematics and from the increasing power of computers. However, the resolution of a transport problem often remains cumbersome, mostly because of the time-dependence of the equations and the numerical stability constraints imposed by their discretization. To avoid these difficulties, another approach is proposed based on the calculation of the temporal moments of a curve of concentration versus time. The transformation into the Laplace domain of the transport equations makes it possible to develop partial derivative equations for the calculation of complete moments or truncated moments between two finite times, and for any point of a bounded domain. The temporal moment equations are stationary equations, independent of time, and with weaker constraints on their stability and diffusion errors compared to the classical advection–dispersion equation, even with simple discrete numerical schemes. Following the complete theoretical development of these equations, they are compared firstly with analytical solutions for simple cases of transport and secondly with a well-performing transport model for advective–dispersive transport in a heterogeneous medium with rate-limited mass transfer between the free water and an immobile phase. Temporal moment equations have a common parametrization with transport equations in terms of their parameters and their spatial distribution on a grid of discretization. Therefore, they can be used to replace the transport equations and thus accelerate the achievement of studies in which a large number of simulations must be carried out, such as the inverse problem conditioned with transport data or for forecasting pollution hazards. 相似文献