Developing a Pedagogical Framework for an Integrated and BIM-Based High-Performance Design Studio: Experimental Case Study

As the current environmental crisis and depletion of our energy resources are pushing the Architecture, Engineering, and Construction (AEC) industry toward the design and construction of High-Performance (HP) buildings, new organizational and technological methods of practice, such as Integrated Design Process (IDP) and Building Information...

Computation of nonlinear ultrasound fields using a linearized contrast source method

Nonlinear ultrasound is important in medical diagnostics because imaging of the higher harmonics improves resolution and reduces scattering artifacts. Second harmonic imaging is currently standard, and higher harmonic imaging is under investigation. The efficient development of novel imaging modalities and equipment requires accurate simulations...

On the Convergence Behavior of IDR(s) and Related Methods

An explanation is given of the convergence behavior of IDR(s) methods. The convergence mechanism of these algorithms has two components. The first consists of damping properties of certain factors in the residual polynomials, which becomes less important for large values of s. The second component depends on the behavior of Lanczos polynomials...

Linearization strategies for the Iterative Nonlinear Contrast Source method for full-wave simulation of nonlinear ultrasound fields

The Iterative Nonlinear Contrast Source (INCS) method is a full-wave method for the accurate computation of wide-angle, pulsed, nonlinear ultrasound fields appearing in, e.g., medical echoscopy. The method is based on the Westervelt equation and considers the occurring nonlinear term as a distributed contrast source that operates in a linear...

The Deflated Preconditioned Conjugate Gradient Method Applied to Composite Materials

Simulations with composite materials often involve large jumps in the coefficients of the underlying stiffness matrix. These jumps can introduce unfavorable eigenvalues in the spectrum of the stiffness matrix. We show that the rigid body modes; the translations and rotations, of the disjunct rigid bodies in the composite material correspond to...

Comparison of fundamental, second harmonic, and superharmonic imaging: A simulation study

In medical ultrasound, fundamental imaging (FI) uses the reflected echoes from the same spectral band as that of the emitted pulse. The transmission frequency determines the trade-off between penetration depth and spatial resolution. Tissue harmonic imaging (THI) employs the second harmonic of the emitted frequency band to construct images....

An experimentally validated and parameterized periodic unit-cell reconstruction of open-cell foams

The physical behavior of open-cell foams depends on their microscopic structure. An open-cell geometrical model is proposed, which can serve as the basis for a future macroscopic analysis. The strut geometry is of particular interest, as it is reported to have substantial influence on the occurring thermo-hydraulic and mechanical phenomena....

Separation of blended data by iterative estimation and subtraction of blending interference noise

Seismic acquisition is a trade-off between economy and quality. In conventional acquisition the time intervals between successive records are large enough to avoid interference in time. To obtain an efficient survey, the spatial source sampling is therefore often (too) large. However, in blending, or simultaneous acquisition, temporal overlap...

Efficient Iterative Solution of Large Linear Systems on Heterogeneous Computing Systems

This dissertation deals mainly with the design, implementation, and analysis of efficient iterative solution methods for large sparse linear systems on distributed and heterogeneous computing systems as found in Grid computing. First, a case study is performed on iteratively solving large symmetric linear systems on both a multi–cluster and a...

A contrast source method for nonlinear acoustic wave fields in media with spatially inhomogeneous attenuation

Experimental data reveals that attenuation is an important phenomenon in medical ultrasound. Attenuation is particularly important for medical applications based on nonlinear acoustics, since higher harmonics experience higher attenuation than the fundamental. Here, a method is presented to accurately solve the wave equation for nonlinear...

Interpreting IDR(s) as a deflation method

In this paper the IDR(s) method is interpreted in the context of deflation methods. It is shown that IDR(s) can be seen as a Richardson iteration preconditioned by a variable deflation–type preconditioner. The main result of this paper is the IDR projection theorem, which relates the spectrum of the deflated system in each IDR(s) cycle to all...

Exploiting the flexibility of IDR(s) for grid computing

The IDR(s) method that is proposed in [26] is an efficient limited memory method for solving large nonsymmetric systems of linear equations. In [11] an IDR(s) variant is described that has a single synchronisation point per iteration step, which makes this variant well-suited for parallel and grid computing. In this paper, we combine this IDR(s)...

An iterative method for the computation of nonlinear, wide-angle, pulsed acoustic fields of medical diagnostic transducers

The development and optimization of medical ultrasound transducers and imaging modalities require a computational method that accurately predicts the nonlinear acoustic pressure field. A prospective method should provide the wide-angle, pulsed field emitted by an arbitrary planar source distribution and propagating in a three-dimensional, large...

A filtered convolution method for the computation of acoustic wave fields in very large spatiotemporal domains

The full-wave computation of transient acoustic fields with sizes in the order of 100x100x100 wavelengths by 100 periods requires a numerical method that is extremely efficient in terms of storage and computation. Iterative integral equation methods offer a good performance on these points, provided that the recurring spatiotemporal convolutions...

Multi-area control of overlapping areas in power systems for FACTS control

As power systems generally are large interconnected systems controlled by several parties, centralized optimal power flow (OPF) control taking the entire grid into account is often not feasible. To use optimal control in power systems nevertheless, the overall system is decomposed into areas with associated subproblems, which are solved in an...

Variable two-stage relaxation based on the principle of aggregation multi-grid

Iterative methods are the preferred solution for solving very large algebraic systems of equations. Well-known schemes are damped Jacobi and Successive Over-Relaxation (SOR). While the optimal choice of the relaxation parameter for these schemes has been established for particular examples of system matrices, it is often by trial and error that...

A new iterative solver for the time-harmonic wave equation

The time-harmonic wave equation, also known as the Helmholtz equation, is obtained if the constant-density acoustic wave equation is transformed from the time domain to the frequency domain. Its discretization results in a large, sparse, linear system of equations. In two dimensions, this system can be solved efficiently by a direct method. In...

Fast solution of nonsymmetric linear systems on Grid computers using parallel variants of IDR(s)

IDR(s) is a family of fast algorithms for iteratively solving large nonsymmetric linear systems [14]. With cluster computing and in particular with Grid computing, the inner product is a bottleneck operation. In this paper, three techniques are combined in order to alleviate this bottleneck. Firstly, the efficient and stable IDR(s) algorithm...