MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION (Ⅱ)——THE DISCRETE-TIME CASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.

Characteristics and Methods of Groundwater Environmental Impact Assessment in Mining Areas of Guizhou Province:A Case Study of Zhijin Area of Zhina Mining Area

In Guizhou Province,karst areas and springs are widely distributed,and hydrogeological conditions are complex in mining areas. Usually there are many hydrogeological units in a mining area,and the hydrogeological conditions are very different from that of most northern mining areas in China. In view of the uniqueness of mining areas in Guizhou Province,taking Zhijin area of Zhina mining area as an example,the planning characteristics of mining areas and characteristics of groundwater environment in Guizhou Province were analyzed firstly,and then the characteristics and key considerations of groundwater environmental impact assessment in mining areas of Guizhou Province were studied. For example,when the height of water flowing fractured zone,impact radius,and the amount of water resources affected are calculated,it is necessary to analyze and evaluate as many typical mines as possible. The impact on springs as the sources of residents＇ drinking water should be analyzed one by one.

NUMERICAL SOLUTION OF FLUID-STRUCTURE INTERACTION IN LIQUID-FILLED PIPES BY METHOD OF CHARACTERISTICS

Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.

Application of material-mesh algebraic collapsing acceleration technique in method of characteristics——based neutron transport code

The algebraic collapsing acceleration(ACA)technique maximizes the use of geometric flexibility of the method of characteristics(MOC).The spatial grids for loworder ACA are the same as the high-order transport,which makes the numerical solution of ACA equations costly,especially for large-size problems.To speed-up the MOC transport iterations effectively for general geometry,a coarse-mesh ACA method that involves selectively merging fine-mesh cells with identical materials,called material-mesh ACA(MMACA),is presented.The energy group batching(EGB)strategy in the tracing process is proposed to increase the parallel efficiency for microscopic crosssection problems.Microscopic and macroscopic crosssection benchmark problems are used to validate and analyse the accuracy and efficiency of the MMACA method.The maximum errors in the multiplication factor and pin power distributions are from the VERA-4 B-2 D case with silver-indium-cadmium(AIC)control rods inserted and are 104 pcm and 1.97%,respectively.Compared with the single-thread ACA solution,the maximum speed-up ratio reached 25 on 12 CPU cores for microscopic cross-section VERA-4-2 D problem.For the C5 G7-2 D and LRA-2 D benchmarks,the MMACA method can reduce the computation time by approximately one half.The present work proposes the MMACA method and demonstrates its ability to effectively accelerate MOC transport iterations.

A hybrid CFD/characteristics method for fast characterization of hypersonic blunt forebody/inlet flow

A hybrid CFD/characteristic method(CCM) was proposed for fast design and evaluation of hypersonic inlet flow with nose bluntness, which targets the combined advantages of CFD and method of characteristics. Both the accuracy and efficiency of the developed CCM were verified reliably, and it was well demonstrated for the external surfaces design of a hypersonic forebody/inlet with nose bluntness. With the help of CCM method, effects of nose bluntness on forebody shock shapes and the flowfield qualities which dominate inlet performance were examined and analyzed on the two-dimensional and axisymmetric configurations. The results showed that blunt effects of a wedge forebody are more substantial than that of related cone cases. For a conical forebody with a properly blunted nose, a recovery of the shock front back to that of corresponding sharp nose is exhibited, accompanied with a gradually fading out of entropy layer effects. Consequently a simplification is thought to be reasonable for an axisymmetric inlet with a proper compression angle, and a blunt nose of limited radius can be idealized as a sharp nose, as the spillage and flow variations at the entrance are negligible, even though the nose scale increases to 10% cowl lip radius. Whereas for two-dimensional inlets, the blunt effects are substantial since not only the inlet capturing/starting capabilities, but also the flow uniformities are obviously degraded.

TEST METHOD OF HIGH-SPEED ON-OFF VALVE DYNAMIC CHARACTERISTICS

According to the valve port features of high speed on-off valve and its actions, the valve port can be simplified into an a-type half bridge construction. A method that tests the dynamic characteristics of the high speed on-off valve by the output pressure signal of the a-type half bridge is proposed. Having analyzed the factors related to the dynamic characteristics of an a-type half bridge, a rule for designing the outlet chamber＇s volume is worked out. According to the rule, a test stand is built to test the self-developed high-speed on-off valve. From the test results, it can be seen that with the outlet chamber＇s volume controlled by the rule the rise time of the pressure signals driven by signals with different frequencies changes very little. The test results conform to the simulation results, which nroves the correctness of the method.

Effects of Varieties and Cooking Methods on Physical and Chemical Characteristics of Cooked Rice

To analyze the effect of different lowland rice varieties and different cooking methods on physical and chemical characteristics of cooked rice. A factorial randomized block design with two factors was used and each combination of the factors was repeated three times. The first factor was rice variety（Ciherang and Ciliwung） and the second factor was the cooking method（stovetop, boiling and steaming, and rice cooker）. Results showed that Ciherang and Ciliwung varieties were classified into slender grain rice type with yellowred color. The amylose content of Ciherang was classified as moderate, while the amylose content of Ciliwung classified as low. The most abundant amino acid contained in Ciherang and Ciliwung varieties was glutamic acid. Statistical analysis showed that cooking method had significant effects on texture, lightness, chroma, hue and moisture content of cooked rice. Rice cooked with liwet method had the lowest texture value, lowest lightness value, highest chroma value, and highest moisture content.

Bearing capacity analysis using the method of characteristics

Using the method of characteristics,the bearing capacity for a strip footing is analyzed.The method of characteristics leads to an exact true limit load when the calculations of the three terms in the bearing capacity formula are consistent with one collapse mechanism and the soil satisfies the associated flow rule.At the same time,the method of characteristics avoids the assumption of arbitrary slip surfaces,and produces zones within which equilibrium and plastic yield are simultaneously satisfied for given boundary stresses.The exact solution without superposition approximation can still be expressed by Terzaghi＇s equation of bearing capacity,in which the bearing capacity factor N γλ is dependent on the dimensionless parameter λ and the friction angle φ.The influence of groundwater on the bearing capacity of the shallow strip footing is considered,which indicates that when the groundwater effect is taken into account,the error induced by the superposition approximation can be reduced as compared with dry soil condition.The results are presented in the form of charts which give the modified value（Nwγλc /Nγλc） of bearing capacity factor.Finally,an approximated analytical expression,which provides results in close agreement with those obtained by numerical analysis in this paper,has been suggested for practical application purposes.

Dynamic Characteristics and Simplified Numerical Methods of An All-Vertical-Piled Wharf in Offshore Deep Water

There has been a growing trend in the development of offshore deep-water ports in China. For such deep sea projects, all-vertical-piled wharves are suitable structures and generally located in open waters, greatly affected by wave action. Currently, no systematic studies or simplified numerical methods are available for deriving the dynamic characteristics and dynamic responses of all-vertical-piled wharves under wave cyclic loads. In this article, we compare the dynamic characteristics of an all-vertical-piled wharf with those of a traditional inshore high-piled wharf through numerical analysis; our research reveals that the vibration period of an all-vertical-piled wharf under cyclic loading is longer than that of an inshore high-piled wharf and is much closer to the period of the loading wave. Therefore, dynamic calculation and analysis should be conducted when designing and calculating the characteristics of an all-vertical-piled wharf. We establish a dynamic finite element model to examine the dynamic response of an all-vertical-piled wharf under wave cyclic loads and compare the results with those under wave equivalent static load; the comparison indicates that dynamic amplification of the structure is evident when the wave dynamic load effect is taken into account. Furthermore, a simplified dynamic numerical method for calculating the dynamic response of an all-vertical-piled wharf is established based on the P-Y curve. Compared with finite element analysis, the simplified method is more convenient to use and applicable to large structural deformation while considering the soil non-linearity. We confirmed that the simplified method has acceptable accuracy and can be used in engineering applications.

DOUBLE METHOD OF CHARACTERISTICS TO ANALYZE HYDRAULIC-THERMAL TRANSIENTS OF PIPELINE FLOW

The hydraulic and thermal transients in pipeline flow were studied. The method of characteristics for hydraulic transient analysis of batch transport of pipeline flow had been improved. The thermal transient equation, in which the term with v 3 was involved, had been inferred, while the corresponding method of characteristics was constructed. The double method of characteristics, which can be used to study the coherent hydraulic-thermal transients of batch transport of pipeline flow, was developed.

LOCAL ERROR ESTIMATES FOR METHODS OF CHARACTERISTICS INCORPORATING STREAMLINE DIFFUSION

Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity axe derived for this scheme outside the boundary layers.

One modified method of characteristics used to analyze the multiconductor transmission lines

To solve the coupling effect of multiconductor transmission lines excited by external electromagnetic wave, the modified method of characteristics is proposed. The modified method of characteristics which can compute the terminal induced voltages excited by the external electromagnetic wave when the terminal networks or intereonnection networks contain the dynamic elements is introduced. The simulation results indicate that the modified method can analyze the terminal induced voltages when the terminal networks or the interconnection networks contain the dynamic elements excited by the external electromagnetic wave. And the results are compared with the results acquired by FDTD method, the two results are completely same. So one effective modified method is implemented to compute the transmission lines.

Effects of Modulation Methods on the Drying Characteristics and Nutritional Quality of Medicago sativa Hay

Two-year-old Medicago sativa at budding initial stage was taken as research materials.Five methods were used to make green hay,including flatting stems ＋ spraying 2.5% K2CO3,flatting stems,sun curing,drying in shade and drying under 105 ℃ condition（CK）.Besides,effects of different green hay making methods on dry characteristics and nutritional quality of M.sativa green hay were studied,and a comprehensive evaluation of M.sativa green hays was conducted.Results showed that,except CK,the drying rates in other making methods were all fast at first,and then slow down.Both of drying under 105 ℃ condition and flatting stems ＋ spraying K2CO3 could speed up drying rate and reduce nutritional losses of green hay.Sun curing could also speed up drying rate,but it could not maintain the quality of green hay.The results of Grey Relational Analysis on five green hay making methods indicated that CK had the best comprehensive performance,followed by green hays made by flatting stems ＋ spraying K2CO3.Therefore,flatting stems ＋ spraying K2CO3 was a quick and easy method to make green hay,and it was worth to be recommended in practical production.

Thermal Characteristics of PVA-PANI-ZnS Nanocomposite Film Synthesized by Gamma Irradiation Method

Gamma irradiation is employed for in situ preparation of PVA-PANI-ZnS nanocomposite. The irradiation dose is varied from 10 to 40 kGy at 10 kGy intervals. The XRD result confirms the formation of crystalline phases corresponding to ZnS nanoparticles, PVA and PANI. Field emission scanning electron microscopy shows the formation of agglomerated PANI along the PVA backbone, within which the ZnS nanoparticles are dispersed.UV-visible spectroscopy is conducted to measure the transmittance spectra of samples revealing the electronic absorption characteristics of ZnS and PANI nanoparticles. Photo-acoustic（PA） setup is installed to investigate the thermal properties of samples. The PA spectroscopy indicates a high value of thermal diffusivity for samples due to the presence of ZnS and PANI nanoparticles. Moreover, at higher doses, the more polymerization and formation of PANI and ZnS nanoparticles result in enhancement of thermal diffusivity.

CHARACTERISTIC GALERKIN METHOD FOR CONVECTION-DIFFUSION EQUATIONS AND IMPLICIT ALGORITHM USING PRECISE INTEGRATION

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.

Identification and analysis of serum samples by surface-enhanced Raman spectroscopy combined with characteristic ratio method and PCA for gastric cancer detection

This study aimed to explore the application of surface-enhanced Raman scattering(SERS)in the rapid diagnosis of gastric cancer.The SERS spectra of 68 serum samples from gastric cancer patients and healthy volunteers were acquired.The characteristic ratio method(CRM)and principal component analysis(PCA)were used to differentiate gastric cancer serum from normal serum.Compared with healthy volunteers,the serum SERS intensity of gastric cancer patients was relatively high at 722 cm^(-1),while it was relatively low at 588,644,861,1008,1235,1397,1445 and 1586 cm^(-1).These results indicated that the relative content of nucleic acids in the serum of gastric cancer patients rises while the relative content of amino acids and carbohydrates decreases.In PCA,the sensitivity and specificity of discriminating gastric cancer were 94.1%and 94.1%,respectively,with the accuracy of 94.1%.Based on the intensity ratios of four characteristic peaks at 722,861,1008 and 1397 cm^(-1),CRM presented the diagnostic sensitivity and specificity of 100%and 97.4%,respectively,and the accuracy of 98.5%.Therefore,the three peak intensity ratios of I_(722)/I_(861),I_(722)/I_(1008)and I_(722)/I_(1397)can be considered as biologicalfingerprint information for gastric cancer diagnosis and can rapidly and directly reflect the physiological and pathological changes associated with gastric cancer development.This study provides an important basis and standards for the early diagnosis of gastric cancer.

PREDICTION OF ACOUSTIC PERFORMANCE IN EXPANSION CHAMBER MUFFLERS WITH MEAN FLOW BY FINITE ELEMENT METHOD

The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect of the mean flow and geometry (length of expansion chamber and expansion ratio)on acoustic attenuation performance is discussed, the predicted values of transmission loss of expansion chamber without and with mean flow are compared with those reported in the literature and they agree well. The accuracy of the prediction of transmission loss implies that finite element approximations are applicable to a lot of practical applications.

FAST COMPUTATION OF WIDEBAND RCS USING CHARACTERISTIC BASIS FUNCTION METHOD AND ASYMPTOTIC WAVEFORM EVALUATION TECHNIQUE

The Asymptotic Waveform Evaluation (AWE) technique is an extrapolation method that provides a reduced-order model of linear system and has already been successfully used to analyze wideband electromagnetic scattering problems. As the number of unknowns increases, the size of Method Of Moments (MOM) impedance matrix grows very rapidly, so it is a prohibitive task for the computation of wideband Radar Cross Section (RCS) from electrically large object or multi-objects using the traditional AWE technique that needs to solve directly matrix inversion. In this paper, an AWE technique based on the Characteristic Basis Function (CBF) method, which can reduce the matrix size to a manageable size for direct matrix inversion, is proposed to analyze electromagnetic scattering from multi-objects over a given frequency band. Numerical examples are presented to il-lustrate the computational accuracy and efficiency of the proposed method.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.

Application of the Characteristic CIP Method to a Shallow Water Model on the Sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the effciency of the characteristic constrained interpolation profile （CIP） method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.