We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach t...
详细信息
We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.
We present some of the design, modeling, and simulation features of a computer-aided engineering tool for microelectromechanical systems (MEMS) called SUGAR. The features include a flexible SPICE-like netlist language...
详细信息
We present some of the design, modeling, and simulation features of a computer-aided engineering tool for microelectromechanical systems (MEMS) called SUGAR. The features include a flexible SPICE-like netlist language for MEMS design, a simple modeling framework for computationally efficient lumped models, an extensible architecture to which users can add features, and the ability to display 3-D circuits together with deflected electromechanical structures. Since SUGAR is programmed in MATLAB, many MATLAB functions and toolboxes may be used with SUGAR. Such attributes facilitate the exploration of design spaces and feature modifications. In this paper, we describe SUGAR's extensible architecture, flexible design methodology, modeling framework, and reduced-order modeling technique. We do not present the many other advances made for SUGAR by other developers. For a test case, we choose an advanced microdevice that is difficult to simulate with conventional MEMS software. We show that the relative errors of our lumped models are less than 3% of the finite-element analysis (FEA), that the computational costs are less than 1% of the FEA, and that simulation of the test case fairly agrees with the experiment.
In this paper, a mixed finite-element time-domain (FETD) method is presented for the simulation of electrically complex materials, including general combinations of linear dispersion, instantaneous nonlinearity, and d...
详细信息
In this paper, a mixed finite-element time-domain (FETD) method is presented for the simulation of electrically complex materials, including general combinations of linear dispersion, instantaneous nonlinearity, and dispersive nonlinearity. Using both edge and face elements, the presented method offers greater geometric flexibility than existing finite-difference time-domain (FDTD) implementations, and in contrast to existing nonlinear FETD methods, also incorporates both linear and nonlinear material dispersions. Dielectric nonlinearity is incorporated into the Crank-Nicolson mixed FETD formulation via a straightforward Newton-Raphson approach, for which the associated Jacobian is derived. Moreover, the dispersion is modeled via the Mobius z transform method, yielding a simpler more general algorithm. The method's accuracy and convergence are verified, and its capability demonstrated via the simulation of several nonlinear phenomena, including temporal and spatial solitons in two spatial dimensions.
The mechanical rotation of a radially positioned permanent-magnet Halbach array above a conducting, nonmagnetic track induces eddy currents in the track that can inductively create suspension and propulsion forces sim...
详细信息
The mechanical rotation of a radially positioned permanent-magnet Halbach array above a conducting, nonmagnetic track induces eddy currents in the track that can inductively create suspension and propulsion forces simultaneously. The parameters that affect the performance of this electrodynamic wheel are studied using a 2-D steady-state finite-element method. Tradeoffs between the lift and thrust force performance are investigated and methods to improve the thrust efficiency are proposed.
An efficient finite-element method for driven time-harmonic wave-propagation problems is proposed. The special properties of tangentially continuous vector finiteelements (TVFEM's) are utilized to formulate an un...
详细信息
An efficient finite-element method for driven time-harmonic wave-propagation problems is proposed. The special properties of tangentially continuous vector finiteelements (TVFEM's) are utilized to formulate an ungauged vector-potential scheme in terms of the field method plus one very sparse "gradient matrix" with two nonzero integer or pointer entries per row. The suggested formalism is intended for use with iterative solvers. It combines the simplicity and modest memory requirements of the field formulation with the superior numerical convergence of the ungauged vector-potential scheme.
We have studied electromagnetic losses of a frequency-converter-fed cage-induction motor by using a numerical machine model that includes eddy-current and hysteresis phenomena in electrical steel sheets. We used the m...
详细信息
We have studied electromagnetic losses of a frequency-converter-fed cage-induction motor by using a numerical machine model that includes eddy-current and hysteresis phenomena in electrical steel sheets. We used the model to solve the two-dimensional (2-D) time-harmonic field and winding equations of a cage-induction machine, utilizing a finite-element method and phasor variables. We used complex reluctivity to couple the hysteresis and eddy currents in the sheets with the 2-D analysis. The model modifies the absolute value of the reluctivity according to a one-dimensional (I-D) eddy-current solution developed in the lamination thickness. To define the argument of the reluctivity, we applied both the 1-D field solution and measured hysteresis data. We compared computations of additional electromagnetic losses in a 37-kW test machine due to the higher harmonics of a frequency-converter supply with experimental results. The agreement is found to be reasonable.
Nano-scale adhesive contact mediated by intermolecular van der Waals forces has become a typical fundamental problem in many areas. Interpretation and control of the strength and efficiency of the nano-scale adhesive ...
详细信息
Nano-scale adhesive contact mediated by intermolecular van der Waals forces has become a typical fundamental problem in many areas. Interpretation and control of the strength and efficiency of the nano-scale adhesive contacts require a proper modeling considering the actual interfacial forces, the varying contact area, and clearance. In this article, the finite-element (FE) method is developed to model the nano-scale adhesive contact of elastic bodies with an adhesive pressure derived from the interatomic interaction Lennard-Jones potential, which permits numerical solutions for a variety of interface geometries. Compared with the analytical results from conventional Hertz, JKR, and DMT models, the validity of the FE model is verified. For nano-scale contact, the assumption of equivalent radius adopted in the Hertz model is initially investigated and proved to be improper for nano-scale adhesive contact due to the distribution variations of interfacial force caused by local contact geometry. Then adhesive contact behaviors of four typical nano-scale contact geometries inspired by tip shapes of bio-adhesive pads are investigated in detail, which are flat punch tip, sphere tip, mushroom tip, and empty cup tip. The simulation results indicate that the nano-scale tip geometry plays a dominant role on the pull-off strength. Within the investigated geometries, cup tip results in a highest adhesion efficiency followed by flat punch tip, sphere tip, and mushroom tip, respectively, which are highly geometry dependent and verified by former experimental results. The dominant effect is found coming from the contact area ratio of the adhesive area to the sticking area or the whole contact area. The FE modeling can serve a useful purpose in revealing the nano-scale geometry-based adhesion contact for surface topography design in MEMS to avoid stiction failure and for the artificial sticky feet in bionics to increase adhesion strength.
Nanofluidic thin-layer cells based on redox cycling allow for extremely sensitive electrochemical detection. Here we establish a physical mass-transfer model for analyte molecules in these transducers which takes into...
详细信息
Nanofluidic thin-layer cells based on redox cycling allow for extremely sensitive electrochemical detection. Here we establish a physical mass-transfer model for analyte molecules in these transducers which takes into account advective and diffusive transport of both oxidized and reduced species as well as reversible dynamic adsorption at the sensor surfaces. We use finite-element modeling to determine the transient response of nanogap sensors;numerically we predict that the response time can be reduced substantially by pressure-driven advection while the faradaic limiting current remains unaffected by this flow for all experimentally accessible flow rates. (C) 2013 Elsevier Ltd. All rights reserved.
A new spherical mesh generator is described. It represents an efficient, deterministic packing of tetrahedra into a solid sphere, a spherical shell, or both. The mesh can be used for finite-element solutions to a wide...
详细信息
A new spherical mesh generator is described. It represents an efficient, deterministic packing of tetrahedra into a solid sphere, a spherical shell, or both. The mesh can be used for finite-element solutions to a wide variety of global numerical modelling problems in the geosciences. The nodes within the mesh are distributed uniformly, and long, thin tetrahedra are avoided. The method proposed here offers several advantages over 3-D Delaunay algorithms for finite-element mesh generation. For the related problem of trivariate scattered data interpolation, which is not considered here, the 3-D Delaunay algorithms are the method of choice.
This article provides an intercomparison of the dispersive and diffusive properties of several standard numerical methods applied to the 1D linearized shallow-water equations without the Coriolis term, including upwin...
详细信息
This article provides an intercomparison of the dispersive and diffusive properties of several standard numerical methods applied to the 1D linearized shallow-water equations without the Coriolis term, including upwind and central finite-volume, spectral finite-volume, discontinuous Galerkin, spectral element, and staggered finite-volume. All methods are studied up to tenth-order accuracy, where possible. A consistent framework is developed which allows for direct intercomparison of the ability of these methods to capture the behaviour of linear gravity waves. The Courant-Friedrichs-Lewy (CFL) condition is also computed, which is important for gauging the stability of these methods, and leads to a measure of approximate equal error cost. The goal of this work is threefold: first, to determine the shortest wavelength which can be considered `resolved' for a particular method;second, to determine the effect of increasing the order of accuracy on the ability of a method to capture wave-like motion;and third, to determine which numerical methods offer the best treatment of wave-like motion.
暂无评论