This paper describes a second-order projection method for the incompressible Navier-Stokes equations on multiply connected domains with a logically rectangular quadrilateral grid. The method uses a second-order fracti...
详细信息
This paper describes a second-order projection method for the incompressible Navier-Stokes equations on multiply connected domains with a logically rectangular quadrilateral grid. The method uses a second-order fractional step scheme in which one first solves diffusion-convection equations to determine intermediate velocities which are then projected onto the space of divergence-free vector fields. The spatial discretizations are accomplished by formally transforming the equations to a computational space with a uniform grid. The diffusion, pressure gradient, and divergence terms are discretized using standard finite difference approximations. The convection terms are discretized using a second-order Godunov method that provides a robust discretization of these terms at high Reynolds number. Numerical results are presented illustrating the performance of the method.
gradient-based homotopy algorithms have previously been developed for synthesizing H-2 optimal reduced-order dynamic compensators. These algorithms are made efficient and avoid high-order singularities along the homot...
详细信息
gradient-based homotopy algorithms have previously been developed for synthesizing H-2 optimal reduced-order dynamic compensators. These algorithms are made efficient and avoid high-order singularities along the homotopy path by constraining the controller realization to a minimal parameter basis. The resultant homotopy algorithms, however, sometimes experience numerical ill conditioning or failure due to the minimal parameterization constraint. A new homotopy algorithm is presented that is based on solving the optimal projection equations, a set of coupled Riccati and Lyapunov equations that characterize the optimal reduced-order dynamic compensator. Path following in the proposed algorithm is accomplished using a predictor/corrector scheme that computes the prediction and correction steps hy efficiently solving a set of four Lyapunov equations coupled by relatively low-rank linear operators. The algorithm does not suffer from ill conditioning because of constraining the controller basis and often exhibits better numerical properties than the gradient-based homotopy algorithms. The performance of the algorithm is illustrated by considering reduced-order control design for the benchmark four disk axial vibration problem and also reduced-order control of the Active Control Technique Evaluation for Spacecraft structure.
It is well known that the spatial frequency spectrum of membrane and thin plate splines exhibit self-affine characteristics and, hence, behave as fractals. This behavior was exploited in generating the constrained fra...
详细信息
It is well known that the spatial frequency spectrum of membrane and thin plate splines exhibit self-affine characteristics and, hence, behave as fractals. This behavior was exploited in generating the constrained fractal surfaces, which were generated by using a Gibbs sampler algorithm in the work of Szeliski and Terzopoulos. The algorithm involves locally perturbing a constrained spline surface with white noise until the spline surface reaches an equilibrium state. In this paper, we introduce a fast generalized Gibbs sampler that combines two novel techniques, namely, a preconditioning technique in a wavelet basis for constraining the splines and a perturbation scheme in which, unlike the traditional Gibbs sampler, all sites (surface nodes) that do not share a common neighbor are updated simultaneously. In addition, we demonstrate the capability to generate arbitrary order fractal surfaces without resorting to blending techniques. Using this fast Gibbs sampler algorithm, we demonstrate the synthesis of realistic terrain models from sparse elevation data.
Lock time and convergence time are the most important challenges in delay-locked loops (DLLs). In this paper we cover French very high frequency band with a novel all-digital fast-lock DLL-based frequency synthesizer....
详细信息
Lock time and convergence time are the most important challenges in delay-locked loops (DLLs). In this paper we cover French very high frequency band with a novel all-digital fast-lock DLL-based frequency synthesizer. Because this new architecture uses a digital signal processing unit instead of using phase frequency detector, charge pump, and loop filter in conventional DLL, therefore, it shows better jitter performance, lock time, and convergence speed than previous related works. Optimization methods are used to make input and output signals of the proposed DLL in phase. The proposed architecture is designed to cover all channels of French very high frequency band by choosing number of delay cells in signal path. Simulation has been done for 22-27 delay cells, and f(REF)=16MHz, which can produce output frequency in range of 176-216MHz. Locking time is approximately 0.3 mu s, which is equal to five clock cycles of reference clock. All of the simulation results show superiority of the proposed structure. Copyright (c) 2013 John Wiley & Sons, Ltd.
We propose a conjugategradient type optimization technique for the computation of the Karcher mean on the set of complex linear subspaces of fixed dimension, modeled by the so-called Grassmannian. The identification ...
详细信息
We propose a conjugategradient type optimization technique for the computation of the Karcher mean on the set of complex linear subspaces of fixed dimension, modeled by the so-called Grassmannian. The identification of the Grassmannian with Hermitian projection matrices allows an accessible introduction of the geometric concepts required for an intrinsic conjugategradient method. In particular, proper definitions of geodesics, parallel transport, and the Riemannian gradient of the Karcher mean function are presented. We provide an efficient step-size selection for the special case of one dimensional complex subspaces and illustrate how the method can be employed for blind identification via numerical experiments. (C) 2012 Elsevier B.V. All rights reserved.
We present a parallel iterative solver for discrete second order elliptic PDEs. It is based on the conjugate gradient algorithm with incomplete factorization preconditioning, using a domain decomposed ordering to allo...
详细信息
We present a parallel iterative solver for discrete second order elliptic PDEs. It is based on the conjugate gradient algorithm with incomplete factorization preconditioning, using a domain decomposed ordering to allow parallelism in the triangular solves, and resorting to some special recently developed parallelization technique to avoid communication bottle-neck for the computation associated to the internal boundary nodes. Numerical results are given for a transputer network with up to 512 processors and a few workstation cluster.
Scattered data approximation refers to the computation of a multi-dimensional function from measurements obtained from scattered spatial locations. For this problem, the class of methods that adopt a roughness minimiz...
详细信息
Scattered data approximation refers to the computation of a multi-dimensional function from measurements obtained from scattered spatial locations. For this problem, the class of methods that adopt a roughness minimization are the best performing ones. These methods are called variational methods and they are capable of handling contrasting levels of sample density. These methods express the required solution as a continuous model containing a weighted sum of thin-plate spline or radial basis functions with centres aligned to the measurement locations, and the weights are specified by a linear system of equations. The main hurdle in this type of method is that the linear system is ill-conditioned. Further, getting the weights that are parameters of the continuous model representing the solution is only a part of the effort. Getting a regular grid image requires re-sampling of the continuous model, which is typically expensive. We develop a computationally efficient and numerically stable method based on roughness minimization. The method leads to an algorithm that uses standard regular grid array operations only, which makes it attractive for parallelization. We demonstrate experimentally that we get these computational advantages only with a little compromise in performance when compared with thin-plate spline methods.
Purpose: The inversion from the magnetic field to the magnetic susceptibility distribution is ill-posed because the dipole kernel, which relates the magnetic susceptibility to the magnetic field, has zeroes at a pair ...
详细信息
Purpose: The inversion from the magnetic field to the magnetic susceptibility distribution is ill-posed because the dipole kernel, which relates the magnetic susceptibility to the magnetic field, has zeroes at a pair of cone surfaces in the k-space, leading to streaking artifacts on the reconstructed quantitative susceptibility maps (QSM). A method to impose consistency on the cone data (CCD) with structural priors is proposed to improve the solutions of k-space methods. Methods: The information in the cone region is recovered by enforcing structural consistency with structural prior, while information in the noncone trust region is enforced to be consistent with the magnetic field measurements in k-space. This CCD method was evaluated by comparing the initial results of existing QSM algorithms to the QSM results after CCD enhancement with respect to the COSMOS results in simulation, phantom, and in vivo human brain. Results: The proposed method demonstrated suppression of streaking artifacts and the resulting QSM showed better agreement with reference standard QSM compared with other k-space based methods. Conclusion: By enforcing consistency with structural priors in the cone region, the missing data in the cone can be recovered and the streaking artifacts in QSM can be suppressed. (C) 2015 Wiley Periodicals, Inc. Key words: quantitative susceptibility mapping;conjugate
A new experimental device has been developed in order to characterize the phase change material (PCM) thermal properties (thermal conductivity k, sensible and latent heat thermal energy storage, c(p) and L-f) in the s...
详细信息
A new experimental device has been developed in order to characterize the phase change material (PCM) thermal properties (thermal conductivity k, sensible and latent heat thermal energy storage, c(p) and L-f) in the solid phase, during the solid-liquid transition and in the liquid phase. It allows to measure cylindrical samples of maximum 60 mm radius and 10 mm thick. A typical measurement consists in imposing a vertical temperature gradient through the PCM sample driven by a heat source, monitoring during the experiment time all the boundary conditions (temperatures and heat fluxes) and measuring temperature evolution in three locations within the PCM sample. In this work, we will focus only on the solid thermal conductivity characterization. These experiment data are used to solve the inverse heat conduction problem by applying the conjugategradient method and finally, to determine the PCM thermal properties. Two types of composite PCM have been thermally characterized: paraffin mixed with synthetic graphite (Timrex SFG75) and paraffin mixed with graphite waste. (C) 2015 Elsevier Ltd. All rights reserved.
A minimax problem is introduced for the terminal control of a generic dynamical system without disturbances. The maximum magnitude of the weighted output of the system is minimized over a finite interval by the contro...
详细信息
A minimax problem is introduced for the terminal control of a generic dynamical system without disturbances. The maximum magnitude of the weighted output of the system is minimized over a finite interval by the control input of a prescribed class. Such important characteristics of the controlled system appear explicitly in the proposed problem as the maximum magnitude and settling property of the output. Two numerical examples are shown to illustrate the problem. A slewing experiment is also presented to demonstrate the application of the minimax optimal control.
暂无评论