Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “implicit”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A NUMERICAL TREATMENT OF ANISOTROPIC RADIATION FIELDS COUPLED WITH RELATIVISTIC RESISTIVE MAGNETOFLUIDS

We develop a numerical scheme for solving fully special relativistic, resistive radiation magnetohydrodynamics. Our code guarantees conservation of total mass, momentum, and energy. The radiation energy density and the radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation field, in contrast to the Eddington approximation, as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly the Harten-Lax-van Leer scheme; the HLLC and HLLD schemes are also solved for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in an optically thick regime or a large magnetic Reynolds number regime. Although we need to invert 4 Multiplication-Sign 4 matrices (for the gas-radiation interaction) and 3 Multiplication-Sign 3 matrices (for the magnetic energy dissipation) at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek-type magnetic reconnection, revealing the reduction of the reconnection rate via radiation drag.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The classification of ambiguity in polarimetric reconstruction of coronal mass ejection

The Thomson scattering theory indicates that there exist explicit and implicit ambiguities in polarimetric analyses of coronal mass ejection (CME) observations. We suggest a classification for these ambiguities in CME reconstruction. Three samples, including double explicit, mixed, and double implicit ambiguity, are shown with the polarimetric analyses of STEREO CME observations. These samples demonstrate that this classification is helpful for improving polarimetric reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

PARTICLE-GAS DYNAMICS WITH ATHENA: METHOD AND CONVERGENCE

The Athena magnetohydrodynamics code has been extended to integrate the motion of particles coupled with the gas via aerodynamic drag in order to study the dynamics of gas and solids in protoplanetary disks (PPDs) and the formation of planetesimals. Our particle-gas hybrid scheme is based on a second-order predictor-corrector method. Careful treatment of the momentum feedback on the gas guarantees exact conservation. The hybrid scheme is stable and convergent in most regimes relevant to PPDs. We describe a semi-implicit integrator generalized from the leap-frog approach. In the absence of drag force, it preserves the geometric properties of a particle orbit. We also present a fully implicit integrator that is unconditionally stable for all regimes of particle-gas coupling. Using our hybrid code, we study the numerical convergence of the nonlinear saturated state of the streaming instability. We find that gas flow properties are well converged with modest grid resolution (128 cells per pressure length {eta}r for dimensionless stopping time {tau} {sub s} = 0.1) and an equal number of particles and grid cells. On the other hand, particle clumping properties converge only at higher resolutions, and finer resolution leads to stronger clumping before convergence is reached. Finally, we find that the measurement of particle transport properties resulted from the streaming instability may be subject to error of about {+-}20%.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2014

AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Fully coupled simulation of cosmic reionization. I. numerical methods and tests

Here, we describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large similar to(100 Mpc)(3) cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. But, we employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use to advance the system of equations. Radiation transport is done in the gray flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. Finally, we illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with 3200(3) Eulerian grid cells and dark matter particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Shape-driven, emergent behavior in active particle mixtures

Abstract Active particle systems can vary greatly from one-component systems of spheres to mixtures of particle shapes at different composition ratios. We investigate computationally the combined effect of anisotropy and stoichiometry on the collective behavior of two-dimensional active colloidal mixtures of polygons. We uncover three emergent phenomena not yet reported in active Brownian particle systems. First, we find that mixtures containing hexagons exhibit micro-phase separation with large grains of hexagonal symmetry. We quantify a measurable, implicit ‘steric attraction’ between the active particles as a result of shape anisotropy and activity. This calculation provides further evidence that implicit interactions in active systems, even without explicit attraction, can lead to an effective preferential attraction between particles. Next, we report stable fluid clusters in mixtures containing one triangle or square component. We attribute the fluidization of the dense cluster to the interplay of cluster destabilizing particles, which introduce grain boundaries and slip planes into the system, causing solid-like clusters to break up into fluid clusters. Third, we show that fluid clusters can coexist with solid clusters within a sparse gas of particles in a steady state of three coexisting phases. Our results highlight the potential for a wide variety of behavior to be accessible to active matter systems and establish a route to control active colloidal systems through simple parameter designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

The different varieties of the Suyama-Yamaguchi consistency relation and its violation as a signal of statistical inhomogeneity

We present the different consistency relations that can be seen as variations of the well known Suyama-Yamaguchi (SY) consistency relation τ{sub NL}≥((6/5)f{sub NL}){sup 2}, the latter involving the levels of non-gaussianity f{sub NL} and τ{sub NL} in the primordial curvature perturbation ζ. It has been (implicitly) claimed that the following variation: τ{sub NL}(k{sub 1},k{sub 3})≥((6/5)){sup 2}f{sub NL}(k{sub 1})f{sub NL}(k{sub 3}), which we call ''the fourth variety'', in the collapsed (for τ{sub NL}) and squeezed (for f{sub NL}) limits is always satisfied independently of any physics; however, the proof depends sensitively on the assumption of scale-invariance (expressing this way the fourth variety of the SY consistency relation as τ{sub NL}≥((6/5)f{sub NL}){sup 2}) which only applies for cosmological models involving Lorentz-invariant scalar fields (at least at tree level), leaving room for a strong violation of this variety of the consistency relation when non-trivial degrees of freedom, for instance vector fields, are in charge of the generation of the primordial curvature perturbation. With this in mind as a motivation, we explicitly state, in the first part of this work, under which conditions the SY consistency relation has been claimed to hold in its different varieties (implicitly) presented in the literature since its inception back in 2008; as a result, we show for the first time that the variety τ{sub NL}(k{sub 1},k{sub 1})≥((6/5)f{sub NL}(k{sub 1})){sup 2}, which we call ''the fifth variety'', is always satisfied even when there is strong scale-dependence and high levels of statistical anisotropy as long as statistical homogeneity holds: thus, an observed violation of this specific variety would prevent the comparison between theory and observation, shaking this way the foundations of cosmology as a science. In the second part, we concern about the existence of non-trivial degrees of freedom, concretely vector fields for which the levels of non-gaussianity have been calculated for very few models; among them, and by making use of the δN formalism at tree level, we study a class of models that includes the vector curvaton scenario, vector inflation, and the hybrid inflation with coupled vector and scalar ''waterfall field'' where ζ is generated at the end of inflation, finding that the fourth variety of the SY consistency relation is indeed strongly violated for some specific wavevector configurations while the fifth variety continues to be well satisfied. Finally, as a byproduct of our investigation, we draw attention to a quite recently demonstrated variety of the SY consistency relation: τ{sup iso}{sub NL}≥((6/5)f{sup iso}{sub NL}){sup 2}, in scenarios where scalar and vector fields contribute to the generation of the primordial curvature perturbation; this variety of the SY consistency relation is satisfied although the isotropic pieces of the non-gaussianity parameters receive contributions from the vector fields. We discuss further implications for observational cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Locally adaptive activation functions with slope recovery for deep and physics-informed neural networks

Here we propose two approaches of locally adaptive activation functions namely, layer-wise and neuron-wise locally adaptive activation functions, which improve the performance of deep and physics-informed neural networks. The local adaptation of activation function is achieved by introducing a scalable parameter in each layer (layer-wise) and for every neuron (neuron-wise) separately, and then optimizing it using a variant of stochastic gradient descent algorithm. In order to further increase the training speed, an activation slope-based slope recovery term is added in the loss function, which further accelerates convergence, thereby reducing the training cost. On the theoretical side, we prove that in the proposed method, the gradient descent algorithms are not attracted to sub-optimal critical points or local minima under practical conditions on the initialization and learning rate, and that the gradient dynamics of the proposed method is not achievable by base methods with any (adaptive) learning rates. We further show that the adaptive activation methods accelerate the convergence by implicitly multiplying conditioning matrices to the gradient of the base method without any explicit computation of the conditioning matrix and the matrix–vector product. The different adaptive activation functions are shown to induce different implicit conditioning matrices. Furthermore, the proposed methods with the slope recovery are shown to accelerate the training process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Electromagnetic and neutral-weak response functions of 4 He and 12 C

A major goal of nuclear theory is to understand the strong interaction in nuclei as it manifests itself in terms of two- and many-body forces among the nuclear constituents, the protons and neutrons, and the interactions of these constituents with external electroweak probes via one- and many-body currents. The objective of the present work is to calculate the quasielastic electroweak response functions in light nuclei within the realistic dynamical framework outlined above. These response functions determine the inclusive cross section as function of the lepton momentum and energy transfers. Their ab initio calculation is a very challenging quantum many-body problem, since it requires summation over the entire excitation spectrum of the nucleus and inclusion in the electroweak currents of one- and many-body terms. Green's functions Monte Carlo methods allow one to circumvent both difficulties by computing the response in imaginary time (the so-called Euclidean response) and hence summing implicitly over the bound and continuum states of the nucleus, and by implementing specific algorithms designed to deal with the complicated spin-isospin structure of nuclear many-body operators. Theoretical predictions for 4 He and 12 C, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasielastic to the dip region and beyond. In conclusion, these results challenge the conventional picture of quasielastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Fast explicit solutions for neutrino-electron scattering: Explicit asymptotic methods

Here, we present results of explicit asymptotic approximations applied to neutrino-electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and time stepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each time step can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Big Data Analytics Framework for HPC Log Data: Three Case Studies Using the Titan Supercomputer Log

Reliability, availability and serviceability (RAS) logs of high performance computing (HPC) resources, when closely investigated in spatial and temporal dimensions, can provide invaluable information regarding system status, performance, and resource utilization. These data are often generated from multiple logging systems and sensors that cover many components of the system. The analysis of these data for finding persistent temporal and spatial insights faces two main difficulties: the volume of RAS logs makes manual inspection difficult and the unstructured nature and unique properties of log data produced by each subsystem adds another dimension of difficulty in identifying implicit correlation among recorded events. To address these issues, we recently developed a multi-user Big Data analytics framework for HPC log data at Oak Ridge National Laboratory (ORNL). This paper introduces three in-progress data analytics projects that leverage this framework to assess system status, mine event patterns, and study correlations between user applications and system events. We describe the motivation of each project and detail their workflows using three years of log data collected from ORNL's Titan supercomputer.

Park, Byung↗
DOE OSTI2015

Grandmaster: Interactive text-based analytics of social media

People use social media resources like Twitter, Facebook, forums etc. to share and discuss various activities or topics. By aggregating topic trends across many individuals using these services, we seek to construct a richer profile of a person’s activities and interests as well as provide a broader context of those activities. This profile may then be used in a variety of ways to understand groups as a collection of interests and affinities and an individual’s participation in those groups. Our approach considers that much of these data will be unstructured, free-form text. By analyzing free-form text directly, we may be able to gain an implicit grouping of individuals with shared interests based on shared conversation, and not on explicit social software linking them. In this paper, we discuss a proof-of-concept application called Grandmaster built to pull short sections of text, a person’s comments or Twitter posts, together by analysis and visualization to allow a gestalt understanding of the full collection of all individuals: how groups are similar and how they differ, based on their text inputs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A substantial role of soil erosion in the land carbon sink and its future changes

Realistic representation of land carbon sink in climate models is vital for predicting carbon climate feedbacks in a changing world. Although soil erosion that removes land organic carbon has increased substantially since the onset of agriculture, it is rarely included in the current generation of climate models. Using an Earth system model with soil erosion represented, we estimated that on average soil erosion displaces 5% of newly fixed land organic carbon downslope annually in the continental US. In the lower Mississippi river basin and the Cascades, the fraction can be as large as 40%. About 12% of the eroded organic carbon is eventually exported to inland waters, which is equal to 14% of the simulated net carbon gain by terrestrial ecosystems. By comparing the eroded organic carbon export to rivers with the particulate organic carbon export to oceans, we demonstrated that a large fraction of the carbon export to rivers could have been mineralized in inland waters. Importantly, with a direct comparison of eroded and exported soil organic carbon and land net carbon uptake, we found that Earth system models that ignore soil erosion likely offset the erosional carbon loss by increasing heterotrophic respiration implicitly. But as soil erosion and heterotrophic respiration respond differently to a warming climate, this unrealistic compensation would lead to biased predictions of future land carbon sink.

Tan, Zeli↗
DOE OSTI2017

Peridynamics using Irregular Domain Discretization with MOOSE-Based Implementation

In this paper, reformulation of classical bond-based peridynamic thermomechanical model for irregular domain decomposition and its MOOSE-based implicit formulation are presented. First, the irregular grid based peridynamic thermomechanical model is formulated and model parameters are derived. Following this, an implicit formulation for the solution of static or quasi-static problems is presented. Some aspects of the MOOSE-based implementation are given. After that, the formulation is verified against benchmark solutions for thermomechanic problems. Crack initiation and propagation in circular (2D) and cylindrical (3D) nuclear fuels at high temperature are studied using irregular grids.

42 ENGINEERING↗
DOE OSTI2013

Experimental and computational results on exciton/free-carrier ratio, hot/thermalized carrier diffusion, and linear/nonlinear rate constants affecting scintillator proportionality

Models of nonproportional response in scintillators have highlighted the importance of parameters such as branching ratios, carrier thermalization times, diffusion, kinetic order of quenching, associated rate constants, and radius of the electron track. For example, the fraction ηeh of excitations that are free carriers versus excitons was shown by Payne and coworkers to have strong correlation with the shape of electron energy response curves from Compton-coincidence studies. Rate constants for nonlinear quenching are implicit in almost all models of nonproportionality, and some assumption about track radius must invariably be made if one is to relate linear energy deposition dE/dx to volume-based excitation density n (eh/cm3) in terms of which the rates are defined. Diffusion, affecting time-dependent track radius and thus density of excitations, has been implicated as an important factor in nonlinear light yield. Several groups have recently highlighted diffusion of hot electrons in addition to thermalized carriers and excitons in scintillators. However, experimental determination of many of these parameters in the insulating crystals used as scintillators has seemed difficult. Subpicosecond laser techniques including interband z scan light yield, fluence-dependent decay time, and transient optical absorption are now yielding experimental values for some of the missing rates and ratios needed for modeling scintillator response. First principles calculations and Monte Carlo simulations can fill in additional parameters still unavailable from experiment. As a result, quantitative modeling of scintillator electron energy response from independently determined material parameters is becoming possible on an increasingly firmer data base. This paper describes recent laser experiments, calculations, and numerical modeling of scintillator response.

Williams, Richard↗
DOE OSTI2014

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-D-16A-07: Photobleaching Predicts Necrosis in Interstitial PDT

Purpose: Dosimetry for PDT has proven to be a challenge thus far, and for prediction of PDT outcome, a singlet oxygen model based on fundamental photophysical parameters has been developed. Previously, the photobleaching effect of photosensitizers was taken into account in the singlet oxygen explicit dosimetry model; here we report of direct measurements of photobleaching in the same model to assess the conditions under which implicit dosimetry using photobleaching can serve as an intermediate surrogate for PDT damage. Methods: Fluorescence spectra were measured interstitially in sensitized mouse tumors prior to after irradiation via a cylindrical diffuser. Photobleaching was determined by the relative decrease in fluorescence amplitude from the initial pre-treatment measurement. Spectra were analyzed by singular value decomposition to determine the photosensitizer concentration. Different photosensitizers were used to see the effect of photobleaching on PDT outcome and the impact of fluence on photobleaching. The drugs used were BPD (at two drug-light intervals), HPPH, and Photofrin. PDT outcome was determined by tumor necrosis radii measured upon sectioning and staining of treated tumors. Results: Post-PDT photosentizer concentrations were compared to initial pre-PDT photosensitizer concentrations, and the decrease was greater with a higher fluence measured during treatment. Furthermore, photobleaching and necrosis radius were found to be positively correlated. The relationship between photobleaching and necrosis radius is sensitizer-dependent, however the differences among sensitizers can be understood in terms of their respective photophysical parameters. Conclusions: Photobleaching is predictive of PDT outcome, but a comprehensive singlet oxygen model, has the potential to further improve the prediction of PDT outcome and the understanding of implicit dosimetry.

60 APPLIED LIFE SCIENCES↗