Recursos de colección

DSpace at MIT (104.280 recursos)

This site is a university repository providing access to the publication output of the institution. Registered users can set up email alerts to notify them of newly added relevant content. A certain level of encryption and security is embedded in the site which may cause some users accessibility problems.

High Performance Computation for Engineered Systems (HPCES)

Mostrando recursos 1 - 20 de 87

  1. Shock Capturing with Discontinuous Galerkin Method

    Nguyen, Vinh Tan; Khoo, Boo Cheong; Peraire, Jaime; Persson, Per-Olof
    Shock capturing has been a challenge for computational fluid dynamicists over the years. This article deals with discontinuous Galerkin method to solve the hyperbolic equations in which solutions may develop discontinuities in finite time. The high order discontinuous Galerkin method combining the basis of finite volume and finite element methods has shown a lot of attractive features for a wide range of applications. Various techniques proposed in the literature to deal with discontinuities basically reduce the order of interpolation in the region around these discontinuities. The accuracy of the scheme therefore may be degraded in the vicinity of the shock....

  2. Real-Time Reliable Prediction of Linear-Elastic Mode-I Stress Intensity Factors for Failure Analysis

    Huynh, Dinh Bao Phuong; Peraire, Jaime; Patera, Anthony T.; Liu, Guirong
    Modern engineering analysis requires accurate, reliable and efficient evaluation of outputs of interest. These outputs are functions of "input" parameter that serve to describe a particular configuration of the system, typical input geometry, material properties, or boundary conditions and loads. In many cases, the input-output relationship is a functional of the field variable - which is the solution to an input-parametrized partial differential equations (PDE). The reduced-basis approximation, adopting off-line/on-line computational procedures, allows us to compute accurate and reliable functional outputs of PDEs with rigorous error estimations. The operation count for the on-line stage depends only on a small number...

  3. A Precorrected-FFT Method for Coupled Electrostatic-Stokes Flow Problem

    Nguyen, Ngoc Son; Lim, Kian-Meng; White, Jacob K.
    We present the application of the boundary integral equation method for solving the motion of biological cell or particle under Stokes flow in the presence of electrostatic field. The huge dense matrix-vector product from the boundary integral method poses a computationally challenging problem for solving the large system of equations generated. In our work, we used the precorrected-FFT (pFFT) method to reduce the computational time and memory usage drastically, so that large scale simulations can be performed quickly on a personal computer. Results on the force field acting on the particle, as well as the behavior of the particle through...

  4. Numerical Study of the Poisson-Boltzmann Equation for Biomolecular Electrostatics

    Tan, Lian Hing; Lim, Kian Meng; White, Jacob K.
    Electrostatics interaction plays a very important role in almost all biomolecular systems. The Poisson-Boltzmann equation is widely used to treat this electrostatic effect in an ionic solution. In this work, a simple mixed discrete-continuum model is considered and boundary element method is used to solve for the solution.

  5. Model Order Reduction for Determining Bubble Parameters to Attain a Desired Fluid Surface Shape

    My-Ha, D.; Lim, K. M.; Khoo, Boo Cheong; Willcox, Karen E.
    In this paper, a new methodology for predicting fluid free surface shape using Model Order Reduction (MOR) is presented. Proper Orthogonal Decomposition combined with a linear interpolation procedure for its coefficient is applied to a problem involving bubble dynamics near to a free surface. A model is developed to accurately and efficiently capture the variation of the free surface shape with different bubble parameters. In addition, a systematic approach is developed within the MOR framework to find the best initial locations and pressures for a set of bubbles beneath the quiescent free surface such that the resultant free surface attained...

  6. Numerical comparison between Maxwell stress method and equivalent multipole approach for calculation of the dielectrophoretic force in octupolar cell traps

    Rosales, C.; Lim, K. M.; Khoo, Boo Cheong
    This work presents detailed numerical calculations of the dielectrophoretic force in octupolar traps designed for single-cell trapping. A trap with eight planar electrodes is studied for spherical and ellipsoidal particles using an indirect implementation of the boundary element method (BEM). Multipolar approximations of orders one to three are compared with the full Maxwell stress tensor (MST) calculation of the electrical force on spherical particles. Ellipsoidal particles are also studied, but in their case only the dipolar approximation is available for comparison with the MST solution. The results show that the full MST calculation is only required in the study of...

  7. Preconditioning and iterative solution of symmetric indefinite linear systems arising from interior point methods for linear programming

    Chai, Joo-Siong; Toh, Kim Chuan
    We study the preconditioning of symmetric indefinite linear systems of equations that arise in interior point solution of linear optimization problems. The preconditioning method that we study exploits the block structure of the augmented matrix to design a similar block structure preconditioner to improve the spectral properties of the resulting preconditioned matrix so as to improve the convergence rate of the iterative solution of the system. We also propose a two-phase algorithm that takes advantage of the spectral properties of the transformed matrix to solve for the Newton directions in the interior-point method. Numerical experiments have been performed on some...

  8. FastAero – A Precorrected FFT – Fast Multipole Tree Steady and Unsteady Potential Flow Solver

    Willis, David; Peraire, Jaime; White, Jacob K.
    In this paper a precorrected FFT-Fast Multipole Tree (pFFT-FMT) method for solving the potential flow around arbitrary three dimensional bodies is presented. The method takes advantage of the efficiency of the pFFT and FMT algorithms to facilitate more demanding computations such as automatic wake generation and hands-off steady and unsteady aerodynamic simulations. The velocity potential on the body surfaces and in the domain is determined using a pFFT Boundary Element Method (BEM) approach based on the Green’s Theorem Boundary Integral Equation. The vorticity trailing all lifting surfaces in the domain is represented using a Fast Multipole Tree, time advected, vortex...

  9. An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains

    Le, Duc-Vinh; Khoo, Boo Cheong; Peraire, Jaime
    We present an immersed interface method for the incompressible Navier Stokes equations capable of handling rigid immersed boundaries. The immersed boundary is represented by a set of Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisfied, singular forces are applied on the fluid at the immersed boundary. The forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines. The strength of singular forces is determined by solving a small system of equations at each time step. The Navier-Stokes...

  10. A Conservative Front Tracking Algorithm

    Nguyen, Vinh Tan; Khoo, Boo Cheong; Peraire, Jaime
    The discontinuities in the solutions of systems of conservation laws are widely considered as one of the difficulties in numerical simulation. A numerical method is proposed for solving these partial differential equations with discontinuities in the solution. The method is able to track these sharp discontinuities or interfaces while still fully maintain the conservation property. The motion of the front is obtained by solving a Riemann problem based on the state values at its both sides which are reconstructed by using weighted essentially non oscillatory (WENO) scheme. The propagation of the front is coupled with the evaluation of "dynamic" numerical...

  11. Certified Rapid Solution of Parametrized Linear Elliptic Equations: Application to Parameter Estimation

    Nguyen, N. C.; Liu, Guirong; Patera, Anthony T.
    We present a technique for the rapid and reliable evaluation of linear-functional output of elliptic partial differential equations with affine parameter dependence. The essential components are (i) rapidly uniformly convergent reduced-basis approximations — Galerkin projection onto a space WN spanned by solutions of the governing partial differential equation at N (optimally) selected points in parameter space; (ii) a posteriori error estimation — relaxations of the residual equation that provide inexpensive yet sharp and rigorous bounds for the error in the outputs; and (iii) offline/online computational procedures — stratagems that exploit affine parameter dependence to de-couple the generation and projection stages of...

  12. Viscoelastic Mobility Problem Using A Boundary Element Method

    Nhan, Phan-Thien; Fan, Xi-Jun
    In this paper, the complete double layer boundary integral equation formulation for Stokes flows is extended to viscoelastic fluids to solve the mobility problem for a system of particles, where the non-linearity is handled by particular solutions of the Stokes inhomogeneous equation. Some techniques of the meshless method are employed and a point-wise solver is used to solve the viscoelastic constitutive equation. Hence volume meshing is avoided. The method is tested against the numerical solution for a sphere settling in the Odroyd-B fluid and some results on a prolate motion in shear flow of the Oldroyd-B fluid are reported and...

  13. A Trajectory Piecewise-Linear Approach to Model Order Reduction and Fast Simulation of Nonlinear Circuits and Micromachined Devices

    Rewieński, Michał
    In this paper we present an approach to the nonlinear model reduction based on representing the nonlinear system with a piecewise-linear system and then reducing each of the pieces with a Krylov projection. However, rather than approximating the individual components to make a system with exponentially many different linear regions, we instead generate a small set of linearizations about the state trajectory which is the response to a 'training input'. Computational results and performance data are presented for a nonlinear circuit and a micromachined fixed-fixed beam example. These examples demonstrate that the macromodels obtained with the proposed reduction algorithm are...

  14. Stochastic Transportation-Inventory Network Design Problem

    Shu, Jia; Teo, Chung Piaw; Shen, Zuo-Jun Max
    In this paper, we study the stochastic transportation-inventory network design problem involving one supplier and multiple retailers. Each retailer faces some uncertain demand. Due to this uncertainty, some amount of safety stock must be maintained to achieve suitable service levels. However, risk-pooling benefits may be achieved by allowing some retailers to serve as distribution centers (and therefore inventory storage locations) for other retailers. The problem is to determine which retailers should serve as distribution centers and how to allocate the other retailers to the distribution centers. Shen et al. (2000) and Daskin et al. (2001) formulated this problem as a...

  15. Statistical Error in Particle Simulations of Low Mach Number Flows

    Hadjiconstantinou, Nicolas G.; Garcia, Alejandro L.
    We present predictions for the statistical error due to finite sampling in the presence of thermal fluctuations in molecular simulation algorithms. Expressions for the fluid velocity, density and temperature are derived using equilibrium statistical mechanics. The results show that the number of samples needed to adequately resolve the flow-field scales as the inverse square of the Mach number. The theoretical results are verified for a dilute gas using direct Monte Carlo simulations. The agreement between theory and simulation verifies that the use of equilibrium theory is justified.

  16. Solving symmetric indefinite systems in an interior-point method for second order cone programming

    Toh, Kim Chuan; Cai, Zhi; Freund, Robert M.
    Many optimization problems can be formulated as second order cone programming (SOCP) problems. Theoretical results show that applying interior-point method (IPM) to SOCP has global polynomial convergence. However, various stability issues arise in the implementation of IPM. The standard normal equation based implementation of IPM encounters stability problems in the computation of search direction. In this paper, an augmented system approach is proposed to overcome the stability problems. Numerical experiments show that the new approach can improve the stability.

  17. Solution Methodologies for the Smallest Enclosing Circle Problem

    Xu, Sheng; Freund, Robert M.; Sun, Jie
    Given a set of circles C = {c₁, ..., cn}on the Euclidean plane with centers {(a₁, b₁), ..., (an, bn)}and radii {r₁..., r

  18. Simulation Study of a Semi-Dynamic AGV-Container Unit Job Deployment Scheme

    Cheng, Yong Leong
    Automated Guided Vehicle (AGV) Container-Job deployment is essentially a vehicle-dispatching problem. In this problem, the impact of vehicle dispatching polices on the ship makespan for discharging and/or loading operations is analyzed. In particular, given a storage location for each container to be discharged from the ship and given the current location of each container to be loaded onto the ship, the problem is to propose an efficient deployment scheme to dispatch vehicles to containers so as to minimize the makespan of the ship so as to increase the throughput. The makespan of the ship refers to the time a ship...

  19. A Simulation Method for Calculating the Path Travel Time in Dynamic Transportation Network

    Lin, G.C.; Peraire, Jaime; Khoo, Boo Cheong; Perakis, Georgia
    The calculation of path travel times is an essential component for the dynamic traffic assignment and equilibrium problems. This paper presents a simulation method for calculating actual path travel times for the traffic network with dynamic demands. The method is based on a path-based macroscopic simulation model of network traffic dynamics. There is no need to explicitly model intersection delays in this method. Discontinuity in the travel time caused by traffic light control can be captured by this method. It's flexible in terms that the model is not limited to a specific velocity-density relationship. Some numerical results for signalized and...

  20. A Simple But Effective Evolutionary Algorithm for Complicated Optimization Problems

    Xu, Y.G.; Liu, Guirong
    A simple but effective evolutionary algorithm is proposed in this paper for solving complicated optimization problems. The new algorithm presents two hybridization operations incorporated with the conventional genetic algorithm. It takes only 4.1% ~ 4.7% number of function evaluations required by the conventional genetic algorithm to obtain global optima for the benchmark functions tested. Application example is also provided to demonstrate its effectiveness.

Aviso de cookies: Usamos cookies propias y de terceros para mejorar nuestros servicios, para análisis estadístico y para mostrarle publicidad. Si continua navegando consideramos que acepta su uso en los términos establecidos en la Política de cookies.