Posicionamento aproximado do estado final para sistemas térmicos descritos pela equação do calor.

Neste trabalho, será considerado um problema de controle ótimo quadrático para a equação do calor em domínios retangulares com condição de fronteira do tipo Dirichlet é nos quais, a função de controle (dependente apenas no tempo) constitui um termo de fonte. Uma caracterização da solução ótima é obtida na forma de uma equação linear em um espaço de funções reais definidas no intervalo de tempo considerado. Em seguida, utiliza-se uma sequência de projeções em subespaços de dimensão finita para obter aproximações para o controle ótimo, o cada uma das quais pode ser gerada por um sistema linear de dimensão finita. A sequência de soluções aproximadas assim obtidas converge para a solução ótima do problema original. Finalmente, são apresentados resultados numéricos para domínios espaciais de dimensão 1.

Uma nova descrição para a transferência de massa em meios porosos com transição saturado-insaturado.

Esse texto trata do problema de um fluido contaminado escoando por um meio poroso, tratando os componentes na mistura como meios contínuos. Na primeira parte, desenvolvemos a teoria de misturas de meios contínuos e discutimos equações da continuidade, momento linear e momento angular. A seguir, descrevemos o problema em detalhe e fazemos hipóteses para simplificar o escoamento. Aplicamos as equações encontradas anteriormente para encontrarmos um sistema de equações diferenciais parciais. Desse ponto em diante, o problema se torna quase puramente matemático. Discutimos o caso insaturado, e depois a saturação do meio poroso. Finalmente, adicionamos um contaminante à mistura e, em seguida, N contaminantes.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Influence of a bow of finite width on bowed string motion: numerical modelling and experimental evidence

An enhanced physical model of the bowed string presented previously [1] is explored. It takes into account: the width of the bow, the angular motion of the string, bow-hair elasticity and string bending stiffness. The results of an analytical investigation of a model system - an infinite string sticking to a bow of finite width and driven on one side of the bow - are compared with experimental results published by Cremer [2] and reinterpreted here. Comparison shows that both the width of the bow and the bow-hair elasticity have a large impact on the reflection and transmission behaviour. In general, bending stiffness plays a minor role. Furthermore, a method of numerical simulation of the stiff string bowed with a bow of finite width is presented along with some preliminary results.

Building flexible user interfaces for solving PDEs

FEniCS is a collection of software tools for the automated solution of differential equations by finite element methods. In this note, we describe how FEniCS can be used to solve a simple nonlinear model problem with varying levels of automation. At one extreme, FEniCS provides tools for the fully automated and adaptive solution of nonlinear partial differential equations. At the other extreme, FEniCS provides a range of tools that allow the computational scientist to experiment with novel solution algorithms. © 2010 American Institute of Physics.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Effects of Lewis number on scalar variance transport in premixed flames

The influences of differential diffusion rates of heat and mass on the transport of the variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been studied using three-dimensional simplified chemistry based Direct Numerical Simulation (DNS) data of statistically planar turbulent premixed flames with global Lewis number ranging from Le = 0.34 to 1.2. The Lewis number effects on the statistical behaviours of the various terms of the transport equations of variances of Favre fluctuations of reaction progress variable and non-dimensional temperature have been analysed in the context of Reynolds Averaged Navier Stokes (RANS) simulations. It has been found that the turbulent fluxes of the progress variable and temperature variances exhibit counter-gradient transport for the flames with Lewis number significantly smaller than unity whereas the extent of this counter-gradient transport is found to decrease with increasing Lewis number. The Lewis number is also shown to have significant influences on the magnitudes of the chemical reaction and scalar dissipation rate contributions to the scalar variance transport. The modelling of the unclosed terms in the scalar variance equations for the non-unity Lewis number flames have been discussed in detail. The performances of the existing models for the unclosed terms are assessed based on a-priori analysis of DNS data. Based on the present analysis, new models for the unclosed terms of the active scalar variance transport equations are proposed, whenever necessary, which are shown to satisfactorily capture the behaviours of unclosed terms for all the flames considered in this study. © 2010 Springer Science+Business Media B.V.

The post-processed Galerkin method applied to non-linear shell vibrations

We present the results of a computational study of the post-processed Galerkin methods put forward by Garcia-Archilla et al. applied to the non-linear von Karman equations governing the dynamic response of a thin cylindrical panel periodically forced by a transverse point load. We spatially discretize the shell using finite differences to produce a large system of ordinary differential equations (ODEs). By analogy with spectral non-linear Galerkin methods we split this large system into a 'slowly' contracting subsystem and a 'quickly' contracting subsystem. We then compare the accuracy and efficiency of (i) ignoring the dynamics of the 'quick' system (analogous to a traditional spectral Galerkin truncation and sometimes referred to as 'subspace dynamics' in the finite element community when applied to numerical eigenvectors), (ii) slaving the dynamics of the quick system to the slow system during numerical integration (analogous to a non-linear Galerkin method), and (iii) ignoring the influence of the dynamics of the quick system on the evolution of the slow system until we require some output, when we 'lift' the variables from the slow system to the quick using the same slaving rule as in (ii). This corresponds to the post-processing of Garcia-Archilla et al. We find that method (iii) produces essentially the same accuracy as method (ii) but requires only the computational power of method (i) and is thus more efficient than either. In contrast with spectral methods, this type of finite-difference technique can be applied to irregularly shaped domains. We feel that post-processing of this form is a valuable method that can be implemented in computational schemes for a wide variety of partial differential equations (PDEs) of practical importance.

Optimal control with stochastic PDE constraints and uncertain controls

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.

Experimental and numerical study on AlGaInAs/AlGaAs distributed feedback lasers with GaInP gratings

Ridge-waveguide AlGaInAs/AlGaAs distributed feedback lasers with lattice-matched GaInP gratings were fabricated and their light-current characteristics, spectrum and far-field characteristics were measured. On the basis of our experimental results we analyze the effect of the electron stopper layer on light-current performance using the commercial laser simulation software PICS3D. The simulator is based on the self-consistent solution of drift diffusion equations, the Schrodinger equation, and the photon rate equation. The simulation results suggest that, with the use of a 80 nm-width p-doped Al0.6GaAs electron stopper layer, the slope efficiency can be increased and the threshold current can be reduced by more than 10 mA.

Mode-coupling analysis of three-dimensional microdisk resonators by the finite-difference time-domain technique

Quality factor enhancement due to mode coupling is observed in a three-dimensional microdisk resonator. The microdisk, which is vertically sandwiched between air and a substrate, with a radius of 1 mu m, a thickness of 0.2 mu m, and a refractive index of 3.4, is considered in a finite-difference time-domain (FDTD) numerical simulation. The mode quality factor of the fundamental mode HE71 decreases with an increase of the refractive index of the substrate, n(sub), from 2.0 to 3.17. However, the mode quality factor of the first-order mode HE72 reaches a peak value at n(sub) = 2.7 because of the mode coupling between the fundamental and the first-order modes. The variation of mode field distributions due to the mode coupling is also observed. This mechanism may be used to realize high-quality-factor modes in microdisks with high-refractive-index substrates. (c) 2006 Optical Society of America.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Hydroelastic dynamic characteristics of a slender axis-symmetric body

The slender axis-symmetric submarine body moving in the vertical plane is the object of our investigation. A coupling model is developed where displacements of a solid body as a Euler beam (consisting of rigid motions and elastic deformations) and fluid pressures are employed as basic independent variables, including the interaction between hydrodynamic forces and structure dynamic forces. Firstly the hydrodynamic forces, depending on and conversely influencing body motions, are taken into account as the governing equations. The expressions of fluid pressure are derived based on the potential theory. The characteristics of fluid pressure, including its components, distribution and effect on structure dynamics, are analyzed. Then the coupling model is solved numerically by means of a finite element method (FEM). This avoids the complicacy, combining CFD (fluid) and FEM (structure), of direct numerical simulation, and allows the body with a non-strict ideal shape so as to be more suitable for practical engineering. An illustrative example is given in which the hydroelastic dynamic characteristics, natural frequencies and modes of a submarine body are analyzed and compared with experimental results. Satisfactory agreement is observed and the model presented in this paper is shown to be valid.