Search results for "Numerical"
showing 10 items of 2002 documents
First measurement of proton's charge form factor at very low $Q^2$ with initial state radiation
2017
We report on a new experimental method based on initial-state radiation (ISR) in e-p scattering, in which the radiative tail of the elastic e-p peak contains information on the proton charge form factor ($G_E^p$) at extremely small $Q^2$. The ISR technique was validated in a dedicated experiment using the spectrometers of the A1-Collaboration at the Mainz Microtron (MAMI). This provided first measurements of $G_E^p$ for $0.001\leq Q^2\leq 0.004 (GeV/c)^2$.
Influence of bed roughness on flow and turbulence structure around a partially-buried, isolated freshwater mussel
2023
The present study uses eddy-resolving numerical simulations to investigate how bed roughness affects flow and turbulence structure around an isolated, partially-buried mussel (Unio elongatulus) aligned with the incoming flow. The rough-bed simulations resolve the flow past the exposed part of a gravel bed, whose surface is obtained from a laboratory experiment that also provides some additional data for validation of the numerical model. Results are also discussed for the limiting case of a horizontal smooth bed. Additionally, the effects of varying the level of burial of the mussel inside the substrate and the discharge through the two mussel siphons are investigated via a set of simulatio…
Finite element approximation for a div-rot system with mixed boundary conditions in non-smooth plane domains
1984
The authors examine a finite element method for the numerical approximation of the solution to a div-rot system with mixed boundary conditions in bounded plane domains with piecewise smooth boundary. The solvability of the system both in an infinite and finite dimensional formulation is proved. Piecewise linear element fields with pointwise boundary conditions are used and their approximation properties are studied. Numerical examples indicating the accuracy of the method are given. peerReviewed
Inverse problems and invisibility cloaking for FEM models and resistor networks
2013
In this paper we consider inverse problems for resistor networks and for models obtained via the finite element method (FEM) for the conductivity equation. These correspond to discrete versions of the inverse conductivity problem of Calderón. We characterize FEM models corresponding to a given triangulation of the domain that are equivalent to certain resistor networks, and apply the results to study nonuniqueness of the discrete inverse problem. It turns out that the degree of nonuniqueness for the discrete problem is larger than the one for the partial differential equation. We also study invisibility cloaking for FEM models, and show how an arbitrary body can be surrounded with a layer …
On FE-grid relocation in solving unilateral boundary value problems by FEM
1992
We consider FE-grid optimization in elliptic unilateral boundary value problems. The criterion used in grid optimization is the total potential energy of the system. It is shown that minimization of this cost functional means a decrease of the discretization error or a better approximation of the unilateral boundary conditions, Design sensitivity analysis is given with respect to the movement of nodal points. Numerical results for the Dirichlet-Signorini problem for the Laplace equation and the plane elasticity problem with unilateral boundary conditions are given. In plane elasticity we consider problems with and without friction. peerReviewed
A fast Fourier transform based direct solver for the Helmholtz problem
2018
This article is devoted to the efficient numerical solution of the Helmholtz equation in a two‐ or three‐dimensional (2D or 3D) rectangular domain with an absorbing boundary condition (ABC). The Helmholtz problem is discretized by standard bilinear and trilinear finite elements on an orthogonal mesh yielding a separable system of linear equations. The main key to high performance is to employ the fast Fourier transform (FFT) within a fast direct solver to solve the large separable systems. The computational complexity of the proposed FFT‐based direct solver is O(N log N) operations. Numerical results for both 2D and 3D problems are presented confirming the efficiency of the method discussed…
The FLO Diffusive 1D-2D Model for Simulation of River Flooding
2016
An integrated 1D-2D model for the solution of the diffusive approximation of the shallow water equations, named FLO, is proposed in the present paper. Governing equations are solved using the MArching in Space and Time (MAST) approach. The 2D floodplain domain is discretized using a triangular mesh, and standard river sections are used for modeling 1D flow inside the section width occurring with low or standard discharges. 1D elements, inside the 1D domain, are quadrilaterals bounded by the trace of two consecutive sections and by the sides connecting their extreme points. The water level is assumed to vary linearly inside each quadrilateral along the flow direction, but to remain constant …
An optimization-based approach for solving a time-harmonic multiphysical wave problem with higher-order schemes
2013
This study considers developing numerical solution techniques for the computer simulations of time-harmonic fluid-structure interaction between acoustic and elastic waves. The focus is on the efficiency of an iterative solution method based on a controllability approach and spectral elements. We concentrate on the model, in which the acoustic waves in the fluid domain are modeled by using the velocity potential and the elastic waves in the structure domain are modeled by using displacement.Traditionally, the complex-valued time-harmonic equations are used for solving the time-harmonic problems. Instead of that, we focus on finding periodic solutions without solving the time-harmonic problem…
Characterizing Cavities in Model Inclusion Fullerenes: A Comparative Study
2001
Abstract: The fullerene-82 cavity is selected as a model system in order to test several methods for characterizing inclusion molecules. The methods are based on different technical foundations such as a square and triangular tessellation of the molecular surface, spherical tessellation of the molecular surface, numerical integration of the atomic volumes and surfaces, triangular tessellation of the molecular surface, and cubic lattice approach to the molecular volume. Accurate measures of the molecular volume and surface area have been performed with the pseudorandom Monte Carlo (MCVS) and uniform Monte Carlo (UMCVS) methods. These calculations serve as a reference for the rest of the meth…
Mass-Metallicity Relation from Cosmological Hydrodynamical Simulations and X-ray Observations of Galaxy Groups and Clusters
2018
Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations ag…