Numerical 3D modeling of heat transfer in human tissues for microwave radiometry monitoring of Brown fat metabolismo

Background: Brown adipose tissue (BAT) plays an important role in whole body metabolism and could potentially mediate weight gain and insulin sensitivity. Although some imaging techniques allow BAT detection, there are currently no viable methods for continuous acquisition of BAT energy expenditure. We present a non-invasive technique for long term monitoring of BAT metabolism using microwave radiometry. Methods: A multilayer 3D computational model was created in HFSS™ with 1.5 mm skin, 3-10 mm subcutaneous fat, 200 mm muscle and a BAT region (2-6 cm3) located between fat and muscle. Based on this model, a log-spiral antenna was designed and optimized to maximize reception of thermal emissions from the target (BAT). The power absorption patterns calculated in HFSS™ were combined with simulated thermal distributions computed in COMSOL® to predict radiometric signal measured from an ultra-low-noise microwave radiometer. The power received by the antenna was characterized as a function of different levels of BAT metabolism under cold and noradrenergic stimulation. Results: The optimized frequency band was 1.5-2.2 GHz, with averaged antenna efficiency of 19%. The simulated power received by the radiometric antenna increased 2-9 mdBm (noradrenergic stimulus) and 4-15 mdBm (cold stimulus) corresponding to increased 15-fold BAT metabolism. Conclusions: Results demonstrated the ability to detect thermal radiation from small volumes (2-6 cm3) of BAT located up to 12 mm deep and to monitor small changes (0.5°C) in BAT metabolism. As such, the developed miniature radiometric antenna sensor appears suitable for non-invasive long term monitoring of BAT metabolism.

3D Surface Modeling From Point Clouds

The goal of this thesis is to implement software for creating 3D models from point clouds. Point clouds are acquired with stereo cameras, monocular systems or laser scanners. The created 3D models are triangular models or NURBS (Non-Uniform Rational B-Splines) models. Triangular models are constructed from selected areas from the point clouds and resulted triangular models are translated into a set of quads. The quads are further translated into an estimated grid structure and used for NURBS surface approximation. Finally, we have a set of NURBS surfaces which represent the whole model. The problem wasn’t so easy to solve. The selected triangular surface reconstruction algorithm did not deal well with noise in point clouds. To handle this problem, a clustering method is introduced for simplificating the model and removing noise. As we had better results with the smaller point clouds produced by clustering, we used points in clusters to better estimate the grids for NURBS models. The overall results were good when the point cloud did not have much noise. The point clouds with small amount of error had good results as the triangular model was solid. NURBS surface reconstruction performed well on solid models.

Hyperresolution global land surface modeling: meeting a grand challenge for monitoring Earth's terrestrial water

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (∼10–100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface‐subsurface interactions due to fine‐scale topography and vegetation; improved representation of land‐atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 109 unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a “grand challenge” to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

Characterization of glycerol kinase from baker's yeast: Response surface modeling of the enzymatic reaction

The present study describes a methodology of dosage of glycerol kinase (GK) from baker's yeast. The standardization of the activity of the glycerol kinase from baker's yeast was accomplished using the diluted enzymatic preparation containing glycerol phosphate oxidase (GPO) and glycerol kinase. The mixture was incubated at 60 degrees C by 15 min and the reaction was stopped by the SDS solution addition. A first set of experiments was carried out in order to investigate the individual effect of temperature (7), pH and substrate concentration (S), on GK activity and stability. The pH and temperature stability tests showed that the enzyme presented a high stability to pH 6.0-8.0 and the thermal stability were completely maintained up to 50 degrees C during 1 h. The K(m) of the enzyme for glycerol was calculated to be 2 mM and V(max) to be 1.15 U/mL. In addition, modeling and optimization of reaction conditions was attempted by response surface methodology (RSM). Higher activity values will be attained at temperatures between 52 and 56 degrees C, pH around 10.2-10.5 and substrate concentrations from 150 to 170 mM.This low cost method for glycerol kinase dosage in a sequence of reactions is of great importance for many industries, like food, sugar and alcohol. RSM showed to be an adequate approach for modeling the reaction and optimization of reaction conditions to maximize glycerol kinase activity. (C) 2007 Elsevier B.V. All rights reserved.

Comparative analysis between experimental characterization results and numerical fdtd modeling of self-assembled photonic crystals

This paper presents a comparative analysis between the experimental characterization and the numerical simulation results for a three-dimensional FCC photonic crystal (PhC) based on a self-assembly synthesis of monodispersive latex spheres. Specifically, experimental optical characterization, by means of reflectance measurements under variable angles over the lattice plane family [1,1, 1], are compared to theoretical calculations based on the Finite Di®erence Time Domain (FDTD) method, in order to investigate the correlation between theoretical predictions and experimental data. The goal is to highlight the influence of crystal defects on the achieved performance.

Numerical optimization,modeling and system evaluation of a thermophotovoltaic hybrid panel

Photovoltaic (PV) solar panels generally produce electricity in the 6% to 16% efficiency range, the rest being dissipated in thermal losses. To recover this amount, hybrid photovoltaic thermal systems (PVT) have been devised. These are devices that simultaneously convert solar energy into electricity and heat. It is thus interesting to study the PVT system globally from different point of views in order to evaluate advantages and disadvantages of this technology and its possible uses. In particular in Chapter II, the development of the PVT absorber numerical optimization by a genetic algorithm has been carried out analyzing different internal channel profiles in order to find a right compromise between performance and technical and economical feasibility. Therefore in Chapter III ,thanks to a mobile structure built into the university lab, it has been compared experimentally electrical and thermal output power from PVT panels with separated photovoltaic and solar thermal productions. Collecting a lot of experimental data based on different seasonal conditions (ambient temperature,irradiation, wind...),the aim of this mobile structure has been to evaluate average both thermal and electrical increasing and decreasing efficiency values obtained respect to separate productions through the year. In Chapter IV , new PVT and solar thermal equation based models in steady state conditions have been developed by software Dymola that uses Modelica language. This permits ,in a simplified way respect to previous system modelling softwares, to model and evaluate different concepts about PVT panel regarding its structure before prototyping and measuring it. Chapter V concerns instead the definition of PVT boundary conditions into a HVAC system . This was made trough year simulations by software Polysun in order to finally assess the best solar assisted integrated structure thanks to F_save(solar saving energy)factor. Finally, Chapter VI presents the conclusion and the perspectives of this PhD work.

Hypnotic and opioid anesthetic drug interactions on the CNS, focus on response surface modeling

This chapter will present the conceptual and applied approaches to capture the interaction of anesthetic hypnotic drugs with opioid drugs, as used in the clinical anesthetic state. The graphic and mathematical approaches used to capture hypnotic/opiate anesthetic drug interactions will be presented. This chapter is not a review article about interaction modeling, but focuses on specific drug interactions within a quite narrow field, anesthesia.

Response surface modeling of the interaction between propofol and sevoflurane

BACKGROUND: Propofol and sevoflurane display additivity for gamma-aminobutyric acid receptor activation, loss of consciousness, and tolerance of skin incision. Information about their interaction regarding electroencephalographic suppression is unavailable. This study examined this interaction as well as the interaction on the probability of tolerance of shake and shout and three noxious stimulations by using a response surface methodology. METHODS: Sixty patients preoperatively received different combined concentrations of propofol (0-12 microg/ml) and sevoflurane (0-3.5 vol.%) according to a crisscross design (274 concentration pairs, 3 to 6 per patient). After having reached pseudo-steady state, the authors recorded bispectral index, state and response entropy and the response to shake and shout, tetanic stimulation, laryngeal mask airway insertion, and laryngoscopy. For the analysis of the probability of tolerance by logistic regression, a Greco interaction model was used. For the separate analysis of bispectral index, state and response entropy suppression, a fractional Emax Greco model was used. All calculations were performed with NONMEM V (GloboMax LLC, Hanover, MD). RESULTS: Additivity was found for all endpoints, the Ce(50, PROP)/Ce(50, SEVO) for bispectral index suppression was 3.68 microg. ml(-1)/ 1.53 vol.%, for tolerance of shake and shout 2.34 microg . ml(-1)/ 1.03 vol.%, tetanic stimulation 5.34 microg . ml(-1)/ 2.11 vol.%, laryngeal mask airway insertion 5.92 microg. ml(-1) / 2.55 vol.%, and laryngoscopy 6.55 microg. ml(-1)/2.83 vol.%. CONCLUSION: For both electroencephalographic suppression and tolerance to stimulation, the interaction of propofol and sevoflurane was identified as additive. The response surface data can be used for more rational dose finding in case of sequential and coadministration of propofol and sevoflurane.

Precisão de plantas planialtimétricas em função do posicionamento de pontos em irradiadas de campo

The present experiment aimed to study the influence of point positioning in an irradiated field to produce a planialtimetric plant. A planialtimetric evaluation was carried out in a 14-acre experimental area with well-defined topographic variations. Planialtimetric maps were designed using manual procedures, Datageosis and Topoesalq. Datageosis built all the curves after numerical surface modeling. Topoesalq provided only height reports and the drawing of curves was done manually. The third method was a manual procedure. Because there were planialtimetry representation differences, longitudinal profiles were used in the sites where there was a great divergence among plants. When obtained profiles and plants were compared, it was verified that the one produced by Datageosis represented the relief plant better. Later, only the irradiated field points were evaluated and each point presented positioned readings before and after each relief map variation. The processing through the three methods resulted in significant plants of the local planialtimetry, according to the control profiles. It was concluded that the planning of the field procedure should be suitable to the posterior treatment method of obtained data in order to make a planialtimetric plant to accord to the evaluated local topography.

Modelagem numérica da superfície como instrumento de avaliação de curvas de nível na preservação de manacial de água

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Semi-automatic modeling of buildings from digital surface models

Semi-automatic building detection and extraction is a topic of growing interest due to its potential application in such areas as cadastral information systems, cartographic revision, and GIS. One of the existing strategies for building extraction is to use a digital surface model (DSM) represented by a cloud of known points on a visible surface, and comprising features such as trees or buildings. Conventional surface modeling using stereo-matching techniques has its drawbacks, the most obvious being the effect of building height on perspective, shadows, and occlusions. The laser scanner, a recently developed technological tool, can collect accurate DSMs with high spatial frequency. This paper presents a methodology for semi-automatic modeling of buildings which combines a region-growing algorithm with line-detection methods applied over the DSM.

Boulder deposition during major tsunami events

A remarkable accumulation of marine boulders located above the present spring tide level has occurred in two coastal lowlands of the Algarve (Portugal). The size-interval of the particles studied here is seldom reported in the literature in association with extreme events of coastal inundation, thus making this study of relevance to many other coasts worldwide. The spreads of boulders extend several hundred meters inland and well beyond the present landward limit of storm activity. The marine origin of the boulders is demonstrated by well-developed macro-bioerosion sculpturing and in situ skeletal remains of endolithic shallow marine bivalves. The good state preservation of the fossils within the boulders indicates that abrasion duringtransport and redeposition was not significant. We envisage boulder deposition as having taken place during the Lisbon tsunami of ad 1755 through the simultaneous landward entrainment of coarse particles from nearshore followed by rapid shoreward suspended-dominated transport and non-graded redeposition that excluded significant sorting by weight or boulder dimensions. We use numerical hydrodynamic modeling of tsunami (and storm) waves to test the observational data on boulder dimensions (density, size, distribution) on the most likely processes of sediment deposition. This work demonstrates the effectiveness of the study of boulder deposits in tsunami reconstruction. Copyright (C) 2011 John Wiley & Sons, Ltd.

What is the shape of an air bubble on a liquid surface?

We have calculated the equilibrium shape of the axially symmetric meniscus along which a spherical bubble contacts a flat liquid surface by analytically integrating the Young-Laplace equation in the presence of gravity, in the limit of large Bond numbers. This method has the advantage that it provides semianalytical expressions for key geometrical properties of the bubble in terms of the Bond number. Results are in good overall agreement with experimental data and are consistent with fully numerical (Surface Evolver) calculations. In particular, we are able to describe how the bubble shape changes from hemispherical, with a flat, shallow bottom, to lenticular, with a deeper, curved bottom, as the Bond number is decreased.

Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.