Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

Modified Augmented Lagrange Multiplier Methods for Large-Scale Chemical Process Optimization

Chemical process optimization can be described as large-scale nonlinear constrained minimization. The modified augmented Lagrange multiplier methods (MALMM) for large-scale nonlinear constrained minimization are studied in this paper. The Lagrange function contains the penalty terms on equality and inequality constraints and the methods can be applied to solve a series of bound constrained sub-problems instead of a series of unconstrained sub-problems. The steps of the methods are examined in full detail. Numerical experiments are made for a variety of problems, from small to very large-scale, which show the stability and effectiveness of the methods in large-scale problems.

Nonlinear constrained optimization using the flexible tolerance method hybridized with different unconstrained methods

This paper proposes the use of the flexible tolerance method(FTM) modified with scaling of variables and hybridized with different unconstrained optimization methods to solve real constrained optimization problems.The benchmark problems used to analyze the performance of the methods were taken from G-Suite functions.The original method(FTM) and other four proposed methods:(i) FTM with scaling of variables(FTMS),(ii) FTMS hybridized with BFGS(FTMS-BFGS),(iii) FTMS hybridized with modified Powell's method(FTMS-Powell)and(iv) FTMS hybridized with PSO(FTMS-PSO), were implemented. The success rates of the methods were 80%,100%, 75%, 95% and 85%, for FTM, FTMS, FTMS-BFGS, FTMS-Powell and FTMS-PSO, respectively. Numerical experiments including real constrained problems indicated that FTMS gave the best performance, followed by FTMSPowell and FTMS-PSO. Despite the inferior performance compared to FTMS and FTMS-Powell, the FTMS-PSO method presented some advantages since good different initial points could be obtained, which allow exploring different routes through the solution space and to escape from local optima. The proposed methods proved to be an effective way of improving the performance of the original FTM.

Propagation of traveling wave solutions to the Vakhnenko-Parkes dynamical equation via modified mathematical methods

In this paper, we investigate some new traveling wave solutions to Vakhnenko-Parkes equation via three modified mathematical methods. The derived solutions have been obtained including periodic and solitons solutions in the form of trigonometric, hyperbolic, and rational function solutions. The graphical representations of some solutions by assigning particular values to the parameters under prescribed conditions in each solutions and comparing of solutions with those gained by other authors indicate that these employed techniques are more effective, efficient and applicable mathematical tools for solving nonlinear problems in applied science.

A Comparative Study on Extraction Methods of Total RNA from Leaves of Acer truncatum ‘Luhong No.1’

Leaves of Acer truncatum ＇ Luhong No. 1 ＇ contain large amounts of polysaccharides and polyphenols, which seriously affect extraction yield and quality of total RNA. In order to explore the appropriate total RNA extraction method, total RNA was extracted from leaves of A. truncatum ＇ Luhong No. 1 ＇ with three methods, including kit method, Trizol method and modified CTAB method. The results showed that ODE6o/OD2so and OD^o/OD2ao ratios of total RNA extracted from leaves of A. truncatum ＇ Luhong No. 1 ＇ with kit method were higher than 1.8, with a general yield and certain level of DNA contamination ; OD^o/OD~ and OD^o/OD^o ratios of total RNA extracted with Trizol method were about 1.5, with the lowest yield; OD260/OD280 and OD260/OD230 ratios of total RNA extracted with modified CTAB method were about 2.0, with the highest yield and distinct eleetrophoresis patterns. The results demonstrated that total RNA extracted with modified CTAB method exhibited high yield and purity, which could meet the demands of subsequent molecular biology research. Thus, modified CTAB method is the appro- oriate method for extracting total RNA from leaves of A. truncatum ＇ Luhong No. 1 ＇

Solving Different Types of Differential Equations Using Modified and New Modified Adomian Decomposition Methods

The Modified Adomian Decomposition Method (MADM) is presented. A number of problems are solved to show the efficiency of the method. Further, a new solution scheme for solving boundary value problems with Neumann conditions is proposed. The scheme is based on the modified Adomian decomposition method and the inverse linear operator theorem. Several differential equations with Neumann boundary conditions are solved to demonstrate the high accuracy and efficiency of the proposed scheme.

Investigation of granite failure precursor under axial load using modified LSTM framework

Granite is usually composed of quartz,biotite,feldspar,and cracks,and the variation characteristics of these components could reflect the deformation and failure process of rock well.Taking granite as an example,the video camera was used to record the deformation and failure process of rock.The distribution of meso-components in video images was then identified.The meso-components of rock failure precursors were also discussed.Moreover,a modified LSTM(long short-term memory method)based on SSA(sparrow search algorithm)was proposed to estimate the change of meso-components of rock failure precursor.It shows that the initiation and expansion of cracks are mainly caused by feldspar and quartz fracture,and when the quartz and feldspar exit the stress framework,rock failure occurs;the second large increase of crack area and the second large decrease of quartz or feldspar area may be used as a precursor of rock failure;the precursor time of rock failure based on meso-scopic components is about 4 s earlier than that observed by the naked eye;the modified LSTM network has the strongest estimation ability for quartz area change,followed by feldspar and biotite,and has the worst estimation ability for cracks;when using the modified LSTM network to predict the precursors of rock instability and failure,quartz and feldspar could be given priority.The results presented herein may provide reference in the investigation of rock failure mechanism.

Travelling wave solutions of nonlinear conformable analytical approaches

The presented study deals with the investigation of nonlinear Bogoyavlenskii equations with conformable time-derivative which has great importance in plasma physics and non-inspectoral scattering problems.Travelling wave solutions of this nonlinear conformable model are constructed by utilizing two powerful analytical approaches,namely,the modified auxiliary equation method and the Sardar sub-equation method.Many novel soliton solutions are extracted using these methods.Furthermore,3D surface graphs,contour plots and parametric graphs are drawn to show dynamical behavior of some obtained solutions with the aid of symbolic software such as Mathematica.The constructed solutions will help to understand the dynamical framework of nonlinear Bogoyavlenskii equations in the related physical phenomena.

Numerical Simulation of Heat Transfer Process and Heat Loss Analysis in Siemens CVD Reduction Furnaces

The modified Siemens method is the dominant process for the production of polysilicon,yet it is characterised by high energy consumption.Two models of laboratory-grade Siemens reduction furnace and 12 pairs of rods industrial-grade Siemens chemical vapor deposition(CVD)reduction furnace were established,and the effects of factors such as the diameter of silicon rods,the surface temperature of silicon rods,the air inlet velocity and temperature on the heat transfer process inside the reduction furnace were investigated by numerical simulation.The results show that the convective and radiant heat losses in the furnace increased with the diameter of the silicon rods.Furthermore,the radiant heat loss of the inner and outer rings of silicon rods was inconsistent for the industrial-grade reduction furnace.As the surface temperature of the silicon rods increases,the convective heat loss in the furnace increases,while the radiative heat loss remains relatively constant.When the inlet temperature and inlet velocity increase,the convective heat loss decreases,while the radiant heat loss remains relatively constant.Furthermore,the furnace wall surface emissivity increases,resulting in a significant increase in the amount of radiant heat loss in the furnace.In practice,this can be mitigated by polishing or adding coatings to reduce the furnace wall surface emissivity.

A Modified Lagrange Method for Solving Convex Quadratic Optimization Problems

In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.

A Method Suitable for Extracting Genomic DNA from Animal and Plant——Modified CTAB Method

[Objective] The study aimed to introduce a rapid and effective method that is suitable for extracting genomic DNA from animal and plant. [ Method ] The genomic DNAs were extracted from tender leaves of 24 peanut cuhivars and from the liver, lung and kidney of white mouse through the specifically modified CTAB method. The DNAs were run on agarose gel, next detected by DNA/Protein analyzer. Finally PCR amplification was conducted to detect the quality of DNAs extracted using the modified CTAB method. [ Result] The clear and orderly bands were observed in gel detection, and the values of OD200/OD200 for DNAs extracted via modified CTAB method were between 1.77 - 1.83. The DNAs performed well in PCR amplification. [ Conclusion] The DNAs extracted by modified CTAB method could satisfy the requirement of PCR amplification.

Modified CTAB Method for Extracting Genomic DNA from Wheat Leaf

ObjectiveThe aim was to seek for a rapid DNA minipreparation method from wheat leaf. MethodThe total DNA of wheat leaf was extracted using CTAB, SDS and boiling water, separately, with some modifications. Integrity and purity of nucleic acids were detected through agarose gel electrophoresis, ultraviolet absorption and PCR. ResultThe DNA extracted by the modified CTAB method had high quality and purity, and was not degraded. Two hundreds of DNA samples could be extracted each workday by per capita using this method; and the PCR detection of wheat transgenic plants showed that amplified bands of target gene were clear, without false-positive, and the test results were satisfactory. The DNA purity and concentration extracted by modified SDS method were not as good as that extracted by modified CTAB method, but it also met the DNA requirements of major molecular research. The DNA quantity extracted by modified boiling method was small and there were a lot of impurities in it, PCR detection of this DNA showed no amplified band. ConclusionModified CTAB method is a simple and rapid method for DNA minipreparation from wheat leaf, and was suitable for PCR amplification and other molecular biology researches.

Modified embedded-atom interatomic potential for Co-W and Al-W systems

A semi-empirical interatomic potential formalism,the second-nearest-neighbor modified embedded-atom method(2NN MEAM),has been applied to obtaining interatomic potentials for the Co-W and Al-W binary system using previously developed MEAM potentials of Co,Al and W.The potential parameters were determined by fitting the experimental data on the enthalpy of formation,lattice parameter,melting point and elastic constants.The present potentials generally reproduce the fundamental physical properties of the Co-W and Al-W systems accurately.The lattice parameters,the enthalpy of formation,the thermal stability and the elastic constants match well with experiment and the first-principles results.The enthalpy of mixing and the enthalpy of formation and mixing of liquid are in good agreement with CALPHAD calculations.The potentials can be easily combined with already-developed MEAM potentials for binary cobalt systems and can be used to describe Co-Al-W-based multicomponent alloys,especially for interfacial properties.

Nonlinear Resonance of the Rotating Circular Plate under Static Loads in Magnetic Field

The rotating circular plate is widely used in mechanical engineering, meanwhile the plates are often in the electromagnetic field in modern industry with complex loads. In order to study the resonance of a rotating circular plate under static loads in magnetic field, the nonlinear vibration equation about the spinning circular plate is derived according to Hamilton principle. The algebraic expression of the initial deflection and the magneto elastic forced disturbance differential equation are obtained through the application of Galerkin integral method. By mean of modified Multiple scale method, the strongly nonlinear amplitude-frequency response equation in steady state is established. The amplitude frequency characteristic curve and the relationship curve of amplitude changing with the static loads and the excitation force of the plate are obtained according to the numerical calculation. The influence of magnetic induction intensity, the speed of rotation and the static loads on the amplitude and the nonlinear characteristics of the spinning plate are analyzed. The proposed research provides the theory reference for the research of nonlinear resonance of rotating plates in engineering.

Grouping response method for equivalent static wind loads based on a modified LRC method

Wind loading is one of the most important loads for controlling the design of large-span roof structures. Equivalent static wind loads, which can generally aim at determining a specific response, are widely used by structural designers. A method for equivalent static wind loads applicable to multi-responses is proposed in this paper. A modified load- response-correlation （LRC） method corresponding to a particular peak response is presented, and the similarity algorithm implemented for the group response is described. The main idea of the algorithm is that two responses can be put into one group if the value of one response is close to that of the other response, when the structure is subjected to equivalent static wind loads aiming at the other response. Based on the modified LRC, the grouping response method is put forward to construct equivalent static wind loading. This technique can simultaneously reproduce peak responses for some grouped responses. To verify its computational accuracy, the method is applied to an actual large-span roof structure. Calculation results show that when the similarity of responses in the same group is high, equivalent static wind loads with high accuracy and reasonable magnitude of equivalent static wind distribution can be achieved.

Carbon-supported ultrafine Pt nanoparticles modified with trace amounts of cobalt as enhanced oxygen reduction reaction catalysts for proton exchange membrane fuel cells

To accelerate the kinetics of the oxygen reduction reaction(ORR)in proton exchange membrane fuel cells,ultrafine Pt nanoparticles modified with trace amounts of cobalt were fabricated and decorated on carbon black through a strategy involving modified glycol reduction and chemical etching.The obtained Pt36Co/C catalyst exhibits a much larger electrochemical surface area(ECSA)and an improved ORR electrocatalytic activity compared to commercial Pt/C.Moreover,an electrode prepared with Pt36Co/C was further evaluated under H2-air single cell test conditions,and exhibited a maximum specific power density of 10.27 W mgPt^-1,which is 1.61 times higher than that of a conventional Pt/C electrode and also competitive with most state-of-the-art Pt-based architectures.In addition,the changes in ECSA,power density,and reacting resistance during the accelerated degradation process further demonstrate the enhanced durability of the Pt36Co/C electrode.The superior performance observed in this work can be attributed to the synergy between the ultrasmall size and homogeneous distribution of catalyst nanoparticles,bimetallic ligand and electronic effects,and the dissolution of unstable Co with the rearrangement of surface structure brought about by acid etching.Furthermore,the accessible raw materials and simplified operating procedures involved in the fabrication process would result in great cost-effectiveness for practical applications of PEMFCs.

MODIFIED LAYER REMOVAL METHOD FOR MEASUREMENT OF RESIDUAL STRESS DISTRIBUTION IN THICK PRE-STRETCHED ALUMINUM PLATE

The integrated structure parts are widely used in aircraft. The distortion caused by residual stresses in thick pre-stretched aluminum plates during machining integrated parts is a common and serious problem. To predict and control the machining distortion, the residual stress distribution in the thick plate must be measured firstly. The modified removal method for measuring residual stress in thick pre-stretched aluminum plates is proposed and the stress-strain relation matrix is deduced by elasticity theory. The residual stress distribution in specimen of 7050T7451 plate is measured by using the method, and measurement results are analyzed and compared with data obtained by other methods. The method is effective to measure the residual stress.

Some new exact solutions of Jacobian elliptic function about the generalized Boussinesq equation and Boussinesq-Burgers equation

By using the modified mapping method, we find some new exact solutions of the generalized Boussinesq equation and the Boussinesq-Burgers equation. The solutions obtained in this paper include Jacobian elliptic function solutions, combined Jacobian elliptic function solutions, soliton solutions, triangular function solutions.

Density of Ni-Al Alloys in Liquid and Solid-Liquid Coexistence State Measured by a Modified Pycnometric Method

The density of Ni-AI alloys in both liquid state and solid-liquid coexistence state was measured with a modified pycnometric method. It was found that the density of Ni-AI alloys decreases with increasing temperature and Al concentration in the alloys. The molar volume of liquid Ni-AI binary alloys increases with the increase of temperature and Al concentration. The partial molar volume of Al in Ni-AI binary alloy was calculated approximately. The molar volume of liquid Ni-AI alloy determined in the present work shows a negative deviation from the ideal linear molar volume.

Calculation of circulation in the South China Sea during summer of 2000 by the modified inverse method

On the basis of hydrographic data obtained in August 2000 cruise, the circulation in the South China Sea (SCS) is computed by the modified inverse method in combination with SSH data from TOPEX/ERS-2 analysis. For study of the dynamical mechanism, which causes the pattern of summer circulation in the SCS, the diagnostic model (Yuan et al. 1982. Acta Oceanologica Sinica,4(1):1-11; Yuan and Su. 1992. Numerical Computation of Physical Oceanography.474-542) is used to simulate numerically the summer circulation in the SCS. The following results have been obtained. (1) The central and southwestern SCSs are dominated mainly by anticy-clonic circulation systems. They are mainly as follows. 1) There is strong anticyclonic eddy southeast of Vietnam (W1). Its horizontal scale is about 300 km, and it extends vertically from the surface to the about 1 000 m level. 2) There are a warm eddy W2 southeast of Zhongsha Islands and the anticyclonic circulation system W3 west off the Luzon Island. 3) There is a stronger cyclonic eddy C1 between the anticyclonic eddies W1 and W2.4) A strong northward coastal jet is present near the coast of Vietnam, and separates from the coast of Vietnam at about 12° N to the northeast.(2)The northern SCS is dominated mainly by a cyclonic circulation system. There is a cyclonic circulation system near and north of Section N2. (3) The southeastern SCS is dominated mainly by the cyclonic circulation system. (4) Comparing the results of circulation in the SCS during the summer of 2000 with those during the summer of 1998, it is found that they agree qualitatively, but there is the some difference between them in quantity.This shows that the circulation in the SCS has obviously seasonal feature. (5) The dynamical mechanism which products the basic pattern of summer circulation is because the following two reasons: 1) the joint effect of the baroclinity and relief (JEBAR) is essential dynamical cause; and 2) it is next important dynamical cause that the interaction between the wind stress and bottom topography under the southerly monsoon. (6) Comparing the hydrographic structure and distribution of stream functions with the SSH data from TOPEX/ERS-2 analysis in the SCS during August of 2000, they agree qualitatively.