Modelling of submarine landslides and generated water waves

The slide of unstable sedimentary bodies and their hydraulic effects are studied by numerical means. A two-dimensional fluid mechanics model based on Navier-Stokes equations has been developed considering the sediments and water as a mixture. Viscoplastic and diffusion laws for the sediments have been introduced into the model. The numerical model is validated with an analytical solution for a Bingham flow. Laboratory experiments consisting in the slide of gravel mass have been carried out. The results of these experiments have shown the importance of the sediment rheology and the diffusion. The model parameters are adjusted by trial and error to match the observed “sandflow”.

Analytical method of predicating the instabilities of a micro arch-shaped beam under electrostatic loading

An arch-shaped beam with different configurations under electrostatic loading experiences either the direct pull-in instability or the snap-through first and then the pull-in instability. When the pull-in instability occurs, the system collides with the electrode and adheres to it, which usually causes the system failure. When the snap-through instability occurs, the system experiences a discontinuous displacement to flip over without colliding with the electrode. The snap-through instability is an ideal actuation mechanism because of the following reasons: (1) after snap-through the system regains the stability and capability of withstanding further loading; (2) the system flips back when the loading is reduced, i.e. the system can be used repetitively; and (3) when approaching snap-through instability the system effective stiffness reduces toward zero, which leads to a fast flipping-over response. To differentiate these two types of instability responses for an arch-shaped beam is vital for the actuator design. For an arch-shaped beam under electrostatic loading, the nonlinear terms of the mid-plane stretching and the electrostatic loading make the analytical solution extremely difficult if not impossible and the related numerical solution is rather complex. Using the one mode expansion approximation and the truncation of the higher-order terms of the Taylor series, we present an analytical solution here. However, the one mode approximation and the truncation error of the Taylor series can cause serious error in the solution. Therefore, an error-compensating mechanism is also proposed. The analytical results are compared with both the experimental data and the numerical multi-mode analysis. The analytical method presented here offers a simple yet efficient solution approach by retaining good accuracy to analyze the instability of an arch-shaped beam under electrostatic loading.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

Hybrid time dependent/independent solution for the He i line ratio temperature and density diagnostic for a thermal helium beam with applications in the scrape-off layer-edge regions in tokamaks

Spectroscopic studies of line emission intensities and ratios offer an attractive option in the\r\ndevelopment of non-invasive plasma diagnostics. Evaluating ratios of selected He I line\r\nemission profiles from the singlet and triplet neutral helium spin systems allows for simultaneous\r\nmeasurement of electron density (ne) and temperature (Te) profiles. Typically, this powerful\r\ndiagnostic tool is limited by the relatively long relaxation times of the 3S metastable term of helium\r\nthat populates the triplet spin system, and on which electron temperature sensitive lines are based.\r\nBy developing a time dependent analytical solution, we model the time evolution of the two spin\r\nsystems. We present a hybrid time dependent/independent line ratio solution that improves the\r\nrange of application of this diagnostic technique in the scrape-off layer (SOL) and edge plasma\r\nregions when comparing it against the current equilibrium line ratio helium model used at\r\nTEXTOR.

Applications of Statistical Modeling in Iterative CT Image Reconstruction

Thesis (Master's)--University of Washington, 2015

Analytical modelling of bacterial migration and accumulation with rate-limited sorption in discrete fractures

An analytical model for bacterial accumulation in a discrete fractllre has been developed. The transport and accumlllation processes incorporate into the model include advection, dispersion, rate-limited adsorption, rate-limited desorption, irreversible adsorption, attachment, detachment, growth and first order decay botl1 in sorbed and aqueous phases. An analytical solution in Laplace space is derived and nlln1erically inverted. The model is implemented in the code BIOFRAC vvhich is written in Fortran 99. The model is derived for two phases, Phase I, where adsorption-desorption are dominant, and Phase II, where attachment-detachment are dominant. Phase I ends yvhen enollgh bacteria to fully cover the substratllm have accllillulated. The model for Phase I vvas verified by comparing to the Ogata-Banks solution and the model for Phase II was verified by comparing to a nonHomogenous version of the Ogata-Banks solution. After verification, a sensitiv"ity analysis on the inpllt parameters was performed. The sensitivity analysis was condllcted by varying one inpllt parameter vvhile all others were fixed and observing the impact on the shape of the clirve describing bacterial concentration verSllS time. Increasing fracture apertllre allovvs more transport and thus more accllffilliation, "Vvhich diminishes the dllration of Phase I. The larger the bacteria size, the faster the sllbstratum will be covered. Increasing adsorption rate, was observed to increase the dllration of Phase I. Contrary to the aSSllmption ofllniform biofilm thickness, the accllffilliation starts frOll1 the inlet, and the bacterial concentration in aqlleous phase moving towards the olitiet declines, sloyving the accumulation at the outlet. Increasing the desorption rate, redllces the dliration of Phase I, speeding IIp the accllmlilation. It was also observed that Phase II is of longer duration than Phase I. Increasing the attachment rate lengthens the accliffililation period. High rates of detachment speeds up the transport. The grovvth and decay rates have no significant effect on transport, althollgh increases the concentrations in both aqueous and sorbed phases are observed. Irreversible adsorption can stop accllillulation completely if the vallIes are high.

Patient-specific anisotropic model of human trunk based on MR data

There are many ways to generate geometrical models for numerical simulation, and most of them start with a segmentation step to extract the boundaries of the regions of interest. This paper presents an algorithm to generate a patient-specific three-dimensional geometric model, based on a tetrahedral mesh, without an initial extraction of contours from the volumetric data. Using the information directly available in the data, such as gray levels, we built a metric to drive a mesh adaptation process. The metric is used to specify the size and orientation of the tetrahedral elements everywhere in the mesh. Our method, which produces anisotropic meshes, gives good results with synthetic and real MRI data. The resulting model quality has been evaluated qualitatively and quantitatively by comparing it with an analytical solution and with a segmentation made by an expert. Results show that our method gives, in 90% of the cases, as good or better meshes as a similar isotropic method, based on the accuracy of the volume reconstruction for a given mesh size. Moreover, a comparison of the Hausdorff distances between adapted meshes of both methods and ground-truth volumes shows that our method decreases reconstruction errors faster. Copyright © 2015 John Wiley & Sons, Ltd.

Exact Solution of the Nonlinear Dynamics of Recurrent Neural Mechanisms for Direction Selectivity

Different theoretical models have tried to investigate the feasibility of recurrent neural mechanisms for achieving direction selectivity in the visual cortex. The mathematical analysis of such models has been restricted so far to the case of purely linear networks. We present an exact analytical solution of the nonlinear dynamics of a class of direction selective recurrent neural models with threshold nonlinearity. Our mathematical analysis shows that such networks have form-stable stimulus-locked traveling pulse solutions that are appropriate for modeling the responses of direction selective cortical neurons. Our analysis shows also that the stability of such solutions can break down giving raise to a different class of solutions ("lurching activity waves") that are characterized by a specific spatio-temporal periodicity. These solutions cannot arise in models for direction selectivity with purely linear spatio-temporal filtering.

Analytical and semi-analytical analysis of an artificial satellite's rotational motion

The rotational motion of an artificial satellite is studied by considering torques produced by gravity gradient and direct solar radiation pressure. A satellite of circular cylinder shape is considered here, and Andoyers variables are used to describe the rotational motion. Expressions for direct solar radiation torque are derived. When the earth's shadow is not considered, an analytical solution is obtained using Lagrange's method of variation of parameters. A semi-analytical procedure is proposed to predict the satellite's attitude under the influence of the earth's shadow. The analytical solution shows that angular variables are linear and periodic functions of time while their conjugates suffer only periodic variations. When compared, numerical and analytical solutions have a good agreement during the time range considered.

Analytical attitude prediction of spin stabilized spacecrafts perturbed by magnetic residual torque

An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques. Assuming an inclined dipole model for the Earth's magnetic field, an analytical averaging method is applied to obtain the mean residual torque every orbital period. The orbit mean anomaly is utilized to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, and non-circular orbits, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Spin-Stabilized Spacecrafts: Analytical Attitude Propagation Using Magnetic Torques

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

RR3D: Uma solução para renderização remota de imagens médicas tridimensionais

The visualization of three-dimensional(3D)images is increasigly being sed in the area of medicine, helping physicians diagnose desease. the advances achived in scaners esed for acquisition of these 3d exames, such as computerized tumography(CT) and Magnetic Resonance imaging (MRI), enable the generation of images with higher resolutions, thus, generating files with much larger sizes. Currently, the images of computationally expensive one, and demanding the use of a righ and computer for such task. The direct remote acess of these images thruogh the internet is not efficient also, since all images have to be trasferred to the user´s equipment before the 3D visualization process ca start. with these problems in mind, this work proposes and analyses a solution for the remote redering of 3D medical images, called Remote Rendering (RR3D). In RR3D, the whole hedering process is pefomed a server or a cluster of servers, with high computational power, and only the resulting image is tranferred to the client, still allowing the client to peform operations such as rotations, zoom, etc. the solution was developed using web services written in java and an architecture that uses the scientific visualization packcage paraview, the framework paraviewWeb and the PACS server DCM4CHEE.The solution was tested with two scenarios where the rendering process was performed by a sever with graphics hadwere (GPU) and by a server without GPUs. In the scenarios without GPUs, the soluction was executed in parallel with several number of cores (processing units)dedicated to it. In order to compare our solution to order medical visualization application, a third scenario was esed in the rendering process, was done locally. In all tree scenarios, the solution was tested for different network speeds. The solution solved satisfactorily the problem with the delay in the transfer of the DICOM files, while alowing the use of low and computers as client for visualizing the exams even, tablets and smart phones

Formation of soliton trains in Bose-Einstein condensates as a nonlinear Fresnel diffraction of matter waves

The problem of generation of atomic soliton trains in elongated Bose-Einstein condensates is considered in framework of Whitham theory of modulations of nonlinear waves. Complete analytical solution is presented for the case when the initial density distribution has sharp enough boundaries. In this case the process of soliton train formation can be viewed as a nonlinear Fresnel diffraction of matter waves. Theoretical predictions are compared with results of numerical simulations of one- and three-dimensional Gross-Pitaevskii equation and with experimental data on formation of Bose-Einstein bright solitons in cigar-shaped traps. (C) 2003 Elsevier B.V. All rights reserved.

Salting kinetics and salt diffusivities in farmed Pantanal caiman muscle

A criação de jacaré do Pantanal (Caiman crocodilus yacare) em cativeiro tem sido estimulada, e entre as técnicas de processamento de sua carne, a salga é um processo de conservação relativamente simples e de baixo custo. O objetivo deste trabalho foi estudar a cinética de difusão de cloreto de sódio em carne de jacaré do Pantanal criado em cativeiro, durante a salga úmida. Foram utilizados volumes limitados de salmoura e os experimentos foram realizados com relações salmoura/músculo de 3, 4 e 5, com concentrações de salmoura de 15%, 20% e 25% em peso e temperaturas de 10, 15 e 20ºC. A solução analítica da segunda lei de Fick, considerando difusão unidimensional em uma placa infinita em contato com uma solução bem agitada de volume limitado, foi utilizada para calcular os coeficientes de difusão efetivos de sal e estimar o conteúdo de cloreto de sódio nos filés. Obteve-se boa concordância entre o modelo analítico considerado e os dados experimentais. As difusividades do sal nos filés ocorreram na faixa de 0,47x10-10 a 9,62x10-10 m²/s.

Thermodynamic properties of water sorption of jackfruit (Artocarpus heterophyllus Lam.) as a function of moisture content

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)