Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

The methane hydrate formation and the resource estimate resulting from free gas migration in seeping seafloor hydrate stability zone

It is a typical multiphase flow process for hydrate formation in seeping seafloor sediments. Free gas can not only be present but also take part in formation of hydrate. The volume fraction of free gas in local pore of hydrate stable zone (HSZ) influences the formation of hydrate in seeping seafloor area, and methane flux determines the abundance and resource of hydrate-bearing reservoirs. In this paper, a multiphase flow model including water (dissolved methane and salt)-free gas hydrate has been established to describe this kind of flow-transfer-reaction process where there exists a large scale of free gas migration and transform in seafloor pore. In the order of three different scenarios, the conversions among permeability, capillary pressure, phase saturations and salinity along with the formation of hydrate have been deducted. Furthermore, the influence of four sorts of free gas saturations and three classes of methane fluxes on hydrate formation and the resource has also been analyzed and compared. Based on the rules drawn from the simulation, and combined information gotten from drills in field, the methane hydrate(MH) formation in Shenhu area of South China Sea has been forecasted. It has been speculated that there may breed a moderate methane flux below this seafloor HSZ. If the flux is about 0.5 kg m-2 a-1, then it will go on to evolve about 2700 ka until the hydrate saturation in pore will arrive its peak (about 75%). Approximately 1.47 109 m3 MH has been reckoned in this marine basin finally, is about 13 times over preliminary estimate.

Parallel Fully-Implicit Computation of Magnetohydrodynamics Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETS,: package. The results of convergence tests arc plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

层析静校正技术研究和应用

Static correction is one of the indispensable steps in the conventional onshore seismic data processing, particularly in the western part of China; it is theoretically and practically significant to resolve the issue of static correction. Conventional refraction static correction is put forward under the assumption that layered medium is horizontal and evenly distributed. The complicated nature of the near surface from western part of China is far from the assumption. Therefore, the essential way to resolve the static correction problem from the complex area is to develop a new theory. In this paper, a high-precision non-linear first arrival tomography is applied to solve the problem, it moved beyond the conventional refraction algorithm based on the layered medium and can be used to modeling the complex near surface. Some of the new and creative work done is as follows: One. In the process of first arrival tomographic image modeling, a fast high-order step algorithm is used to calculate the travel time for first arrival and ray path and various factors concerning the fast step ray tracing algorithm is analyzed. Then the second-order and third-order differential format is applied to the step algorithm which greatly increased the calculation precision of the ray tracing and there is no constraint to the velocity distribution from the complex areas. This method has very strong adaptability and it can meet the needs of great velocity variation from the complicated areas. Based on the numerical calculation, a fast high-order step is a fast, non-conditional and stable high-precision tomographic modeling algorithm. Two, in the tomographic inversion, due to the uneven fold coverage and insufficient information, the inversion result is unstable and less reliable. In the paper, wavelet transform is applied to the tomographic inversion which has achieved a good result. Based on the result of the inversion from the real data, wavelet tomographic inversion has increased the reliability and stability of the inversion. Three. Apply the constrained high-precision wavelet tomographic image to the static correction processing from the complex area. During tomographic imaging, by using uphole survey, refraction shooting or other weathering layer method, weathering layer can be identified before the image. Because the group interval for the shot first arrival is relatively big, there is a lack of precision for the near surface inversion. In this paper, an inversion method of the layer constraint and well constraint is put forward, which can be used to compensate the shallow velocity of the inversion for the shot first arrival and increase the precision of the tomographic inversion. Key words: Tomography ，Fast marching method，Wavelet transform, Static corrections, First break

地震岩性与物性联合反演技术研究和应用

This paper studies how to more effectively invert seismic data and predict reservoir under complicated sedimentary environment, complex rock physical relationships and fewer drills in offshore areas of China. Based on rock physical and seismic amplitude-preserving process, and according to depositional system and laws of hydrocarbon reservoir, in the light of feature of seismic inversion methods present applied, series methods were studied. A joint inversion technology for complex geological condition had been presented, at the same time the process and method system for reservoir prediction had been established. This method consists four key parts. 1)We presented the new conception called generalized wave impedance, established corresponding inversion process, and provided technical means for joint inversion lithology and petrophysical on complex geological condition. 2)At the aspect of high-resolution nonlinear seismic wave impedance joint inversion, this method used a multistage nonlinear seismic convolution model rather than conventional primary structure Robinson seismic convolution model, and used Caianiello neural network implement inversion. Based on the definition of multistage positive and negative wavelet, it adopted both deterministic and statistical physical mechanism, direct inversion and indirect inversion. It integrated geological knowledge, rock physical theory, well data, and seismic data, and improved the resolution and anti-noise ability of wave impedence inversion. 3)At the aspect of high-resolution nonlinear reservoir physical property joint inversion, this method used nonlinear rock physical model which introduced convolution model into the relationship between wave impedance and porosity/clay. Through multistage decomposition, it handles separately the large- and small-scale components of the impedance-porosity/clay relationships to achieve more accurate rock physical relationships. By means of bidirectional edge detection with wavelets, it uses the Caianiello neural network to finish statistical inversion with combined applications of model-based and deconvolution-based methods. The resulted joint inversion scheme can integrate seismic data, well data, rock physical theory, and geological knowledge for estimation of high-resolution petrophysical parameters. 4)At the aspect of risk assessment of lateral reservoir prediction, this method integrated the seismic lithology identification, petrophysical prediction, multi-scale decomposition of petrophysical parameters, P- and H-spectra, and the match relationship of data got from seismics, well logging and geology. It could describe the complexity of medium preferably. Through applications of the joint inversion of seismic data for lithologic and petrophysical parameters in several selected target areas, the resulted high-resolution lithologic and petrophysical sections(impedance, porosity, clay) show that the joint inversion can significantly improve the spatial description of reservoirs in data sets involving complex deposits. It proved the validity and practicality of this method adequately.

内蒙古中部道郎和都格铜多金属矿的地质特征矿床类型及找矿模式研究

Daolangheduge copper polymetallic deposit is located on east edge of Ondor Sum-Bainaimiao metallogenic belt, which is a prospective area of porphyry copper deposit, in Xianghuangqi of central Inner Mongolia. Geotectonically, it occurred in the continental margin accretion belt along the north margin of North China Plate, south of the suture zone between North China Plate and Siberian Plate. The intrusive rocks in this area mainly consist of intermediate-acid magmatic rocks, and the quartz veins, tourmaline veins and the transitional phase are comparatively developed. According to our research, the ore-bearing rock body is mainly quartz diorite while the surrounding rock is mainly biotite granite. Besides, the wall rock alteration are mainly propylitization, pyritization and silicification, which consist of epidotization, actinolitization, chloritzation and so on. The metallic minerals are mainly chalcopyrite and pyrite. In addition, the primary ore is mainly of quartz-chalcopyrite-pyrite type. Above all, Daolangheduge copper polymetallic deposit is suggested to be categorized in the porphyry copper type. With isotopic dating and geochemical research on quartz diorite of ore-bearing rock body, the zircon LA-ICP-MS U-Pb dating of two samples yields an age of 266±2 Ma, falling into the range of late Permian Epoch. It is the first accurate age data in Xianghuangqi area, so it should play a key role in the research of deposit and magmatic rocks in this area. With the major elements and trace elements analysis of 14 samples, the quartz diorite should be among the calc-alkaline series, the geochemical characteristics show higher large-ion lithophile elements of Rb, Sr and LREE, low high-field strength elements of Nb, Ta and high transition elements of Cu, Cr . Also, the REE patterns have negative Eu anomalies. With the same analysis of 4 sample for the biotite granite, the geochemical characteristics show higher Rb, Th,, Zr, Hf and LREE, low Nb, Sm and HREE and Eu has no anomaly. It should be among the calc-alkaline series, over aluminum quality and has characteristics of Adakites. According to isotopic dating and geochemical characteristics of ore-bearing rock body, it is suggested that its materials mainly derived from upper mantle that had fractional crystallization and its magma source region may be affected by fluid metasomatism of paleo-asian ocean. It should be an extensional process of post-orogeny according to regional tectonic evolution. Consequently, because of the decrease of temperature and pressure, the ore forming fluid was raised to surface and mineralized accompanied by magmatic activity which might occur in south of the suture zone. By geological survey, further geophysical and geochemical work is needed. In this area, we have accomplished high precision magnetic prospecting, high density electrical survey, gravity prospecting, soil geochemical prospecting, X-ray fluorescence analyzer prospecting and so on. According to geophysical and geochemical abnormal and surface occurrence, 11 drills are arranged to verification. The type of ores are mainly quartz-chalcopyrite-pyrite ores within 3 drills by drill core logging. Although the grade as well as the scale of already-found Cu deposits are insufficient for industrial exploitation, the mineralization prospect in this region is supposed to be great and the potential in mineral exploration at depth is excellent.

塔里木盆地塔中地区寒武-奥陶系烃源岩古地温及成烃演化史研究

The present maturity of Cambrian and Ordovician source rocks in Tazhong area, Tarim basin, is studied using several organic petrology methods and conodont CAI method. The highest palaeotemperature that the Cambrian-Ordovician undergone is revealed by Laser Roman Microprobe (LRM) analysis and by simulating experiment of the kerogen chemical kinetics. In according to all above study, the thermal history of Cambrian and Ordovician is reconstructed based on numerical simulating approaches. The characteristics of secondary hydrocarbon generating are studied by inclusions analysis. The reflectances of the samples in the drills located in Tazhong area show that the maturities of Cambrian source rocks are in the stages of condense oil-dry gas, and that of Ordovician source rocks range from peak of oil generating to wet gas stage. The palaeotemperature data of Cambrian-Orovician source rocks from well Tacan 1, based on LRM analysis, are in coincidence with that from other methods. Also are the palaeotemperature data of Cambrian-Orovician source rocks in well Tacan 1 based on the simulating experiments of kerogen pyrolysis, similar to the homogenization temperatures of inclusions in the source rocks. Aaccording to the vitrinite inflectance data of the TZ12 well and Tacan 1 well, the paleotemperature gradients are analysized and reconstructed. These data show that the paleotemperature gradient in Tazhong area was the highest during Cambrian-Ordovician period, it was up to 3.5°C/100m. Following, the temperature gradient descended gradually and it reached to the lowest at present (2.2°C/100m). The histories of maturation and hydrocarbon generation of Cambrian and Ordovician source rocks in Tazhong area are researched systematically and quantitatively, the results show that periods of oil generation from Cambrian and Ordovician source rocks lasted for a long time from Ordovician to Carbonferious periods because the central Cambrian stratum in the north slope of Tazhong area is buried differently in depth. The top of the Cambrian entered into the peak of oil generation in middle-late Ordovician, and most area of the north slope of Tazhong area entered into the peak of oil generation in Carbonferious period, and on the uplift belt some of source rocks entered into the peak of oil generation in Permian period. In early Devonian, the central of the Lower Ordovician source rocks near the Manjiaer depression reached the peak of oil generation and near the top of the Tazhong uplift did not reached the peak of oil generation until early Cretaceous. The middle-upper Ordovician entered into the peak of oil generation in early-middle Jurassic. The time of the middle-upper Ordovician in the top of the uplift belt entering into the peak of oil generation was delayed, because the source rock was buried shallowly, and it did not reached the peak of oil generation until middle Cretaceous. Middle-upper Ordovician in the top of the north slope has been in the peak of oil generation now, it is consistent with the maturity (1.0-1.2%Ro) of the source rocks. The characteristics of the inclusions formed by kerogens are different from that by crystal-enclosed organic matters(OM) during secondary hydrocarbon generation of Cambrian and Ordovician source rocks. The secondary hydrocarbon generation mainly occurred in Mesozoic-Cenozoic period, in an area of about 9000km2 in the north slope. The intensity of the secondary hydrocarbon generation of Cambrian and Ordovician is up to 21kg/torg and 36kg/torg) respectively. Using the staged gas chromatography, the high-over maturated carbonate source rocks are analysized to release the adsorbed OM, inclusions OM and crystal-enclosed OM, respectively, and to evaluate their relative contributions to secondary hydrocarbon generation. The three periods of oil and gas migration and petroleum pools formation in Tazhong area are determined according to organic inclusions and solid bitumen.

地质体科学计算可视化研究及其应用

Describing visually space-time properties of geological phenomena consists of one of the most important parts in geology research. Such visual images are of usually helpful for analyzing geological phenomena and for discovering the regulations behind geological phenomena. This report studies mainly three application problems of scientific visualization in geology: (Dvisualizing geological body A new geometric modeling technique with trimmed surface patches has been eveloped to visualize geological body. Constructional surfaces are represented as trimmed surfaces and a constructional solid is represented by the upper and lower surface composed of trimmed surface patches from constructional surfaces. The technique can completely and definitely represent the structure of geological body. It has been applied in visualization for the coal deposit in Huolinhe, the aquifer thermal energy storage in Tianjin and the structure of meteorite impact in Cangshan et al. (2)visualizing geological space field Efficient visualization methods have been discussed. Marching-Cube algorithm used has been improved and is used to extract iso~surface from 3D data set, iso-line from 2D data set and iso-point from ID data set. The improved method has been used to visualize distribution and evolution of the abnormal pressures in Zhungaer Basin. (3)visualizing porous space a novel way was proposed to define distance from any point to a convex set. Thus a convex set skeleton-based implicit surface modeling technique is developed and used to construct a simplified porous space model. A Buoyancy Percolation numerical simulation platform has been developed to simulate the process of migration of oil in the porous media saturated with water.