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.

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.

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.

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.

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.

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。