16,163 research outputs found
Modeling seismic wave propagation and amplification in 1D/2D/3D linear and nonlinear unbounded media
To analyze seismic wave propagation in geological structures, it is possible
to consider various numerical approaches: the finite difference method, the
spectral element method, the boundary element method, the finite element
method, the finite volume method, etc. All these methods have various
advantages and drawbacks. The amplification of seismic waves in surface soil
layers is mainly due to the velocity contrast between these layers and,
possibly, to topographic effects around crests and hills. The influence of the
geometry of alluvial basins on the amplification process is also know to be
large. Nevertheless, strong heterogeneities and complex geometries are not easy
to take into account with all numerical methods. 2D/3D models are needed in
many situations and the efficiency/accuracy of the numerical methods in such
cases is in question. Furthermore, the radiation conditions at infinity are not
easy to handle with finite differences or finite/spectral elements whereas it
is explicitely accounted in the Boundary Element Method. Various absorbing
layer methods (e.g. F-PML, M-PML) were recently proposed to attenuate the
spurious wave reflections especially in some difficult cases such as shallow
numerical models or grazing incidences. Finally, strong earthquakes involve
nonlinear effects in surficial soil layers. To model strong ground motion, it
is thus necessary to consider the nonlinear dynamic behaviour of soils and
simultaneously investigate seismic wave propagation in complex 2D/3D geological
structures! Recent advances in numerical formulations and constitutive models
in such complex situations are presented and discussed in this paper. A crucial
issue is the availability of the field/laboratory data to feed and validate
such models.Comment: of International Journal Geomechanics (2010) 1-1
A numerical method for junctions in networks of shallow-water channels
There is growing interest in developing mathematical models and appropriate
numerical methods for problems involving networks formed by, essentially,
one-dimensional (1D) domains joined by junctions. Examples include hyperbolic
equations in networks of gas tubes, water channels and vessel networks for
blood and lymph in the human circulatory system. A key point in designing
numerical methods for such applications is the treatment of junctions, i.e.
points at which two or more 1D domains converge and where the flow exhibits
multidimensional behaviour. This paper focuses on the design of methods for
networks of water channels. Our methods adopt the finite volume approach to
make full use of the two-dimensional shallow water equations on the true
physical domain, locally at junctions, while solving the usual one-dimensional
shallow water equations away from the junctions. In addition to mass
conservation, our methods enforce conservation of momentum at junctions; the
latter seems to be the missing element in methods currently available. Apart
from simplicity and robustness, the salient feature of the proposed methods is
their ability to successfully deal with transcritical and supercritical flows
at junctions, a property not enjoyed by existing published methodologies.
Systematic assessment of the proposed methods for a variety of flow
configurations is carried out. The methods are directly applicable to other
systems, provided the multidimensional versions of the 1D equations are
available
On the stability of projection methods for the incompressible Navier-Stokes equations based on high-order discontinuous Galerkin discretizations
The present paper deals with the numerical solution of the incompressible
Navier-Stokes equations using high-order discontinuous Galerkin (DG) methods
for discretization in space. For DG methods applied to the dual splitting
projection method, instabilities have recently been reported that occur for
coarse spatial resolutions and small time step sizes. By means of numerical
investigation we give evidence that these instabilities are related to the
discontinuous Galerkin formulation of the velocity divergence term and the
pressure gradient term that couple velocity and pressure. Integration by parts
of these terms with a suitable definition of boundary conditions is required in
order to obtain a stable and robust method. Since the intermediate velocity
field does not fulfill the boundary conditions prescribed for the velocity, a
consistent boundary condition is derived from the convective step of the dual
splitting scheme to ensure high-order accuracy with respect to the temporal
discretization. This new formulation is stable in the limit of small time steps
for both equal-order and mixed-order polynomial approximations. Although the
dual splitting scheme itself includes inf-sup stabilizing contributions, we
demonstrate that spurious pressure oscillations appear for equal-order
polynomials and small time steps highlighting the necessity to consider inf-sup
stability explicitly.Comment: 31 page
Renormalization-group investigation of a superconducting -phase transition using five loops calculations
We have studied a Fermi system with attractive -symmetric interaction
at the finite temperatures by the quantum field renormalization group (RG)
method. The RG functions have been calculated in the framework of dimensional
regularization and minimal subtraction scheme up to five loops. It has been
found that for the RG flux leaves the system's stability region --
the system undergoes a first order phase transition. To estimate the
temperature of the transition to superconducting or superfluid phase the RG
analysis for composite operators has been performed using three-loops
approximation. As the result this analysis shows that for systems
estimated phase transition temperature is higher then well known theoretical
estimations based on continuous phase transition formalism
Multi-Dimensional Astrophysical Structural and Dynamical Analysis I. Development of a Nonlinear Finite Element Approach
A new field of numerical astrophysics is introduced which addresses the
solution of large, multidimensional structural or slowly-evolving problems
(rotating stars, interacting binaries, thick advective accretion disks, four
dimensional spacetimes, etc.). The technique employed is the Finite Element
Method (FEM), commonly used to solve engineering structural problems. The
approach developed herein has the following key features:
1. The computational mesh can extend into the time dimension, as well as
space, perhaps only a few cells, or throughout spacetime.
2. Virtually all equations describing the astrophysics of continuous media,
including the field equations, can be written in a compact form similar to that
routinely solved by most engineering finite element codes.
3. The transformations that occur naturally in the four-dimensional FEM
possess both coordinate and boost features, such that
(a) although the computational mesh may have a complex, non-analytic,
curvilinear structure, the physical equations still can be written in a simple
coordinate system independent of the mesh geometry.
(b) if the mesh has a complex flow velocity with respect to coordinate space,
the transformations will form the proper arbitrary Lagrangian- Eulerian
advective derivatives automatically.
4. The complex difference equations on the arbitrary curvilinear grid are
generated automatically from encoded differential equations.
This first paper concentrates on developing a robust and widely-applicable
set of techniques using the nonlinear FEM and presents some examples.Comment: 28 pages, 9 figures; added integral boundary conditions, allowing
very rapidly-rotating stars; accepted for publication in Ap.
- …