Abstract: We propose a new approach to the numerical solution of radiative transfer equations with certified a posteriori error bounds. A key role is played by stable Petrov--Galerkin type variational formulations of parametric transport equations and corresponding radiative transfer equations. This allows us to formulate an iteration in a suitable, infinite dimensional function space that is guaranteed to converge with a fixed error reduction per step. The numerical scheme is then based on approximately realizing this iteration within dynamically updated accuracy tolerances that still ensure convergence to the exact solution. To advance this iteration two operations need to be performed within suitably tightened accuracy tolerances. First, the global scattering operator needs to be approximately applied to the current iterate within a tolerance comparable to the current accuracy level. Second, parameter dependent linear transport equations need to be solved, again at the required accuracy of the iteration. To ensure that the stage dependent error tolerances are met, one has to employ rigorous a posteriori error bounds which, in our case, rest on a Discontinuous Petrov--Galerkin (DPG) scheme. These a posteriori bounds are not only crucial for guaranteeing the convergence of the perturbed iteration but are also used to generate adapted parameter dependent spatial meshes. This turns out to significantly reduce overall computational complexity. Since the global operator is only applied, we avoid the need to solve linear systems with densely populated matrices. Moreover, the approximate application of the global scatterer accelerated through low-rank approximation and matrix compression techniques. The theoretical findings are illustrated and complemented by numerical experiments with non-trivial scattering kernels.
Abstract: We prove the large-time asymptotic orbital stability of strictly entropic Riemann shock solutions of first order scalar hyperbolic balance laws, under piecewise regular perturbations provided that the source term is dissipative about endstates of the shock. Moreover the convergence towards a shifted reference state is exponential with a rate predicted by the linearized equations about constant endstates.
Abstract: A fluid model is developed for multicomponent two-temperature magnetized plasmas in chemical non-equilibrium from the partially- to fully-ionized collisional regimes. We focus on transport phenomena aiming at representing the chromosphere of the Sun. Graille et al. [M3AS 19(04):527-599, 2009] have derived an asymptotic fluid model for multicomponent plamas from kinetic theory, yielding a rigorous description of the dissipative effects. The governing equations and consistent transport properties are obtained using a multiscale Chapman-Enskog perturbative solution to the Boltzmann equation based on a non-dimensional analysis. The mass disparity between the electrons and heavy particles is accounted for, as well as the influence of the electromagnetic field. We couple this model to the Maxwell equations for the electromagnetic field and derive the generalized Ohm's law for multicomponent plasmas. The model inherits a well-identified mathematical structure leading to an extended range of validity for the Sun chromosphere conditions. We compute consistent transport properties by means of a spectral Galerkin method using the Laguerre-Sonine polynomial approximation. Two non-vanishing polynomial terms are used when deriving the transport systems for electrons, whereas only one term is retained for heavy particles. In a simplified framework where the plasma is fully ionized, we compare the transport properties for the Sun chromosphere to conventional expressions for magnetized plasmas due to Braginskii, showing a good agreement between both results. For more general partially ionized conditions, representative of the Sun chromosphere, we compute the muticomponent transport properties corresponding to the species diffusion velocities, heavy-particle and electron heat fluxes, and viscous stress tensor of the model, for a Helium-Hydrogen mixture in local thermodynamic equilibrium. The model is assessed for the 3D radiative magnetohydrodynamic simulation of a pore, in the highly turbulent upper layer of the solar convective zone. The resistive term is found to dominate mainly the dynamics of the electric field at the pore location. The battery term for heavy particles appears to be higher at the pore location and at some intergranulation boundaries.
Abstract: Motivation: Genome-Wide Association Studies (GWAS) seek to identify causal genomic variants associated with rare human diseases. The classical statistical approach for detecting these variants is based on univariate hypothesis testing, with healthy individuals being tested against affected individuals at each locus. Given that an individual's genotype is characterized by up to one million SNPs, this approach lacks precision, since it may yield a large number of false positives that can lead to erroneous conclusions about genetic associations with the disease. One way to improve the detection of true genetic associations is to reduce the number of hypotheses to be tested by grouping SNPs. Results: We propose a dimension-reduction approach which can be applied in the context of GWAS by making use of the haplotype structure of the human genome. We compare our method with standard univariate and multivariate approaches on both synthetic and real GWAS data, and we show that reducing the dimension of the predictor matrix by aggregating SNPs gives a greater precision in the detection of associations between the phenotype and genomic regions.
Ege Tıp Dergisi, Volume 56, pp 145-147; doi:10.19161/etd.392821
Published: 18 October 2018
Hacettepe Journal of Mathematics and Statistics, Volume 47; doi:10.15672/hjms.2017.518
World Journal of Social Science Research, Volume 4; doi:10.22158/wjssr.v4n4p287
Abstract: ADPN (Asynchronous Dynamic Pushdown Networks) are an abstract model for concurrent programs with recursive procedures and dynamic thread creation. Usually, asynchronous dynamic pushdown networks are described with interleaving semantics, in which the backward analysis is not effective. In order to improve interleaving semantics, tree semantics approach was introduced. This paper extends the tree semantics to ADPN. Because the reachability problem of ADPN is also undecidable, we address the context-bounded reachability problem and provide an algorithm for backward reachability analysis with tree-based semantics Approach.
Published: 18 October 2018
The publisher has not yet granted permission to display this abstract.
Published: 18 October 2018
Atatürk Üniversitesi Diş Hekimliği Fakültesi Dergisi, Volume 26; doi:10.17567/ataunidfd.285136