A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

A finite element-spherical harmonics model for radiative transfer in inhomogeneous clouds. Part II. some applications

EVENT has been used to examine the effects of 3D cloud structure, distribution, and inhomogeneity on the scattering of visible solar radiation and the resulting 3D radiation field. Large eddy simulation and aircraft measurements are used to create realistic cloud fields which are continuous or broken with smooth or uneven tops. The values, patterns and variance in the resulting downwelling and upwelling radiation from incident visible solar radiation at different angles are then examined and compared to measurements. The results from EVENT confirm that 3D cloud structure is important in determining the visible radiation field, and that these results are strongly influenced by the solar zenith angle. The results match those from other models using visible solar radiation, and are supported by aircraft measurements of visible radiation, providing confidence in the new model.

Repressor element 1 silencing transcription factor couples loss of pluripotency with neural induction and neural differentiation

Neural differentiation of embryonic stem cells (ESCs) requires coordinated repression of the pluripotency regulatory program and reciprocal activation of the neurogenic regulatory program. Upon neural induction, ESCs rapidly repress expression of pluripotency genes followed by staged activation of neural progenitor and differentiated neuronal and glial genes. The transcriptional factors that underlie maintenance of pluripotency are partially characterized whereas those underlying neural induction are much less explored, and the factors that coordinate these two developmental programs are completely unknown. One transcription factor, REST (repressor element 1 silencing transcription factor), has been linked with terminal differentiation of neural progenitors and more recently, and controversially, with control of pluripotency. Here, we show that in the absence of REST, coordination of pluripotency and neural induction is lost and there is a resultant delay in repression of pluripotency genes and a precocious activation of both neural progenitor and differentiated neuronal and glial genes. Furthermore, we show that REST is not required for production of radial glia-like progenitors but is required for their subsequent maintenance and differentiation into neurons, oligodendrocytes, and astrocytes. We propose that REST acts as a regulatory hub that coordinates timely repression of pluripotency with neural induction and neural differentiation.

A mathematical model of the sterol regulatory element binding protein 2 cholesterol biosynthesis pathway

Cholesterol is one of the key constituents for maintaining the cellular membrane and thus the integrity of the cell itself. In contrast high levels of cholesterol in the blood are known to be a major risk factor in the development of cardiovascular disease. We formulate a deterministic nonlinear ordinary differential equation model of the sterol regulatory element binding protein 2 (SREBP-2) cholesterol genetic regulatory pathway in an hepatocyte. The mathematical model includes a description of genetic transcription by SREBP-2 which is subsequently translated to mRNA leading to the formation of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR), a main precursor of cholesterol synthesis. Cholesterol synthesis subsequently leads to the regulation of SREBP-2 via a negative feedback formulation. Parameterised with data from the literature, the model is used to understand how SREBP-2 transcription and regulation affects cellular cholesterol concentration. Model stability analysis shows that the only positive steady-state of the system exhibits purely oscillatory, damped oscillatory or monotic behaviour under certain parameter conditions. In light of our findings we postulate how cholesterol homestasis is maintained within the cell and the advantages of our model formulation are discussed with respect to other models of genetic regulation within the literature.

A high frequency boundary element method for scattering by a class of nonconvex obstacles

In this paper we propose and analyse a hybrid numerical-asymptotic boundary element method for the solution of problems of high frequency acoustic scattering by a class of sound-soft nonconvex polygons. The approximation space is enriched with carefully chosen oscillatory basis functions; these are selected via a study of the high frequency asymptotic behaviour of the solution. We demonstrate via a rigorous error analysis, supported by numerical examples, that to achieve any desired accuracy it is sufficient for the number of degrees of freedom to grow only in proportion to the logarithm of the frequency as the frequency increases, in contrast to the at least linear growth required by conventional methods. This appears to be the first such numerical analysis result for any problem of scattering by a nonconvex obstacle. Our analysis is based on new frequency-explicit bounds on the normal derivative of the solution on the boundary and on its analytic continuation into the complex plane.

A frequency-independent boundary element method for scattering by two-dimensional screens and apertures

We propose and analyse a hybrid numerical–asymptotic hp boundary element method (BEM) for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation space enriched with oscillatory basis functions, chosen to capture the high-frequency asymptotics of the solution. We provide a rigorous frequency-explicit error analysis which proves that the method converges exponentially as the number of degrees of freedom N increases, and that to achieve any desired accuracy it is sufficient to increase N in proportion to the square of the logarithm of the frequency as the frequency increases (standard BEMs require N to increase at least linearly with frequency to retain accuracy). Our numerical results suggest that fixed accuracy can in fact be achieved at arbitrarily high frequencies with a frequency-independent computational cost, when the oscillatory integrals required for implementation are computed using Filon quadrature. We also show how our method can be applied to the complementary ‘breakwater’ problem of propagation through an aperture in an infinite sound-hard screen.

Acoustic scattering : high frequency boundary element methods and unified transform methods

We describe some recent advances in the numerical solution of acoustic scattering problems. A major focus of the paper is the efficient solution of high frequency scattering problems via hybrid numerical-asymptotic boundary element methods. We also make connections to the unified transform method due to A. S. Fokas and co-authors, analysing particular instances of this method, proposed by J. A. De-Santo and co-authors, for problems of acoustic scattering by diffraction gratings.

Genome-wide mapping of 5-hydroxymethylcytosine in three rice cultivars reveals its preferential localization in transcriptionally silent transposable element genes

5-Hydroxymethylcytosine (5hmC), a modified form of cytosine that is considered the sixth nucleobase in DNA, has been detected in mammals and is believed to play an important role in gene regulation. In this study, 5hmC modification was detected in rice by employing a dot-blot assay, and its levels was further quantified in DNA from different rice tissues using liquid chromatography-multistage mass spectrometry (LC-MS/MS/MS). The results showed large intertissue variation in 5hmC levels. The genome-wide profiles of 5hmC modification in three different rice cultivars were also obtained using a sensitive chemical labelling followed by a next-generation sequencing method. Thousands of 5hmC peaks were identified, and a comparison of the distributions of 5hmC among different rice cultivars revealed the specificity and conservation of 5hmC modification. The identified 5hmC peaks were significantly enriched in heterochromatin regions,and mainly located in transposable element (TE) genes, especially around retrotransposons. The correlation analysis of 5hmC and gene expression data revealed a close association between 5hmC and silent TEs. These findings provide a resource for plant DNA 5hmC epigenetic studies and expand our knowledge of 5hmC modification.

Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

Noether-type discrete conserved quantities arising from a finite element approximation of a variational problem

In this work, we prove a weak Noether-type Theorem for a class of variational problems that admit broken extremals. We use this result to prove discrete Noether-type conservation laws for a conforming finite element discretisation of a model elliptic problem. In addition, we study how well the finite element scheme satisfies the continuous conservation laws arising from the application of Noether’s first theorem (1918). We summarise extensive numerical tests, illustrating the conservation of the discrete Noether law using the p-Laplacian as an example and derive a geometric-based adaptive algorithm where an appropriate Noether quantity is the goal functional.

Thyroid hormone can increase estrogen-mediated transcription from a consensus estrogen response element in neuroblastoma cells

Thyroid hormones (T) and estrogens (E) are nuclear receptor ligands with at least two molecular mechanisms of action: (i) relatively slow genomic effects, such as the regulation of transcription by cognate T receptors (TR) and E receptors (ER); and (ii) relatively rapid nongenomic effects, such as kinase activation and calcium release initiated at the membrane by putative membrane receptors. Genomic and nongenomic effects were thought to be disparate and independent. However, in a previous study using a two-pulse paradigm in neuroblastoma cells, we showed that E acting at the membrane could potentiate transcription from an E-driven reporter gene in the nucleus. Because both T and E can have important effects on mood and cognition, it is possible that the two hormones can act synergistically. In this study, we demonstrate that early actions of T via TRalpha1 and TRbeta1 can potentiate E-mediated transcription (genomic effects) from a consensus E response element (ERE)-driven reporter gene in transiently transfected neuroblastoma cells. Such potentiation was reduced by inhibition of mitogen-activated protein kinase. Using phosphomutants of ERalpha, we also show that probable mitogen-activated protein kinase phosphorylation sites on the ERalpha, the serines at position 167 and 118, are important in TRbeta1-mediated potentiation of ERalpha-induced transactivation. We suggest that crosstalk between T and E includes potential interactions through both nuclear and membrane-initiated molecular mechanisms of hormone signaling.

NEONATAL HANDLING AND THE MATERNAL ODOR PREFERENCE IN RAT PUPS: INVOLVEMENT OF MONOAMINES AND CYCLIC AMP RESPONSE ELEMENT-BINDING PROTEIN PATHWAY IN THE OLFACTORY BULB

Early-life environmental events, such as the handling procedure, can induce long-lasting alterations upon several behavioral and neuroendocrine systems. However, the changes within the pups that could be causally related to the effects in adulthood are still poorly understood. In the present study, we analyzed the effects of neonatal handling on behavioral (maternal odor preference) and biochemical (cyclic AMP response element-binding protein (CREB) phosphorylation, noradrenaline (NA), and serotonin (5-HT) levels in the olfactory bulb (OB)) parameters in 7-day-old male and female rat pups. Repeated handling (RH) abolished preference for the maternal odor in female pups compared with nonhandled (NH) and the single-handled (SH) ones, while in RH males the preference was not different than NH and SH groups. In both male and female pups, RH decreased NA activity in the OB, but 5-HT activity increased only in males. Since preference for the maternal odor involves the synergic action of NA and 5-HT in the OB, the maintenance of the behavior in RH males could be related to the increased 5-HT activity, in spite of reduction in the NA activity in the OB. RH did not alter CREB phosphorylation in the OB of both male and females compared with NH pups. The repeated handling procedure can affect the behavior of rat pups in response to the maternal odor and biochemical parameters related to the olfactory learning mechanism. Sex differences were already detected in 7-day-old pups. Although the responsiveness of the hypothalamic-pituitary-adrenal axis to stressors is reduced in the neonatal period, environmental interventions may impact behavioral and biochemical mechanisms relevant to the animal at that early age. (C) 2009 IBRO. Published by Elsevier Ltd. All rights reserved.

Characterization of Metal and Trace Element Contents of Particulate Matter (PM10) Emitted by Vehicles Running on Brazilian FuelsHydrated Ethanol and Gasoline with 22% of Anhydrous Ethanol

Emission of fine particles by mobile sources has been a matter of great concern due to its potential risk both to human health and the environment. Although there is no evidence that one sole component may be responsible for the adverse health outcomes, it is postulated that the metal particle content is one of the most important factors, mainly in relation to oxidative stress. Data concerning the amount and type of metal particles emitted by automotive vehicles using Brazilian fuels are limited. The aim of this study was to identify inhalable particles (PM10) and their trace metal content in two light-duty vehicles where one was fueled with ethanol while the other was fueled with gasoline mixed with 22% of anhydrous ethanol (gasohol); these engines were tested on a chassis dynamometer. The elementary composition of the samples was evaluated by the particle-induced x-ray emission technique. The experiment showed that total emission factors ranged from 2.5 to 11.8 mg/km in the gasohol vehicle, and from 1.2 to 3 mg/km in the ethanol vehicle. The majority of particles emitted were in the fine fraction (PM2.5), in which Al, Si, Ca, and Fe corresponded to 80% of the total weight. PM10 emissions from the ethanol vehicle were about threefold lower than those of gasohol. The elevated amount of fine particulate matter is an aggravating factor, considering that these particles, and consequently associated metals, readily penetrate deeply into the respiratory tract, producing damage to lungs and other tissues.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.