53 resultados para discontinuous Galerkin
Resumo:
The optimal control of problems that are constrained by partial differential equations with uncertainties and with uncertain controls is addressed. The Lagrangian that defines the problem is postulated in terms of stochastic functions, with the control function possibly decomposed into an unknown deterministic component and a known zero-mean stochastic component. The extra freedom provided by the stochastic dimension in defining cost functionals is explored, demonstrating the scope for controlling statistical aspects of the system response. One-shot stochastic finite element methods are used to find approximate solutions to control problems. It is shown that applying the stochastic collocation finite element method to the formulated problem leads to a coupling between stochastic collocation points when a deterministic optimal control is considered or when moments are included in the cost functional, thereby forgoing the primary advantage of the collocation method over the stochastic Galerkin method for the considered problem. The application of the presented methods is demonstrated through a number of numerical examples. The presented framework is sufficiently general to also consider a class of inverse problems, and numerical examples of this type are also presented. © 2011 Elsevier B.V.
Resumo:
The innately highly efficient light-powered separation of charge that underpins natural photosynthesis can be exploited for applications in photoelectrochemistry by coupling nanoscale protein photoreaction centers to man-made electrodes. Planar photoelectrochemical cells employing purple bacterial reaction centers have been constructed that produce a direct current under continuous illumination and an alternating current in response to discontinuous illumination. The present work explored the basis of the open-circuit voltage (V(OC)) produced by such cells with reaction center/antenna (RC-LH1) proteins as the photovoltaic component. It was established that an up to ~30-fold increase in V(OC) could be achieved by simple manipulation of the electrolyte connecting the protein to the counter electrode, with an approximately linear relationship being observed between the vacuum potential of the electrolyte and the resulting V(OC). We conclude that the V(OC) of such a cell is dependent on the potential difference between the electrolyte and the photo-oxidized bacteriochlorophylls in the reaction center. The steady-state short-circuit current (J(SC)) obtained under continuous illumination also varied with different electrolytes by a factor of ~6-fold. The findings demonstrate a simple way to boost the voltage output of such protein-based cells into the hundreds of millivolts range typical of dye-sensitized and polymer-blend solar cells, while maintaining or improving the J(SC). Possible strategies for further increasing the V(OC) of such protein-based photoelectrochemical cells through protein engineering are discussed.
Resumo:
Over the last 50 years, the city of Venice, Italy, has observed a significant increase in the frequency of flooding. Numerous engineering solutions have been proposed, including the use of movable gates located at the three lagoon inlets. A key element in the prediction of performance is the estimation of settlements of the foundation system of the gates. The soils of Venice Lagoon are characterized by very erratic depositional patterns of clayey silts, resulting in an extremely heterogeneous stratigraphy with discontinuous layering. The soils are also characterized by varying contents of coarse and fine-grained particles. In contrast, the mineralogical composition of these deposits is quite uniform, which allows us to separate the influence of mineralogy from that of grain size distribution. A comprehensive geotechnical testing program was performed to assess the one-dimensional compression of Venice soils and examine the factors affecting the response in the transition from one material type to another. The compressibility of these natural silty clayey soils can be described by a single set of constitutive laws incorporating the relative fraction of granular to cohesive material. © 2007 ASCE.
Resumo:
This paper studies the dynamical response of a rotary drilling system with a drag bit, using a lumped parameter model that takes into consideration the axial and torsional vibration modes of the bit. These vibrations are coupled through a bit-rock interaction law. At the bit-rock interface, the cutting process introduces a state-dependent delay, while the frictional process is responsible for discontinuous right-hand sides in the equations governing the motion of the bit. This complex system is characterized by a fast axial dynamics compared to the slow torsional dynamics. A dimensionless formulation exhibits a large parameter in the axial equation, enabling a two-time-scales analysis that uses a combination of averaging methods and a singular perturbation approach. An approximate model of the decoupled axial dynamics permits us to derive a pseudoanalytical expression of the solution of the axial equation. Its averaged behavior influences the slow torsional dynamics by generating an apparent velocity weakening friction law that has been proposed empirically in earlier work. The analytical expression of the solution of the axial dynamics is used to derive an approximate analytical expression of the velocity weakening friction law related to the physical parameters of the system. This expression can be used to provide recommendations on the operating parameters and the drillstring or the bit design in order to reduce the amplitude of the torsional vibrations. Moreover, it is an appropriate candidate model to replace empirical friction laws encountered in torsional models used for control. © 2009 Society for Industrial and Applied Mathematics.
Resumo:
Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.
Resumo:
The aerodynamic design of turbomachinery presents the design optimisation community with a number of exquisite challenges. Chief among these are the size of the design space and the extent of discontinuity therein. This discontinuity can serve to limit the full exploitation of high-fidelity computational fluid dynamics (CFD): such codes require detailed geometric information often available only sometime after the basic configuration of the machine has been set by other means. The premise of this paper is that it should be possible to produce higher performing designs in less time by exploiting multi-fidelity techniques to effectively harness CFD earlier in the design process, specifically by facilitating its participation in configuration selection. The adopted strategy of local multi-fidelity correction, generated on demand, combined with a global search algorithm via an adaptive trust region is first tested on a modest, smooth external aerodynamic problem. Speed-up of an order of magnitude is demonstrated, comparable to established techniques applied to smooth problems. A number of enhancements aimed principally at effectively evaluating a wide range of configurations quickly is then applied to the basic strategy, and the emerging technique is tested on a generic aeroengine core compression system. A similar order of magnitude speed-up is achieved on this relatively large and highly discontinuous problem. A five-fold increase in the number of configurations assessed with CFD is observed. As the technique places constraints neither on the underlying physical modelling of the constituent analysis codes nor on first-order agreement between those codes, it has potential applicability to a range of multidisciplinary design challenges. © 2012 by Jerome Jarrett and Tiziano Ghisu.
Resumo:
The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.
Resumo:
The construction of protein-based photoelectrochemical cells that produce a variety of alternating currents in response to discontinuous illumination is reported. The photovoltaic component is a protein complex from the purple photosynthetic bacterium Rhodobacter sphaeroides which catalyses photochemical charge separation with a high quantum yield. Photoelectrochemical cells formed from this protein, a mobile redox mediator and a counter electrode formed from cobalt disilicide, titanium nitride, platinum, or multi-walled carbon nanotubes (MWCNT) generate a direct current during continuous illumination and an alternating current with different characteristics during discontinuous illumination. In particular, the use of superhydrophobic MWCNT as the back electrode results in a near symmetrical forward and reverse current upon light on and light off, respectively. The symmetry of the AC output of these cells is correlated with the wettability of the counter electrode. Potential applications of a hybrid biological/synthetic solar cell capable of generating an approximately symmetrical alternating current are discussed. Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.