Codomain Rigidity of the Dirichlet to Neumann Operator for the Riemannian Wave Equation

We study the Dirichlet to Neumann operator for the Riemannian wave equation on a compact Riemannian manifold. If the Riemannian manifold is modelled as an elastic medium, this operator represents the data available to an observer on the boundary of the manifold when the manifold is set into motion through boundary vibrations. We study the Dirichlet to Neumann operator when vibrations are imposed and data recorded on disjoint sets, a useful setting for applications. We prove that this operator determines the Dirichlet to Neumann operator where sources and observations are on the same set, provided a spectral condition on the Laplace-Beltrami operator for the manifold is satisfied. We prove this by providing an implementable procedure for determining a portion of the Riemannian manifold near the area where sources are applied. Drawing on established results, an immediate corollary is that a compact Riemannian manifold can be reconstructed from the Dirichlet to Neumann operator where sources and observations are on disjoint sets.

Codomain Rigidity of the Dirichlet to Neumann Operator for the Riemannian Wave Equation

We study the Dirichlet to Neumann operator for the Riemannian wave equation on a compact Riemannian manifold. If the Riemannian manifold is modelled as an elastic medium, this operator represents the data available to an observer on the boundary of the manifold when the manifold is set into motion through boundary vibrations. We study the Dirichlet to Neumann operator when vibrations are imposed and data recorded on disjoint sets, a useful setting for applications. We prove that this operator determines the Dirichlet to Neumann operator where sources and observations are on the same set, provided a spectral condition on the Laplace-Beltrami operator for the manifold is satisfied. We prove this by providing an implementable procedure for determining a portion of the Riemannian manifold near the area where sources are applied. Drawing on established results, an immediate corollary is that a compact Riemannian manifold can be reconstructed from the Dirichlet to Neumann operator where sources and observations are on disjoint sets.

An impulsive differential equation model for Marek's disease

Many dynamical processes are subject to abrupt changes in state. Often these perturbations can be periodic and of short duration relative to the evolving process. These types of phenomena are described well by what are referred to as impulsive differential equations, systems of differential equations coupled with discrete mappings in state space. In this thesis we employ impulsive differential equations to model disease transmission within an industrial livestock barn. In particular we focus on the poultry industry and a viral disease of poultry called Marek's disease. This system lends itself well to impulsive differential equations. Entire cohorts of poultry are introduced and removed from a barn concurrently. Additionally, Marek's disease is transmitted indirectly and the viral particles can survive outside the host for weeks. Therefore, depopulating, cleaning, and restocking of the barn are integral factors in modelling disease transmission and can be completely captured by the impulsive component of the model. Our model allows us to investigate how modern broiler farm practices can make disease elimination difficult or impossible to achieve. It also enables us to investigate factors that may contribute to virulence evolution. Our model suggests that by decrease the cohort duration or by decreasing the flock density, Marek's disease can be eliminated from a barn with no increase in cleaning effort. Unfortunately our model also suggests that these practices will lead to disease evolution towards greater virulence. Additionally, our model suggests that if intensive cleaning between cohorts does not rid the barn of disease, it may drive evolution and cause the disease to become more virulent.

Influencing effect of heat-treatment on radon emanation and exhalation characteristic of red mud

The reuse of industrial by-products is important for members of numerous industrial sectors. However, though the benefits of reuse are evident from an economical point of view, some compounds in these materials can have a negative effect on users' health.In this study, the radon emanation and exhalation features of red mud were surveyed using heat-treatment (100-1200 °C). As a result of the 1200°C-treated samples, massic radon exhalation capacity reduced from 75 ± 10 mBq kg^-1 h^-1 to 7 ± 4 mBq kg^-1 h^-1, approximately 10% of the initial exhalation rate.To find an explanation for internal structural changes, the porosity features of the heat-treated samples were also investigated. It was found that the cumulative pore volume reduced significantly in less than 100 nm, which can explain the reduced massic exhalation capacity in the high temperature treated range mentioned above.SEM snapshots were taken of the surfaces of the samples as visual evidence for superficial morphological changes. It was found that the surface of the high temperature treated samples had changed, proving the decrement of open pores on the surface.

Electron acoustic soliton energy of the Kadomtsev-Petviashvili equation in the Earth’s magnetotail region at critical ion density

The nonlinear properties of small amplitude electron-acoustic solitary waves (EAWs) in a homogeneous system of unmagnetized collisionless plasma consisted of a cold electron fluid and isothermal ions with two different temperatures obeying Boltzmann type distributions have been investigated. A reductive perturbation method was employed to obtain the Kadomstev-Petviashvili (KP) equation. At the critical ion density, the KP equation is not appropriate for describing the system. Hence, a new set of stretched coordinates
is considered to derive the modified KP equation. Moreover, the solitary solution, soliton energy and the associated electric field at the critical ion density were computed. The present investigation can be of relevance to the electrostatic solitary structures observed in various space plasma environments, such as Earth’s magnetotail region.

Symmetries and exact solutions of KP equation with an arbitrary nonlinear term

The generalized KP (GKP) equations with an arbitrary nonlinear term model and characterize many nonlinear physical phenomena. The symmetries of GKP equation with an arbitrary nonlinear term are obtained. The condition that must satisfy for existence the symmetries group of GKP is derived and also the obtained symmetries are classified according to different forms of the nonlinear term. The resulting similarity reductions are studied by performing the bifurcation and the phase portrait of GKP and also the corresponding solitary wave solutions of GKP
equation are constructed.

A Simultaneous Equation Approach to Estimating HIV Prevalence with Non-Ignorable Missing Responses

Estimates of HIV prevalence are important for policy in order to establish the health status of a country's population and to evaluate the effectiveness of population-based interventions and campaigns. However, participation rates in testing for surveillance conducted as part of household surveys, on which many of these estimates are based, can be low. HIV positive individuals may be less likely to participate because they fear disclosure, in which case estimates obtained using conventional approaches to deal with missing data, such as imputation-based methods, will be biased. We develop a Heckman-type simultaneous equation approach which accounts for non-ignorable selection, but unlike previous implementations, allows for spatial dependence and does not impose a homogeneous selection process on all respondents. In addition, our framework addresses the issue of separation, where for instance some factors are severely unbalanced and highly predictive of the response, which would ordinarily prevent model convergence. Estimation is carried out within a penalized likelihood framework where smoothing is achieved using a parametrization of the smoothing criterion which makes estimation more stable and efficient. We provide the software for straightforward implementation of the proposed approach, and apply our methodology to estimating national and sub-national HIV prevalence in Swaziland, Zimbabwe and Zambia.

Meningeal alveolar soft part sarcoma confirmed by characteristic ASPCR1-TFE3 fusion

Sarcoma metastatic to the brain is uncommon and rarely occurs as the initial manifestation of tumor. Alveolar soft part sarcoma (ASPS) is a rare but well-studied subtype of sarcoma. A 39-year-old man presented with seizures due to a left temporal meningeal-enhancing lesion with striking brain edema on MRI. The patient underwent neurosurgical resection for suspected meningioma. Histology showed large tumor cells clustering and forming small nests, in places with pseudoalveolar pattern. Diastase-resistant periodic acid-Schiff revealed very rare granular and rod-like cytoplasmic inclusions. Immunohistochemistry showed convincing positivity only with vimentin and smooth muscle actin. The histological features were strongly suggestive of ASPS. At the molecular level RT-PCR and sequencing analysis demonstrated ASPCR1-TFE3 fusion confirming the histological diagnosis of ASPS. There was no evidence of primary extracranial tumor by physical examination and on chest and abdominal CT scan 11 months after presentation. ASPS typically arise from the soft tissues of the extremities and develop multiple metastatic deposits usually with a long clinical course. This case may represent primary meningeal ASPS although metastatic deposit from an undiscovered primary site cannot be entirely excluded.

Automated Equation Formulation for Causal Loop Diagrams

The annotation of Business Dynamics models with parameters and equations, to simulate the system under study and further evaluate its simulation output, typically involves a lot of manual work. In this paper we present an approach for automated equation formulation of a given Causal Loop Diagram (CLD) and a set of associated time series with the help of neural network evolution (NEvo). NEvo enables the automated retrieval of surrogate equations for each quantity in the given CLD, hence it produces a fully annotated CLD that can be used for later simulations to predict future KPI development. In the end of the paper, we provide a detailed evaluation of NEvo on a business use-case to demonstrate its single step prediction capabilities.

Convection-diffusion in MHD: Comparison of primitive variable and vector potential solutions of the magnetic induction equation

Abstract not available

The Helmholtzian Equation: Stabilizing Mechanisms for Quadratic Nonlinear Fields

Hexagonal Resonant Triad patterns are shown to exist as stable solutions of a particular type of nonlinear field where no cubic field nonlinearity is present. The zero ‘dc’ Fourier mode is shown to stabilize these patterns produced by a pure quadratic field nonlinearity. Closed form solutions and stability results are obtained near the critical point, complimented by numerical studies far from the critical point. These results are obtained using a neural field based on the Helmholtzian operator. Constraints on structure and parameters for a general pure quadratic neural field which supports hexagonal patterns are obtained.

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

Thesis (Ph.D.)--University of Washington, 2016-08

Computer simulations of soft particles in flow

Understanding the dynamics of blood cells is a crucial element to discover biological mechanisms, to develop new efficient drugs, design sophisticated microfluidic devices, for diagnostics. In this work, we focus on the dynamics of red blood cells in microvascular flow. Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. We investigate the development of blood flow and its resistance starting from a dispersed configuration of red blood cells in simulations for different hematocrits, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow, which allows an estimation of the length of a vessel required for full flow development, $l_c \approx 25D$, with vessel diameter $D$. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length $l_c$. The presence of aggregation interactions between red blood cells lead in general to a reduction of blood flow resistance. The development of the cell-free layer thickness looks similar for both cases with and without aggregation interactions. Although, attractive interactions result in a larger cell-free layer plateau values. However, because the aggregation forces are short-ranged at high enough shear rates ($\bar{\dot{\gamma}} \gtrsim 50~\text{s}^{-1}$) aggregation of red blood cells does not bring a significant change to the blood flow properties. Also, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness with respect to flow rate assuming steady-state flow. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells. Finally, we investigate segregation phenomena in blood as a two-component suspension under Poiseuille flow, consisting of red blood cells and target cells. The spatial distribution of particles in blood flow is very important. For example, in case of nanoparticle drug delivery, the particles need to come closer to microvessel walls, in order to adhere and bring the drug to a target position within the microvasculature. Here we consider that segregation can be described as a competition between shear-induced diffusion and the lift force that pushes every soft particle in a flow away from the wall. In order to investigate the segregation, on one hand, we have 2D DPD simulations of red blood cells and target cell of different sizes, on the other hand the Fokker-Planck equation for steady state. For the equation we measure force profile, particle distribution and diffusion constant across the channel. We compare simulation results with those from the Fokker-Planck equation and find a very good correspondence between the two approaches. Moreover, we investigate the diffusion behavior of target particles for different hematocrit values and shear rates. Our simulation results indicate that diffusion constant increases with increasing hematocrit and depends linearly on shear rate. The third part of the study describes development of a simulation model of complex vascular geometries. The development of the model is important to reproduce vascular systems of small pieces of tissues which might be gotten from MRI or microscope images. The simulation model of the complex vascular systems might be divided into three parts: modeling the geometry, developing in- and outflow boundary conditions, and simulation domain decomposition for an efficient computation. We have found that for the in- and outflow boundary conditions it is better to use the SDPD fluid than DPD one because of the density fluctuations along the channel of the latter. During the flow in a straight channel, it is difficult to control the density of the DPD fluid. However, the SDPD fluid has not that shortcoming even in more complex channels with many branches and in- and outflows because the force acting on particles is calculated also depending on the local density of the fluid.