A sub-grid scale model for Burgers turbulence based on the artificial neural network method

The present study proposes a sub-grid scale model for the one-dimensional Burgers turbulence based on the neuralnetwork and deep learning method.The filtered data of the direct numerical simulation is used to establish thetraining data set,the validation data set,and the test data set.The artificial neural network(ANN)methodand Back Propagation method are employed to train parameters in the ANN.The developed ANN is applied toconstruct the sub-grid scale model for the large eddy simulation of the Burgers turbulence in the one-dimensionalspace.The proposed model well predicts the time correlation and the space correlation of the Burgers turbulence.

A SUPERLINEARLY CONVERGENT SPLITTING FEASIBLE SEQUENTIAL QUADRATIC OPTIMIZATION METHOD FOR TWO-BLOCK LARGE-SCALE SMOOTH OPTIMIZATION

This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.

A LIMITED MEMORY QUASI-NEWTON METHOD FOR LARGE SCALE PROBLEM

We study how to use the SR1 update to realize minimization methods for problems where the storage is critical. We give an update formula which generates matrices using information from the last m iterations. The numerical tests show that the method is efficent.

Hazard degree identification of goafs based on scale effect of structure by RS-TOPSIS method

In order to precisely predict the hazard degree of goaf(HDG), the RS-TOPSIS model was built based on the results of expert investigation. To evaluate the HDG in the underground mine, five structure size factors, i.e. goaf span, exposed area, goaf height, goaf depth, and pillar width, were selected as the evaluation indexes. And based on rough dependability in rough set(RS)theory, the weights of evaluation indexes were identified by calculating rough dependability between evaluation indexes and evaluation results. Fourty goafs in some mines of western China, whose indexes parameters were measured by cavity monitoring system(CMS), were taken as evaluation objects. In addition, the characteristic parameters of five grades' typical goafs were built according to the interval limits value of single index evaluation. Then, using the technique for order preference by similarity to ideal solution(TOPSIS), five-category classification of HDG was realized based on closeness degree, and the HDG was also identified.Results show that the five-category identification of mine goafs could be realized by RS-TOPSIS method, based on the structure-scale-effect. The classification results are consistent with those of numerical simulation based on stress and displacement,while the coincidence rate is up to 92.5%. Furthermore, the results are more conservative to safety evaluation than numerical simulation, thus demonstrating that the proposed method is more easier, reasonable and more definite for HDG identification.

Research on Analytical Method of Fatigue Characteristics of Soft Yoke Mooring System Based on Full-Scale Measurement

By focusing on the vulnerability of the structure of marine equipments, together with considering the randomness of meta-ocean load in statistics, a kind of analytical method of fatigue characteristics of marine structure based on full- scale and actual measurement of prototype is proposed. On the basis of short-term field measurement results of structural response, research is carried out on the fatigue analysis of hinge joints at the upper part of the Soft Yoke single point Mooring System （SYMS） by simultaneously monitoring the environmental load and considering the design criteria of offshore structure. Through analysis of finite element modeling, the time-histories of typical stress response are obtained, and then the assessment of fatigue damage at key hinge joints is conducted. The simulation results indicate that the proposed method can accurately analyze the fatigue damage of offshore engineering structure caused by the effect of wave load. The present analytical method of fatigue characteristics can be extended on other offshore engineering structures subjected to meta-ocean load.

Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.

Modified slow-fast analysis method for slow-fast dynamical systems with two scales in frequency domain

A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term as a slow-varying parameter,a generalized autonomous fast subsystem can be defined,the equilibrium branches as well as the bifurcations of which can be employed to account for the mechanism of the bursting oscillations by combining the transformed phase portrait introduced.As an example,a typical periodically excited Hartley model is used to demonstrate the validness of the method,in which the exciting frequency is far less than the natural frequency.The equilibrium branches and their bifurcations of the fast subsystem with the variation of the slow-varying parameter are presented.Bursting oscillations for two typical cases are considered,which reveals that,fold bifurcation may cause the the trajectory to jump between different equilibrium branches,while Hopf bifurcation may cause the trajectory to oscillate around the stable limit cycle.

A scaled boundary node method applied to two-dimensional crack problems

A boundary-type meshless method called the scaled boundary node method （SBNM） is developed to directly evaluate mixed mode stress intensity factors （SIFs） without extra post-processing. The SBNM combines the scaled boundary equations with the moving Kriging （MK） interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter. As a result, the SBNM requires only a set of scattered nodes on the boundary, and the displacement field is approximated by using the MK interpolation technique, which possesses the 5 function property. This makes the developed method efficient and straightforward in imposing the essential boundary conditions, and no special treatment techniques are required. Besides, the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction. Therefore, the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip. Numerical examples using the SBNM for computing the SIFs are presented. Good agreements with available results in the literature are obtained.

Determination of the natural frequencies of axially moving beams by the method of multiple scales

The natural frequencies of an axially moving beam were determined by using the method of multiple scales. The method of second-order multiple scales could be directly applied to the governing equation if the axial motion of the beam is assumed to be small. It can be concluded that the natural frequencies affected by the axial motion are proportional to the square of the velocity of the axially moving beam. The results obtained by the perturbation method were compared with those given with a numerical method and the comparison shows the correctness of the multiple-scale method if the velocity is rather small.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

Prediction of grain scale plasticity of NiTi shape memory alloy based on crystal plasticity finite element method

Grain scale plasticity of NiTi shape memory alloy(SMA)during uniaxial compression deformation at 400℃was investigated through two-dimensional crystal plasticity finite element simulation and corresponding analysis based on the obtained orientation data.Stress and strain distributions of the deformed NiTi SMA samples confirm that there exhibits a heterogeneous plastic deformation at grain scale.Statistically stored dislocation(SSD)density and geometrically necessary dislocation(GND)density were further used in order to illuminate the microstructure evolution during uniaxial compression.SSD is responsible for sustaining plastic deformation and it increases along with the increase of plastic strain.GND plays an important role in accommodating compatible deformation between individual grains and thus it is correlated with the misorientation between neighboring grains,namely,a high GND density corresponds to large misorientation between grains and a low GND density corresponds to small misorientation between grains.

Numerical manifold method modeling of coupled processes in fractured geological media at multiple scales

The greatest challenges of rigorously modeling coupled hydro-mechanical(HM)processes in fractured geological media at different scales are associated with computational geometry.These challenges include dynamic shearing and opening of intersecting fractures at discrete fracture scales as a result of coupled processes,and contact alteration along rough fracture surfaces that triggers structural and physical changes of fractures at micro-asperity scale.In this paper,these challenges are tackled by developing a comprehensive modeling approach for coupled processes in fractured geological media based on numerical manifold method(NMM)at multiple scales.Based on their distinct geometric features,fractures are categorized into three different scales:dominant fracture,discrete fracture,and discontinuum asperity scales.Here the scale is relative,that of the fracture relative to that of the research interest or domain.Different geometric representations of fractures at different scales are used,and different governing equations and constitutive relationships are applied.For dominant fractures,a finite thickness zone model is developed to treat a fracture as a porous nonlinear domain.Nonlinear fracture mechanical behavior is accurately modeled with an implicit approach based on strain energy.For discrete fractures,a zero-dimensional model was developed for analyzing fluid flow and mechanics in fractures that are geometrically treated as boundaries of the rock matrix.With the zero-dimensional model,these fractures can be modeled with arbitrary orientations and intersections.They can be fluid conduits or seals,and can be open,bonded or sliding.For the discontinuum asperity scale,the geometry of rough fracture surfaces is explicitly represented and contacts involving dynamic alteration of contacts among asperities are rigorously calculated.Using this approach,fracture alteration caused by deformation,re-arrangement and sliding of rough surfaces can be captured.Our comprehensive model is able to handle the computational challenges with accurate representation of intersections and shearing of fractures at the discrete fracture scale and rigorously treats contacts along rough fracture surfaces at the discontinuum asperity scale.With future development of three-dimensional(3D)geometric representation of discrete fracture networks in porous rock and contacts among multi-body systems,this model is promising as a basis of 3D fully coupled analysis of fractures at multiple scales,for advancing understanding and optimizing energy recovery and storage in fractured geological media.

The SA Method on Large Scale Goal Programming with Independent Subsystems

ＴｈｅＳＡＭｅｔｈｏｄｏｎＬａｒｇｅＳｃａｌｅＧｏａｌＰｒｏｇｒａｍｍｉｎｇｗｉｔｈＩｎｄｅｐｅｎｄｅｎｔＳｕｂｓｙｓｔｅｍｓＺＨＡＮＧＪｉｅ；ＦＥＮＧＹｉｎｇｊｕｎ张杰，冯英浚（Ｄｅｐｔ．ｏｆＭａｔｈｅｍａｔｉｃｓ，ＨａｒｂｉｎＩｎｓｔｉｔｕｔｅｏ...

Investigation on the biodegradation levels of super heavy oils by parameter-striping method and refined Manco scale: a case study from the Chepaizi Uplift of Junggar Basin

The Carboniferous volcanic reservoir in the Chepaizi Uplift became an exploration hot target in recent years for its substantial amount of oils discovered. However, most of the Carboniferous heavy oils were biodegraded to PM7 or higher with orders of magnitude variation in oil viscosities. Two oil groups (I and II) exactly corresponding to the western and eastern Chepaizi Uplift were distinguished according to their source diagnose. Furthermore, three oil families (II1, II2 and II3), with the biodegradation level of PM7, PM8–8+, PM9+, respectively, were classified based on molecular compositions and parameter-stripping method of strongly bioresistant parameters. Allowing for this extremely high biodegradation case, more biodegradation refractory compound class were added to establish a refined Manco scale to quantitatively evaluate the biodegradation extent. Refined Manco number (RMN2) positively correlated with the oil density, NSO contents, and absolute concentrations of diasteranes and gammacerane, negatively correlated with the absolute concentrations of diahopane, summed tricyclic terpanes and pentacyclic terpanes. This refined scale showed higher resolution than the PM one to differentiate the biodegradation extent of Carboniferous heavy oils from the Chepaizi Uplift, especially those with same PM values but different oil viscosities.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

Determination of an Unknown Source in the Heat Equation by the Method of Tikhonov Regularization in Hilbert Scales

In this paper, we consider the problem for determining an unknown source in the heat equation. The Tikhonov regularization method in Hilbert scales is presented to deal with ill-posedness of the problem and error estimates are obtained with a posteriori choice rule to find the regularization parameter. The smoothness parameter and the a priori bound of exact solution are not needed for the choice rule. Numerical tests show that the proposed method is effective and stable.

Dynamic Crack Propagation Analysis Using Scaled Boundary Finite Element Method

The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.

A two-scale method for identifying mechanical parameters of composite materials with periodic configuration

In this paper, a two-scale method （TSM） is presented for identifying the mechanics parameters such as stiffness and strength of composite materials with small periodic configuration. Firstly, a formulation is briefly given for two-scale analysis （TSA） of the composite materials. And then a two-scale computation formulation of strains and stresses is developed by displacement solution with orthotropic material coefficients for three kinds of such composites structures, i.e., the tension column with a square cross section, the bending cantilever with a rectangular cross section and the torsion column with a circle cross section. The strength formulas for the three kinds of structures are derived and the TSM procedure is discussed. Finally the numerical results of stiffness and strength are presented and compared with experimental data. It shows that the TSM method in this paper is feasible and valid for predicting both the stiffness and the strength of the composite materials with periodic configuration.

Two-scale finite element method for piezoelectric problem in periodic structure

The prediction of the mechanical and electric properties of piezoelectric fibre composites has become an active research area in recent years. By means of introducing a boundary layer problem, some new kinds of two-scale finite element methods for solutions to the electric potential and the displacement for composite material in periodic struc- ture under the coupled piezoelectricity are derived. The coupled two-scale relation of the electric potential and the displacement is set up, and some finite element approximate estimates and numerical examples which show the effectiveness of the method are presented.

Fabrication of micro-scale gratings for moiré method with a femtosecond laser

Fabrication of micro gratings using a femtosecond laser exposure system is experimentally investigated for the electron moire method. Micro holes and lines are firstly etched for parameter study. Grating profile is theoretically optimized to form high quality moire patterns. For a demonstration, a parallel grating is fabricated on a specimen of quartz glass. The minimum line width and the distance between two adjacent lines are both set to be 1 μm, and the frequency of grating is 500 lines/ram. The experimental results indicate that the quality of gratings is good and the relative error of the gratings pitch is about 1.5%. Based on molte method, scanning electron microscope （SEM） moire patterns are observed clearly, which manifests that gratings fabricated with the femtosecond laser exposure is suitable for micro scale deformation measurement.