A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Short-term temporal changes of soil carbon losses after tillage described by a first-order decay model

Tillage stimulates soil carbon (C) losses by increasing aeration, changing temperature and moisture conditions, and thus favoring microbial decomposition. In addition, soil aggregate disruption by tillage exposes once protected organic matter to decomposition. We propose a model to explain carbon dioxide (CO2) emission after tillage as a function of the no-till emission plus a correction due to the tillage disturbance. The model assumes that C in the readily decomposable organic matter follows a first-order reaction kinetics equation as: dC(sail)(t)/dt = -kC(soil)(t) and that soil C-CO2 emission is proportional to the C decay rate in soil, where C-soil(t) is the available labile soil C (g m(-2)) at any time (t). Emissions are modeled in terms soil C available to decomposition in the tilled and non-tilled plots, and a relationship is derived between no-till (F-NT) and tilled (F-Gamma) fluxes, which is: F-T = a1F(NT)e(-a2t), where t is time after tillage. Predicted and observed fluxes showed good agreement based on determination coefficient (R-2), index of agreement and model efficiency, with R-2 as high as 0.97. The two parameters included in the model are related to the difference between the decay constant (k factor) of tilled and no-till plots (a(2)) and also to the amount of labile carbon added to the readily decomposable soil organic matter due to tillage (a,). These two parameters were estimated in the model ranging from 1.27 and 2.60 (a(1)) and - 1.52 x 10(-2) and 2.2 x 10(-2) day(-1) (a(2)). The advantage is that temporal variability of tillage-induced emissions can be described by only one analytical function that includes the no-till emission plus an exponential term modulated by tillage and environmentally dependent parameters. (C) 2008 Elsevier B.V. All rights reserved.

Short-term temporal changes of bare soil CO2 fluxes after tillage described by first-order decay models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Otimização topológica de estruturas termoelásticas tridimensionais

This work presents an optimization technique based on structural topology optimization methods, TOM, designed to solve problems of thermoelasticity 3D. The presented approach is based on the adjoint method of sensitivity analysis unified design and is intended to loosely coupled thermomechanical problems. The technique makes use of analytical expressions of sensitivities, enabling a reduction in the computational cost through the use of a coupled field adjoint equation, defined in terms the of temperature and displacement fields. The TOM used is based on the material aproach. Thus, to make the domain is composed of a continuous distribution of material, enabling the use of classical models in nonlinear programming optimization problem, the microstructure is considered as a porous medium and its constitutive equation is a function only of the homogenized relative density of the material. In this approach, the actual properties of materials with intermediate densities are penalized based on an artificial microstructure model based on the SIMP (Solid Isotropic Material with Penalty). To circumvent problems chessboard and reduce dependence on layout in relation to the final optimal initial mesh, caused by problems of numerical instability, restrictions on components of the gradient of relative densities were applied. The optimization problem is solved by applying the augmented Lagrangian method, the solution being obtained by applying the finite element method of Galerkin, the process of approximation using the finite element Tetra4. This element has the ability to interpolate both the relative density and the displacement components and temperature. As for the definition of the problem, the heat load is assumed in steady state, i.e., the effects of conduction and convection of heat does not vary with time. The mechanical load is assumed static and distributed

A bounded upwinding scheme for computing convection-dominated transport problems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Processos difusivos generalizados

We investigate several diffusion equations which extend the usual one by considering the presence of nonlinear terms or a memory effect on the diffusive term. We also considered a spatial time dependent diffusion coefficient. For these equations we have obtained a new classes of solutions and studied the connection of them with the anomalous diffusion process. We start by considering a nonlinear diffusion equation with a spatial time dependent diffusion coefficient. The solutions obtained for this case generalize the usual one and can be expressed in terms of the q-exponential and q-logarithm functions present in the generalized thermostatistics context (Tsallis formalism). After, a nonlinear external force is considered. For this case the solutions can be also expressed in terms of the q-exponential and q-logarithm functions. However, by a suitable choice of the nonlinear external force, we may have an exponential behavior, suggesting a connection with standard thermostatistics. This fact reveals that these solutions may present an anomalous relaxation process and then, reach an equilibrium state of the kind Boltzmann- Gibbs. Next, we investigate a nonmarkovian linear diffusion equation that presents a kernel leading to the anomalous diffusive process. Particularly, our first choice leads to both a the usual behavior and anomalous behavior obtained through a fractionalderivative equation. The results obtained, within this context, correspond to a change in the waiting-time distribution for jumps in the formalism of random walks. These modifications had direct influence in the solutions, that turned out to be expressed in terms of the Mittag-Leffler or H of Fox functions. In this way, the second moment associated to these distributions led to an anomalous spread of the distribution, in contrast to the usual situation where one finds a linear increase with time

Controlling dispersal dynamics of Aedes aegypti

The dengue virus is transmitted in regions previously infested with the mosquito Aedes aegypti. To assess the spreading and establishment of the dengue disease vector, a mathematical model is developed that takes into account the diffusion and advection phenomena. A discrete model based on the cellular automata approach, which is a good framework to deal with small populations, is also developed to be compared with the continuous modeling.

Traveling waves in the Lethargic Crab Disease

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modal analysis of anisotropic diffused-channel waveguide by a scalar finite element method

A procedure to model optical diffused-channel waveguides is presented in this work. The dielectric waveguides present anisotropic refractive indexes which are calculated from the proton concentration. The proton concentration inside the channel is calculated by the anisotropic 2D-linear diffusion equation and converted to the refractive indexes using mathematical relations obtained from experimental data, the arbitrary refractive index profile is modeled by a. nodal expansion in the base functions. The TE and TM-like propagation properties (effective index) and the electromagnetic fields for well-annealed proton-exchanged (APE) LiNbO3 waveguides are computed by the finite element method.

Periodic waves and solitons in a nonlinear fibre with resonant impurities

We shall consider a coupled nonlinear Schrodinger equation- Bloch system of equations describing the propagation of a single pulse through a nonlinear dispersive waveguide in the presence of resonances; this could be, for example, a doped optical fibre. By making use of the integrability of the dynamic equations, we shall apply the finite-gap integration method to obtain periodic solutions for this system. Next, we consider the problem of the formation of solitons at a sharp front pulse and, by means of the Whitham modulational theory, we derive the amplitude and velocity of the largest soliton.

Spontaneous symmetry breaking of Bose-Fermi mixtures in double-well potentials

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

First-order decay models to describe soil C-CO2 Loss after rotary tillage

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Drying of persimmon: Mathematical model for diffusivity as a simultaneous function of moisture content and shrinkage

Fresh persimmon has a high moisture content (about 85% wet basis) making it highly perishable and requiring adequate drying conditions to obtain an acceptable dehydrated product. Drying kinetics of persimmon cv. Rama Forte was studied in a fixed bed dryer at temperatures ranging from 50 to 80 degreesC and air velocity of 0.8 m/s. Shrinkage during drying was described by a linear correlation with respect to water content. Evaluation of effective diffusivity as a function of moisture content, with undergoing shrinkage during drying was based on Fourier series solution of Fick's diffusion equation. Effective diffusivity values at moisture contents between 0.09 - 4.23 kg water/kg dry matter were found to be in the range of 2.6 x 10(-10) m(2)/s to 5.4 x 10(-10) m(2)/s, and its dependence on air drying temperature was represented by an Arrhenius type equation. Activation energy increased with decreasing water content in persimmons.

Pressure field along the axis of a high-power klystron amplifier

The pressure field in a high-power klystron amplifier is investigated to scale the ionic vacuum pump used to maintain the ultra high-vacuum in the device in order to increase its life-time. The investigation is conducted using an 1.3 GHz, 100 A - 240 keV high-power klystron with five reentrant coaxial cavities, assembled in a cylindrical drift tube 1.2 m long. The diffusion equation is solved to the regime molecular flow to obtain the pressure profile along the axis of the klystron drift tube. The model, solved by both analytical and numerical procedures, is able to determine the pressure values in steady-state case. This work considers the specific conductance and all important gas sources, as in the degassing of the drift tube and cavities walls, cathode, and collector. For the drift tube degassing rate equals to q(deg) = 2x10(-12) (-)mbar.L.s(-1) cm(-2) (degassing rate per unit area), to cavities q(cavity) = 3x10(-13) mbar.L.s(-1)cm(-2), to the cathode q(cathode) = 6x10(-9)_mbar.L.s(-1) and to the collector q(collector) = 6x10(-9) mbar.L.s(-1), it was found that a 10 L.s(-1) ionic vacuum pump connected in the output waveguide wall is suitable. In this case, the pressure obtained in the cathode is p(cathode) = 6.3x10(-9) mbar, in the collector p(collector) = 2.7x10(-9) mbar, and in the output waveguide p = 2.1x10(-9) mbar. Although only the steady-state case is analyzed, some aspects that may be relevant in a transient situation, for instance, when the beam hits the drift tube walls, producing a gas burst, is also commented.

Cephalosporin C production by immobilized Cephalosporium acremonium cells in a repeated batch tower bioreactor

The industrial production of antibiotics with filamentous fungi is usually carried out in conventional aerated and agitated tank fermentors. Highly viscous non-Newtonian broths are produced and a compromise must be found between convenient shear stress and adequate oxygen transfer. In this work, cephalosporin C production by bioparticles of immobilized cells of Cephalosporium acremonium ATCC 48272 was studied in a repeated batch tower bioreactor as an alternative to the conventional process. Also, gas-liquid oxygen transfer volumetric coefficients, k(L)a, were determined at various air flow-rates and alumina contents in the bioparticle. The bioparticles were composed of calcium alginate (2.0% w/w), alumina (<44 micra), cells, and water. A model describing the cell growth, cephalosporin C production, oxygen, glucose, and sucrose consumption was proposed. To describe the radial variation of oxygen concentration within the pellet, the reaction-diffusion model forecasting a dead core bioparticle was adopted. The k(L)a measurements with gel beads prepared with 0.0, 1.0, 1.5, and 2.0% alumina showed that a higher k(L)a value is attained with 1.5 and 2.0%. An expression relating this coefficient to particle density, liquid density, and air velocity was obtained and further utilized in the simulation of the proposed model. Batch, followed by repeated batch experiments, were accomplished by draining the spent medium, washing with saline solution, and pouring fresh medium into the bioreactor. Results showed that glucose is consumed very quickly, within 24 h, followed by sucrose consumption and cephalosporin C production. Higher productivities were attained during the second batch, as cell concentration was already high, resulting in rapid glucose consumption and an early derepression of cephalosporin C synthesizing enzymes. The model incorporated this improvement predicting higher cephalosporin C productivity. (C) 2004 Wiley Periodicals, Inc.