Numerical simulation of undersea cable dynamics

A fully three-dimensional code has been written to compute the motion of a towed cable. The code is based on a robust and stable finite difference approximation to the differential equations derived from basic dynamics. A 3500-ft (1.07 km) cable pulled at 18.5 knots (34.3 km hr⁻¹) through a circular turn of 700 yd (0.64 km) radius has been computed in about half of the real time of the maneuver. The computed displacements are close to the measured ones; the changes in depth are within 2%.     

High-efficiency calculation method for watershed rainfall-runoff routing using one-dimensional dynamic wave equations

Efficiency in solving the Saint-Venant equations for watershed rainfall-runoff routing is important in flood hydrology. This paper presents a high-efficiency numerical solution of one-dimensional dynamic wave equations (HEDWE) for watershed rainfall-runoff routing, in which the full momentum equation is written as a quadratic equation with only one unknown variable Q, water depth is derived from the continuity equation using the two-step predictor-corrector method, and the discrete scheme is the explicit upwind scheme. The results of numerical tests showed the HEDWE approach has several major advantages. 1) It is a stable numerical method, even for an initially dry area. 2) Its computational efficiency is higher than 4.76E+05 times/s. 3) It can be used for overland flow, river flow, and combinations thereof. The primary disadvantages of the HEDWE approach are its unsuitability for rapidly varying flow, such as dam-break floods.     

Identification of Gleissberg Cycles and a Rising Trend in a 315-Year-Long Series of Sunspot Numbers

We show in this short note that the method of singular spectrum analysis (SSA) is able to clearly extract a strong, clean, and clear component from the longest available sunspot (International Sunspot Number, ISN) time series (1700?–?2015) that cannot be an artifact of the method and that can be safely identified as the Gleissberg cycle. This is not a small component, as it accounts for 13% of the total variance of the total original signal. Almost three and a half clear Gleissberg cycles are identified in the sunspot number series. Four extended solar minima (XSM) are determined by SSA, the latest around 2000 (Cycle 23/24 minimum). Several authors have argued in favor of a double-peaked structure for the Gleissberg cycle, with one peak between 55 and 59 years and another between 88 and 97 years. We find no evidence of the former: solar activity contains an important component that has undergone clear oscillations of \(\approx90\) years over the past three centuries, with some small but systematic longer-term evolution of “instantaneous” period and amplitude. Half of the variance of solar activity on these time scales can be satisfactorily reproduced as the sum of a monotonous multi-secular increase, a \(\approx90\)-year Gleissberg cycle, and a double-peaked (\(\approx10.0\) and 11.0 years) Schwabe cycle (the sum amounts to 46% of the total variance of the signal). The Gleissberg-cycle component definitely needs to be addressed when attempting to build dynamo models of solar activity. The first SSA component offers evidence of an increasing long-term trend in sunspot numbers, which is compatible with the existence of the modern grand maximum.     

Simulation of diagenesis and permeability variation in two-dimensional rock structure

A model is suggested to simulate the physical aspect of diagenesis in porous rocks. A bidisperse ballistic deposition model with relaxation of deposited grains is used to generate the porous structure. Sedimentation and erosion are allowed to restructure the pore space as a fluid flows through the rock. The effect of this restructuring of the pore space on permeability is studied. The Navier–Stokes equation is solved numerically by the finite difference method to determine the pressure and velocity distributions in the pore space. We find that though deposition is the dominant process in our model of diagenesis, reducing the porosity, the permeability may increase dramatically in some cases. These are when the erosion takes place at a single narrow constriction in the pore channel.     

Automatic Determination of Stellar Atmospheric Parameters Using Neural Networks and Instance-Based Machine Learning

In this article we show how machine learning methods can beeffectively applied to the problem of automatically predictingstellar atmospheric parameters from spectral information, a veryimportant problem in stellar astronomy. We apply feedforwardneural networks, Kohonen's self-organizing maps andlocally-weighted regression to predict the stellar atmosphericparameters effective temperature, surface gravity and metallicityfrom spectral indices. Our experimental results show that thethree methods are capable of predicting the parameters with verygood accuracy. Locally weighted regression gives slightly betterresults than the other methods using the original dataset asinput, while self-organizing maps outperform the other methods when significant amounts of noise are added. We also implemented a heterogeneous ensemble of predictors, combining the results given by the three algorithms. This ensemble yields better results than any of the three algorithms alone, using both the original and the noisy data.     

Long-wavelength MHD instability in the prefront of collisionless shocks with accelerated particles

Collisionless shocks in turbulent space plasmas accelerate particles by the Fermi mechanism to ultrarelativistic energies. The interaction of accelerated particles with the plasma inflow produces extended supersonic MHD flows of multicomponent plasma. We investigate the instabilities of a flow of three-component turbulent plasma with relativistic particles against long-wavelength perturbations with scales larger than the accelerated particle transport mean free path and the initial turbulence scales. The presence of turbulence allows us to formulate the system of single-fluid equations, the equation of motion for the medium as a whole, and the induction equation for the magnetic field with turbulent magnetic and kinematic viscosities. The current of accelerated particles enters into the induction equation with an effective magnetic diffusion coefficient. We have calculated the local growth rates of the perturbations related to the nonresonant long-wavelength instability of the current of accelerated particles for MHD perturbations in the WKB approximation. The amplification of long-wavelength magnetic field perturbations in the flow upstream of the shock front can affect significantly the maximum energies of the particles accelerated by a collisionless shock and can lead to the observed peculiarities of the synchrotron X-ray radiation in supernova remnants.     

Analysis of mutations in the p53 tumor suppressor gene and Ki- and Ha-ras proto-oncogenes in hepatic tumors of European flounder (Platichthys flesus).

DNA from five tumors, three other-injured livers and two normal liver tissue samples from the European flounder were analyzed for mutations in exons 5-8 of the tumor suppressor gene p53 and at codons 12, 13 and 61 of Ha- and Ki-ras proto-oncogenes. No tumor-specific mutations were identified by direct sequencing and single-strand conformation polymorphism of these genes. A number of silent polymorphisms were noted in p53. In addition to a need for more extensive analyses of flounder liver tumor samples for ras or p53 mutations, other cancer-related genes should be investigated.     

Fe and O isotope composition of meteorite fusion crusts: Possible natural analogues to chondrule formation?

Meteorite fusion crust formation is a brief event in a high‐temperature (2000–12,000 K) and high‐pressure (2–5 MPa) regime. We studied fusion crusts and bulk samples of 10 ordinary chondrite falls and 10 ordinary chondrite finds. The fusion crusts show a typical layering and most contain vesicles. All fusion crusts are enriched in heavy Fe isotopes, with δ⁵⁶Fe values up to +0.35‰ relative to the solar system mean. On average, the δ⁵⁶Fe of fusion crusts from finds is +0.23‰, which is 0.08‰ higher than the average from falls (+0.15‰). Higher δ⁵⁶Fe in fusion crusts of finds correlate with bulk chondrite enrichments in mobile elements such as Ba and Sr. The δ⁵⁶Fe signature of meteorite fusion crusts was produced by two processes (1) evaporation during atmospheric entry and (2) terrestrial weathering. Fusion crusts have either the same or higher δ¹⁸O (0.9–1.5‰) than their host chondrites, and the same is true for Δ¹⁷O. The differences in bulk chondrite and fusion crust oxygen isotope composition are explained by exchange of oxygen between the molten surface of the meteorites with the atmosphere and weathering. Meteorite fusion crust formation is qualitatively similar to conditions of chondrule formation. Therefore, fusion crusts may, at least to some extent, serve as a natural analogue to chondrule formation processes. Meteorite fusion crust and chondrules exhibit a similar extent of Fe isotope fractionation, supporting the idea that the Fe isotope signature of chondrules was established in a high‐pressure environment that prevented large isotope fractionations. The exchange of O between a chondrule melt and an ¹⁶O‐poor nebula as the cause for the observed nonmass dependent O isotope compositions in chondrules is supported by the same process, although to a much lower extent, in meteorite fusion crusts.