Inclusion of nonlinear demand-supply relationships within large-scale partial equilibrium linear programming models

This paper presents a new method for the inclusion of nonlinear demand and supply relationships within a linear programming model. An existing method for this purpose is described first and its shortcomings are pointed out before showing how the new approach overcomes those difficulties and how it provides a more accurate and 'smooth' (rather than a kinked) approximation of the nonlinear functions as well as dealing with equilibrium under perfect competition instead of handling just the monopolistic situation. The workings of the proposed method are illustrated by extending a previously available sectoral model for the UK agriculture.

Assays for enhanced activity of low efficacy partial agonists at the D-2 dopamine receptor

Background and purpose: Low efficacy partial agonists at the D-2 dopamine receptor may be useful for treating schizophrenia. In this report we describe a method for assessing the efficacy of these compounds based on stimulation of [S-35]GTP gamma S binding. Experimental approach: Agonist efficacy was assessed from [S-35]GTP gamma S binding to membranes of CHO cells expressing D2 dopamine receptors in buffers with and without Na+. Effects of Na+ on receptor/G protein coupling were assessed using agonist/[H-3] spiperone competition binding assays. Key results: When [S-35]GTP gamma S binding assays were performed in buffers containing Na+, some agonists (aripiprazole, AJ-76, UH-232) exhibited very low efficacy whereas other agonists exhibited measurable efficacy. When Na+ was substituted by N-methyl D-glucamine, the efficacy of all agonists increased (relative to that of dopamine) but particularly for aripiprazole, aplindore, AJ-76, (-)-3-PPP and UH-232. In ligand binding assays, substitution of Na+ by N-methyl D-glucamine increased receptor/G protein coupling for some agonists -. aplindore, dopamine and (-)-3-PPP-but for aripiprazole, AJ-76 and UH-232 there was little effect on receptor/G protein coupling. Conclusions and implications: Substitution of Na+ by NMDG increases sensitivity in [S-35] GTPgS binding assays so that very low efficacy agonists were detected clearly. For some agonists the effect seems to be mediated via enhanced receptor/G protein coupling whereas for others the effect is mediated at another point in the G protein activation cycle. AJ-76, aripiprazole and UH-232 seem particularly sensitive to this change in assay conditions. This work provides a new method to discover these very low efficacy agonists.

Velocity-based moving mesh methods for nonlinear partial differential equations

This article describes a number of velocity-based moving mesh numerical methods formultidimensional nonlinear time-dependent partial differential equations (PDEs). It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation. Finite element algorithms are derived for both mass-conserving and non mass-conserving problems, and results shown for a number of multidimensional nonlinear test problems, including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem. Further applications and extensions are referenced.

GSHMC: an efficient method for molecular simulation

The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.

A modified PNN algorithm with optimal PD modeling using the orthogonal least squares method

In this paper a modified algorithm is suggested for developing polynomial neural network (PNN) models. Optimal partial description (PD) modeling is introduced at each layer of the PNN expansion, a task accomplished using the orthogonal least squares (OLS) method. Based on the initial PD models determined by the polynomial order and the number of PD inputs, OLS selects the most significant regressor terms reducing the output error variance. The method produces PNN models exhibiting a high level of accuracy and superior generalization capabilities. Additionally, parsimonious models are obtained comprising a considerably smaller number of parameters compared to the ones generated by means of the conventional PNN algorithm. Three benchmark examples are elaborated, including modeling of the gas furnace process as well as the iris and wine classification problems. Extensive simulation results and comparison with other methods in the literature, demonstrate the effectiveness of the suggested modeling approach.

Identification of partial resetting using De as a function of illumination time

Modern age samples from various depositional environments were examined for signal resetting. For 19 modern aeolian/beach samples all De values obtained were View the MathML source, with ∼70% having View the MathML source. For 21 fluvial/colluvial samples, all De values were View the MathML source with ∼80% being View the MathML source. De as a function of illumination (OSL measurement) time (De(t)) plots were examined for all samples. Based on previous laboratory experiments, increases in De(t) were expected for partially reset samples, and constant De(t) for fully reset samples. All aeolian samples, both modern age and additional ‘young’ samples (<1000 years), showed constant (flat) De(t) while all modern, non-zero De, fluvial/colluvial samples showed increasing De(t). ‘Replacement plots’, where a regenerated signal is substituted for the natural, yielded constant (flat) De(t). These findings support strongly the use of De(t) as a method of identifying incomplete resetting in fluvial samples. Potential complicating factors, such as illumination (bleaching) spectrum, thermal instability and component composition are discussed and a series of internal checks on the applicability of the De(t) for each individual aliquot/grain level are outlined.

Modified section method for laser-welding of ill-fitting cp Ti and Ni-Cr alloy one-piece cast implant-supported frameworks

P>This study aimed to verify the effect of modified section method and laser-welding on the accuracy of fit of ill-fitting commercially pure titanium (cp Ti) and Ni-Cr alloy one-piece cast frameworks. Two sets of similar implant-supported frameworks were constructed. Both groups of six 3-unit implant-supported fixed partial dentures were cast as one-piece [I: Ni-Cr (control) and II: cp Ti] and evaluated for passive fitting in an optical microscope with both screws tightened and with only one screw tightened. All frameworks were then sectioned in the diagonal axis at the pontic region (III: Ni-Cr and IV: cp Ti). Sectioned frameworks were positioned in the matrix (10-Ncm torque) and laser-welded. Passive fitting was evaluated for the second time. Data were submitted to anova and Tukey-Kramer honestly significant difference tests (P < 0 center dot 05). With both screws tightened, one-piece cp Ti group II showed significantly higher misfit values (27 center dot 57 +/- 5 center dot 06 mu m) than other groups (I: 11 center dot 19 +/- 2 center dot 54 mu m, III: 12 center dot 88 +/- 2 center dot 93 mu m, IV: 13 center dot 77 +/- 1 center dot 51 mu m) (P < 0 center dot 05). In the single-screw-tightened test, with readings on the opposite side to the tightened side, Ni-Cr cast as one-piece (I: 58 center dot 66 +/- 14 center dot 30 mu m) was significantly different from cp Ti group after diagonal section (IV: 27 center dot 51 +/- 8 center dot 28 mu m) (P < 0 center dot 05). On the tightened side, no significant differences were found between groups (P > 0 center dot 05). Results showed that diagonally sectioning ill-fitting cp Ti frameworks lowers misfit levels of prosthetic implant-supported frameworks and also improves passivity levels of the same frameworks when compared to one-piece cast structures.

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.

On the statistical estimation of diffusion processes - a partial survey

Data available on continuous-time diffusions are always sampled discretely in time. In most cases, the likelihood function of the observations is not directly computable. This survey covers a sample of the statistical methods that have been developed to solve this problem. We concentrate on some recent contributions to the literature based on three di§erent approaches to the problem: an improvement of the Euler-Maruyama discretization scheme, the employment of Martingale Estimating Functions, and the application of Generalized Method of Moments (GMM).

A Numerical method for the semilinear stochastic transport equation

Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014.

Estimation of genetic parameters for partial egg production periods by means of random regression models

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

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Influence of Ridge Inclination and Implant Localization on the Association of Mandibular Kennedy Class I Removable Partial Denture

The aim of this study was to evaluate the tendency of displacement of the supporting structures of the distal extension removable partial denture (DERPD) associated to the implant with different inclinations of alveolar ridge and implant localizations through a two-dimensional finite-element method. Sixteen mandibular models were fabricated, presenting horizontal, distally descending, distally ascending, or descending-ascending ridges. All models presented the left canine and were rehabilitated with conventional DERPD or implant-retained prosthesis with the ERA system. The models were obtained by the AutoCAD software and transferred to the finite-element software ANSYS 9.0 for analysis. A force of 50 N was applied on the cusp tips of the teeth, with 5 points of loading of 10 N. The results were visualized by displacement maps. For all ridge inclinations, the assembly of the DERPD with distal plate retained by an anterior implant exhibited the lowest requisition of the supporting structures. The highest tendency of displacement occurred in the model with distally ascending ridge with incisal rest. It was concluded that the association of the implant decreased the displacement of the DERPD, and the anterior positioning of the implant associated to the DERPD with the distal plate preserved the supporting structures for all ridges.

Influence of Implant Inclination Associated With Mandibular Class I Removable Partial Denture

The aim of this study was to use two-dimensional finite element method to evaluate the displacement and stress distribution transmitted by a distal extension removable partial denture (DERPD) associated with an implant placed at different inclinations (0, 5, 15, and 30 degrees) in the second molar region of the edentulous mandible ridge. Six hemimandibular models were created: model A, only with the presence of the natural tooth 33; model B, similar to model A, with the presence of a conventional DERPD replacing the missing teeth; model C, similar to the previous model, with a straight implant (0 degrees) in the distal region of the ridge, under the denture base; model D, similar to model C, with the implant angled at 5 degrees in the mesial direction; model E, similar to model C, with the implant angled at 15 degrees in the mesial direction; and model F, similar to ME, with the implant angled at 30 degrees in the mesial direction. The models were created with the use of the AutoCAD 2000 program (Autodesk, Inc, San Rafael, CA) and processed for finite element analysis by the ANSYS 8.0 program (Swanson Analysis Systems, Houston, PA). The force applied was vertical of 50 N on each cusp tip. The results showed that the introduction of the RPD overloaded the supporting structures of the RPD and that the introduction of the implant helped to relieve the stresses of the mucosa alveolar, cortical bone, and trabecular bone. The best stress distribution occurred in model D with the implant angled at 5 degrees. The use of an implant as a support decreased the displacement of alveolar mucosa for all inclinations simulated. The stress distribution transmitted by the DERPD to the supporting structures was improved by the use of straight or slightly inclined implants. According to the displacement analysis and von Mises stress, it could be expected that straight or slightly inclined implants do not represent biomechanical risks to use.

Evaluation of the influence of location of osseointegrated implants associated with mandibular removable partial dentures

Purpose: The aim of this research was to assess, by means of, the bi-dimensional finite element method, the best implant location in the alveolar edge, through stress distribution and support structure displacement of a distal extension removable partial denture associated with an osseointegrated implant of 10.0 x .75 mm, acting as abutment for the denture base.Methods and Materials: Five models in sagittal cut were used to represent: model A-hemi arch containing natural tooth 33 and the distal alveolar edge; model B-similar to model A, but with a conventional removable partial denture to replace the absent teeth; model C (MC)-similar to the previous one, with an implant in the distal region of the edge under the denture base; model D-similar to MC, with the implant in the central region of the edge; model E-similar to MC, with an implant in the mesial region of the edge. With the aid of the finite element program ANSYS 8.0, the models were loaded with strictly vertical forces of 50 N on each cusp tip. Displacement and von Mises Maps were plotted for visualization of results.Results: The introduction of implant diminished the tendency of intrusion of the removable partial denture in all situations. The maximum stress was observed on implant in all situations. Approximating implant in direction of support teeth was benefit for stress distribution.Conclusion: Model D presented the lowest value for maximum tendency to displacement when compared with those found in the other models; model E demonstrated better relief with regard to demand from the abutment tooth; locating the implant near of the abutment tooth influenced positively the distribution of stresses on the analyzed structures.