8 research outputs found
Domain decomposition methods for domain composition purpose: Chimera, overset, gluing and sliding mesh methods
Domain composition methods (DCM) consist in
obtaining a solution to a problem, from the formulations of the same problem expressed on various subdomains. These methods have therefore the opposite objective of domain
decomposition methods (DDM). Indeed, in contrast to
DCM, these last techniques are usually applied to matching
meshes as their purpose consists mainly in distributing the
work in parallel environments. However, they are sometimes
based on the same methodology as after decomposing,
DDM have to recompose. As a consequence, in the
literature, the term DDM has many times substituted DCM.
DCM are powerful techniques that can be used for different
purposes: to simplify the meshing of a complex geometry
by decomposing it into different meshable pieces; to perform
local refinement to adapt to local mesh requirements;
to treat subdomains in relative motion (Chimera, sliding
mesh); to solve multiphysics or multiscale problems, etc.
The term DCM is generic and does not give any clue about
how the fragmented solutions on the different subdomains
are composed into a global one. In the literature, many
methodologies have been proposed: they are mesh-based,
equation-based, or algebraic-based. In mesh-based formulations,
the coupling is achieved at the mesh level, before the governing equations are assembled into an algebraic
system (mesh conforming, Shear-Slip Mesh Update,
HERMESH). The equation-based counterpart recomposes
the solution from the strong or weak formulation itself, and
are implemented during the assembly of the algebraic
system on the subdomain meshes. The different coupling
techniques can be formulated for the strong formulation at
the continuous level, for the weak formulation either at the
continuous or at the discrete level (iteration-by-subdomains,
mortar element, mesh free interpolation). Although
the different methods usually lead to the same solutions at
the continuous level, which usually coincide with the
solution of the problem on the original domain, they have
very different behaviors at the discrete level and can be
implemented in many different ways. Eventually, algebraic-
based formulations treat the composition of the
solutions directly on the matrix and right-hand side of the
individual subdomain algebraic systems. The present work
introduces mesh-based, equation-based and algebraicbased
DCM. It however focusses on algebraic-based
domain composition methods, which have many advantages
with respect to the others: they are relatively problem
independent; their implicit implementation can be hidden
in the iterative solver operations, which enables one to
avoid intensive code rewriting; they can be implemented in
a multi-code environment
A Chimera method for the incompressible Navier-Stokes equations
The Chimera method was developed three decades ago as a meshing simplification tool. Different components are meshed independently and then glued together using a domain decomposition technique to couple the equations solved on each component. This coupling is achieved via transmission conditions (in the finite element context) or by imposing the continuity of fluxes (in the finite volume context). Historically, the method has then been used extensively to treat moving objects, as the independent meshes are free to move with respect to the others. At each time step, the main task consists in recomputing the interpolation of the transmission conditions or fluxes. This paper presents a Chimera method applied to the Navier-Stokes equations. After an introduction on the Chimera method, we describe in two different sections the two independent steps of the method: the hole cutting to create the interfaces of the subdomains and the coupling of the subdomains. Then, we present the Navier-Stokes solver considered in this work. Implementation aspects are then detailed in order to apply efficiently the method to this specific parallel Navier-Stokes solver. We conclude with some examples to demonstrate the reliability and application of the proposed method.Peer Reviewe
蛙皮素, 胃泌素釋放胜之拮抗劑調節受體信使核醣核酸含量及受體密度
Deposition of polydisperse particles representing nasal spray application in a human nasal cavity was performed under transient breathing profiles of sniffing, constant flow, and breath hold. The LES turbulence model was used to describe the fluid phase. Particles were introduced into the flow field with initial spray conditions, including spray cone angle, insertion angle, and initial velocity. Since nasal spray atomizer design determines the particle conditions, fifteen particle size distributions were used, each defined by a log-normal distribution with a different volume mean diameter (Dv50). Particle deposition in the anterior region was approximately 80% when Dv50 > 50μm, and this decreased to 45% as Dv50 decreased to 10μ m for constant and sniff breathing conditions. The decrease in anterior deposition was countered with increased deposition in the middle and posterior regions. The significance of increased deposition in the middle region for drug delivery shows there is potential for nasal delivered drugs to reach the highly vascularised mucosal walls in the main nasal passages. For multiple targeted deposition sites, an optimisation equation was introduced where deposition results of any two targeted sites could be combined and a weighting between 0 to 1 was applied to each targeted site, representing the relative importance of each deposition site
Alya: computational solid mechanics for supercomputers
While solid mechanics codes are now conventional tools both in industry and research, the increasingly more exigent requirements of both sectors are fuelling the need for more computational power and more advanced algorithms. For obvious reasons, commercial codes are lagging behind academic codes often dedicated either to the implementation of one new technique, or the upscaling of current conventional codes to tackle massively large scale computational problems. Only in a few cases, both approaches have been followed simultaneously. In this article, a solid mechanics simulation strategy for parallel supercomputers based on a hybrid approach is presented. Hybrid parallelization exploits the thread-level parallelism of multicore architectures, combining MPI tasks with OpenMP threads. This paper describes the proposed strategy, programmed in Alya, a parallel multi-physics code. Hybrid parallelization is specially well suited for the current trend of supercomputers, namely large clusters of multicores. The strategy is assessed through transient non-linear solid mechanics problems, both for explicit and implicit schemes, running on thousands of cores. In order to demonstrate the flexibility of the proposed strategy under advance algorithmic evolution of computational mechanics, a non-local parallel overset meshes method (Chimera-like) is implemented and the conservation of the scalability is demonstrated
Domain Decomposition Methods for Domain Composition Purpose: Chimera, Overset, Gluing and Sliding Mesh Methods
The final publication is available at link.springer.com via http://dx.doi.org/10.1007/s11831-016-9198-8Domain composition methods (DCM) consist in
obtaining a solution to a problem, from the formulations of the same problem expressed on various subdomains. These methods have therefore the opposite objective of domain
decomposition methods (DDM). Indeed, in contrast to
DCM, these last techniques are usually applied to matching
meshes as their purpose consists mainly in distributing the
work in parallel environments. However, they are sometimes
based on the same methodology as after decomposing,
DDM have to recompose. As a consequence, in the
literature, the term DDM has many times substituted DCM.
DCM are powerful techniques that can be used for different
purposes: to simplify the meshing of a complex geometry
by decomposing it into different meshable pieces; to perform
local refinement to adapt to local mesh requirements;
to treat subdomains in relative motion (Chimera, sliding
mesh); to solve multiphysics or multiscale problems, etc.
The term DCM is generic and does not give any clue about
how the fragmented solutions on the different subdomains
are composed into a global one. In the literature, many
methodologies have been proposed: they are mesh-based,
equation-based, or algebraic-based. In mesh-based formulations,
the coupling is achieved at the mesh level, before the governing equations are assembled into an algebraic
system (mesh conforming, Shear-Slip Mesh Update,
HERMESH). The equation-based counterpart recomposes
the solution from the strong or weak formulation itself, and
are implemented during the assembly of the algebraic
system on the subdomain meshes. The different coupling
techniques can be formulated for the strong formulation at
the continuous level, for the weak formulation either at the
continuous or at the discrete level (iteration-by-subdomains,
mortar element, mesh free interpolation). Although
the different methods usually lead to the same solutions at
the continuous level, which usually coincide with the
solution of the problem on the original domain, they have
very different behaviors at the discrete level and can be
implemented in many different ways. Eventually, algebraic-
based formulations treat the composition of the
solutions directly on the matrix and right-hand side of the
individual subdomain algebraic systems. The present work
introduces mesh-based, equation-based and algebraicbased
DCM. It however focusses on algebraic-based
domain composition methods, which have many advantages
with respect to the others: they are relatively problem
independent; their implicit implementation can be hidden
in the iterative solver operations, which enables one to
avoid intensive code rewriting; they can be implemented in
a multi-code environment