Continuation-conjugate gradient methods for the least squares solution of nonlinear boundary value problems

We discuss in this paper a new combination of methods for solving nonlinear boundary value problems containing a parameter. Methods of the continuation type are combined with least squares formulations, preconditioned conjugate gradient algorithms and finite element approximations. We can compute branches of solutions with limit points, bifurcation points, etc. Several numerical tests illustrate the possibilities of the methods discussed in the present paper; these include the Bratu problem in one and two dimensions, one-dimensional bifurcation and perturbed bifurcation problems, the driven cavity problem for the Navier–Stokes equations

Parametric Manifolds I: Extrinsic Approach

A parametric manifold can be viewed as the manifold of orbits of a (regular) foliation of a manifold by means of a family of curves. If the foliation is hypersurface orthogonal, the parametric manifold is equivalent to the 1-parameter family of hypersurfaces orthogonal to the curves, each of which inherits a metric and connection from the original manifold via orthogonal projections; this is the well-known Gauss-Codazzi formalism. We generalize this formalism to the case where the foliation is not hypersurface orthogonal. Crucial to this generalization is the notion of deficiency, which measures the failure of the orthogonal tangent spaces to be surface-forming, and which behaves very much like torsion. Some applications to initial value problems in general relativity will be briefly discussed.Comment: Plain TeX, 21 pages, no figure

Simulasi Penanganan Potensi Aliran Debris di Gunung Sago (Studi Kasus di Batang Lakin, Kecamatan Lareh Sago Halaban, Kabupaten Lima Puluh Kota)

The regions in foothills of Sago mountain are flood-prone area due to debris flow. As occurred on March 22, 2010, there has been a catastrophic overflow of debris flow from Sago mountain. The disaster resulted in severe damage around the rivers downstream Sago mountains, including Batang Lakin river. This research study debris flow potential and how to mitigate it in Batang Lakin river, West Sumatra. Analysis of potential debris flow hazard of Batang Lakin river and alternative debris mitigation is simulated using the debris flow simulator Kanako 2D version 2.051. Simulation is important for verifying effect of controlling flow of debris prior to construction work carried out. Rain data input was calculated based on fifty years time period and one hundred years time period Research findings show that at Batang Lakin river, debris flow occurred and overflowing river channel. Alternative countermeasure chosen is sabo dam. For fifty years period when debris flow peak discharge of 59.50 m3/second required 2 units of sabo dams (closed type) with positions at Sta 0 +200 (Sabo height 6 m) and at Sta 0 +450 (Sabo height 4 m). For one hundred years period when debris flow peak discharge of 62.66 m3/second required 2 units of sabo dams (closed type) with positions at Sta 0 +200 (Sabo height 6 m) and at Sta 0 +450 (Sabo height 5 m) to prevent overflow of debris flow to the settlement. Thus, the right efforts to control debris flow on Batang Lakin is the sabo dam

On the classification of type D spacetimes

We give a classification of the type D spacetimes based on the invariant differential properties of the Weyl principal structure. Our classification is established using tensorial invariants of the Weyl tensor and, consequently, besides its intrinsic nature, it is valid for the whole set of the type D metrics and it applies on both, vacuum and non-vacuum solutions. We consider the Cotton-zero type D metrics and we study the classes that are compatible with this condition. The subfamily of spacetimes with constant argument of the Weyl eigenvalue is analyzed in more detail by offering a canonical expression for the metric tensor and by giving a generalization of some results about the non-existence of purely magnetic solutions. The usefulness of these results is illustrated in characterizing and classifying a family of Einstein-Maxwell solutions. Our approach permits us to give intrinsic and explicit conditions that label every metric, obtaining in this way an operational algorithm to detect them. In particular a characterization of the Reissner-Nordstr\"{o}m metric is accomplished.Comment: 29 pages, 0 figure

Long term stable integration of a maximally sliced Schwarzschild black hole using a smooth lattice method

We will present results of a numerical integration of a maximally sliced Schwarzschild black hole using a smooth lattice method. The results show no signs of any instability forming during the evolutions to t=1000m. The principle features of our method are i) the use of a lattice to record the geometry, ii) the use of local Riemann normal coordinates to apply the 1+1 ADM equations to the lattice and iii) the use of the Bianchi identities to assist in the computation of the curvatures. No other special techniques are used. The evolution is unconstrained and the ADM equations are used in their standard form.Comment: 47 pages including 26 figures, plain TeX, also available at http://www.maths.monash.edu.au/~leo/preprint

The Topology of Branching Universes

The purpose of this paper is to survey the possible topologies of branching space-times, and, in particular, to refute the popular notion in the literature that a branching space-time requires a non-Hausdorff topology

When Micro Prudence Increases Macro Risk: The Destabilizing Effects of Financial Innovation, Leverage, and Diversification

By exploiting basic common practice accounting and risk-management rules, we propose a simple analytical dynamical model to investigate the effects of microprudential changes on macroprudential outcomes. Specifically, we study the consequence of the introduction of a financial innovation that allows reducing the cost of portfolio diversification in a financial system populated by financial institutions having capital requirements in the form of Value at Risk (VaR) constraint and following standard mark-to-market and risk-management rules. We provide a full analytical quantification of the multivariate feedback effects between investment prices and bank behavior induced by portfolio rebalancing in presence of asset illiquidity and show how changes in the constraints of the bank portfolio optimization endogenously drive the dynamics of the balance sheet aggregate of financial institutions and, thereby, the availability of bank liquidity to the economic system and systemic risk. The model shows that when financial innovation reduces the cost of diversification below a given threshold, the strength (because of higher leverage) and coordination (because of similarity of bank portfolios) of feedback effects increase, triggering a transition from a stationary dynamics of price returns to a nonstationary one characterized by steep growths (bubbles) and plunges (bursts) of market prices

Cohomological tautness for Riemannian foliations

In this paper we present some new results on the tautness of Riemannian foliations in their historical context. The first part of the paper gives a short history of the problem. For a closed manifold, the tautness of a Riemannian foliation can be characterized cohomologically. We extend this cohomological characterization to a class of foliations which includes the foliated strata of any singular Riemannian foliation of a closed manifold

Late time behaviour of the maximal slicing of the Schwarzschild black hole

A time-symmetric Cauchy slice of the extended Schwarzschild spacetime can be evolved into a foliation of the $r>3m/2$-region of the spacetime by maximal surfaces with the requirement that time runs equally fast at both spatial ends of the manifold. This paper studies the behaviour of these slices in the limit as proper time-at-infinity becomes arbitrarily large and gives an analytic expression for the collapse of the lapse.Comment: 18 pages, Latex, no figure