Home About us Contact
|
Home About us Contact | ||
|
|
||
Computational Domain (computational + domain)
Selected AbstractsComputational cardiac atlases: from patient to population and backEXPERIMENTAL PHYSIOLOGY, Issue 5 2009Alistair A. Young Integrative models of cardiac physiology are important for understanding disease and planning intervention. Multimodal cardiovascular imaging plays an important role in defining the computational domain, the boundary/initial conditions, and tissue function and properties. Computational models can then be personalized through information derived from in vivo and, when possible, non-invasive images. Efforts are now established to provide Web-accessible structural and functional atlases of the normal and pathological heart for clinical, research and educational purposes. Efficient and robust statistical representations of cardiac morphology and morphodynamics can thereby be obtained, enabling quantitative analysis of images based on such representations. Statistical models of shape and appearance can be built automatically from large populations of image datasets by minimizing manual intervention and data collection. These methods facilitate statistical analysis of regional heart shape and wall motion characteristics across population groups, via the application of parametric mathematical modelling tools. These parametric modelling tools and associated ontological schema also facilitate data fusion between different imaging protocols and modalities as well as other data sources. Statistical priors can also be used to support cardiac image analysis with applications to advanced quantification and subject-specific simulations of computational physiology. [source] 2-D/3-D multiply transmitted, converted and reflected arrivals in complex layered media with the modified shortest path methodGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2009Chao-Ying Bai SUMMARY Grid-cell based schemes for tracing seismic arrivals, such as the finite difference eikonal equation solver or the shortest path method (SPM), are conventionally confined to locating first arrivals only. However, later arrivals are numerous and sometimes of greater amplitude than the first arrivals, making them valuable information, with the potential to be used for precise earthquake location, high-resolution seismic tomography, real-time automatic onset picking and identification of multiple events on seismic exploration data. The purpose of this study is to introduce a modified SPM (MSPM) for tracking multiple arrivals comprising any kind of combination of transmissions, conversions and reflections in complex 2-D/3-D layered media. A practical approach known as the multistage scheme is incorporated into the MSPM to propagate seismic wave fronts from one interface (or subsurface structure for 3-D application) to the next. By treating each layer that the wave front enters as an independent computational domain, one obtains a transmitted and/or converted branch of later arrivals by reinitializing it in the adjacent layer, and a reflected and/or converted branch of later arrivals by reinitializing it in the incident layer. A simple local grid refinement scheme at the layer interface is used to maintain the same accuracy as in the one-stage MSPM application in tracing first arrivals. Benchmark tests against the multistage fast marching method are undertaken to assess the solution accuracy and the computational efficiency. Several examples are presented that demonstrate the viability of the multistage MSPM in highly complex layered media. Even in the presence of velocity variations, such as the Marmousi model, or interfaces exhibiting a relatively high curvature, later arrivals composed of any combination of the transmitted, converted and reflected events are tracked accurately. This is because the multistage MSPM retains the desirable properties of a single-stage MSPM: high computational efficiency and a high accuracy compared with the multistage FMM scheme. [source] A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous mediaGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2006M. De Kool SUMMARY We present a practical grid-based method in 3-D spherical coordinates for computing multiple phases comprising any number of reflection and transmission branches in heterogeneous layered media. The new scheme is based on a multistage approach which treats each layer that the wave front enters as a separate computational domain. A finite-difference eikonal solver known as the fast-marching method (FMM) is reinitialized at each interface to track the evolving wave front as either a reflection back into the incident layer or a transmission through to the adjacent layer. Unlike the standard FMM, which only finds first arrivals, this multistage approach can track those later arriving phases explicitly caused by the presence of discontinuities. Notably, the method does not require an irregular mesh to be constructed in order to connect interface nodes to neighbouring velocity nodes which lie on a regular grid. To improve accuracy, local grid refinement is used in the neighbourhood of a source point where wave front curvature is high. The method also provides a way to trace reflections from an interface that are not the first arrival (e.g. the global PP phase). These are computed by initializing the multistage FMM from both the source and receiver, propagating the two wave fronts to the reflecting interface, and finding stationary points of the sum of the two traveltime fields on the reflecting interface. A series of examples are presented to test the efficiency, accuracy and robustness of the new scheme. As well as efficiently computing various global phases to an acceptable accuracy through the ak135 model, we also demonstrate the ability of the scheme to track complex crustal phases that may be encountered in coincident reflection, wide-angle reflection/refraction or local earthquake surveys. In one example, a variety of phases are computed in the presence of a realistic subduction zone, which includes several layer pinch-outs and a subducting slab. Our numerical tests show that the new scheme is a practical and robust alternative to conventional ray tracing for finding various phases in layered media at a variety of scales. [source] He's homotopy perturbation method for two-dimensional heat conduction equation: Comparison with finite element methodHEAT TRANSFER - ASIAN RESEARCH (FORMERLY HEAT TRANSFER-JAPANESE RESEARCH), Issue 4 2010M. Jalaal Abstract Heat conduction appears in almost all natural and industrial processes. In the current study, a two-dimensional heat conduction equation with different complex Dirichlet boundary conditions has been studied. An analytical solution for the temperature distribution and gradient is derived using the homotopy perturbation method (HPM). Unlike most of previous studies in the field of analytical solution with homotopy-based methods which investigate the ODEs, we focus on the partial differential equation (PDE). Employing the Taylor series, the gained series has been converted to an exact expression describing the temperature distribution in the computational domain. Problems were also solved numerically employing the finite element method (FEM). Analytical and numerical results were compared with each other and excellent agreement was obtained. The present investigation shows the effectiveness of the HPM for the solution of PDEs and represents an exact solution for a practical problem. The mathematical procedure proves that the present mathematical method is much simpler than other analytical techniques due to using a combination of homotopy analysis and classic perturbation method. The current mathematical solution can be used in further analytical and numerical surveys as well as related natural and industrial applications even with complex boundary conditions as a simple accurate technique. © 2010 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/htj.20292 [source] Assessment of flooding in urbanized ungauged basins: a case study in the Upper Tiber area, ItalyHYDROLOGICAL PROCESSES, Issue 10 2005T. Moramarco Abstract The reliability of a procedure for investigation of flooding into an ungauged river reach close to an urban area is investigated. The approach is based on the application of a semi-distributed rainfall,runoff model for a gauged basin, including the flood-prone area, and that furnishes the inlet flow conditions for a two-dimensional hydraulic model, whose computational domain is the urban area. The flood event, which occurred in October 1998 in the Upper Tiber river basin and caused significant damage in the town of Pieve S. Stefano, was used to test the approach. The built-up area, often inundated, is included in the gauged basin of the Montedoglio dam (275 km2), for which the rainfall,runoff model was adapted and calibrated through three flood events without over-bank flow. With the selected set of parameters, the hydrological model was found reasonably accurate in simulating the discharge hydrograph of the three events, whereas the flood event of October 1998 was simulated poorly, with an error in peak discharge and time to peak of ,58% and 20%, respectively. This discrepancy was ascribed to the combined effect of the rainfall spatial variability and a partial obstruction of the bridge located in Pieve S. Stefano. In fact, taking account of the last hypothesis, the hydraulic model reproduced with a fair accuracy the observed flooded urban area. Moreover, incorporating into the hydrological model the flow resulting from a sudden cleaning of the obstruction, which was simulated by a ,shock-capturing' one-dimensional hydraulic model, the discharge hydrograph at the basin outlet was well represented if the rainfall was supposed to have occurred in the region near the main channel. This was simulated by reducing considerably the dynamic parameter, the lag time, of the instantaneous unit hydrograph for each homogeneous element into which the basin is divided. The error in peak discharge and time to peak decreased by a few percent. A sensitivity analysis of both the flooding volume involved in the shock wave and the lag time showed that this latter parameter requires a careful evaluation. Moreover, the analysis of the hydrograph peak prediction due to error in rainfall input showed that the error in peak discharge was lower than that of the same input error quantity. Therefore, the obtained results allowed us to support the hypothesis on the causes which triggered the complex event occurring in October 1998, and pointed out that the proposed procedure can be conveniently adopted for flood risk evaluation in ungauged river basins where a built-up area is located. The need for a more detailed analysis regarding the processes of runoff generation and flood routing is also highlighted. Copyright © 2005 John Wiley & Sons, Ltd. [source] Simulating the hydraulic characteristics of the lower Yellow River by the finite-volume techniqueHYDROLOGICAL PROCESSES, Issue 14 2002Qing Wan Abstract The finite-volume technique is used to solve the two-dimensional shallow-water equations on unstructured mesh consisting of quadrilateral elements. In this paper the algorithm of the finite-volume method is discussed in detail and particular attention is paid to accurately representing the complex irregular computational domain. The lower Yellow River reach from Huayuankou to Jiahetan is a typical meandering river. The generation of the computational mesh, which is used to simulate the flood, is affected by the distribution of water works in the river channel. The spatial information about the two Yellow River levee, the protecting dykes, and those roads that are obviously higher than the ground, need to be used to generate the computational mesh. As a result these dykes and roads locate the element interfaces of the computational mesh. In the model the finite-volume method is used to solve the shallow-wave equations, and the Osher scheme of the empirical function is used to calculate the flux through the interface between the neighbouring elements. The finite-volume method has the advantage of using computational domain with complex geometry, and the Osher scheme is a method based on characteristic theory and is a monotone upwind numerical scheme with high resolution. The flood event with peak discharge of 15 300 m3/s, occurring in the period from 30 July to 10 August 1982, is simulated. The estimated result indicates that the simulation method is good for routing the flood in a region with complex geometry. Copyright © 2002 John Wiley & Sons, Ltd. [source] Matched interface and boundary (MIB) method for the vibration analysis of platesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 9 2009S. N. Yu Abstract This paper proposes a novel approach, the matched interface and boundary (MIB) method, for the vibration analysis of rectangular plates with simply supported, clamped and free edges, and their arbitrary combinations. In previous work, the MIB method was developed for three-dimensional elliptic equations with arbitrarily complex material interfaces and geometric shapes. The present work generalizes the MIB method for eigenvalue problems in structural analysis with complex boundary conditions. The MIB method utilizes both uniform and non-uniform Cartesian grids. Fictitious values are utilized to facilitate the central finite difference schemes throughout the entire computational domain. Boundary conditions are enforced with fictitious values,a common practice used in the previous discrete singular convolution algorithm. An essential idea of the MIB method is to repeatedly use the boundary conditions to achieve arbitrarily high-order accuracy. A new feature in the proposed approach is the implementation of the cross derivatives in the free boundary conditions. The proposed method has a banded matrix. Nine different plates, particularly those with free edges and free corners, are employed to validate the proposed method. The performance of the proposed method is compared with that of other established methods. Convergence and comparison studies indicate that the proposed MIB method works very well for the vibration analysis of plates. In particular, modal bending moments and shear forces predicted by the proposed method vanish at boundaries for free edges. Copyright © 2008 John Wiley & Sons, Ltd. [source] Discrete singular convolution methodology for free vibration and stability analyses of arbitrary straight-sided quadrilateral platesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 11 2008Ömer Civalek Abstract A new discrete singular convolution (DSC) method is developed for vibration, buckling and static analyses of arbitrary straight-sided quadrilateral plates. The straight-sided quadrilateral domain is mapped into a square domain in the computational space using a four-node element. By using the geometric transformation, the governing equations and boundary conditions of the plate are transformed from the physical domain into a square computational domain. Numerical examples illustrating the accuracy and convergence of the DSC method for straight-sided quadrilateral thin plates such as rectangular, skew, trapezoidal and rhombic plates are presented. The results obtained by the DSC method were compared with those obtained by the other numerical and analytical methods. Copyright © 2007 John Wiley & Sons, Ltd. [source] Simulation technique for wave generationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 5 2003S. Aliabadi Abstract In this paper, we present a new finite element technique for simulation of water waves impacting on floating structures. The emphasis will be on the numerical methods for water wave generation and propagation. In our approach, the governing equations are the Navier,Stokes equations written for two incompressible fluids. An interface function with two distinct values serves as a marker identifying the location of the free-surface. This function is transported throughout the computational domain with a time-dependent advection equation. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian,Eulerian domain. This allows us to handle the motion of the physical boundaries, such as the wave generator surface by moving the computational nodes. In the mesh-moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The numerical examples include 3D wave generation and wave breaking as they approach the coast, and the waves impacting on near-shore support columns. Copyright © 2003 John Wiley & Sons, Ltd. [source] Approximate imposition of boundary conditions in immersed boundary methodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 11 2009Ramon Codina Abstract We analyze several possibilities to prescribe boundary conditions in the context of immersed boundary methods. As basic approximation technique we consider the finite element method with a mesh that does not match the boundary of the computational domain, and therefore Dirichlet boundary conditions need to be prescribed in an approximate way. As starting variational approach we consider Nitsche's methods, and we then move to two options that yield non-symmetric problems but that turned out to be robust and efficient. The essential idea is to use the degrees of freedom of certain nodes of the finite element mesh to minimize the difference between the exact and the approximated boundary condition. Copyright © 2009 John Wiley & Sons, Ltd. [source] A parallel implicit/explicit hybrid time domain method for computational electromagneticsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2009Z. Q. Xie Abstract The numerical solution of Maxwell's curl equations in the time domain is achieved by combining an unstructured mesh finite element algorithm with a cartesian finite difference method. The practical problem area selected to illustrate the application of the approach is the simulation of three-dimensional electromagnetic wave scattering. The scattering obstacle and the free space region immediately adjacent to it are discretized using an unstructured mesh of linear tetrahedral elements. The remainder of the computational domain is filled with a regular cartesian mesh. These two meshes are overlapped to create a hybrid mesh for the numerical solution. On the cartesian mesh, an explicit finite difference method is adopted and an implicit/explicit finite element formulation is employed on the unstructured mesh. This approach ensures that computational efficiency is maintained if, for any reason, the generated unstructured mesh contains elements of a size much smaller than that required for accurate wave propagation. A perfectly matched layer is added at the artificial far field boundary, created by the truncation of the physical domain prior to the numerical solution. The complete solution approach is parallelized, to enable large-scale simulations to be effectively performed. Examples are included to demonstrate the numerical performance that can be achieved. Copyright © 2009 John Wiley & Sons, Ltd. [source] Parallel DSMC method using dynamic domain decompositionINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 1 2005J.-S. Wu Abstract A general parallel direct simulation Monte Carlo method using unstructured mesh is introduced, which incorporates a multi-level graph-partitioning technique to dynamically decompose the computational domain. The current DSMC method is implemented on an unstructured mesh using particle ray-tracing technique, which takes the advantages of the cell connectivity information. In addition, various strategies applying the stop at rise (SAR) (IEEE Trans Comput 1988; 39:1073,1087) scheme is studied to determine how frequent the domain should be re-decomposed. A high-speed, bottom-driven cavity flow, including small, medium and large problems, based on the number of particles and cells, are simulated. Corresponding analysis of parallel performance is reported on IBM-SP2 parallel machine up to 64 processors. Analysis shows that degree of imbalance among processors with dynamic load balancing is about ,,½ of that without dynamic load balancing. Detailed time analysis shows that degree of imbalance levels off very rapidly at a relatively low value with increasing number of processors when applying dynamic load balancing, which makes the large problem size fairly scalable for processors more than 64. In general, optimal frequency of activating SAR scheme decreases with problem size. At the end, the method is applied to compute two two-dimensional hypersonic flows, a three-dimensional hypersonic flow and a three-dimensional near-continuum twin-jet gas flow to demonstrate its superior computational capability and compare with experimental data and previous simulation data wherever available. Copyright © 2005 John Wiley & Sons, Ltd. [source] A non-reflecting layer method for non-linear wave-type equations on unbounded domains with applications to shape memory alloy rodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 15 2005M. Newman Abstract In this paper a new technique is introduced and applied in solving one-dimensional linear and non-linear wave-type equations on an unbounded spatial domain. This new technique referred to as the non-reflecting layer method (NRLM) extends the computational domain with an artificial layer on which a one-way wave equation is solved. The method will be applied to compute stress waves in long rods consisting of NiTi shape memory alloy material subjected to impact loading and undergoing detwinning and pseudo-elastic material responses. The NRLM has been tested on model problems and it has been found that the computed solutions agree well with the exact solutions, i.e. normalized error levels are in ranges acceptable for engineering computations. Copyright © 2005 John Wiley & Sons, Ltd. [source] A computational model for impact failure with shear-induced dilatancyINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 14 2003Z. Chen Abstract It has been observed in plate impact experiments that some brittle solids may undergo elastic deformation at the shock wave front, and fail catastrophically at a later time when they are shocked near but below the apparent Hugoniot elastic limit. Because this phenomenon appears to have features different from those of usual inelastic waves, it has been interpreted as the failure wave. To design an effective numerical procedure for simulating impact failure responses, a three-dimensional computational damage model is developed in this paper. The propagation of the failure wave behind the elastic shock wave is described by a non-linear diffusion equation. Macroscopic shear-induced dilatancy is assumed and treated as a one-to-one measure of the mean intensity of microcracking. The damage evolution in time is determined based on the assumption that the deviatoric strain energy in the elastically compressed material (undamaged) is converted, through the damaging process, into the volumetric potential energy in the comminuted and dilated material. For the ease in large-scale simulations, the coupled damage diffusion equation and the stress wave equation are solved via a staggered manner in a single computational domain. Numerical solutions by using both the finite element method and the material point method, i.e. with and without a rigid mesh connectivity, are presented and compared with the experimental data available. It is shown that the model simulations capture the essential features of the failure wave phenomenon observed in shock glasses, and that the numerical solutions for localized failure are not mesh-dependent. Copyright © 2003 John Wiley & Sons, Ltd. [source] Application of a new differential quadrature methodology for free vibration analysis of platesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 6 2003G. Karami Abstract A new methodology is introduced in the differential quadrature (DQ) analysis of plate problems. The proposed approach is distinct from other DQ methods by employing the multiple boundary conditions in a different manner. For structural and plate problems, the methodology employs the displacement within the domain as the only degree of freedom, whereas along the boundaries the displacements as well as the second derivatives of the displacements with respect to the co-ordinate variable normal to the boundary in the computational domain are considered as the degrees of freedom for the problem. Employing such a procedure would facilitate the boundary conditions to be implemented exactly and conveniently. In order to demonstrate the capability of the new methodology, all cases of free vibration analysis of rectangular isotropic plates, in which the conventional DQ methods have had some sort of difficulty to arrive at a converged or accurate solution, are carried out. Excellent convergence behaviour and accuracy in comparison with exact results and/or results obtained by other approximate methods were obtained. The analogous DQ formulation for a general rectangular plate is derived and for each individual boundary condition the general format for imposing the given conditions is devised. It must be emphasized that the computational efforts of this new methodology are not more than for the conventional differential quadrature methods. Copyright © 2002 John Wiley & Sons, Ltd. [source] A box scheme for transcritical flowINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2002T. C. Johnson Abstract The accurate computer simulation of river and pipe flow is of great importance in the design of urban drainage networks. The use of implicit numerical schemes allows the time step to be chosen on the basis of accuracy rather than stability, offering a potential computational saving over explicit methods. The highly successful Box Scheme is an implicit method which can be used to model a wide range of subcritical and supercritical flows. However, care must be taken over the modelling of transcritical flows since, unless the correct internal boundary conditions are imposed, the scheme becomes unstable. The necessity of accurately tracking all the critical interfaces and treating them accordingly can be algorithmically complex and in practice the underlying mathematical model is often modified to ensure that the flow remains essentially subcritical. Such a modification however inevitably leads to additional errors and incorrect qualitative behaviour can be observed. In this paper we show how the technique of ,residual distribution' can be successfully implemented in order to accurately model unsteady transcritical flow without the need to know a priori which regions of the computational domain correspond to subcritical and supercritical flow. When used in conjunction with a form of artificial smoothing, the resulting method generates very high resolution results even for transcritical problems involving shocks, as can be seen in the numerical results. Copyright © 2002 John Wiley & Sons, Ltd. [source] Coupled solution of the species conservation equations using unstructured finite-volume methodINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2010Ankan Kumar Abstract A coupled solver was developed to solve the species conservation equations on an unstructured mesh with implicit spatial as well as species-to-species coupling. First, the computational domain was decomposed into sub-domains comprised of geometrically contiguous cells,a process similar to additive Schwarz decomposition. This was done using the binary spatial partitioning algorithm. Following this step, for each sub-domain, the discretized equations were developed using the finite-volume method, and solved using an iterative solver based on Krylov sub-space iterations, that is, the pre-conditioned generalized minimum residual solver. Overall (outer) iterations were then performed to treat explicitness at sub-domain interfaces and nonlinearities in the governing equations. The solver is demonstrated for both two-dimensional and three-dimensional geometries for laminar methane,air flame calculations with 6 species and 2 reaction steps, and for catalytic methane,air combustion with 19 species and 24 reaction steps. It was found that the best performance is manifested for sub-domain size of 2000 cells or more, the exact number depending on the problem at hand. The overall gain in computational efficiency was found to be a factor of 2,5 over the block (coupled) Gauss,Seidel procedure. All calculations were performed on a single processor machine. The largest calculations were performed for about 355 000 cells (4.6 million unknowns) and required 900,MB of peak runtime memory and 19,h of CPU on a single processor. Copyright © 2009 John Wiley & Sons, Ltd. [source] Nonreflecting boundary conditions based on nonlinear multidimensional characteristicsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2010Qianlong Liu Abstract Nonlinear characteristic boundary conditions based on nonlinear multidimensional characteristics are proposed for 2- and 3-D compressible Navier,Stokes equations with/without scalar transport equations. This approach is consistent with the flow physics and transport properties. Based on the theory of characteristics, which is a rigorous mathematical technique, multidimensional flows can be decomposed into acoustic, entropy, and vorticity waves. Nonreflecting boundary conditions are derived by setting corresponding characteristic variables of incoming waves to zero and by partially damping the source terms of the incoming acoustic waves. In order to obtain the resulting optimal damping coefficient, analysis is performed for problems of pure acoustic plane wave propagation and arbitrary flows. The proposed boundary conditions are tested on two benchmark problems: cylindrical acoustic wave propagation and the wake flow behind a cylinder with strong periodic vortex convected out of the computational domain. This new approach substantially minimizes the spurious wave reflections of pressure, density, temperature, and velocity as well as vorticity from the artificial boundaries, where strong multidimensional flow effects exist. The numerical simulations yield accurate results, confirm the optimal damping coefficient obtained from analysis, and verify that the method substantially improves the 1-D characteristics-based nonreflecting boundary conditions for complex multidimensional flows. Copyright © 2009 John Wiley & Sons, Ltd. [source] Effect of blockage on free vibration of a circular cylinder at low ReINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2008T. K. Prasanth Abstract The effect of the blockage on vortex-induced vibrations of a circular cylinder of low non-dimensional mass (m*=10) in the laminar flow regime is investigated in detail. A stabilized space,time finite element formulation is utilized to solve the incompressible flow equations in primitive variables form in two dimensions. The transverse response of the cylinder is found to be hysteretic at both ends of synchronization/lock-in region for 5% blockage. However, for the 1% blockage hysteresis occurs only at the higher Re end of synchronization/lock-in region. Computations are carried out at other blockages to understand its effect on the hysteretic behavior. The hysteresis loop at the lower Re end of the synchronization decreases with decrease in blockage and is completely eliminated for blockage of 2.5% and less. On the other hand, hysteresis persists for all values of blockage at the higher Re end of synchronization/lock-in. Although the peak transverse oscillation amplitude is found to be same for all blockage (,0.6D), the peak value of the aerodynamic coefficients vary significantly with blockage. The r.m.s. values show lesser variation with blockage. The effect of streamwise extent of computational domain on hysteretic behavior is also studied. The phase between the lift force and transverse displacement shows a jump of almost 180° at, approximately, the middle of the synchronization region. This jump is not hysteretic and is independent of blockage. Copyright © 2008 John Wiley & Sons, Ltd. [source] The dispersal processes within the tide-modulated Changjiang River plume, ChinaINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2007Li-Feng Lu Abstract The dispersal processes of the tide-modulated Changjiang River plume, China, are studied by using a three-dimensional hydrodynamical module of the COHERENS (A COupled Hydrodynamical,Ecological model for REgional and Shelf Seas). The model is driven by the river discharge and the M2 tidal constituent. Modelled results show: (1) the fresh water, which forms the Changjiang River plume expanding southeastwards, is discharged mostly into the North Channel, the North Passage, and the South Passage; (2) the larger horizontal gradient outside the North Channel and the North Passage forms a strong plume front; (3) the Changjiang River plume is homogeneous vertically, and dispersing gradually within the computational domain, with an averaged propagating rate of 3.38,km/day, while the plume front is surface-to-bottom type, and trapped between ,10 and ,18m isobaths; and (4) both the plume length and the plume front intensity vary periodically. The maximum plume length occurs about 2,h after low slack water and the minimum plume length during high slack water. The maximum plume front intensity occurs during high slack water and the minimum plume front intensity during low slack water. Copyright © 2007 John Wiley & Sons, Ltd. [source] A parallel cell-based DSMC method on unstructured adaptive meshesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2004Min Gyu Kim Abstract A parallel DSMC method based on a cell-based data structure is developed for the efficient simulation of rarefied gas flows on PC-clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two-dimensional supersonic leading-edge flow, the axi-symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC-based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd. [source] Numerical simulation of vortical ideal fluid flow through curved channelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2003N. P. Moshkin Abstract A numerical algorithm to study the boundary-value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co-ordinate system. The convergence of the finite-difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka,Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two-dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd. [source] Evaluation of Smagorinsky-based subgrid-scale models in a finite-volume computationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2002Petri Majander Abstract Smagorinsky-based models are assessed in a turbulent channel flow simulation at Reb=2800 and Reb=12500. The Navier,Stokes equations are solved with three different grid resolutions by using a co-located finite-volume method. Computations are repeated with Smagorinsky-based subgrid-scale models. A traditional Smagorinsky model is implemented with a van Driest damping function. A dynamic model assumes a similarity of the subgrid and the subtest Reynolds stresses and an explicit filtering operation is required. A top-hat test filter is implemented with a trapezoidal and a Simpson rule. At the low Reynolds number computation none of the tested models improves the results at any grid level compared to the calculations with no model. The effect of the subgrid-scale model is reduced as the grid is refined. The numerical implementation of the test filter influences on the result. At the higher Reynolds number the subgrid-scale models stabilize the computation. An analysis of an accurately resolved flow field reveals that the discretization error overwhelms the subgrid term at Reb=2800 in the most part of the computational domain. Copyright © 2002 John Wiley & Sons, Ltd. [source] A ,-coordinate three-dimensional numerical model for surface wave propagationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2002Pengzhi Lin Abstract A three-dimensional numerical model based on the full Navier,Stokes equations (NSE) in , -coordinate is developed in this study. The , -coordinate transformation is first introduced to map the irregular physical domain with the wavy free surface and uneven bottom to the regular computational domain with the shape of a rectangular prism. Using the chain rule of partial differentiation, a new set of governing equations is derived in the , -coordinate from the original NSE defined in the Cartesian coordinate. The operator splitting method (Li and Yu, Int. J. Num. Meth. Fluids 1996; 23: 485,501), which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified NSE. The model is first tested for mass and energy conservation as well as mesh convergence by using an example of water sloshing in a confined tank. Excellent agreements between numerical results and analytical solutions are obtained. The model is then used to simulate two- and three-dimensional solitary waves propagating in constant depth. Very good agreements between numerical results and analytical solutions are obtained for both free surface displacements and velocities. Finally, a more realistic case of periodic wave train passing through a submerged breakwater is simulated. Comparisons between numerical results and experimental data are promising. The model is proven to be an accurate tool for consequent studies of wave-structure interaction. Copyright © 2002 John Wiley & Sons, Ltd. [source] Numerical simulation of high-Reynolds number flow around circular cylinders by a three-step FEM,BEM modelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2001D. L. Young Abstract An innovative computational model, developed to simulate high-Reynolds number flow past circular cylinders in two-dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier,Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three-step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high-Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd. [source] Assessment of the performances of first- and second-order time-domain ABC's for the truncation of finite element gridsMICROWAVE AND OPTICAL TECHNOLOGY LETTERS, Issue 1 2003Salvatore Caorsi Abstract In this paper we investigate the performances of first- and second-order time-domain absorbing boundary conditions (ABCs) when introduced in a finite-element algorithm to solve electromagnetic scattering problems. Attention is focused on the analysis of the ABC's absorbing characteristics when different geometries are considered for the truncation of the computational domain. Numerical results will be given by considering, as a first analysis, two-dimensional scattering problems. © 2003 Wiley Periodicals, Inc. Microwave Opt Technol Lett 38: 11,16, 2003 [source] Applications of transformed-space non-uniform PSTD (TSNU-PSTD) in scattering analysis without the use of the non-uniform FFTMICROWAVE AND OPTICAL TECHNOLOGY LETTERS, Issue 1 2003Xiaoping Liu Abstract In this work, we extend the transformed-space, non-uniform pseudo-spectral time domain (TSNU-PSTD) Maxwell solver for a 2D scattering analysis. Prior to implementing the PSTD in this analysis, we first transform the non-uniform grids {xi} and {yj} sampled in the real space for describing complex geometries to uniform ones {ui} and {vj}, in order to fit the dimensions of practical structures and utilize the standard fast Fourier transform (FFT). Next, we use a uniform-sampled, standard FFT to represent spatial derivatives in the space domain of (u, v). It is found that this scheme is as efficient as the conventional uniform PSTD with the computational complexity of O(N log N), since the difference is only the factors of du/dx and dv/dy between the conventional PSTD and the TSNU-PSTD technique. Additionally, we apply an anisotropic version of the Berenger's perfectly matched layers (APML) to suppress the wraparound effect at the open boundaries of the computational domain, which is caused by the periodicity of the FFT. We also employ the pure scattered-field formulation and develop a near-to-far-zone field transformation in order to calculate scattered far fields. © 2003 Wiley Periodicals, Inc. Microwave Opt Technol Lett 38: 16,21, 2003 [source] A global Eta model on quasi-uniform gridsTHE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 623 2007Hai Zhang Abstract The application of quasi-uniform grids in global models of the atmosphere is an attempt to increase the computational efficiency by a more cost-effective exploitation of the computing infrastructure. This paper describes the development of a global version of NCEP's regional, step-coordinate, Eta model on two quasi-uniform grids: cubic and octagonal. The governing equations are expressed in a general curvilinear form, so that the cubic and the octagonal versions of the model share the same code in spite of different mapping of the computational domain. The dynamical core of the derived global Eta model is successfully tested in the benchmark test of Held and Suarez. The model with the step-wise formulation of the terrain and full physics is integrated in a series of tests with real data, and the results are compared both with the analysis and the results of the regional Eta model. Copyright © 2007 Royal Meteorological Society [source] Multiscale resolution of shortwave-longwave interactionCOMMUNICATIONS ON PURE & APPLIED MATHEMATICS, Issue 1 2009Chris Stucchio In the study of time-dependent waves, it is computationally expensive to solve a problem in which high frequencies (shortwaves, with wavenumber k = kmax) and low frequencies (longwaves, near k = kmin) mix. Consider a problem in which low frequencies scatter off a sharp impurity. The impurity generates high frequencies that propagate and spread throughout the computational domain, while the domain must be large enough to contain several longwaves. Conventional spectral methods have a computational cost that is proportional to O(kmax/kmin log(kmax/kmin)). We present here a multiscale algorithm (implemented for the Schrödinger equation but generally applicable) that solves the problem with cost (in space and time) O(kmaxL log(kmax/kmin) log(kmaxL)). Here, L is the width of the region in which the algorithm resolves all frequencies and is independent of kmin. © 2008 Wiley Periodicals, Inc. [source] An a posteriori error estimator for the mimetic finite difference approximation of elliptic problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 11 2008Lourenço Beirão da Veiga Abstract We present an a posteriori error indicator for the mimetic finite difference approximation of elliptic problems in the mixed form. We show that this estimator is reliable and efficient with respect to an energy-type error comprising both flux and pressure. Its performance is investigated by numerically solving the diffusion equation on computational domains with different shapes, different permeability tensors, and different types of computational meshes. Copyright © 2008 John Wiley & Sons, Ltd. [source] | ||||||||||||||||