Developing discontinuous Galerkin methods for the resolution of the incompressible Navier-Stokes equations

Informe de investigación elaborado a partir de una estancia en el Laboratorio de Diseño Computacional en Aeroespacial en el Massachusetts Institute of Technology (MIT), Estados Unidos, entre noviembre de 2006 y agosto de 2007. La aerodinámica es una rama de la dinámica de fluidos referida al estudio de los movimientos de los líquidos o gases, cuya meta principal es predecir las fuerzas aerodinámicas en un avión o cualquier tipo de vehículo, incluyendo los automóviles. Las ecuaciones de Navier-Stokes representan un estado dinámico del equilibrio de las fuerzas que actúan en cualquier región dada del fluido. Son uno de los sistemas de ecuaciones más útiles porque describen la física de una gran cantidad de fenómenos como corrientes del océano, flujos alrededor de una superficie de sustentación, etc. En el contexto de una tesis doctoral, se está estudiando un flujo viscoso e incompresible, solucionando las ecuaciones de Navier- Stokes incompresibles de una manera eficiente. Durante la estancia en el MIT, se ha utilizado un método de Galerkin discontinuo para solucionar las ecuaciones de Navier-Stokes incompresibles usando, o bien un parámetro de penalti para asegurar la continuidad de los flujos entre elementos, o bien un método de Galerkin discontinuo compacto. Ambos métodos han dado buenos resultados y varios ejemplos numéricos se han simulado para validar el buen comportamiento de los métodos desarrollados. También se han estudiado elementos particulares, los elementos de Raviart y Thomas, que se podrían utilizar en una formulación mixta para obtener un algoritmo eficiente para solucionar problemas numéricos complejos.

On separation of a degenerate differential operator in Hilbert space

A coercive estimate for a solution of a degenerate second order di fferential equation is installed, and its applications to spectral problems for the corresponding dif ferential operator is demonstrated. The suffi cient conditions for existence of the solutions of one class of the nonlinear second order diff erential equations on the real axis are obtained.

Uniform estimates in the Poincare-Aronszajn theorem on the separation of singularities of analytic functions

We study the possibility of splitting any bounded analytic function $f$ with singularities in a closed set $E\cup F$ as a sum of two bounded analytic functions with singularities in $E$ and $F$ respectively. We obtain some results under geometric restrictions on the sets $E$ and $F$ and we provide some examples showing the sharpness of the positive results.

A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [13, 12], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [13] to further simplify the kinetic theory.

Temperamental Dimension and Anxiety Problems in a Clinical Sample of Three- to Six-year old Children: A Study of Variables

In the last few years, many researchers have studied the presence of common dimensions of temperament in subjects with symptoms of anxiety. The aim of this study is to examine the association between temperamental dimensions (high negative affect and activity level) and anxiety problems in clinicalpreschool children. A total of 38 children, ages 3 to 6 years, from the Infant and Adolescent Mental Health Center of Girona and the Center of Diagnosis and Early Attention of Sabadell and Olot were evaluated by parents and psychologists. Their parents completed several screening scales and, subsequently, clinical child psychopathology professionals carried out diagnostic interviews with children from the sample who presented signs of anxiety. Findings showed that children with high levels of negative affect and low activity level have pronounced symptoms of anxiety. However, children with anxiety disorders do not present different temperament styles from their peers without these pathologies

General clas of inhomogeneous perfect-fluid solutions

We present a general class of solutions to Einstein's field equations with two spacelike commuting Killing vectors by assuming the separation of variables of the metric components. The solutions can be interpreted as inhomogeneous cosmological models. We show that the singularity structure of the solutions varies depending on the different particular choices of the parameters and metric functions. There exist solutions with a universal big-bang singularity, solutions with timelike singularities in the Weyl tensor only, solutions with singularities in both the Ricci and the Weyl tensors, and also singularity-free solutions. We prove that the singularity-free solutions have a well-defined cylindrical symmetry and that they are generalizations of other singularity-free solutions obtained recently.

Simulated Annealing, High-Order Statistics and Mutual Information for Separation of Sources

In this article, the fusion of a stochastic metaheuristic as Simulated Annealing (SA) with classical criteria for convergence of Blind Separation of Sources (BSS), is shown. Although the topic of BSS, by means of various techniques, including ICA, PCA, and neural networks, has been amply discussed in the literature, to date the possibility of using simulated annealing algorithms has not been seriously explored. From experimental results, this paper demonstrates the possible benefits offered by SA in combination with high order statistical and mutual information criteria for BSS, such as robustness against local minima and a high degree of flexibility in the energy function.

Localized prostate cancer in Norway, the United States, and Spain: between-country differences of variables before treatment among patients eligible for curative treatment

Between-country differences in medical and sociodemographic variables, and patient-related outcomes (PROs) before treatment might explain published variations of side effects after radical prostatecomy (RP) or radiotherapy (RAD) for prostate cancer (PCa). This hypothesis was tested among 1908 patients from the United States, Spain, and Norway. Significant between-country differences were observed for most factors investigated before treatment. The observations should be considered in comparison of the frequency and severity of internationally published studies. Background: In men with PCa, large variations of PROs after RP or high-dose RAD might be related to betweencountry differences of medical and sociodemographic variables, and differences in PROs before treatment in the sexual and urinary domains. Patients and Methods: In 1908 patients with localized PCa from Norway, the United States, or Spain, the relation between medical (prostate-specific antigen, Gleason score, cT-category) and sociodemographic variables (age, education, marital status) before treatment was investigated. Using the Expanded Prostate Cancer Index Composite questionnaire, PROs before treatment within the sexual and urinary domains were also considered. Results: Compared with the European patients, American patients were younger, fewer had comorbid conditions, and more had a high education level. Fifty-three percent of the US men eligible for RP had low-risk tumors compared with 42% and 31% among the Norwegian and the Spanish patients, respectively. Among the Spanish RAD patients, 54% had had low-risk tumors compared with 34% of the American and 21% of the Norwegian men planned for RAD, respectively. Compared with the European patients, significantly fewer US patients reported moderate or severe sexual dysfunction and related problems. In most subgroups, the number of patients with sexual or urinary dysfunction exceeded that of patients with bother related to the reported dysfunction. Conclusion: Statistically significant between-country differences were observed in medical and sociodemographic variables, and in PROs before treatment within the sexual and urinary domains. Large differences between reported dysfunction and related problems within the sexual and urinary domains indicate that dysfunction and bother should be reported separately in addition to calculation of summary scores. The documented differences, not at least regarding PROs, might in part explain the large variation of side effects after treatment evident in the medical literature

Transient shear banding in time-dependent fluids

We study the dynamics of shear-band formation and evolution using a simple rheological model. The description couples the local structure and viscosity to the applied shear stress. We consider in detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation to obtain the time evolution of the local velocity and viscosity fields. It is found that the underlying reason for dynamic effects is the nonhomogeneous shear distribution, which is amplified due to a positive feedback between the flow field and the viscosity response of the shear thinning fluid. This offers a simple explanation for the recent observations of transient shear banding in time-dependent fluids. Extensions to more complicated rheological systems are considered.

From lattice-gas models to nonlinear diffusion models: A derivation of macroscopic equations and fluctuations

We derive nonlinear diffusion equations and equations containing corrections due to fluctuations for a coarse-grained concentration field. To deal with diffusion coefficients with an explicit dependence on the concentration values, we generalize the Van Kampen method of expansion of the master equation to field variables. We apply these results to the derivation of equations of phase-separation dynamics and interfacial growth instabilities.

Interaction of Spiral Waves in the Complex Ginzburg-Landau Equation

Solutions of the general cubic complex Ginzburg-Landau equation comprising multiple spiral waves are considered, and laws of motion for the centers are derived. The direction of the motion changes from along the line of centers to perpendicular to the line of centers as the separation increases, with the strength of the interaction algebraic at small separations and exponentially small at large separations. The corresponding asymptotic wave number and frequency are also determined, which evolve slowly as the spirals move

On moduli of smoothness of fractional order

In this paper we consider the properties of moduli of smoothness of fractional order. The main result of the paper describes the equivalence of the modulus of smoothness and a function from some class.

Contact melting of a three-dimensional phase change material on a flat substrate

In this paper a model is developed to describe the three dimensional contact melting process of a cuboid on a heated surface. The mathematical description involves two heat equations (one in the solid and one in the melt), the Navier-Stokes equations for the flow in the melt, a Stefan condition at the phase change interface and a force balance between the weight of the solid and the countering pressure in the melt. In the solid an optimised heat balance integral method is used to approximate the temperature. In the liquid the small aspect ratio allows the Navier-Stokes and heat equations to be simplified considerably so that the liquid pressure may be determined using an igenfunction expansion and finally the problem is reduced to solving three first order ordinary differential equations. Results are presented showing the evolution of the melting process. Further reductions to the system are made to provide simple guidelines concerning the process. Comparison of the solutions with experimental data on the melting of n-octadecane shows excellent agreement.

Country effects in ISSP-1993 environmental data: Comparison of SEM approaches

Structural equation models (SEM) are commonly used to analyze the relationship between variables some of which may be latent, such as individual ``attitude'' to and ``behavior'' concerning specific issues. A number of difficulties arise when we want to compare a large number of groups, each with large sample size, and the manifest variables are distinctly non-normally distributed. Using an specific data set, we evaluate the appropriateness of the following alternative SEM approaches: multiple group versus MIMIC models, continuous versus ordinal variables estimation methods, and normal theory versus non-normal estimation methods. The approaches are applied to the ISSP-1993 Environmental data set, with the purpose of exploring variation in the mean level of variables of ``attitude'' to and ``behavior''concerning environmental issues and their mutual relationship across countries. Issues of both theoretical and practical relevance arise in the course of this application.

Why do they separate it or not? Attitudes and behaviors towards organic waste separation

This article shows the results of an exploratory study related to the separation of organic waste in order to offer suggestions for the improvement of waste disposal communication campaigns. The overall objective is to analyze attitude and behavior of those who do and those who do not separate organic waste, related to a specific promotional campaign carried out in two neighborhoods, in the municipality of Badalona (Spain), within the framework of the study of proenvironmental attitudes and behaviors and based on the Psychosocial Four Spheres Model. 1,010 interviews were conducted and data was analyzed using Chi-Squared Automatic Interaction Detector (CHAID). Waste separation behavior was used as a dependent variable. The reasons given to explain why people do or do not separate organic waste and sociodemographic variables, have been introduced as independent variables. In accordance with the Four Spheres Model, results show significant differences in waste separation. Based on the profiles obtained, we find some predictive variables that facilitate the development of communication campaigns according to the requirements of each community.