Goethite on Mars - A laboratory study of physically and chemically bound water in ferric oxides

Thermogravimetric study of physically and chemically bound water in ferric oxides of limonite with application to goethite on Mar

Eigen-analysis of Inviscid Fluid Structure Interaction (FSI) Systems with Complex Boundary Conditions

A method for extracting the eigenvalues and eigenmodes from complex coupled fluid-structure interaction (FSI) systems is presented. The FSI system under consideration in this case is a one-sided, inviscid flow over a finite-length compliant surface with complex boundary conditions, although the method could be applied to any FSI system. The flow is solved for the inviscid case using a boundary-element method solution of Laplace’s equation, while the finite compliant surface is solved through a finite-difference solution of the one-dimensional beam equation. The crux of the method lies in reducing the coupled fluid and structural equations down to a set of coupled linear differential equations. Standard Krylov subspace projection methods may then be used to determine the eigenvalues of the large system of linear equations. This method is applied to the analysis of hydroelastic FSI systems with complex boundary conditions that would be difficult or otherwise impossible to analyse using standard Galerkin methods. Specifically, the complex cases of inhomogeneous and discontinuous compliant wall properties and arbitrary hinge-joint conditions along the compliant surface are considered

Eigen-analysis of a Fully Viscous Boundary-Layer flow Interacting with a Finite Compliant Surface

A method and preliminary results are presented for the determination of eigenvalues and eigenmodes from fully viscous boundary layer flow interacting with a finite length one-sided compliant wall. This is an extension to the analysis of inviscid flow-structure systems which has been established in previous work. A combination of spectral and finite-difference methods are applied to a linear perturbation form of the full Navier-Stokes equations and one-dimensional beam equation. This yields a system of coupled linear equations that accurately define the spatio-temporal development of linear perturbations to a boundary layer flow over a finite-length compliant surface. Standard Krylov subspace projection methods are used to extract the eigenvalues from this complex system of equations. To date, the analysis of the development of Tollmien-Schlichting (TS) instabilities over a finite compliant surface have relied upon DNS-type results across a narrow (or even singular) spectrum of TS waves. The results from this method have the potential to describe conclusively the role that a finite length compliant surface has in the development of two-dimensional TS instabilities and other FSI instabilities across a broad spectrum

The development and technology transfer of software engineering technology at NASA. Johnson Space Center

The United State's big space projects of the next decades, such as Space Station and the Human Exploration Initiative, will need the development of many millions of lines of mission critical software. NASA-Johnson (JSC) is identifying and developing some of the Computer Aided Software Engineering (CASE) technology that NASA will need to build these future software systems. The goal is to improve the quality and the productivity of large software development projects. New trends are outlined in CASE technology and how the Software Technology Branch (STB) at JSC is endeavoring to provide some of these CASE solutions for NASA is described. Key software technology components include knowledge-based systems, software reusability, user interface technology, reengineering environments, management systems for the software development process, software cost models, repository technology, and open, integrated CASE environment frameworks. The paper presents the status and long-term expectations for CASE products. The STB's Reengineering Application Project (REAP), Advanced Software Development Workstation (ASDW) project, and software development cost model (COSTMODL) project are then discussed. Some of the general difficulties of technology transfer are introduced, and a process developed by STB for CASE technology insertion is described

Hydraulic flow through a channel contraction: multiple steady states

We have investigated shallow water flows through a channel with a contraction by experimental and theoretical means. The horizontal channel consists of a sluice gate and an upstream channel of constant width $b_0$ ending in a linear contraction of minimum width $b_c$. Experimentally, we observe upstream steady and moving bores/shocks, and oblique waves in the contraction, as single and multiple steady states, as well as a steady reservoir with a complex hydraulic jump in the contraction occurring in a small section of the $b_c/b_0$ and Froude number parameter plane. One-dimensional hydraulic theory provides a comprehensive leading-order approximation, in which a turbulent frictional parametrization is used to achieve quantitative agreement. An analytical and numerical analysis is given for two-dimensional supercritical shallow water flows. It shows that the one-dimensional hydraulic analysis for inviscid flows away from hydraulic jumps holds surprisingly well, even though the two-dimensional oblique hydraulic jump patterns can show large variations across the contraction channel

Bessel bridges decomposition with varying dimension. Applications to finance

We consider a class of stochastic processes containing the classical and well-studied class of Squared Bessel processes. Our model, however, allows the dimension be a function of the time. We first give some classical results in a larger context where a time-varying drift term can be added. Then in the non-drifted case we extend many results already proven in the case of classical Bessel processes to our context. Our deepest result is a decomposition of the Bridge process associated to this generalized squared Bessel process, much similar to the much celebrated result of J. Pitman and M. Yor. On a more practical point of view, we give a methodology to compute the Laplace transform of additive functionals of our process and the associated bridge. This permits in particular to get directly access to the joint distribution of the value at t of the process and its integral. We finally give some financial applications to illustrate the panel of applications of our results

Solution of the Fokker-Planck equation with a logarithmic potential and mixed eigenvalue spectrum

Motivated by a problem in climate dynamics, we investigate the solution of a Bessel-like process with negative constant drift, described by a Fokker-Planck equation with a potential V(x) = - [b \ln(x) + a\, x], for b>0 and a<0. The problem belongs to a family of Fokker-Planck equations with logarithmic potentials closely related to the Bessel process, that has been extensively studied for its applications in physics, biology and finance. The Bessel-like process we consider can be solved by seeking solutions through an expansion into a complete set of eigenfunctions. The associated imaginary-time Schroedinger equation exhibits a mix of discrete and continuous eigenvalue spectra, corresponding to the quantum Coulomb potential describing the bound states of the hydrogen atom. We present a technique to evaluate the normalization factor of the continuous spectrum of eigenfunctions that relies solely upon their asymptotic behavior. We demonstrate the technique by solving the Brownian motion problem and the Bessel process both with a negative constant drift. We conclude with a comparison with other analytical methods and with numerical solutions.Comment: 21 pages, 8 figure

How long does it take to pull an ideal polymer into a small hole?

We present scaling estimates for characteristic times $\tau_{\rm lin}$ and $\tau_{\rm br}$ of pulling ideal linear and randomly branched polymers of $N$ monomers into a small hole by a force $f$. We show that the absorbtion process develops as sequential straightening of folds of the initial polymer configuration. By estimating the typical size of the fold involved into the motion, we arrive at the following predictions: $\tau_{\rm lin}(N) \sim N^{3/2}/f$ and $\tau_{\rm br}(N) \sim N^{5/4}/f$, and we also confirm them by the molecular dynamics experiment.Comment: 4 pages, 3 figure

Windings of the 2D free Rouse chain

We study long time dynamical properties of a chain of harmonically bound Brownian particles. This chain is allowed to wander everywhere in the plane. We show that the scaling variables for the occupation times T_j, areas A_j and winding angles \theta_j (j=1,...,n labels the particles) take the same general form as in the usual Brownian motion. We also compute the asymptotic joint laws P({T_j}), P({A_j}), P({\theta_j}) and discuss the correlations occuring in those distributions.Comment: Latex, 17 pages, submitted to J. Phys.

Levy-Student Distributions for Halos in Accelerator Beams

We describe the transverse beam distribution in particle accelerators within the controlled, stochastic dynamical scheme of the Stochastic Mechanics (SM) which produces time reversal invariant diffusion processes. This leads to a linearized theory summarized in a Shchr\"odinger--like (\Sl) equation. The space charge effects have been introduced in a recent paper~\cite{prstab} by coupling this \Sl equation with the Maxwell equations. We analyze the space charge effects to understand how the dynamics produces the actual beam distributions, and in particular we show how the stationary, self--consistent solutions are related to the (external, and space--charge) potentials both when we suppose that the external field is harmonic (\emph{constant focusing}), and when we \emph{a priori} prescribe the shape of the stationary solution. We then proceed to discuss a few new ideas~\cite{epac04} by introducing the generalized Student distributions, namely non--Gaussian, L\'evy \emph{infinitely divisible} (but not \emph{stable}) distributions. We will discuss this idea from two different standpoints: (a) first by supposing that the stationary distribution of our (Wiener powered) SM model is a Student distribution; (b) by supposing that our model is based on a (non--Gaussian) L\'evy process whose increments are Student distributed. We show that in the case (a) the longer tails of the power decay of the Student laws, and in the case (b) the discontinuities of the L\'evy--Student process can well account for the rare escape of particles from the beam core, and hence for the formation of a halo in intense beams.Comment: revtex4, 18 pages, 12 figure