Quantitative unique continuation for the heat equations with inverse square potential.
Ontology highlight
ABSTRACT:
SUBMITTER: Zheng G
PROVIDER: S-EPMC6244741 | biostudies-other | 2018
REPOSITORIES: biostudies-other
Similar Datasets
Project description:A classical scenario for tipping is that a dynamical system experiences a slow parameter drift across a fold tipping point, caused by a run-away positive feedback loop. We study what happens if one turns around after one has crossed the threshold. We derive a simple criterion that relates how far the parameter exceeds the tipping threshold maximally and how long the parameter stays above the threshold to avoid tipping in an inverse-square law to observable properties of the dynamical system near the fold. For the case when the dynamical system is subject to stochastic forcing we give an approximation to the probability of tipping if a parameter changing in time reverses near the tipping point. The derived approximations are valid if the parameter change in time is sufficiently slow. We demonstrate for a higher-dimensional system, a model for the Indian summer monsoon, how numerically observed escape from the equilibrium converge to our asymptotic expressions. The inverse-square law between peak of the parameter forcing and the time the parameter spends above a given threshold is also visible in the level curves of equal probability when the system is subject to random disturbances.
| S-EPMC6405451 | biostudies-literature
Project description:Experimental study on free convection heat transfer was carried out inside shallow square enclosures filled with water. Two enclosures were used with size L × L × H (m3), where L and H are the side length and the inside thickness of the enclosure, respectively. Two different aspect ratios, ? = L/H, 7.143 and 12.0 were used. Constant heat flux boundary condition was applied at the bottom surface and stream of ambient air was applied at the top surface. Different values of constant heat fluxes were used as boundary conditions. For each aspect ratio of the enclosure, average Nusselt numbers were developed and correlated with the modified Rayleigh number. The results show that the heat transfer coefficient increases with the modified Rayleigh number with observed overlapping region between the two aspect ratios. A general overall correlation was developed using the aspect ratio as a parameter. The results also show that at the overlapping zone, the Nusselt number decreases as the aspect ratio increases for the fixed modified Rayleigh number.
| S-EPMC6209140 | biostudies-literature
Project description:We investigate the inverse problem of identifying a conditional probability measure in measure-dependent evolution equations arising in size-structured population modeling. We formulate the inverse problem as a least squares problem for the probability measure estimation. Using the Prohorov metric framework, we prove existence and consistency of the least squares estimates and outline a discretization scheme for approximating a conditional probability measure. For this scheme, we prove general method stability. The work is motivated by Partial Differential Equation (PDE) models of flocculation for which the shape of the post-fragmentation conditional probability measure greatly impacts the solution dynamics. To illustrate our methodology, we apply the theory to a particular PDE model that arises in the study of population dynamics for flocculating bacterial aggregates in suspension, and provide numerical evidence for the utility of the approach.
| S-EPMC5352987 | biostudies-literature
Project description:In this report, we quantify the divergence from the inverse square law (ISL) of the beam output as a function of distance (standoff) from closed-ended applicators for a modern clinical orthovoltage unit. The divergence is clinically significant exceeding 3% at a 1.2 cm distance for 4 × 4 and 10 × 10 cm2 closed-ended applicators. For all investigated cases, the measured dose falloff is more rapid than that predicted by the ISL and, therefore, causes a systematic underdose when using the ISL for dose calculations at extended SSD. The observed divergence from the ISL in closed-ended applicators can be explained by the end-plate scattering contribution not accounted for in the ISL calculation. The standoff measurements were also compared to the predictions from a home-built kV dose computation algorithm, kVDoseCalc. The kVDoseCalc computation predicted a more rapid falloff with distance than observed experimentally. The computation and measurements agree to within 1.1% for standoff distances of 3 cm or less for 4 × 4 cm2 and 10 × 10 cm2 field sizes. The overall agreement is within 2.3% for all field sizes and standoff distances measured. No significant deviation from the ISL was observed for open-ended applicators for standoff distances up to 10 cm.
| S-EPMC5875524 | biostudies-literature
Project description:Many models in Systems Biology are described as a system of Ordinary Differential Equations, which allows for transient, steady-state or bifurcation analysis when kinetic information is available. Complementary structure-related qualitative analysis techniques have become increasingly popular in recent years, like qualitative model checking or pathway analysis (elementary modes, invariants, flux balance analysis, graph-based analyses, chemical organization theory, etc.). They do not rely on kinetic information but require a well-defined structure as stochastic analysis techniques equally do. In this article, we look into the structure inference problem for a model described by a system of Ordinary Differential Equations and provide conditions for the uniqueness of its solution. We describe a method to extract a structured reaction network model, represented as a bipartite multigraph, for example, a continuous Petri net (CPN), from a system of Ordinary Differential Equations (ODEs). A CPN uniquely defines an ODE, and each ODE can be transformed into a CPN. However, it is not obvious under which conditions the transformation of an ODE into a CPN is unique, that is, when a given ODE defines exactly one CPN. We provide biochemically relevant sufficient conditions under which the derived structure is unique and counterexamples showing the necessity of each condition. Our method is implemented and available; we illustrate it on some signal transduction models from the BioModels database. A prototype implementation of the method is made available to modellers at http://contraintes.inria.fr/~soliman/ode2pn.html, and the data mentioned in the "Results" section at http://contraintes.inria.fr/~soliman/ode2pn_data/. Our results yield a new recommendation for the import/export feature of tools supporting the SBML exchange format.
| S-EPMC3008708 | biostudies-literature
Project description:Transient mass-transfer phenomena occurring in natural and engineered systems consist of convection, diffusion, and reaction processes. The coupled phenomena can be described by using the unsteady convection-diffusion-reaction (CDR) equation, which is classified in mathematics as a linear, parabolic partial-differential equation. The availability of analytic solutions is limited to simple cases, e.g., unsteady diffusion and steady convective diffusion. The CDR equation has been considered analytically intractable, depending on the initial and boundary conditions. If spatial adsorption and desorption of matter are super-positioned in the CDR equation as sink and source functions, respectively, then the governing equation becomes an unsteady convection-diffusion-reaction-source (CDRS) equation, of which general solutions are unknown. In this study, a general 1D analytic solution of the CDRS equation is obtained by using a one-sided Laplace transform, by assuming constant diffusivity, velocity, and reactivity. This paper also provides a general formalism to derive 1D analytic solutions for Dirichlet/Dirichlet and Dirichlet/Neumann boundary conditions. Derivations of the analytic solutions are found to be straightforward if a combination of the source function and the initial concentration provide a non-zero singularity pole of inverse Laplace transform.
| S-EPMC7228951 | biostudies-literature
Project description:Conjugate natural convection-conduction heat transfer in a square enclosure with a finite wall thickness is studied numerically in the present paper. The governing parameters considered are the Rayleigh number (5 × 10(3) ≤ Ra ≤ 10(6)), the wall-to-fluid thermal conductivity ratio (0.5 ≤ Kr ≤ 10), and the ratio of wall thickness to its height (0.2 ≤ D ≤ 0.4). The staggered grid arrangement together with MAC method was employed to solve the governing equations. It is found that the fluid flow and the heat transfer can be controlled by the thickness of the bottom wall, the thermal conductivity ratio, and the Rayleigh number.
| S-EPMC4058157 | biostudies-other
Project description:We describe a novel 3-D ultrasound technology, the quantitative transmission ultrasound system and algorithm to image a pendent breast in a water bath. Quantitative accuracy is verified using phantoms. Morphological accuracy is verified using cadaveric breast and in vivo images, and spatial resolution is estimated. This paper generalizes an earlier 2-D algorithm to a full 3-D inversion algorithm and shows the importance of such a 3-D algorithm for artifact suppression as compared with the 2-D algorithm. The resultant high-resolution ultrasound images, along with quantitative information regarding tissue speed of sound/stiffness, provide a more accurate depiction of the breast anatomy and lesions, contributing to improved breast care.
| S-EPMC6214813 | biostudies-literature
Project description:The improved capacity to acquire quantitative data in a clinical setting has generally failed to improve outcomes in acutely ill patients, suggesting a need for advances in computer-supported data interpretation and decision making. In particular, the application of mathematical models of experimentally elucidated physiological mechanisms could augment the interpretation of quantitative, patient-specific information and help to better target therapy. Yet, such models are typically complex and nonlinear, a reality that often precludes the identification of unique parameters and states of the model that best represent available data. Hypothesizing that this non-uniqueness can convey useful information, we implemented a simplified simulation of a common differential diagnostic process (hypotension in an acute care setting), using a combination of a mathematical model of the cardiovascular system, a stochastic measurement model, and Bayesian inference techniques to quantify parameter and state uncertainty. The output of this procedure is a probability density function on the space of model parameters and initial conditions for a particular patient, based on prior population information together with patient-specific clinical observations. We show that multimodal posterior probability density functions arise naturally, even when unimodal and uninformative priors are used. The peaks of these densities correspond to clinically relevant differential diagnoses and can, in the simplified simulation setting, be constrained to a single diagnosis by assimilating additional observations from dynamical interventions (e.g., fluid challenge). We conclude that the ill-posedness of the inverse problem in quantitative physiology is not merely a technical obstacle, but rather reflects clinical reality and, when addressed adequately in the solution process, provides a novel link between mathematically described physiological knowledge and the clinical concept of differential diagnoses. We outline possible steps toward translating this computational approach to the bedside, to supplement today's evidence-based medicine with a quantitatively founded model-based medicine that integrates mechanistic knowledge with patient-specific information.
| S-EPMC2065888 | biostudies-literature
Project description:Recent published results in inverse scattering generally show the difficulty in dealing with moderate to high contrast inhomogeneities when employing linearized or iteratively linearized algorithms (e.g., distorted Born iterative method). This paper presents a fully nonlinear algorithm utilizing full wave field data, that results in ultrasound computed tomographic images from a laboratory breast scanner, and shows several such unique images from volunteer subjects. The forward problem, data collection process and inverse scattering algorithm used are discussed. A functional that represents the "best fit" between predicted and measured data is minimized, and therefore requires a very fast forward problem solver, Jacobian calculation, and gradient estimation, all of which are described. The data collection device is described. The algorithm and device yield quantitative estimates of human breast tissue in vivo. Several high resolution images, measuring ∼150 by 150 wavelengths, obtained from the 2D inverse scattering algorithms, using data collected from a first prototype, are shown and discussed. The quantitative values are compared with previous published work.
| S-EPMC3356315 | biostudies-other