Study of storm surge due to Typhoon Linda (1997) in the Gulf of Thailand using a three dimensional ocean model

Numerical integrations using the three dimensional ocean model based on the princeton ocean model (POM) were applied for the study of both sea level elevation and ocean circulation patterns forced by the wind fields during typhoons that moved over the Gulf of Thailand (GoT). The simulation concerned a case of Typhoon Linda which occurred during November 1-4, 1997. Typhoon Linda was one of the worst storms that passed the Gulf of Thailand and hit the southern coastal provinces of Thailand on November 3, 1997. It caused flooding and a strong wind covering large areas of agriculture and fisheries, which destroyed households, utilities and even human lives. The model is the time-dependent, primitive equation, Cartesian coordinates in a horizontal and sigma coordinate in the vertical. The model grid has 37 × 97 orthogonal curvilinear grid points in the horizontal, with variable spacing from 2 km near the head of the GoT to 55 km at the eastern boundary, with 10 sigma levels in the vertical conforming to a realistic bottom topography. Open boundary conditions are determined by using radiation conditions, and the sea surface elevation is prescribed from the archiving, validation and interpretation of satellite oceanographic data (AVISO). The initial condition is determined from the spin up phase of the first model run, which was executed by using wind stress calculated from climatological monthly mean wind, restoring-type surface heat and salt and climatological monthly mean freshwater flux. The model was run in spin up phase until an ocean model reached an equilibrium state under the applied force. A spatially variable wind field taken from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used to compute the wind stress directly from the velocity fluctuations. Comparison of tendency between the sea surface elevations from model and the observed significant wave heights of moored buoys in the Gulf of Thailand under Seawatch project is investigated. The model predicts the sea level elevation up to 68.5 cm at the Cha-Am area located in the north of where the typhoon strands to the shore. Results of sea level elevation show that there is an area of peak set-up in the upper gulf, particularly in the western coast, and the effects of the storm surge are small at the lower gulf. During the entire period of this study, the surge in the gulf was induced by the northeasterly wind blowing over it.     

Numerical modeling of hydrodynamic and hydrothermal characteristics in subtropical alpine lake

A three-dimensional, time-dependent hydrodynamic and hydrothermal model was performed and applied to the subtropical alpine Yuan-Yang Lake (YYL) in northeastern region of Taiwan. The model was driven with discharge inflow, heat, and wind stress to simulate the hydrodynamic and hydrothermal in the lake. The model was validated with measured water surface elevation, current, and temperature in 2008. The overall model simulation results are in quantitative agreement with the available field data. The validated model was then used to investigate wind-driven current, mean circulation, and residence time in the YYL. The modeling results reveal that the velocity field along the wind axis present the variations over depth with return current where the velocity at the surface layer is along the wind direction while it is opposite near 1 m below water surface. The simulated mean current indicates that the surface currents flow towards the southwest direction and form a clock-wise rotation. The calculated residence time is strongly dependent on the inflows and wind effects. Regression analysis of model results reveals that an exponential regression equation can be employed to correlate the residence time to change of discharge input. The residence time without wind stress is higher than that with wind effect, indicating that wind plays an important role in lake mixing. The calculated residence time is approximately 2-2.5 days under low inflow with wind effect.     

The study of structure optimization of blast furnace cast steel cooling stave based on heat transfer analysis

The three-dimensional mathematical model of temperature and thermal stress field of cast steel cooling stave in a blast furnace has been modeled. Kinds of the parameters optimization of cast steel cooling stave in a blast furnace are proposed based on the heat transfer analysis. The results indicate that the values of the parameters optimization for a cast steel cooling stave are 200 mm for cooling channels interdistance, 25 mm for inner radius of the water channel, 180 mm for thickness of the cooling stave body, 70 mm for thickness of inlaid brick and 1.5 m/s for speed of cooling water. Reducing the water temperature would be uneconomical. The water temperature can be chosen according to the local conditions. The best choice for lining material is silicon nitrogen bond silicon carbide brick or silicon carbide brick.     

Dependence and the dimensionality reduction principle

Stone’s dimensionality reduction principle has been confirmed on several occasions for independent observations. When dependence is expressed with ϕ-mixing, a minimum distance estimate is proposed for a smooth projection pursuit regression-type function θ∈Я, that is either additive or multiplicative, in the presence of or without interactions. Upper bounds on theL ₁-risk and theL ₁-error of are obtained, under restrictions on the order of decay of the mixing coefficient. The bounds show explicitly the addive effect of ϕ-mixing on the error, and confirm the dimensionality reduction principle.     

To the theory of optimal splines

On June 18, 2008 at the Plenary Meeting of the International Conference “Differential Equations and Topology” dedicated to the 100th anniversary of Pontryagin, the report [1] was submitted by Isaev and Leitmann. This report in a summary form included a section dedicated to the research of scientists of TsAGI in the field of automation of full life-cycle (i.e. engineering-design-manufacturing, or CAE/CAD/CAM, or CALS-technologies) of wind tunnel models [2]. Within this framework, methods of geometric modeling [3] and [4] were intensively developed, new classes of optimal splines have been built, including the Pontryagin splines and the Chebyshev splines [5], [6], [7] and [8]. This paper reviews some results on the Pontryagin splines. We also give some results on the Lurie splines, that arise in the problem of interpolation of a cylindrical type surface given by the family of table coplanar planes.     

Numerical bifurcation and stability analysis for steady-states of reaction diffusion equations

Solutions of a semilinear elliptic boundary value problem, (with bounded below) can be put into a one-to-one correspondence with zeros of a function . Often d is small. The function is called the bifurcation function. It can also be shown that the eigenvalues of the matrix characterize the stability properties of the solutions of the elliptic problem as rest points of . A finite element method that can be used for computing B and B _c has recently been proposed. An overview of these results and the finite element method is given. This revised version was published online in June 2006 with corrections to the Cover Date.     

Double stage estimation of population variance

Summary Consider a normal population with mean μ and variance σ². We are interested in the estimation of population variance with the help of guess value σ ₀ ² and a sample of observations. In this paper, a double stage shrinkage estimator based on the shrinkage estimatorks ₁ ² +(1-k)σ ₀ ² ifs ₁ ² ∈R and the usual estimator ifs ₁ ² ∋R, whereR is some specified region, have been proposed. The expressions for bias and mean squared error have been obtained. Comparison with the usual estimators ² have been made. It was found that though the largest gain is obtained fork=0, we can use with 0≦k≦1/2 even when σ² is very close to σ ₀ ²     

Source mechanism model for ground motion simulation

A simple analytical model for computing ground motion in a layered half-space due to a buried seismic source is presented in this paper. The buried earthquake source is represented as a distribution of double couples varying in time as a ramp function on the fault plane. The analysis is simplified by first decoupling the governing equations into P-SV and SH problem by a coordinate transformation in the frequency-wave number domain. These two problems are solved separately and the final solution is obtained by the sum of solutions of these individual problems. Explicit expressions for ground motion in a layered half-space due to an impulsive double couple are derived. In the sequel, Green’s function for the displacement field in an infinite medium is also presented. The developed source mechanism model is also demonstrated by simulating ground motion for the Kucth earthquake (M_w = 7.7) of 26th January 2001.     

Scheduling of nonresumable jobs and flexible maintenance activities on a single machine to minimize makespan

This study addresses a single machine scheduling problem with periodic maintenance, where the machine is assumed to be stopped periodically for maintenance for a constant time w during the scheduling period. Meanwhile, the maintenance period [u, v] is assumed to have been previously arranged and the time w is assumed not to exceed the available maintenance period [u, v] (i.e. w ? v − u). The time u(v) is the earliest (latest) time at which the machine starts (stops) its maintenance. The objective is to minimize the makespan. Two mixed binary integer programming (BIP) models are provided for deriving the optimal solution. Additionally, an efficient heuristic is proposed for finding the near-optimal solution for large-sized problems. Finally, computational results are provided to demonstrate the efficiency of the models and the effectiveness of the heuristics. The mixed BIP model can optimally solve up to 100-job instances, while the average percentage error of the heuristic is below 1%.     

Verification of model for adsorption and reduction of chromium(VI) by Escherichia coli 33456 using chitosan bead as a supporting medium

A non-steady-state mathematical model system for the kinetics of adsorption and reduction of Cr(VI) by Escherichia coli (E. coli) biofilm on chitosan bead (EBCB) process was derived. The mechanisms in the model system included Cr(VI) adsorption by chitosan beads, Cr(VI) bioreduction by E. coli cells and Cr(VI) mass transport diffusion. Batch kinetic tests were performed to determine surface diffusivity of Cr(VI), adsorption parameters for Cr(VI) and biokinetic parameters of E. coli 33456. Experiments were conducted using an EBCB reactor system with high recycled rate to approximate a completely-mixed flow reactor for model verification. The experimental results indicated that E. coli biofilm bioregenerated the chitosan beads after E. coli biofilm has grown significantly. Cr(VI) reducing efficiency by E. coli was about 84% when Cr(VI) concentration in the influent was 5 mg/L at a steady-state condition. The concentration of suspended E. coli cells reached up to 10 mg/L while the thickness of attached E. coli cells was estimated to be 150 μm at a steady-state condition by model prediction. The comparisons of experimental data and model simulation show that EBCB model system for Cr(VI) adsorption and reduction can predict the experimental results well.     

Closed polynomials and saturated subalgebras of polynomial algebras

The behavior of closed polynomials, i.e., polynomials such that the subalgebra is integrally closed in , is studied under extensions of the ground field. Using some properties of closed polynomials, we prove that, after shifting by constants, every polynomial can be factorized into a product of irreducible polynomials of the same degree. We consider some types of saturated subalgebras , i.e., subalgebras such that, for any , a generative polynomial of f is contained in A. Published in Ukrains’kyi Matematychnyi Zhurnal, Vol. 59, No. 12, pp. 1587–1593, December, 2007.     

Some trace formulae involving the split sequences of a Leonard pair

Let K denote a field, and let V denote a vector space over K with finite positive dimension. We consider a pair of linear transformations A : V → V and A^∗ : V → V that satisfy (i) and (ii) below:

(i)

There exists a basis for V with respect to which the matrix representing A is irreducible tridiagonal and the matrix representing A^∗ is diagonal.

(ii)

There exists a basis for V with respect to which the matrix representing A^∗ is irreducible tridiagonal and the matrix representing A is diagonal.

We call such a pair a Leonard pair on V. Let diag(θ₀, θ₁, … , θ_d) denote the diagonal matrix referred to in (ii) above and let denote the diagonal matrix referred to in (i) above. It is known that there exists a basis u₀, u₁, … , u_d for V and there exist scalars ?₁, ?₂, … , ?_d in K such that Au_i = θ_iu_i + u_i+1 (0 ? i ? d − 1), Au_d = θ_du_d, , . The sequence ?₁, ?₂, … , ?_d is called the first split sequence of the Leonard pair. It is known that there exists a basis v₀, v₁, … , v_d for V and there exist scalars ?₁, ?₂, … , ?_d in K such that Av_i = θ_d−iv_i + v_i+1 (0 ? i ? d − 1),Av_d = θ₀v_d, , . The sequence ?₁, ?₂, … , ?_d is called the second split sequence of the Leonard pair. We display some attractive formulae for the first and second split sequence that involve the trace function.     

Compact differences of weighted composition operators from weighted Bergman spaces to weighted-type spaces

We characterize the compactness of differences of weighted composition operators from the weighted Bergman space , 0 < p < ∞, α > −1, to the weighted-type space of analytic functions on the open unit disk D in terms of inducing symbols and . For the case 1 < p < ∞ we find an asymptotically equivalent expression to the essential norm of these operators.