Numerical solutions of elliptic inverse problems via the equation error method

To estimate a parameter in an elliptic boundary value problem, the method of equation error chooses the value that minimizes the error in the PDE and boundary condition (the solution of the BVP having been replaced by a measurement). The estimated parameter converges to the exact value as the measured data converge to the exact value, provided Tikhonov regularization is used to control the instability inherent in the problem. The error in the estimated solution can be bounded in an appropriate quotient norm; estimates can be derived for both the underlying (infinite-dimensional) problem and a finite-element discretization that can be implemented in a practical algorithm. Numerical experiments demonstrate the efficacy and limitations of the method.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Finite element study on fracture patterns in human scaphoid wrist bone due to free fall

Scaphoid is one of the 8 carpal bones found adjacent to the thumb supported proximally by Radius bone. During the free fall, on outstretched hand, the impact load gets transferred to the scaphoid at its free anterior end. Unique arrangement of other carpal bones in the palm is also one of the reasons for the load to get transferred to scaphoid. About half of the total load acting upon carpal bone gets transferred to scaphoid at its distal pole. There are about 10 to 12 clinically observed fracture pattern in the scaphoid due to free fall. The aim of the study is to determine the orientation of the load, magnitude of the load and the corresponding fracture pattern. This study includes both static and dynamic finite element models validated by experiments. The scaphoid model has been prepared from CT scans of a 27 year old person. The 2D slices of the CT scans have been converted to 3D model by using MIMICS software. There are four cases of loading studied which are considered to occur clinically more frequently. In case (i) the load is applied at the posterior end at distal pole whereas in case (ii), (iii) and (iv), the load is applied at anterior end at different directions. The model is given a fixed boundary condition at the region which is supported by Radius bone during the impact. Same loading and boundary conditions have been used in both static and dynamic explicit finite element analysis. The site of fracture initiation and path of fracture propagation have been identified by using max principal stress / gradient and max principal strain / gradient criterion respectively in static and dynamic explicit finite element analysis. Static and dynamic impact experiments were performed on the polyurethane foam specimens to validate the finite element results. Experimental results such as load at fracture, site of fracture initiation and path of fracture propagation have been compared with the results of finite element analysis. Four different types of fracture patterns observed in clinical studies have been identified in this study.

Experimental characterization of cement–bentonite interaction using core infiltration techniques and 4D computed tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel drums, and as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predict significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this project is to characterize and quantify the cement/bentonite skin effects spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used, which allows performing X-ray computed tomography (CT) periodically without interrupting a running experiment. A pre-saturated cylindrical MX-80 bentonite sample (1920 kg/m3 average wet density) is subjected to a confining pressure as a constant total pressure boundary condition. The infiltration of a hyperalkaline (pH 13.4), artificial OPC (ordinary Portland cement) pore water into the bentonite plug alters the mineral assemblage over time as an advancing reaction front. The related changes in X-ray attenuation values are related to changes in phase densities, porosity and local bulk density and are tracked over time periodically by non-destructive CT scans.

Three-dimensional Vortex Structures on Heated Micro-lattice in a Gas

This paper addresses the microscale heat transfer problem from heated lattice to the gas. A micro-device for enhanced heat transfer is presented and numerically investigated. Thermal creep induces 3-D vortex structures in the vicinity of the lattice. The gas flow is in the slip flow regime (Knudsen number Kn⩽0.1Kn⩽0.1). The simulations are performed using slip flow Navier–Stokes equations with boundary condition formulations proposed by Maxwell and Smoluchowski. In this study the wire thicknesses and distances of the heated lattice are varied. The surface geometrical properties alter significantly heat flux through the surface.

Next-to-Minimal SOFTSUSY

We describe an extension to the SOFTSUSY program that provides for the calculation of the sparticle spectrum in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), where a chiral superfield that is a singlet of the Standard Model gauge group is added to the Minimal Supersymmetric Standard Model (MSSM) fields. Often, a Z3 symmetry is imposed upon the model. SOFTSUSY can calculate the spectrum in this case as well as the case where general Z3 violating (denoted as ) terms are added to the soft supersymmetry breaking terms and the superpotential. The user provides a theoretical boundary condition for the couplings and mass terms of the singlet. Radiative electroweak symmetry breaking data along with electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. This paper serves as a manual to the NMSSM mode of the program, detailing the approximations and conventions used.

Fluid-structure interaction of an elastic pipe immersed inside a coaxial rigid pipe: a model for the Tullio phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo).

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.

Fate of accidental symmetries of the relativistic hydrogen atom in a spherical cavity

The non-relativistic hydrogen atom enjoys an accidental SO(4) symmetry, that enlarges the rotational SO(3) symmetry, by extending the angular momentum algebra with the Runge–Lenz vector. In the relativistic hydrogen atom the accidental symmetry is partially lifted. Due to the Johnson–Lippmann operator, which commutes with the Dirac Hamiltonian, some degeneracy remains. When the non-relativistic hydrogen atom is put in a spherical cavity of radius R with perfectly reflecting Robin boundary conditions, characterized by a self-adjoint extension parameter γ, in general the accidental SO(4) symmetry is lifted. However, for R=(l+1)(l+2)a (where a is the Bohr radius and l is the orbital angular momentum) some degeneracy remains when γ=∞ or γ = 2/R. In the relativistic case, we consider the most general spherically and parity invariant boundary condition, which is characterized by a self-adjoint extension parameter. In this case, the remnant accidental symmetry is always lifted in a finite volume. We also investigate the accidental symmetry in the context of the Pauli equation, which sheds light on the proper non-relativistic treatment including spin. In that case, again some degeneracy remains for specific values of R and γ.

Pressure from the vacuum of confined spinor matter

Charged spinor matter field is quantized in a spatial region bounded by two parallel neutral plates. The most general set of boundary conditions ensuring the confinement of matter within the plates is considered. We study a response of the vacuum of the confined matter to the background uniform magnetic field which is directed orthogonally to the plates. It is proven that, in the case of a sufficiently strong magnetic field, the vacuum pressure onto the plates is positive and independent of the boundary condition, as well as of the distance between the plates.

Influence of quantized massive matter fields on the Casimir effect

Charged massive matter fields of spin-0 and spin- 1/2 are quantized in the presence of an external uniform magnetic field in a spatial region bounded by two parallel plates. The most general set of boundary conditions at the plates, that is required by mathematical consistency and the self-adjointness of the Hamiltonian operator, is employed. The vacuum fluctuations of the matter field in the case of the magnetic field orthogonal to the plates are analyzed, and it is shown that the pressure from the vacuum onto the plates is positive and independent of the boundary condition, as well as of the distance between the plates. Possibilities of the detection of this new-type Casimir effect are discussed. Read More: http://www.worldscientific.com/doi/10.1142/S0217732315500996

Casimir effect with quantized charged spinor matter in background magnetic field

We study the influence of a background uniform magnetic field and boundary conditions on the vacuum of a quantized charged spinor matter field confined between two parallel neutral plates; the magnetic field is directed orthogonally to the plates. The admissible set of boundary conditions at the plates is determined by the requirement that the Dirac Hamiltonian operator be self-adjoint. It is shown that, in the case of a sufficiently strong magnetic field and a sufficiently large separation of the plates, the generalized Casimir force is repulsive, being independent of the choice of a boundary condition, as well as of the distance between the plates. The detection of this effect seems to be feasible in the foreseeable future.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Sediment geochemical evidence for an early-middle Gilbert (early Pliocene) productivity peak in the North Pacific Red Clay Province

Current attempts to understand climatic variability during the early to middle Pliocene require paleoceanographic information from the Pacific and Indian Oceans that may serve to test and/or constrain future circulation models. Ocean Drilling Program (ODP) Sites 885/886 are located in the central subarctic North Pacific at water depths exceeding 5700 m. Recent studies of rock magnetic properties suggest that the fine-grained Fe oxide component in sediment at Sites 885/886 experienced reductive dissolution during the early-middle Gilbert. Because such an interval in the North Pacific Red Clay Province suggests a maximum in the sedimentary flux of organic carbon and/or a minimum in bottom water dissolved O2 concentrations (and hence, a peak change in North Pacific oceanographic conditions), a geochemical investigation was conducted to test the hypothesis. Quaternary sediment at Hole 886B was subjected to an oxyhydroxide removal procedure, and chemical analyses indicate that bulk sediment concentrations of Fe and the Fe/Sc ratio decrease significantly upon reductive dissolution. Downcore chemical analyses of untreated sediment at Hole 886B demonstrate that similar depletions also occur across the proposed interval of reduced sediment. Downcore chemical analyses also indicate that a pronounced increase in the Ba/Sc ratio occurs across the interval. These results are consistent with an interpretation that abyssal sediment of the North Pacific experienced a decrease in redox conditions during the early-middle Gilbert, and that this change in oxidation state was related to a peak in paleoproductivity. If the zenith of late Miocene to middle Pliocene enhanced productivity observed at other Indo-Pacific divergence regions similarly can be constrained to the early-middle Gilbert, there exists an oceanographic boundary condition in which to test future models concerning Pliocene warmth.

Mass balance measurements on King George Island ice cap in 2007 and 2008

The Antarctic Peninsula has been identified as a region of rapid on-going climate change with impacts on the cryosphere. The knowledge of glacial changes and freshwater budgets resulting from intensified glacier melt is an important boundary condition for many biological and integrated earth system science approaches. We provide a case study on glacier and mass balance changes for the ice cap of King George Island. The area loss between 2000 and 2008 amounted to about 20 km**2 (about 1.6% of the island area) and compares to glacier retreat rates observed in previous years. Measured net accumulation rates for two years (2007 and 2008) show a strong interannual variability with maximum net accumulation rates of 4950 mm w.e./a and 3184 mm w.e./a, respectively. These net accumulation rates are at least 4 times higher than reported mean values (1926-95) from an ice core. An elevation dependent precipitation rate of 343 mm w.e./a (2007) and 432 mm w.e./a (2008) per 100 m elevation increase was observed. Despite these rather high net accumulation rates on the main ice cap, consistent surface lowering was observed at elevations below 270 m above ellipsoid over an 11-year period. These DGPS records reveal a linear dependence of surface lowering with altitude with a maximum annual surface lowering rate of 1.44 m/a at 40 m and -0.20 m/a at 270 m above ellipsoid. These results fit well to observations by other authors and surface lowering rates derived from the ICESat laser altimeter. Assuming that climate conditions of the past 11 years continue, the small ice cap of Bellingshausen Dome will disappear in about 285 years.