A stress relaxation model for the viscoelastic solids based on the steady-state creep equation

Creep and stress relaxation are inherent mechanical behaviors of viscoelastic materials. It is considered that both are different performances of one identical physical phenomenon. The relationship between the decay stress and time during stress relaxation has been derived from the power law equation of the steady-state creep. The model was used to analyse the stress relaxation curves of various different viscoelastic materials (such as pure polycrystalline ice, polymers, foods, bones, metal, animal tissues, etc.). The calculated results using the theoretical model agree with the experimental data very well. Here we show that the new mathematical formula is not only simple but its parameters have the clear physical meanings. It is suitable to materials with a very broad scope and has a strong predictive ability.

Data assimilation for state and parameter estimation: application to morphodynamic modelling

Data assimilation is predominantly used for state estimation; combining observational data with model predictions to produce an updated model state that most accurately approximates the true system state whilst keeping the model parameters fixed. This updated model state is then used to initiate the next model forecast. Even with perfect initial data, inaccurate model parameters will lead to the growth of prediction errors. To generate reliable forecasts we need good estimates of both the current system state and the model parameters. This paper presents research into data assimilation methods for morphodynamic model state and parameter estimation. First, we focus on state estimation and describe implementation of a three dimensional variational(3D-Var) data assimilation scheme in a simple 2D morphodynamic model of Morecambe Bay, UK. The assimilation of observations of bathymetry derived from SAR satellite imagery and a ship-borne survey is shown to significantly improve the predictive capability of the model over a 2 year run. Here, the model parameters are set by manual calibration; this is laborious and is found to produce different parameter values depending on the type and coverage of the validation dataset. The second part of this paper considers the problem of model parameter estimation in more detail. We explain how, by employing the technique of state augmentation, it is possible to use data assimilation to estimate uncertain model parameters concurrently with the model state. This approach removes inefficiencies associated with manual calibration and enables more effective use of observational data. We outline the development of a novel hybrid sequential 3D-Var data assimilation algorithm for joint state-parameter estimation and demonstrate its efficacy using an idealised 1D sediment transport model. The results of this study are extremely positive and suggest that there is great potential for the use of data assimilation-based state-parameter estimation in coastal morphodynamic modelling.

Ensemble forecasts of a flood-producing storm: comparison of the influence of model-state perturbations and parameter modifications

High-resolution ensemble simulations (Δx = 1 km) are performed with the Met Office Unified Model for the Boscastle (Cornwall, UK) flash-flooding event of 16 August 2004. Forecast uncertainties arising from imperfections in the forecast model are analysed by comparing the simulation results produced by two types of perturbation strategy. Motivated by the meteorology of the event, one type of perturbation alters relevant physics choices or parameter settings in the model's parametrization schemes. The other type of perturbation is designed to account for representativity error in the boundary-layer parametrization. It makes direct changes to the model state and provides a lower bound against which to judge the spread produced by other uncertainties. The Boscastle has genuine skill at scales of approximately 60 km and an ensemble spread which can be estimated to within ∼ 10% with only eight members. Differences between the model-state perturbation and physics modification strategies are discussed, the former being more important for triggering and the latter for subsequent cell development, including the average internal structure of convective cells. Despite such differences, the spread in rainfall evaluated at skilful scales is shown to be only weakly sensitive to the perturbation strategy. This suggests that relatively simple strategies for treating model uncertainty may be sufficient for practical, convective-scale ensemble forecasting.

Atmospheric blocking in a high resolution climate model: influences of mean state, orography and eddy forcing

An underestimate of atmospheric blocking occurrence is a well-known limitation of many climate models. This article presents an analysis of Northern Hemisphere winter blocking in an atmospheric model with increased horizontal resolution. European blocking frequency increases with model resolution, and this results from an improvement in the atmospheric patterns of variability as well as a simple improvement in the mean state. There is some evidence that the transient eddy momentum forcing of European blocks is increased at high resolution, which could account for this. However, it is also shown that the increase in resolution of the orography is needed to realise the improvement in blocking, consistent with the increase in height of the Rocky Mountains acting to increase the tilt of the Atlantic jet stream and giving higher mean geopotential heights over northern Europe. Blocking frequencies in the Pacific sector are also increased with atmospheric resolution, but in this case the improvement in orography actually leads to a decrease in blocking

Changes in reactive stratospheric gases due to a change in Brewer–Dobson circulation: results from a simple model

Amounts of source gases with stratospheric sinks (CFCs, N2O, CH4) are affected by changes in Brewer–Dobson circulation. Source gases and their degradation products are important for atmospheric chemistry and climate. With a simple model, we examine how amounts and lifetimes of source gases and products depend on speed of the circulation. Transient results differ from steady-state and stratospheric results differ from those for stratosphere plus troposphere. Increases in speed increase the stratospheric burden of source gases, but reduce products and reduce total burdens and lifetimes of source gases

Present state of global wetland extent and wetland methane modelling: conclusions from a model inter-comparison project (WETCHIMP)

Global wetlands are believed to be climate sensitive, and are the largest natural emitters of methane (CH4). Increased wetland CH4 emissions could act as a positive feedback to future warming. The Wetland and Wetland CH4 Inter-comparison of Models Project (WETCHIMP) investigated our present ability to simulate large-scale wetland characteristics and corresponding CH4 emissions. To ensure inter-comparability, we used a common experimental protocol driving all models with the same climate and carbon dioxide (CO2) forcing datasets. The WETCHIMP experiments were conducted for model equilibrium states as well as transient simulations covering the last century. Sensitivity experiments investigated model response to changes in selected forcing inputs (precipitation, temperature, and atmospheric CO2 concentration). Ten models participated, covering the spectrum from simple to relatively complex, including models tailored either for regional or global simulations. The models also varied in methods to calculate wetland size and location, with some models simulating wetland area prognostically, while other models relied on remotely sensed inundation datasets, or an approach intermediate between the two. Four major conclusions emerged from the project. First, the suite of models demonstrate extensive disagreement in their simulations of wetland areal extent and CH4 emissions, in both space and time. Simple metrics of wetland area, such as the latitudinal gradient, show large variability, principally between models that use inundation dataset information and those that independently determine wetland area. Agreement between the models improves for zonally summed CH4 emissions, but large variation between the models remains. For annual global CH4 emissions, the models vary by ±40% of the all-model mean (190 Tg CH4 yr−1). Second, all models show a strong positive response to increased atmospheric CO2 concentrations (857 ppm) in both CH4 emissions and wetland area. In response to increasing global temperatures (+3.4 °C globally spatially uniform), on average, the models decreased wetland area and CH4 fluxes, primarily in the tropics, but the magnitude and sign of the response varied greatly. Models were least sensitive to increased global precipitation (+3.9 % globally spatially uniform) with a consistent small positive response in CH4 fluxes and wetland area. Results from the 20th century transient simulation show that interactions between climate forcings could have strong non-linear effects. Third, we presently do not have sufficient wetland methane observation datasets adequate to evaluate model fluxes at a spatial scale comparable to model grid cells (commonly 0.5°). This limitation severely restricts our ability to model global wetland CH4 emissions with confidence. Our simulated wetland extents are also difficult to evaluate due to extensive disagreements between wetland mapping and remotely sensed inundation datasets. Fourth, the large range in predicted CH4 emission rates leads to the conclusion that there is both substantial parameter and structural uncertainty in large-scale CH4 emission models, even after uncertainties in wetland areas are accounted for.

A simple method for integrating a complex model into an ensemble data assimilation system using MPI

This paper details a strategy for modifying the source code of a complex model so that the model may be used in a data assimilation context, {and gives the standards for implementing a data assimilation code to use such a model}. The strategy relies on keeping the model separate from any data assimilation code, and coupling the two through the use of Message Passing Interface (MPI) {functionality}. This strategy limits the changes necessary to the model and as such is rapid to program, at the expense of ultimate performance. The implementation technique is applied in different models with state dimension up to $2.7 \times 10^8$. The overheads added by using this implementation strategy in a coupled ocean-atmosphere climate model are shown to be an order of magnitude smaller than the addition of correlated stochastic random errors necessary for some nonlinear data assimilation techniques.

A joint state and parameter estimation scheme for nonlinear dynamical systems

We present a novel algorithm for concurrent model state and parameter estimation in nonlinear dynamical systems. The new scheme uses ideas from three dimensional variational data assimilation (3D-Var) and the extended Kalman filter (EKF) together with the technique of state augmentation to estimate uncertain model parameters alongside the model state variables in a sequential filtering system. The method is relatively simple to implement and computationally inexpensive to run for large systems with relatively few parameters. We demonstrate the efficacy of the method via a series of identical twin experiments with three simple dynamical system models. The scheme is able to recover the parameter values to a good level of accuracy, even when observational data are noisy. We expect this new technique to be easily transferable to much larger models.

Conditioning and preconditioning of the optimal state estimation problem

Optimal state estimation is a method that requires minimising a weighted, nonlinear, least-squares objective function in order to obtain the best estimate of the current state of a dynamical system. Often the minimisation is non-trivial due to the large scale of the problem, the relative sparsity of the observations and the nonlinearity of the objective function. To simplify the problem the solution is often found via a sequence of linearised objective functions. The condition number of the Hessian of the linearised problem is an important indicator of the convergence rate of the minimisation and the expected accuracy of the solution. In the standard formulation the convergence is slow, indicating an ill-conditioned objective function. A transformation to different variables is often used to ameliorate the conditioning of the Hessian by changing, or preconditioning, the Hessian. There is only sparse information in the literature for describing the causes of ill-conditioning of the optimal state estimation problem and explaining the effect of preconditioning on the condition number. This paper derives descriptive theoretical bounds on the condition number of both the unpreconditioned and preconditioned system in order to better understand the conditioning of the problem. We use these bounds to explain why the standard objective function is often ill-conditioned and why a standard preconditioning reduces the condition number. We also use the bounds on the preconditioned Hessian to understand the main factors that affect the conditioning of the system. We illustrate the results with simple numerical experiments.

The dynamics of cargo driven by molecular motors in the context of asymmetric simple exclusion processes

We consider the dynamics of cargo driven by a collection of interacting molecular motors in the context of ail asymmetric simple exclusion process (ASEP). The model is formulated to account for (i) excluded-volume interactions, (ii) the observed asymmetry of the stochastic movement of individual motors and (iii) interactions between motors and cargo. Items (i) and (ii) form the basis of ASEP models and have already been considered to study the behavior of motor density profile [A. Parmeggiani. T. Franosch, E. Frey, Phase Coexistence in driven one-dimensional transport, Phys. Rev. Lett. 90 (2003) 086601-1-086601-4]. Item (iii) is new. It is introduced here as an attempt to describe explicitly the dependence of cargo movement on the dynamics of motors in this context. The steady-state Solutions Of the model indicate that the system undergoes a phase transition of condensation type as the motor density varies. We study the consequences of this transition to the behavior of the average cargo velocity. (C) 2009 Elsevier B.V. All rights reserved.

Time-dependent or state-dependent pricing? : evidence from a large currency devaluation episode

Duas classes de modelos buscam explicar o padrão de ajustamento de preço das firmas: modelos tempo-dependente e estado-dependente. O objetivo deste trabalho é levantar algumas evidencias empíricas de modo a distinguir os modelos, ou seja, identificar de que maneira as firmas realmente precificam. Para isso, escolheu-se a grande desvalorização cambial de 1999 como principal ferramenta e ambiente de análise. A hipótese fundamental é que o choque cambial impacta significativamente o custo de algumas indústrias, em alguns casos induzindo-as a alterarem seus preço após o choque. A partir de uma imensa base de micro dados formada por preços que compõem o CPI, algumas estimações importantes como a probabilidade e a magnitude média das trocas foram levantadas. A magnitude é dada por uma média simples, enquanto a probabilidade é estimada pelo método da máxima verossimilhança. Os resultados indicam um comportamento de precificação similar ao proposto por modelos estado-dependente.

Separating interface state response from parasitic effects in conductance measurements on organic metal-insulator-semiconductor capacitors

A simple model is developed for the admittance of a metal-insulator-semiconductor (MIS) capacitor which includes the effect of a guard ring surrounding the Ohmic contact to the semiconductor. The model predicts most of the features observed in a MIS capacitor fabricated using regioregular poly(3-hexylthiophene) as the active semiconductor and polysilsesquioxane as the gate insulator. In particular, it shows that when the capacitor is driven into accumulation, the parasitic transistor formed by the guard ring and Ohmic contact can give rise to an additional feature in the admittance-voltage plot that could be mistaken for interface states. When this artifact and underlying losses in the bulk semiconductor are accounted for, the remaining experimental feature, a peak in the loss-voltage plot when the capacitor is in depletion, is identified as an interface (or near interface) state of density of similar to 4 x 10(10) cm(-2) eV(-1). Application of the model shows that exposure of a vacuum-annealed device to laboratory air produces a rapid change in the doping density in the channel region of the parasitic transistor but only slow changes in the bulk semiconductor covered by the gold Ohmic contact. (C) 2008 American Institute of Physics.

Prevalence of intestinal parasites in dogs from São Paulo State, Brazil

The prevalence of gastrointestinal parasites in stray dogs, and dogs with owners was investigated by fecal examinations from 271 dogs employing sedimentation, simple flotation and centrifugation-flotation methods. The centrifugation-flotation method, when compared to simple flotation or sedimentation methods was generally more accurate in the diagnosis of all intestinal parasites, but statistical differences were detected only in relation to Giardia spp. and Cystoisospora spp. (synonym Isospora spp.). The following parasites, with their respective prevalence, were diagnosed in the fecal samples: Ancylostoma spp. (23.6%); Toxocara canis (5.5%); Trichuris vulpis (4.8%); Spirocerca lupi (1.9%); Dipylidium caninum (0.7%); Giardia spp. (12.2%); Hammondia heydorni (2.6%); Cystoisospora spp. (8.5%); and Sarcocystis spp. (2.2%). The prevalence of most parasites was similar for dogs of mixed-breed and for dogs of a defined-breed, except for Cystoisospora spp. and T canis which showed a significantly higher prevalence in mixed-breed dogs. The prevalence of Ancylostoma spp. (17.1%) was significantly lower in stray dogs than in those with an owner (31.9%) and the prevalence of Giardia spp. and Cystoisospora spp. was higher in stray dogs (P < 0.05). No effect of season on the occurrence of the different parasite genera could be observed, except for Ancylostoma spp., for which an increase in the percentage of dogs shedding eggs was observed at the beginning of Summer with a peak occurrence during April and May (Autumn). The prevalence of Ancylostoma spp., T canis, T vulpis, Giardia spp. and Cystoisospora spp. was higher in adult males than in adult females, but significant differences between the two groups occurred only with Giardia spp. Young animals were found to more frequently shed Nematode eggs in feces than adult animals. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Production of xylanase and CMCase on solid state fermentation in different residues by Thermoascus aurantiacus miehe

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

Water quality index as a simple indicator of aquaculture effects on aquatic bodies

This paper proposes a water quality index (WQI) to subsidize management actions in the Medio Paranapanema Watershed in São Paulo State, Brazil, as a simple pollution indicator for aquaculture activity. Water quality of the Macuco and Queixada rivers was investigated for 2 years (from May 2003 to May 2005). The index proposed in this work is composed of three measurable environmental parameters-turbidity, total phosphorus and dissolved oxygen. Concentrations of these three variables were normalized on a scale from 0 to 100 and translated into statements of water quality (excellent, good, regular, fair and poor). The index was applied to seventeen monitoring points in the aquatic bodies described above and compared to others, one being that used by the Environmental Protection Agency of United States and proposed for the National Sanitation Foundation, other employing minimal index and the last one considering the minimum operator concept. The results show that the degradation in this watershed from aquaculture activity can be easily inferred with this index, which is more restricted than the others routinely used to infer water quality. (C) 2007 Elsevier Ltd. All rights reserved.