A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

GENERICITY OF NONDEGENERATE GEODESICS WITH GENERAL BOUNDARY CONDITIONS

Let M be a possibly noncompact manifold. We prove, generically in the C(k)-topology (2 <= k <= infinity), that semi-Riemannian metrics of a given index on M do not possess any degenerate geodesics satisfying suitable boundary conditions. This extends a result of L. Biliotti, M. A. Javaloyes and P. Piccione [6] for geodesics with fixed endpoints to the case where endpoints lie on a compact submanifold P subset of M x M that satisfies an admissibility condition. Such condition holds, for example, when P is transversal to the diagonal Delta subset of M x M. Further aspects of these boundary conditions are discussed and general conditions under which metrics without degenerate geodesics are C(k)-generic are given.

The impact of different glacial boundary conditions on atmospheric dynamics and precipitation in the North Atlantic region

Using a highly resolved atmospheric general circulation model, the impact of different glacial boundary conditions on precipitation and atmospheric dynamics in the North Atlantic region is investigated. Six 30-yr time slice experiments of the Last Glacial Maximum at 21 thousand years before the present (ka BP) and of a less pronounced glacial state – the Middle Weichselian (65 ka BP) – are compared to analyse the sensitivity to changes in the ice sheet distribution, in the radiative forcing and in the prescribed time-varying sea surface temperature and sea ice, which are taken from a lower-resolved, but fully coupled atmosphere-ocean general circulation model. The strongest differences are found for simulations with different heights of the Laurentide ice sheet. A high surface elevation of the Laurentide ice sheet leads to a southward displacement of the jet stream and the storm track in the North Atlantic region. These changes in the atmospheric dynamics generate a band of increased precipitation in the mid-latitudes across the Atlantic to southern Europe in winter, while the precipitation pattern in summer is only marginally affected. The impact of the radiative forcing differences between the two glacial periods and of the prescribed time-varying sea surface temperatures and sea ice are of second order importance compared to the one of the Laurentide ice sheet. They affect the atmospheric dynamics and precipitation in a similar but less pronounced manner compared with the topographic changes.

Effect of boundary conditions on yield properties of human femoral trabecular bone

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7–72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient R2≥0.93. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Entropy generation for microscale forced convection: effects of different thermal boundary conditions, velocity slip, temperature jump, viscous dissipation, and duct geometry

This work presents closed form solutions for fully developed temperature distribution and entropy generation due to forced convection in microelectromechanical systems (MEMS) in the Slip-flow regime, for which the Knudsen number lies within the range 0.001

Constant wall heat flux boundary conditions in micro-channels filled with porous material with internal heat generation under local thermal non-equilibrium conditions

Forced convection heat transfer in a micro-channel filled with a porous material saturated with rarefied gas with internal heat generation is studied analytically in this work. The study is performed by analysing the boundary conditions for constant wall heat flux under local thermal non-equilibrium (LTNE) conditions. Invoking the velocity slip and temperature jump, the thermal behaviour of the porous-fluid system is studied by considering thermally and hydrodynamically fully-developed conditions. The flow inside the porous material is modelled by the Darcy–Brinkman equation. Exact solutions are obtained for both the fluid and solid temperature distributions for two primary approaches models A and B using constant wall heat flux boundary conditions. The temperature distributions and Nusselt numbers for models A and B are compared, and the limiting cases resulting in the convergence or divergence of the two models are also discussed. The effects of pertinent parameters such as fluid to solid effective thermal conductivity ratio, Biot number, Darcy number, velocity slip and temperature jump coefficients, and fluid and solid internal heat generations are also discussed. The results indicate that the Nusselt number decreases with the increase of thermal conductivity ratio for both models. This contrasts results from previous studies which for model A reported that the Nusselt number increases with the increase of thermal conductivity ratio. The Biot number and thermal conductivity ratio are found to have substantial effects on the role of temperature jump coefficient in controlling the Nusselt number for models A and B. The Nusselt numbers calculated using model A change drastically with the variation of solid internal heat generation. In contrast, the Nusselt numbers obtained for model B show a weak dependency on the variation of internal heat generation. The velocity slip coefficient has no noticeable effect on the Nusselt numbers for both models. The difference between the Nusselt numbers calculated using the two models decreases with an increase of the temperature jump coefficient.

Steps towards fully nonlinear simulations of arrays of OWSC

Abandonment of farming systems on upland areas in southwest Britain during the Late Bronze Age – some 3000 years ago – is widely considered a ‘classic’ demonstration of the impact of deteriorating climate on the vulnerability of populations in such marginal environments. Here we test the hypothesis that climate change drove the abandonment of upland areas by developing new chronologies for human activity on upland areas during the Bronze Age across southwest Britain (Dartmoor, Exmoor and Bodmin Moor). We find Bronze Age activity in these areas spanned 3900–2950 calendar years ago with abandonment by 2900 calendar years ago. Holocene Irish bog and lake oak tree populations provide evidence of major shifts in hydroclimate across western Britain and Ireland, coincident with ice rafted debris layers recognized in North Atlantic marine sediments, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. We observe abandonment of upland areas in southwest Britain coinciding with a sustained period of extreme wet conditions that commenced 3100 calendar years ago. Our results are consistent with the view that climate change increased the vulnerability of these early farming communities and led to a less intensive use of such marginal environments across Britain.

Dynamical Casimir effect for a massless scalar field between two concentric spherical shells with mixed boundary conditions

In this work we analyze the dynamical Casimir effect for a massless scalar field confined between two concentric spherical shells considering mixed boundary conditions. We thus generalize a previous result in literature [Phys. Rev. A 78, 032521 (2008)], where the same problem is approached for the field constrained to the Dirichlet-Dirichlet boundary conditions. A general expression for the average number of particle creation is deduced considering an arbitrary law of radial motion of the spherical shells. This expression is then applied to harmonic oscillations of the shells, and the number of particle production is analyzed and compared with the results previously obtained under Dirichlet-Dirichlet boundary conditions.

Effects of viscous dissipation and boundary conditions on forced convection in a channel occupied by a saturated porous medium

Forced convection with viscous dissipation in a parallel plate channel filled by a saturated porous medium is investigated numerically. Three different viscous dissipation models are examined. Two different sets of wall conditions are considered: isothermal and isoflux. Analytical expressions are also presented for the asymptotic temperature profile and the asymptotic Nusselt number. With isothermal walls, the Brinkman number significantly influences the developing Nusselt number but not the asymptotic one. At constant wall heat flux, both the developing and the asymptotic Nusselt numbers are affected by the value of the Brinkman number. The Nusselt number is sensitive to the porous medium shape factor under all conditions considered.

Effects of viscous dissipation and boundary conditions on forced convection in a channel occupied by a saturated porous medium

Forced convection with viscous dissipation in a parallel plate channel filled by a saturated porous medium is investigated numerically. Three different viscous dissipation models are examined. Two different sets of wall conditions are considered: isothermal and isoflux. Analytical expressions are also presented for the asymptotic temperature profile and the asymptotic Nusselt number. With isothermal walls, the Brinkman number significantly influences the developing Nusselt number but not the asymptotic one. At constant wall heat flux, both the developing and the asymptotic Nusselt numbers are affected by the value of the Brinkman number. The Nusselt number is sensitive to the porous medium shape factor under all conditions considered.

An improved novel method for solving nonlinear boundary value problems

A modified formula for the integral transform of a nonlinear function is proposed for a class of nonlinear boundary value problems. The technique presented in this paper results in analytical solutions. Iterations and initial guess, which are needed in other techniques, are not required in this novel technique. The analytical solutions are found to agree surprisingly well with the numerically exact solutions for two examples of power law reaction and Langmuir-Hinshelwood reaction in a catalyst pellet.

Integrable open-boundary conditions for the q-deformed supersymmetric U model of strongly correlated electrons

A general graded reflection equation algebra is proposed and the corresponding boundary quantum inverse scattering method is formulated. The formalism is applicable to all boundary lattice systems where an invertible R-matrix exists. As an application, the integrable open-boundary conditions for the q-deformed supersymmetric U model of strongly correlated electrons are investigated. The diagonal boundary K-matrices are found and a class of integrable boundary terms are determined. The boundary system is solved by means of the coordinate space Bethe ansatz technique and the Bethe ansatz equations are derived. As a sideline, it is shown that all R-matrices associated with a quantum affine superalgebra enjoy the crossing-unitarity property. (C) 1998 Elsevier Science B.V.

New integrable boundary conditions for the q-deformed supersymmetric U model and Bethe ansatz equations

New classes of integrable boundary conditions for the q-deformed (or two-parameter) supersymmetric U model are presented. The boundary systems are solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.

Nine classes of integrable boundary conditions for the eight-state supersymmetric fermion model

Nine classes of integrable boundary conditions for the eight-state supersymmetric model of strongly correlated fermions are presented. The boundary systems are solved by using the coordinate Bethe ansatz method and the Bethe ansatz equations for all nine cases are given.