Search results for "Solver"
showing 10 items of 157 documents
Systematic Approach for Calculating the Concentrations of Chemical Species in Multiequilibrium Problems: Inclusion of the Ionic Strength Effects
2012
A general systematic approach including ionic strength effects is proposed for the numerical calculation of concentrations of chemical species in multiequilibrium problems. This approach extends the versatility of the approach presented in a previous article and is applied using the Solver option of the Excel spreadsheet to solve real problems such as the calculation of the pH of buffer solutions at any ionic strength. It is useful for undergraduate programs, in post-graduate programs, and in professional laboratories to predict experimental conditions.
A Meshfree Solver for the MEG Forward Problem
2015
Noninvasive estimation of brain activity via magnetoencephalography (MEG) involves an inverse problem whose solution requires an accurate and fast forward solver. To this end, we propose the Method of Fundamental Solutions (MFS) as a meshfree alternative to the Boundary Element Method (BEM). The solution of the MEG forward problem is obtained, via the Method of Particular Solutions (MPS), by numerically solving a boundary value problem for the electric scalar potential, derived from the quasi-stationary approximation of Maxwell’s equations. The magnetic field is then computed by the Biot-Savart law. Numerical experiments have been carried out in a realistic single-shell head geometry. The p…
Molecular dynamics simulations in hybrid particle-continuum schemes: Pitfalls and caveats
2017
Heterogeneous multiscale methods (HMM) combine molecular accuracy of particle-based simulations with the computational efficiency of continuum descriptions to model flow in soft matter liquids. In these schemes, molecular simulations typically pose a computational bottleneck, which we investigate in detail in this study. We find that it is preferable to simulate many small systems as opposed to a few large systems, and that a choice of a simple isokinetic thermostat is typically sufficient while thermostats such as Lowe-Andersen allow for simulations at elevated viscosity. We discuss suitable choices for time steps and finite-size effects which arise in the limit of very small simulation bo…
Total Variation Regularization in Digital Breast Tomosynthesis
2013
We developed an iterative algebraic algorithm for the reconstruction of 3D volumes from limited-angle breast projection images. Algebraic reconstruction is accelerated using the graphics processing unit. We varied a total variation (TV)-norm parameter in order to verify the influence of TV regularization on the representation of small structures in the reconstructions. The Barzilai-Borwein algorithm is used to solve the inverse reconstruction problem. The quality of our reconstructions was evaluated with the Quart Mam/Digi Phantom, which features so-called Landolt ring structures to verify perceptibility limits. The evaluation of the reconstructions was done with an automatic LR detection a…
Decentralized Subspace Projection for Asymmetric Sensor Networks
2020
A large number of applications in Wireless Sensor Networks include projecting a vector of noisy observations onto a subspace dictated by prior information about the field being monitored. In general, accomplishing such a task in a centralized fashion, entails a large power consumption, congestion at certain nodes and suffers from robustness issues against possible node failures. Computing such projections in a decentralized fashion is an alternative solution that solves these issues. Recent works have shown that this task can be done via the so-called graph filters where only local inter-node communication is performed in a distributed manner using a graph shift operator. Most of the existi…
An optimization approach for communal home meal delivery service
2009
Abstract: This paper is the first to discuss the communal home meal delivery problem. The problem can be modelled as a multiple travelling salesman problem with time windows, that is closely related to the well-studied vehicle routing problem with time windows. Experimental results are reported for a real-life case study from Central Finland over several alternative scenarios using the SPIDER commercial solver. The comparison with current practice reveals that a significant savings potential can be obtained using off-the-shelf optimization tools. As such, the potential for supporting real-life communal routing problems can be considered to be important for VRP practitioners.
Numerical experiments with a parallel fast direct elliptic solver on Cray T3E
1997
A parallel fast direct O(N log N) solver is shortly described for linear systems with separable block tridiagonal matrices. A good parallel scalability of the proposed method is demonstrated on a Cray T3E parallel computer using MPI in communication. Also, the sequential performance is compared with the well-known BLKTRI-implementation of the generalized. cyclic reduction method using a single processor of Cray T3E.
Parallelization of a Lattice Boltzmann Suspension Flow Solver
2002
We have applied a parallel Lattice Boltzmann method to solve the behaviour of the suspension flow. The complex behaviour of the suspension flow cannot be solved by analytical methods, so simulations are the only way to study it. Usually the size of an interesting problem is so big that calculation time on one processor is too long, and this can be solved by parallel program. We have written a parallel suspension flow solver and tested it on massive parallel computers. The measured performance of our program show that the parallelization of suspension particles was successful. We also show that over one million particles can be simulated.
Chebyshev’s Method on Projective Fluids
2020
We demonstrate the acceleration potential of the Chebyshev semi-iterative approach for fluid simulations in Projective Dynamics. The Chebyshev approach has been successfully tested for deformable bodies, where the dynamical system behaves relatively linearly, even though Projective Dynamics, in general, is fundamentally nonlinear. The results for more complex constraints, like fluids, with a particular nonlinear dynamical system, remained unknown so far. We follow a method describing particle-based fluids in Projective Dynamics while replacing the Conjugate Gradient solver with Chebyshev’s method. Our results show that Chebyshev’s method can be successfully applied to fluids and potentially…
Fine-Mesh Numerical Simulations for 2D Riemann Problems with a Multilevel Scheme
2001
The numerical simulation of physical problems modeled by systems of conservation laws can be difficult due to the occurrence of discontinuities and other non-smooth features in the solution.