NUMERICAL SIMULATION STUDY OF THE INNER-CORE STRUCTURES AND THE MECHANISM FOR INSHORE STRENGTHENING OF SOUTH CHINA SEA TYPHOON VONGFONG (0214) DURING LANDFALL

An explicit simulation with a fine mesh at intervals of 6 km is used to explore the inner-core structures of Vongfong (0214). The dynamic mechanism for the inshore strengthening of Vongfong is examined. It is found as follows. (1) The radius of maximum wind of the axisymmetric structures of the typhoon decreased with height during its mature stage. When Vongfong was inshore, the strongest low-layer inflow located in front of it and the outflow was to the rear of it, which was just reversed from the Atlantic hurricanes and other Pacific typhoons. (2) The dynamic and thermodynamic fields were highly asymmetric in structure. Convection was stronger in the northwest quadrant of the typhoon than in the southeast; the strongest convective cloud bands were consistent with the maximum wind region. During its strengthening stage, it was cold west of and warm east of the eye in the lower layer but warm in the west and cold in the east of the mid-upper layer. During its mature stage, a warm-core structure was evident in the lower and middle-upper layer. (3) The interactions between a mid-latitude cold low in the middle-upper troposphere and the typhoon were responsible for the latter to strengthen inshore. Firstly, the outer circulation of the cold low entered the typhoon from the middle troposphere when an outer cold airflow from the cold low flowed into the northwest quadrant of the typhoon so that geopotentially instable energy increased and convection developed. Secondly, the downdraft in the cold low was just the corresponding branch of the secondary circulation of the typhoon system; when the cold low weakened while moving south, the typhoon strengthened inshore. Due to the CISK mechanism, these two phenomena might be realized.

A study on the numerical prediction method for the vertical thermal structure in the Bohai Sea and the Huanghai Sea-I.One-dimensional numerical prediction model

In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( Tt ) and the thickness of the upper homogeneous layer ( h ) is developed in terms of the dimensionless temperature θT and depth η and self-simulation function θT - f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T, Th, h , the thickness and gradient of thermocline are satisfactory to some extent.

Numerical study on three-dimensional flow field of continuously rotating detonation in a toroidal chamber

Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.

Comparison of structure and physical fields in 400 kA aluminum reduction cells

To investigate the differences and the development trends of the 400 kA aluminum reduction cell, four representative cells were deeply analyzed. By using numerical simulation methods in ANSYS software, the structure parameters were firstly compared, and then three-dimensional models of electric-magnetic-flow field were built and solved with finite element method(FEM). The comparison of the structures reveals that the cell bodies are similar while the current flow path and distribution ratio of bus bars are different. It appears that most of the current(70%-80%) in side A are used as the magnetic field compensation current and flow through two ends. The numerical simulation results indicate that the distributions of magnetic fields are different but all satisfy with the magnetohydrodynamics(MHD) stabilization, and the flow patterns are all two or multi vortexes with appropriate velocities. The comparison shows that all studied cells can satisfy with the physical field requirement, and the commercial applications also verify that the 400 kA cells have become the product of the mature and world's leading technology.

Numerical simulation of vortex-induced drag of elastic swimmer models

We present numerical simulations of simplified models for swimming organisms or robots, using chordwise flexible elastic plates. We focus on the tip vortices originating from three-dimensional effects due to the finite span of the plate. These effects play an important role when predicting the swimmer＇s cruising velocity, since they contribute significantly to the drag force. First we simulate swimmers with rectangular plates of different aspect ratios and compare the results with a recent experimental study. Then we consider plates with expanding and contracting shapes. We find the cruising velocity of the contracting swimmer to be higher than the rectangular one, which in turn is higher than the expanding one. We provide some evidence that this result is due to the tip vortices interacting differently with the swimmer.

Catastrophe mechanism and disaster countermeasure for soft rock roadway surrounding rock in Meihe mine

The soft rock's heterogeneity and nonlinear mechanical behavior cause extremely difficult maintenance on the soft rock roadway. Aiming at the asymmetric deformation and destruction phenomenon appearing after excavating and supporting the 7101 air return way in Meihe mine, this paper comprehensively adopted a variety of methods to analyze the roadway surrounding rock deformation rule, obtaining the roadway surrounding rock stress and plastic zone distribution rule under no supporting condition and the roadway surrounding rock deformation features under original symmetric supporting condition.Furthermore, this paper revealed the catastrophe mechanism, and proposed the concept of ‘‘weak structure'' and the disaster countermeasure of ‘‘overall stabilizing the roadway and strengthening the support of weak structure''. The industrial test shows that the disaster control technology can realize the coordination deformation of the supporting structure and roadway surrounding rock, thus significantly controlling the deformation of roadway surrounding rock.

Modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies

We present an efficient, robust computational method for modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies. With appropriate mathematical treatments, the apparent numerical difficulties associated with singularities in computing elliptic integrals are completely removed. Using a boundary element discretization procedure, the governing equations are transformed into a linear algebra matrix equation that can be solved by straightforward Gauss elimination in one step without further iterations. The numerical code implemented according to our algorithm can accurately determine the surface mass density distribution in a disk galaxy from a measured rotation curve （or vice versa）. For a disk galaxy with a typical flat rotation curve, our modeling results show that the surface mass density monotonically decreases from the galactic center toward the periphery, according to Newtonian dynamics. In a large portion of the galaxy, the surface mass density follows an approximately exponential law of decay with respect to the galactic radial coordinate. Yet the radial scale length for the surface mass density seems to be generally larger than that of the measured brightness distribution, suggesting an increasing mass-tolight ratio with the radial distance in a disk galaxy. In a nondimensionalized form, our mathematical system contains a dimensionless parameter which we call the ＂galactic rotation number＂ that represents the gross ratio of centrifugal force and gravitational force. The value of this galactic rotation number is determined as part of the numerial solution. Through a systematic computational analysis, we have illustrated that the galactic rotation number remains within 4-10% of 1.70 for a wide variety of rotation curves. This implies that the total mass in a disk galaxy is proportional to V02 Rg, with V0 denoting the characteristic rotation velocity （such as the ＂flat＂ value in a typical ro- tation curve） and Rg the radius of the galactic disk. The predicted total galactic mass of the Milky Way is in good agreement with the star-count data.

Reduced technique for modeling electromagnetic immunity on braid shielding cable bundles

In this paper, an efficient multi-conductor simplification technique is proposed to model the electromagnetic immunity on cable bundles within a braid shielding structure over a large frequency range. By grouping together the conductors based on the knowledge of Z-Smith chart, the required computation time is markedly reduced and the complexity of modeling the completely shielding cable bundles is significantly simplified with a good accuracy. After a brief description of the immunity problems in shielding structure, a six-phase procedure is detailed to generate the geometrical characteristics of the reduced cable bundles. Numerical simulation is carried out by using a commercial software CST to validate the efficiency and advantages of the proposed approach. The research addressed in this paper is considered as a simplified modeling technique for the electromagnetic immunity within a shielding structure.

Thermal analysis and optimization of the EAST ICRH antenna

The ion cyclotron resonance of frequency heating（ICRH） plays an important role in plasma heating.Two ICRH antennas were designed and applied on the EAST tokamak.In order to meet the requirement imposed by high-power and long-pulse operation of EAST in the future,an active cooling system is mandatory to be designed to remove the heat load deposited on the components.Thermal analyses for high heat-load components have been carried out,which presented clear temperature distribution on each component and provided the reference data to do the optimization.Meanwhile,heat pipes were designed to satisfy the high requirement imposed by a Faraday shield and lateral limiter.

Made-to-measure galaxy modelling utilising absorption line strength data

We enhance the Syer ＆ Tremaine made-to-measure （M2M） particle method of stellar dynamical modelling to model simultaneously both kinematic data and absorption line strength data, thus creating a ＇chemo-M2M＇ modelling scheme. We apply the enhanced method to four galaxies （NGC 1248, NGC 3838, NGC 4452, NGC 4551） observed using the SAURON integral-field spectrograph as part of the ATLAS3D programme. We are able to reproduce successfully the 2D line strength data achieving mean X2 per bin values of ≈ 1 with 〉 95% of particles having converged weights. Because M2M uses a 3D particle system, we are also able to examine the underlying 3D line strength distributions. The extent to which these dis- tributions are plausible representations of real galaxies requires further consideration. Overall, we consider the modelling exercise to be a promising first step in developing a ＇chemo-M2M＇ modelling system and in understanding some of the issues to be addressed. While the made-to-measure techniques developed have been applied to absorption line strength data, they are in fact general and may be of value in modelling other aspects of galaxies.

Recovering 3D structural properties of galaxies from SDSS-like photometry

Because of the 3D nature of galaxies, an algorithm for constructing spatial density distribution models of galaxies on the basis of galaxy images has many advan- tages over approximations of the surface density distribution. We present a method for deriving the spatial structure and overall parameters of galaxies from images and estimate its accuracy and derived parameter degeneracies on a sample of idealised model galaxies. The test galaxies consist component with varying proportions and of a disc-like component and a spheroidal properties. Both components are assumed to be axially symmetric and coplanar. We simulate these test galaxies as if they had been observed in the SDSS project through ugriz filters, thus gaining a set of realis- tically imperfect images of galaxies with known intrinsic properties. These artificial SDSS galaxies were thereafter remodelled by approximating the surface brightness distribution with a 2D projection of a bulge＋disc spatial distribution model and the restored parameters were compared to the initial ones. Down to the r-band limiting magnitude of 18, errors in the restored integral luminosities and colour indices re- main within 0.05 mag and errors in the luminosities of individual components within 0.2 mag. Accuracy of the restored bulge-to-disc luminosity ratio （B/D） is within 40% in most cases, and becomes worse for galaxies with low B/D, but the general balance between bulges and discs is not shifted systematically. Assuming that the intrinsic disc axial ratio is ≤ 0.3, then the inclination angles can be estimated with errors 〈 5° for most of the galaxies with B/D 〈 2 and with errors 〈 15° up to B/D = 6. Errors in the recovered sizes of the galactic components are below 10% in most cases. The axial ratios and the shape parameter N of Einasto＇s distribution （similar to the Sersic index） are relatively inaccurate, but can provide statistical estimates for large samples. In general, models of disc components are more accurate than models of spheroidal components for geometrical reasons.

Numerical study on three-dimensional CJ detonation waves interacting with isotropic turbulence

The three-dimensional structures of a cellular detonation wave interacting with different turbulent flows were investigated using a one-step irreversible Arrhenius kinetics model. High-resolution bandwidth-optimized WENO scheme of spatial discretization and total variation diminishing temporal integration are used to solve the three dimensional chemically reactive Navier-Stokes equations. The turbulent vertical and entropic forcing effects on the three dimensional detonation wave structures and dynam- ics are analyzed, as well as the detonation effects on tur- bulent vortex structures. It has been found that the turbulence field imposed has created small scale wrinkles embedded in the detonation front, apart from the large scale features of detonation without turbulence. The deto- nation propagating velocity over the leading shock front varies from 0.8 to 1.6 times of CJ velocity and its proba- bility density function （pdf） skews towards sub-CJ velocity and peaks at about 0.9. The recorded detonation velocity always preferentially decays with time, with very rapid accelerations through triple point interactions. Its pdf also skews to sub-CJ velocity, while its overall shape agrees well with W3. The reaction zone is greatly influenced by the vortex, much more irregular and elongated for the turbulent cases. Distributed burning pockets are more likely to be found there. The turbulent kinetic energy is amplified across the detonation, and periodically oscillates downstream the detonation. The off-diagonal components of Reynolds stress also show a rapid rise across the deto- nation and present to be non-zero downstream of detona- tion. Vortex structures are compound results of the convected vortex and the generated vortex by the collision of triple points. The convection term and baroclinic gen- eration term in the transport equation of enstrophy are compared in detail.

Simulation of 3D parachute fluid–structure interaction based on nonlinear finite element method and preconditioning finite volume method

A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual（GMRES） method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor（HHT） time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC（Harten–Lax–van Leer contact） scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel（LUSGS） approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.

Light-weighting in aerospace component and system design

Light-weighting involves the use of advanced materials and engineering methods to enable structural elements to deliver the same,or enhanced,technical performance while using less material.The concept has been extensively explored and utilised in many industries from automotive applications to fashion and packaging and offers significant potential in the aviation sector.Typical implementations of light-weighting have involved use of high performance materials such as composites and optimisation of structures using computational aided engineering approaches with production enabled by advanced manufacturing methods such as additive manufacture,foam metals and hot forming.This paper reviews the principal approaches used in light-weighting,along with the scope for application of light-weighting in aviation applications from power-plants to airframe components.A particular area identified as warranting attention and amenable to the use of lightweighting approaches is the design of solar powered aircraft wings.The high aspect ratio typically used for these can be associated with insufficient stiffness,giving rise to non-linear deformation,aileron reversal,flutter and rigid-elastic coupling.Additional applications considered include ultralight aviation components and sub-systems,UAVs,and rockets.Advanced optimisation approaches can be applied to optimise the layout of structural elements,as well as geometrical parameters in order to maximise structural stiffness,minimise mass and enable incorporation of energy storage features.The use of additive manufacturing technologies,some capable of producing composite or multi-material components is an enabler for light-weighting,as features formally associated with one principal function can be designed to fulfil multiple functionalities。