80,648 research outputs found
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
Deep D-Bar: Real-Time Electrical Impedance Tomography Imaging With Deep Neural Networks
The mathematical problem for electrical impedance tomography (EIT) is a highly nonlinear ill-posed inverse problem requiring carefully designed reconstruction procedures to ensure reliable image generation. D-bar methods are based on a rigorous mathematical analysis and provide robust direct reconstructions by using a low-pass filtering of the associated nonlinear Fourier data. Similarly to low-pass filtering of linear Fourier data, only using low frequencies in the image recovery process results in blurred images lacking sharp features, such as clear organ boundaries. Convolutional neural networks provide a powerful framework for post-processing such convolved direct reconstructions. In this paper, we demonstrate that these CNN techniques lead to sharp and reliable reconstructions even for the highly nonlinear inverse problem of EIT. The network is trained on data sets of simulated examples and then applied to experimental data without the need to perform an additional transfer training. Results for absolute EIT images are presented using experimental EIT data from the ACT4 and KIT4 EIT systems
Solving ill-posed inverse problems using iterative deep neural networks
We propose a partially learned approach for the solution of ill posed inverse
problems with not necessarily linear forward operators. The method builds on
ideas from classical regularization theory and recent advances in deep learning
to perform learning while making use of prior information about the inverse
problem encoded in the forward operator, noise model and a regularizing
functional. The method results in a gradient-like iterative scheme, where the
"gradient" component is learned using a convolutional network that includes the
gradients of the data discrepancy and regularizer as input in each iteration.
We present results of such a partially learned gradient scheme on a non-linear
tomographic inversion problem with simulated data from both the Sheep-Logan
phantom as well as a head CT. The outcome is compared against FBP and TV
reconstruction and the proposed method provides a 5.4 dB PSNR improvement over
the TV reconstruction while being significantly faster, giving reconstructions
of 512 x 512 volumes in about 0.4 seconds using a single GPU
Statistical thinking: From Tukey to Vardi and beyond
Data miners (minors?) and neural networkers tend to eschew modelling, misled
perhaps by misinterpretation of strongly expressed views of John Tukey. I
discuss Vardi's views of these issues as well as other aspects of Vardi's work
in emision tomography and in sampling bias.Comment: Published at http://dx.doi.org/10.1214/074921707000000210 in the IMS
Lecture Notes Monograph Series
(http://www.imstat.org/publications/lecnotes.htm) by the Institute of
Mathematical Statistics (http://www.imstat.org
Conditions for duality between fluxes and concentrations in biochemical networks
Mathematical and computational modelling of biochemical networks is often
done in terms of either the concentrations of molecular species or the fluxes
of biochemical reactions. When is mathematical modelling from either
perspective equivalent to the other? Mathematical duality translates concepts,
theorems or mathematical structures into other concepts, theorems or
structures, in a one-to-one manner. We present a novel stoichiometric condition
that is necessary and sufficient for duality between unidirectional fluxes and
concentrations. Our numerical experiments, with computational models derived
from a range of genome-scale biochemical networks, suggest that this
flux-concentration duality is a pervasive property of biochemical networks. We
also provide a combinatorial characterisation that is sufficient to ensure
flux-concentration duality. That is, for every two disjoint sets of molecular
species, there is at least one reaction complex that involves species from only
one of the two sets. When unidirectional fluxes and molecular species
concentrations are dual vectors, this implies that the behaviour of the
corresponding biochemical network can be described entirely in terms of either
concentrations or unidirectional fluxes
- …