Gravitational actions upon a tether in a non-uniform gravity field with arbitrary number of zonal harmonics

We develop general closed-form expressions for the mutual gravitational potential, resultant and torque acting upon a rigid tethered system moving in a non-uniform gravity field produced by an attracting body with revolution symmetry, such that an arbitrary number of zonal harmonics is considered. The final expressions are series expansion in two small parameters related to the reference radius of the primary and the length of the tether, respectively, each of which are scaled by the mutual distance between their centers of mass. A few numerical experiments are performed to study the convergence behavior of the final expressions, and conclude that for high precision applications it might be necessary to take into account additional perturbation terms, which come from the mutual Two-Body interaction.

Some aspects of lubrication in heavy regimes: thermal effects, stability and turbulence

In this work, a combination of numerical methods applied to thermohydrodynamic lubrication problems with cavitation is presented. It should be emphasized the difficulty of the nonlinear mathematical coupled model involving a free boundary problem, but also the simplicity of the algorithms employed to solve it. So, finite element discretizations for the hydrodynamic and thermal equations combined with upwind techniques for the convection terms and duality methods for nonlinear features are proposed. Additionally, a model describing the movement of the shaft is provided. Considering the shaft as a rigid body this model will consist of an ODE system relating acceleration of the center of gravity and external and pressure loads. The numerical experiments of mechanical stability try to clarify the position of the neutral stability curve. Finally, a rotating machine for ship propulsion involving both axial and radial bearings operating with nonconventional lubricants (seawater to avoid environmental pollution) is analyzed by using laminar and turbulent inertial flows.

A subgrid viscosity Lagrance-Galerkin method for convection-diffusion problems

We present and analyze a subgrid viscosity Lagrange-Galerk in method that combines the subgrid eddy viscosity method proposed in W. Layton, A connection between subgrid scale eddy viscosity and mixed methods. Appl. Math. Comp., 133: 14 7-157, 2002, and a conventional Lagrange-Galerkin method in the framework of P1⊕ cubic bubble finite elements. This results in an efficient and easy to implement stabilized method for convection dominated convection diffusion reaction problems. Numerical experiments support the numerical analysis results and show that the new method is more accurate than the conventional Lagrange-Galerkin one.

Spectral bifurcations in dispersive wave turbulence

Dispersive wave turbulence is studied numerically for a class of one-dimensional nonlinear wave equations. Both deterministic and random (white noise in time) forcings are studied. Four distinct stable spectra are observed—the direct and inverse cascades of weak turbulence (WT) theory, thermal equilibrium, and a fourth spectrum (MMT; Majda, McLaughlin, Tabak). Each spectrum can describe long-time behavior, and each can be only metastable (with quite diverse lifetimes)—depending on details of nonlinearity, forcing, and dissipation. Cases of a long-live MMT transient state dcaying to a state with WT spectra, and vice-versa, are displayed. In the case of freely decaying turbulence, without forcing, both cascades of weak turbulence are observed. These WT states constitute the clearest and most striking numerical observations of WT spectra to date—over four decades of energy, and three decades of spatial, scales. Numerical experiments that study details of the composition, coexistence, and transition between spectra are then discussed, including: (i) for deterministic forcing, sharp distinctions between focusing and defocusing nonlinearities, including the role of long wavelength instabilities, localized coherent structures, and chaotic behavior; (ii) the role of energy growth in time to monitor the selection of MMT or WT spectra; (iii) a second manifestation of the MMT spectrum as it describes a self-similar evolution of the wave, without temporal averaging; (iv) coherent structures and the evolution of the direct and inverse cascades; and (v) nonlocality (in k-space) in the transferral process.

Self-similar intermediate asymptotics for a degenerate parabolic filtration-absorption equation

The equation ∂tu = u∂xx2u − (c − 1)(∂xu)2 is known in literature as a qualitative mathematical model of some biological phenomena. Here this equation is derived as a model of the groundwater flow in a water-absorbing fissurized porous rock; therefore, we refer to this equation as a filtration-absorption equation. A family of self-similar solutions to this equation is constructed. Numerical investigation of the evolution of non-self-similar solutions to the Cauchy problems having compactly supported initial conditions is performed. Numerical experiments indicate that the self-similar solutions obtained represent intermediate asymptotics of a wider class of solutions when the influence of details of the initial conditions disappears but the solution is still far from the ultimate state: identical zero. An open problem caused by the nonuniqueness of the solution of the Cauchy problem is discussed.

Stability and sensitivity analysis of Be-CoDiS, an epidemiological model to predict the spread of human diseases between countries. Validation with data from the 2014-16 West African Ebola Virus Disease epidemic.

Ebola virus disease is a lethal human and primate disease that requires a particular attention from the international health authorities due to important recent outbreaks in some Western African countries and isolated cases in European and North-America continents. Regarding the emergency of this situation, various decision tools, such as mathematical models, were developed to assist the authorities to focus their efforts in important factors to eradicate Ebola. In a previous work, we have proposed an original deterministic spatial-temporal model, called Be-CoDiS (Between-Countries Disease Spread), to study the evolution of human diseases within and between countries by taking into consideration the movement of people between geographical areas. This model was validated by considering numerical experiments regarding the 2014-16 West African Ebola Virus Disease epidemic. In this article, we propose to perform a stability analysis of Be-CoDiS. Our first objective is to study the equilibrium states of simplified versions of this model, limited to the cases of one an two countries, and to determine their basic reproduction ratios. Then, in order to give some recommendations for the allocation of resources used to control the disease, we perform a sensitivity analysis of those basic reproduction ratios regarding the model parameters. Finally, we validate the obtained results by considering numerical experiments based on data from the 2014-16 West African Ebola Virus Disease epidemic.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

On the Complex Dimensions of Nonlattice Fractal Strings in Connection with Dirichlet Polynomials

In this paper we give a new characterization of the closure of the set of the real parts of the zeros of a particular class of Dirichlet polynomials that is associated with the set of dimensions of fractality of certain fractal strings. We show, for some representative cases of nonlattice Dirichlet polynomials, that the real parts of their zeros are dense in their associated critical intervals, confirming the conjecture and the numerical experiments made by M. Lapidus and M. van Frankenhuysen in several papers.

Optimization of the Modular Educational Program Structure

Optimized structure of the educational program consisting of a set of the interconnected educational objects is offered by means of problem solution of optimum partition of the acyclic weighed graph. The condition of acyclicity preservation for subgraphs is formulated and the quantitative assessment of decision options is executed. The original algorithm of search of quasioptimum partition using the genetic algorithm scheme with coding chromosomes by permutation is offered. Object-oriented realization of algorithm in language C++ is described and results of numerical experiments are presented.

Composition of basalts from the north part of the Bougault anomaly, Mid-Atlantic Ridge 15°N

A quantitative model of development of magmatic and ore-magmatic systems under crests of mid-ocean ridges is constructed. Correct physical models of melting zone formation in approximation to active spreading, non-stationary dynamics of magma intrusion from a center of generation, filling of magma chambers of various shapes, feeding of fissure-type volcanoes, and retrograde boiling of melts during solidification of intrusive bodies beneath axial zones of spreading in crests of ridges are proposed. Physicochemical and mathematical theories of disintegration of multi-component solutions, growth of liquational drops of ore melts, and sublimation of components from magmatic gases are elaborated. Methods for constructing physically correct models of heat and mass transfer in heterophase media are devised. Modeling of development of magmatic and ore-magmatic systems on the basis of the Usov-Kuznetsov facies method and the Pospelov system approach are advanced. For quantitative models numerical circuits are developed and numerical experiments are carried out.

Dynamic, convex, and robust optimization with Bayesian learning for response-guided dosing

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

Testing Independence in High Dimensions & Identifiability of Graphical Models

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

Application of the cross-entropy method to the buffer allocation problem in a simulation-based environment

The buffer allocation problem (BAP) is a well-known difficult problem in the design of production lines. We present a stochastic algorithm for solving the BAP, based on the cross-entropy method, a new paradigm for stochastic optimization. The algorithm involves the following iterative steps: (a) the generation of buffer allocations according to a certain random mechanism, followed by (b) the modification of this mechanism on the basis of cross-entropy minimization. Through various numerical experiments we demonstrate the efficiency of the proposed algorithm and show that the method can quickly generate (near-)optimal buffer allocations for fairly large production lines.

Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.

Tidal impacts on riparian salinities near estuaries

Groundwater-dependent riparian biota is known to be sensitive to changes in soil and groundwater salinity in estuarine systems. The groundwater flow and salinity behaviour in a phreatic aquifer adjoining a partially penetrating, tidal. estuary is investigated through two-dimensional numerical experiments for a lateral cross-section, which explore the influence of factors, such as aquifer and soil materials, tidal amptitudes, and regional groundwater hydraulic gradients. The density contrast between estuarine water and the fresh groundwater drives saltwater penetration of the aquifer even in the case of a marked groundwater hydraulic gradient towards the estuary. We show that tidal fluctuations in estuaries can significantly affect the groundwater salinity distribution in adjacent density-stratified phreatic aquifers. This has consequences for the expected distribution of salinity-sensitive biota in the hyporheic zone as well as vegetation and fauna dependent on water in the riparian soil and aquifer. The shape of the dense saltwater wedge propagating into the adjacent groundwater system is also modified by the estuarine tidal signal, although this effect appears to have only minor influence on the maximum distance penetrated into the aquifer (i.e., location of the 'toe' of the wedge). Tide-induced changes to riparian groundwater salinity are advection-driven, as evidenced by the modified time-averaged groundwater flow dynamics. (c) 2006 Elsevier B.V. All rights reserved.