The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Improved detection in CDMA for biased sources

We consider the detection of biased information sources in the ubiquitous code-division multiple-access (CDMA) scheme. We propose a simple modification to both the popular single-user matched-filter detector and a recently introduced near-optimal message-passing-based multiuser detector. This modification allows for detecting modulated biased sources directly with no need for source coding. Analytical results and simulations with excellent agreement are provided, demonstrating substantial improvement in bit error rate in comparison with the unmodified detectors and the alternative of source compression. The robustness of error-performance improvement is shown under practical model settings, including bias estimation mismatch and finite-length spreading codes. © 2007 IOP Publishing Ltd.

A variational method and approximations of a Cauchy problem for elliptic equations

A Cauchy problem for general elliptic second-order linear partial differential equations in which the Dirichlet data in H½(?1 ? ?3) is assumed available on a larger part of the boundary ? of the bounded domain O than the boundary portion ?1 on which the Neumann data is prescribed, is investigated using a conjugate gradient method. We obtain an approximation to the solution of the Cauchy problem by minimizing a certain discrete functional and interpolating using the finite diference or boundary element method. The minimization involves solving equations obtained by discretising mixed boundary value problems for the same operator and its adjoint. It is proved that the solution of the discretised optimization problem converges to the continuous one, as the mesh size tends to zero. Numerical results are presented and discussed.

High performance shift invariant motion estimation and compensation in wavelet domain video compression

The contributions of this dissertation are in the development of two new interrelated approaches to video data compression: (1) A level-refined motion estimation and subband compensation method for the effective motion estimation and motion compensation. (2) A shift-invariant sub-decimation decomposition method in order to overcome the deficiency of the decimation process in estimating motion due to its shift-invariant property of wavelet transform. ^ The enormous data generated by digital videos call for an intense need of efficient video compression techniques to conserve storage space and minimize bandwidth utilization. The main idea of video compression is to reduce the interpixel redundancies inside and between the video frames by applying motion estimation and motion compensation (MEMO) in combination with spatial transform coding. To locate the global minimum of the matching criterion function reasonably, hierarchical motion estimation by coarse to fine resolution refinements using discrete wavelet transform is applied due to its intrinsic multiresolution and scalability natures. ^ Due to the fact that most of the energies are concentrated in the low resolution subbands while decreased in the high resolution subbands, a new approach called level-refined motion estimation and subband compensation (LRSC) method is proposed. It realizes the possible intrablocks in the subbands for lower entropy coding while keeping the low computational loads of motion estimation as the level-refined method, thus to achieve both temporal compression quality and computational simplicity. ^ Since circular convolution is applied in wavelet transform to obtain the decomposed subframes without coefficient expansion, symmetric-extended wavelet transform is designed on the finite length frame signals for more accurate motion estimation without discontinuous boundary distortions. ^ Although wavelet transformed coefficients still contain spatial domain information, motion estimation in wavelet domain is not as straightforward as in spatial domain due to the shift variance property of the decimation process of the wavelet transform. A new approach called sub-decimation decomposition method is proposed, which maintains the motion consistency between the original frame and the decomposed subframes, improving as a consequence the wavelet domain video compressions by shift invariant motion estimation and compensation. ^

Seismic performance of hybrid fiber reinforced polymer-concrete pier columns

As part of a multi-university research program funded by NSF, a comprehensive experimental and analytical study of seismic behavior of hybrid fiber reinforced polymer (FRP)-concrete column is presented in this dissertation. Experimental investigation includes cyclic tests of six large-scale concrete-filled FRP tube (CFFT) and RC columns followed by monotonic flexural tests, a nondestructive evaluation of damage using ultrasonic pulse velocity in between the two test sets and tension tests of sixty-five FRP coupons. Two analytical models using ANSYS and OpenSees were developed and favorably verified against both cyclic and monotonic flexural tests. The results of the two methods were compared. A parametric study was also carried out to investigate the effect of three main parameters on primary seismic response measures. The responses of typical CFFT columns to three representative earthquake records were also investigated. The study shows that only specimens with carbon FRP cracked, whereas specimens with glass or hybrid FRP did not show any visible cracks throughout cyclic tests. Further monotonic flexural tests showed that carbon specimens both experienced flexural cracks in tension and crumpling in compression. Glass or hybrid specimens, on the other hand, all showed local buckling of FRP tubes. Compared with conventional RC columns, CFFT column possesses higher flexural strength and energy dissipation with an extended plastic hinge region. Among all CFFT columns, the hybrid lay-up demonstrated the highest flexural strength and initial stiffness, mainly because of its high reinforcement index and FRP/concrete stiffness ratio, respectively. Moreover, at the same drift ratio, the hybrid lay-up was also considered as the best in term of energy dissipation. Specimens with glassfiber tubes, on the other hand, exhibited the highest ductility due to better flexibility of glass FRP composites. Furthermore, ductility of CFFTs showed a strong correlation with the rupture strain of FRP. Parametric study further showed that different FRP architecture and rebar types may lead to different failure modes for CFFT columns. Transient analysis of strong ground motions showed that the column with off-axis nonlinear filament-wound glass FRP tube exhibited a superior seismic performance to all other CFFTs. Moreover, higher FRP reinforcement ratios may lead to a brittle system failure, while a well-engineered FRP reinforcement configuration may significantly enhance the seismic performance of CFFT columns.

Seismic Performance of Hybrid Fiber Reinforced Polymer-Concrete Pier Columns

As part of a multi-university research program funded by NSF, a comprehensive experimental and analytical study of seismic behavior of hybrid fiber reinforced polymer (FRP)-concrete column is presented in this dissertation. Experimental investigation includes cyclic tests of six large-scale concrete-filled FRP tube (CFFT) and RC columns followed by monotonic flexural tests, a nondestructive evaluation of damage using ultrasonic pulse velocity in between the two test sets and tension tests of sixty-five FRP coupons. Two analytical models using ANSYS and OpenSees were developed and favorably verified against both cyclic and monotonic flexural tests. The results of the two methods were compared. A parametric study was also carried out to investigate the effect of three main parameters on primary seismic response measures. The responses of typical CFFT columns to three representative earthquake records were also investigated. The study shows that only specimens with carbon FRP cracked, whereas specimens with glass or hybrid FRP did not show any visible cracks throughout cyclic tests. Further monotonic flexural tests showed that carbon specimens both experienced flexural cracks in tension and crumpling in compression. Glass or hybrid specimens, on the other hand, all showed local buckling of FRP tubes. Compared with conventional RC columns, CFFT column possesses higher flexural strength and energy dissipation with an extended plastic hinge region. Among all CFFT columns, the hybrid lay-up demonstrated the highest flexural strength and initial stiffness, mainly because of its high reinforcement index and FRP/concrete stiffness ratio, respectively. Moreover, at the same drift ratio, the hybrid lay-up was also considered as the best in term of energy dissipation. Specimens with glassfiber tubes, on the other hand, exhibited the highest ductility due to better flexibility of glass FRP composites. Furthermore, ductility of CFFTs showed a strong correlation with the rupture strain of FRP. Parametric study further showed that different FRP architecture and rebar types may lead to different failure modes for CFFT columns. Transient analysis of strong ground motions showed that the column with off-axis nonlinear filament-wound glass FRP tube exhibited a superior seismic performance to all other CFFTs. Moreover, higher FRP reinforcement ratios may lead to a brittle system failure, while a well-engineered FRP reinforcement configuration may significantly enhance the seismic performance of CFFT columns.

Performance and Assessment of Rehabilitated Steel Culverts

This thesis explores the effects of rehabilitation on the structural performance of corrugated steel culverts. A full-scale laboratory experiment investigated the effects of grouted slip-liners on the performance of two buried circular corrugated steel culverts. One culvert was slip-lined and grouted using low strength grout, while the other was slip-lined and grouted using high strength grout. The performances of the culverts were measured before and after rehabilitation under service loads using single wheel pair loading at 0.45m of cover. Then, the rehabilitated culverts were loaded to their ultimate limit states. Results showed that the low and high strength grouted slip-liners provided strength well beyond requirements, with the low strength specimen failing at a load 2.4 times the fully factored service load, while the high strength specimen did not reach an ultimate limit state before bearing failure of the soil stopped testing. Results also showed that the low strength specimen behaved rigidly under service loads and flexibly under higher loads, while the high strength specimen behaved rigidly under all loads. A second full-scale experiment investigated the effect of a paved invert rehabilitation procedure on the performance of a deteriorated horizontal ellipse culvert. The performance of the culvert before and after rehabilitation was examined under service loads using tandem axle loading at 0.45m of cover. The rehabilitated culvert was then loaded up to its ultimate limit state. The culvert failed due to the formation of a plastic hinge at the West shoulder, while the paved invert cracked at the invert. Results showed that the rehabilitation increased the structural performance of the culvert, increasing the system stiffness and reducing average strains and local bending at critical locations in the culvert under service loads. A sustainability rating tool specifically for the evaluation of deteriorated culvert replacement or rehabilitation projects was also developed. A module for an existing tool, called GoldSET, was created and tested using two case studies, each comparing the replacement of a culvert using a traditional open-cut method with two trenchless rehabilitation techniques. In each case, the analyses showed that the trenchless techniques were the better alternatives in terms of sustainability.

Aspects of Chiral Symmetry Breaking in Lattice QCD

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

Sloshing and slamming oscillations in collapsible channel flow

We consider laminar high-Reynolds-number flow through a finite-length planar channel, where a portion of one wall is replaced by a thin massless elastic membrane that is held under longitudinal tension T and subject to an external pressure distribution. The flow is driven by a fixed pressure drop along the full length of the channel. We investigate the global stability of two-dimensional Poiseuille flow using a method of matched local eigenfunction expansions, which is compared to direct numerical simulations. We trace the neutral stability curve of the primary oscillatory instability of the system, illustrating a transition from high-frequency ‘sloshing’ oscillations at high T to vigorous ‘slamming’ motion at low T . Small-amplitude sloshing at high T can be captured using a low-order eigenmode truncation involving four surface-based modes in the compliant segment of the channel coupled to Womersley flow in the rigid segments. At lower tensions, we show that hydrodynamic modes contribute increasingly to the global instability and we demonstrate a change in the mechanism of energy transfer from the mean flow, with viscous effects being destabilising. Simulations of finite-amplitude oscillations at low T reveal a generic slamming motion, in which the the flexible membrane is drawn close to the opposite rigid wall before rapidly recovering. A simple model is used to demonstrate how fluid inertia in the downstream rigid channel segment, coupled to membrane curvature downstream of the moving constriction, together control slamming dynamics.

Ordered mesoporous titanium oxide for thin film microbatteries with enhanced lithium storage

A 3D mesoporous TiO2 material with well-developed mesostructure is prepared in the form of a binder-free thin (100 nm) film and studied as potential candidate for the negative electrode in lithium microbatteries. By appropriate thermal treatments, the selected crystal structure (anatase, rutile, or amorphous), and micro-/mesostructure of the materials was obtained. The effects of voltage window and prelithiation treatment improved first cycle reversibility up to 86% and capacity retention of 90% over 100 cycles. After a prolonged intercalation of lithium ions in ordered mesoporous TiO2 appeared small particles assigned to Li2Ti2O4 with cubic structure as observed from ex-situ TEM micrographs. This study highlights the flexibility of the potential window to which the electrode can operate. Maximum capacity values over 100 cycles of 470 μA h cm−2 μm−1 and 177 μA h cm−2 μm−1 are obtained for voltage ranges of 0.1–2.6 V and 1.0–2.6 V, respectively. The observed values are between 6 and 2 times higher than those obtained for films with 600 nm (80 μA h cm−2 μm−1) and 900 nm (92 μA h cm−2 μm−1) lengths. This indicates that 100 nm thin TiO2 films with high accessibility show finite-length type diffusion which is interesting for this particular application.

Root canal filling: fracture strength of fiber-reinforced composite-restored roots and finite element analysis

The aims of this study were to evaluate the effect of root canal filling techniques on root fracture resistance and to analyze, by finite element analysis (FEA), the expansion of the endodontic sealer in two different root canal techniques. Thirty single-rooted human teeth were instrumented with rotary files to a standardized working length of 14 mm. The specimens were embedded in acrylic resin using plastic cylinders as molds, and allocated into 3 groups (n=10): G(lateral) - lateral condensation; G(single-cone) - single cone; G(tagger) - Tagger's hybrid technique. The root canals were prepared to a length of 11 mm with the #3 preparation bur of a tapered glass fiber-reinforced composite post system. All roots received glass fiber posts, which were adhesively cemented and a composite resin core was built. All groups were subjected to a fracture strength test (1 mm/min, 45°). Data were analyzed statistically by one-way ANOVA with a significance level of 5%. FEA was performed using two models: one simulated lateral condensation and Tagger's hybrid technique, and the other one simulated the single-cone technique. The second model was designed with an amount of gutta-percha two times smaller and a sealer layer two times thicker than the first model. The results were analyzed using von Mises stress criteria. One-way ANOVA indicated that the root canal filling technique affected the fracture strength (p=0.004). The G(lateral) and G(tagger) produced similar fracture strength values, while G(single-cone) showed the lowest values. The FEA showed that the single-cone model generated higher stress in the root canal walls. Sealer thickness seems to influence the fracture strength of restored endodontically treated teeth.

Three-Dimensional Finite Element Analysis of Stress Distribution on Different Bony Ridges With Different Lengths of Morse Taper Implants and Prosthesis Dimensions

This finite element analysis (FEA) compared stress distribution on different bony ridges rehabilitated with different lengths of morse taper implants, varying dimensions of metal-ceramic crowns to maintain the occlusal alignment. Three-dimensional FE models were designed representing a posterior left side segment of the mandible: group control, 3 implants of 11 mm length; group 1, implants of 13 mm, 11 mm and 5 mm length; group 2, 1 implant of 11 mm and 2 implants of 5 mm length; and group 3, 3 implants of 5 mm length. The abutments heights were 3.5 mm for 13- and 11-mm implants (regular), and 0.8 mm for 5-mm implants (short). Evaluation was performed on Ansys software, oblique loads of 365N for molars and 200N for premolars. There was 50% higher stress on cortical bone for the short implants than regular implants. There was 80% higher stress on trabecular bone for the short implants than regular implants. There was higher stress concentration on the bone region of the short implants neck. However, these implants were capable of dissipating the stress to the bones, given the applied loads, but achieving near the threshold between elastic and plastic deformation to the trabecular bone. Distal implants and/or with biggest occlusal table generated greatest stress regions on the surrounding bone. It was concluded that patients requiring short implants associated with increased proportions implant prostheses need careful evaluation and occlusal adjustment, as a possible overload in these short implants, and even in regular ones, can generate stress beyond the physiological threshold of the surrounding bone, compromising the whole system.

Finite element modeling of straightening of thin-walled seamless tubes of austenitic stainless steel

During this thesis work a coupled thermo-mechanical finite element model (FEM) was builtto simulate hot rolling in the blooming mill at Sandvik Materials Technology (SMT) inSandviken. The blooming mill is the first in a long line of processes that continuously or ingotcast ingots are subjected to before becoming finished products. The aim of this thesis work was twofold. The first was to create a parameterized finiteelement (FE) model of the blooming mill. The commercial FE software package MSCMarc/Mentat was used to create this model and the programing language Python was used toparameterize it. Second, two different pass schedules (A and B) were studied and comparedusing the model. The two pass series were evaluated with focus on their ability to healcentreline porosity, i.e. to close voids in the centre of the ingot. This evaluation was made by studying the hydrostatic stress (σm), the von Mises stress (σeq)and the plastic strain (εp) in the centre of the ingot. From these parameters the stress triaxiality(Tx) and the hydrostatic integration parameter (Gm) were calculated for each pass in bothseries using two different transportation times (30 and 150 s) from the furnace. The relationbetween Gm and an analytical parameter (Δ) was also studied. This parameter is the ratiobetween the mean height of the ingot and the contact length between the rolls and the ingot,which is useful as a rule of thumb to determine the homogeneity or penetration of strain for aspecific pass. The pass series designed with fewer passes (B), many with greater reduction, was shown toachieve better void closure theoretically. It was also shown that a temperature gradient, whichis the result of a longer holding time between the furnace and the blooming mill leads toimproved void closure.

Finite element analysis of hot rolling in the blooming mill

During this thesis work a coupled thermo-mechanical finite element model (FEM) was builtto simulate hot rolling in the blooming mill at Sandvik Materials Technology (SMT) inSandviken. The blooming mill is the first in a long line of processes that continuously or ingotcast ingots are subjected to before becoming finished products. The aim of this thesis work was twofold. The first was to create a parameterized finiteelement (FE) model of the blooming mill. The commercial FE software package MSCMarc/Mentat was used to create this model and the programing language Python was used toparameterize it. Second, two different pass schedules (A and B) were studied and comparedusing the model. The two pass series were evaluated with focus on their ability to healcentreline porosity, i.e. to close voids in the centre of the ingot. This evaluation was made by studying the hydrostatic stress (σm), the von Mises stress (σeq)and the plastic strain (εp) in the centre of the ingot. From these parameters the stress triaxiality(Tx) and the hydrostatic integration parameter (Gm) were calculated for each pass in bothseries using two different transportation times (30 and 150 s) from the furnace. The relationbetween Gm and an analytical parameter (Δ) was also studied. This parameter is the ratiobetween the mean height of the ingot and the contact length between the rolls and the ingot,which is useful as a rule of thumb to determine the homogeneity or penetration of strain for aspecific pass. The pass series designed with fewer passes (B), many with greater reduction, was shown toachieve better void closure theoretically. It was also shown that a temperature gradient, whichis the result of a longer holding time between the furnace and the blooming mill leads toimproved void closure.