Cone Effect in Astronomical Adaptive Optics System Investigated by a Pure Numerical Simulation

It is well-known that cone effect or focus anisoplanatism is produced by the limited distance of a laser guide star (LGS) which is created within the Earth atmosphere and consequently located at a finite distance from the observer. In this paper, the cone effect of the LGS for different vertical profiles of the refractive index structure constant Cn2 is numerically investigated by using a revised computer program of atmospheric propagation of optical wave and an adaptive optics (AO) system including dynamic control process. According to the practice, the overall tilt for the tilt-correction mirror is obtained from a natural star and the aberrated wavefront for phase correction of the deformable mirror is obtained from a LGS in our numerical simulation. It is surprisingly found that the effect of altitude of the LGS on the AO phase compensation effectiveness by using the commonly-available vertical profiles of Cn2 and the lateral wind speed in the atmosphere is relatively weak, and the cone effect for some Cn2 profiles is even negligible. It is found that the cone effect does not have obvious relationship with the turbulence strength, however, it depends on the vertical distribution profile of Cn 2 apparently. On the other hand, the cone effect depends on the vertical distribution of the lateral wind speed as well. In comparison to a longer wavelength, the cone effect becomes more obvious in the case of a shorter wavelength. In all cases concerned in this paper, an AO system by using a sodium guide star has almost same phase compensation effectiveness as that by using the astronomical target itself as a beacon. Effect of dynamic control process in an AO system on the cone effect is studied in this paper for the first time within our knowledge.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

Coalition and distribution in fuzzy dynamic supply chain systems

In this paper, two models of coalition and income's distribution in FSCS (fuzzy supply chain systems) are proposed based on the fuzzy set theory and fuzzy cooperative game theory. The fuzzy dynamic coalition choice's recursive equations are constructed in terms of sup-t composition of fuzzy relations, where t is a triangular norm. The existence of the fuzzy relations in FSCS is also proved. On the other hand, the approaches to ascertain the fuzzy coalition through the choice's recursive equations and distribute the fuzzy income in FSCS by the fuzzy Shapley values are also given. These models are discussed in two parts: the fuzzy dynamic coalition choice of different units in FSCS; the fuzzy income's distribution model among different participators in the same coalition. Furthermore, numerical examples are given aiming at illustrating these models., and the results show that these models are feasible and validity in FSCS.

Compatibility relationship among reaction equilibria, equivalence principle of reaction approaches, and silicon contamination in semiconductors

This paper discovers some shortcomings in the algorithm for the incorporation of Si into GaAs in the GaAs VPE process. These faults arise from neglecting a link, the compatibility relationship, in chemical thermodynamics. The meaning of said relationship is as follows: In an equilibrium complex system, each species can only contribute one and the same quantity (its equilibrium quantity) to the different equilibria of the various reactions involving it; yet even under this restriction, every equilibrium constant is satisfied, and all the reaction equilibria coexist compatibly in the system. Only by adding the relationship can the equilibrium theory for the complex system be complete. This paper also tells its position in chemical thermodynamics. Such a compatibility concept directly leads to an equivalence principle: In a complex system, a certain species can usually be simultaneously formed by many chemical reactions; when the system has reached equilibrium under fixed environmental conditions, the equilibrium quantity of said species calculated according to each chemical equation of these reactions will be equal and the various reaction approaches will be equivalent, provided that for all the reactants and all the other products of these reactions their equilibrium quantities in the system are respectively taken as corresponding knowns for the calculations, which is extremely useful for seeking a functional relation among the species' equilibrium quantities in a system (Si contamination is one of the examples). Under the guidance of those arguments, the various schools' algorithms for the Si contamination can be uniformized and simplified, and the contamination quantity relation between Si and O, two very important impurities, is found.

Numerical study on onset of oscillatory thermocapillary flow in rectangular liquid pool

Thermocapillary flow in a rectangular liquid pool of large Prandtl fluid (Pr = 105.6) is numerically studied in microgravity. Oscillatory thermocapillary flow arises when the imposed temperature difference between the sidewalls exceeds a critical value. The fluctuations of the oscillatory flow, accompanied by the propagation of the hydrothermal wave from the cold sidewall to the hot one, are much smaller than the time-averaged velocity and temperature fields. The corresponding disturbance cells arise in the centre of the liquid pool initially, and extend to the whole region with the increasing imposed temperature difference. The present study reveals the different characteristics of the oscillatory themocapillary flow in the rectangular liquid pool as compared to the cases in other configurations.

Numerical Study on the Spatial and Temporal Characteristics of Water Quality in the Taihu Lake of China

Taihu Lake is the third largest fresh water lake in China. With the fast economic development, abundant industrial and agricultural waste water has been discharged into Taihu Lake, causing the eutrophication of the water quality, which greatly affected the water utility. In the past decades, the treatment of Taihu Lake has witnessed limited success. Therefore, it is practically and theoretically significant to study the eutrophication of Taihu Lake. This research has focused on the issue of water quality including the characteristics of spatial and temporal distributions, and the rules of nutrient diffusion in the Taihu lake area. Based on the monitoring data, the basis distribution characteristics of water quality in Taihu Lake are analyzed. Comparing Taihu Lake with other Lakes shows that one important reason for Taihu eutrophication is the long period of water retention. A transporting and diffusing model of Taihu nutrient is developed by combining with the hydrodynamics model. Using the model, the concentration field of the total phosphorus (TP) and the influence of wind-driven current are numerically investigated, which leads to the conclusion that the flow field has a great influence on the spatial and temporal distributions of TP in Taihu Lake. Furthermore, the effect for improving the water quality by the project of water diversion from the Yangtze River to Taihu Lake was analyzed by simulation. The results demonstrate that short-term water diversion cannot improve the water quality of the heavily-polluted Meiliang Bay and the western bank areas of Taihu Lake.

Numerical simulation of thermal loading produced by shaped high power laser onto engine parts

Recently a new method for simulating the thermal loading on pistons of diesel engines was reported. The spatially shaped high power laser is employed as the heat source, and some preliminary experimental and numerical work was carried out. In this paper, a further effort was made to extend this simulation method to some other important engine parts such as cylinder heads. The incident Gaussian beam was transformed into concentric multi-circular patterns of specific intensity distributions, with the aid of diffractive optical elements (DOEs). By incorporating the appropriate repetitive laser pulses, the designed transient temperature fields and thermal loadings in the engine parts could be simulated. Thermal-structural numerical models for pistons and cylinder heads were built to predict the transient temperature and thermal stress. The models were also employed to find the optimal intensity distributions of the transformed laser beam that could produce the target transient temperature fields. Comparison of experimental and numerical results demonstrated that this systematic approach is effective in simulating the thermal loading on the engine parts. (C) 2009 Elsevier Ltd. All rights reserved.

COMPARATIVE-STUDY OF KALMAN FILTERING, SYNCHRONOUS EXCITATION AND NUMERICAL DERIVATIVE TECHNIQUES IN FLUOROMETRY

The performance of Kalman filtering, synchronous excitation and numerical derivative techniques for the resolution of overlapping emission spectra in spectrofluorimetry was studied. The extent of spectrum overlap was quantitatively described by the separation degree D(s), defined as the ratio of the peak separation to the full width at half-maximum of the emission spectrum of the interferent. For the system of Rhodamine B and Rhodamine 6G with a large D(s) of about 0.4, both Kalman filtering and synchronous techniques are able to resolve the overlapping spectra well and to give satisfactory results while the derivative spectra are still overlapped with each other. Moreover, the sensitivities are greatly decreased in derivative techniques. For more closely spaced spectra emitted by the complexes of Al and Zn with 7-iodo-8-hydroxyquinoline-5-sulphonic acid (Ferron)-hexadecyltrimethylammonium bromide, the synchronous excitation technique cannot completely separate the overlapping peaks, although it increases the separation degree from 0.25 in the conventional spectra to 0.37 in the synchronous spectra. On the other hand, Kalman filtering is capable of resolving this system. When the Al/Zn intensity ratio at the central wavelength of Al was > 1, however, the accuracy and precision of the estimates for Zn concentration produced by the Kalman filter became worse. In this event, the combination of synchronous excitation and Kalman filtering can much improve the analytical results.

IMPROVEMENT OF SELECTIVITY WITH NUMERICAL DERIVATIVE TECHNIQUES IN INDUCTIVELY COUPLED PLASMA ATOMIC EMISSION-SPECTROMETRY

One of the most attractive features of derivative spectrometry is its higher resolving power. In the present power, numerical derivative techniques are evaluated from the viewpoint of increase in selectivity, the latter being expressed in terms of the interferent equivalent concentration (IEC). Typical spectral interferences are covered, including flat background, sloped background, simple curved background and various types of line overlap with different overlapping degrees, which were defined as the ratio of the net interfering signal at the analysis wavelength to the peak signal of the interfering line. the IECs in the derivative spectra are decreased by one to two order of magnitudes compared to those in the original spectra, and in the most cases, assume values below the conventional detection limits. The overlapping degree is the dominant factor that determines whether an analysis line can be resolved from an interfering line with the derivative techniques. Generally, the second derivative technique is effective only for line overlap with an overlapping degree of less than 0.8. The effects of other factors such as line shape, data smoothing, step size and the intensity ratio of analyte to interferent on the performance of the derivative techniques are also discussed. All results are illustrated with practical examples.

Shallow water tidal constituents in the Bohai Sea and the Yellow Sea from a numerical adjoint model with TOPEX/POSEIDON altimeter data

A numerical adjoint model with TOPEX/POSEIDON (T/P) altimeter data was set up to investigate the shallow water tidal constituents in the Bohai Sea and the Yellow Sea. Shallow water tidal constituents W-4, MS4 and M-6) in the Bohai Sea and the Yellow Sea were first extracted from nearly 10 years of T/P data and then assimilated into a nonlinear barotropic tidal model by using adjoint method in order to fully describe the tides in this area. The general patterns of M-4 and MS4 solutions were in good agreement with those of Kang et al. (Cont. Shelf. Res. IS (1998) 739.) and Lefevre et al., (J. Geophys. Res. 105 (2000) 8707.). The RMS values for the principal constituents and coastal constituents were obviously less than those calculated by Kang et al. (1998) and Lefevre et al. (2000). It was shown that the calculated tidal constituents charts obtained in the present study were more accurate than those in other models. In the future the model will be applied to other coastal areas and some semi-enclosed seas. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical simulation of principal tidal constituents in the South China Sea, Gulf of Tonkin and Gulf of Thailand

The principal tidal constituents M-2, S-2, K-1 and O-1 in the South China Sea, Gulf of Tonkin and Gulf of Thailand are simulated simultaneously using the numerical scheme of Kwok et al. (1995 Proceedings of the 1st Asian Computational Fluid Dynamics Conference, pp. 16-19). The average differences between the computed and observed harmonic constants are mostly within 5 cm and 10 degrees for amplitudes and phase-lags, respectively. The simulated tidal regimes in the present model are believed to be more accurate than the previous numerical results. Our studies confirm that a clockwise rotating M-2 amphidromic system lies in the southeast of the Gulf of Thailand and an S-2 amphidromic system at the near-shore area of the northeast South China Sea. The linear tidal energy equation developed by Garrett (1975 Deep-Sea Research 22, 23-35) is generalized to the nonlinear case. Based on the numerical results, the energy budgets in the South China Sea and its subareas, namely the Taiwan Strait, the Gulf of Tonkin, the Gulf of Thailand and the remaining area are investigated. The tidal motion in the Taiwan Strait is maintained mainly by the energy fluxes from the East China Sea for both semidiurnal and diurnal species and partially from the Luzon Strait for semidiurnal species. For the other parts of the South China Sea, the tidal motion is mainly maintained by the energy fluxes through the Luzon Strait. The energy inputs from the tide-generating force are negative for semidiurnal species and positive for diurnal species. (C) 1999 Elsevier Science Ltd. All rights reserved.

Numerical investigation of nutrient limitations in the Bohai Sea

River discharges are the important freshwater and nutrient sources for Bohai Sea (BS), and have a profound impact on the local marine environment. In this paper, the annual cycles of nutrient and phytoplankton dynamics in 1980s were reproduced using a coupled biogeochemical-physical model. Based on the validated simulations, the nutrient limitation characters were further investigated by running the model with the riverine nutrient altered, first enriching nitrogen and then phosphorus. It was found that although the riverine N:P ratios in Yellow and Haihe Rivers were much higher than the Redfield number, the nitrogen enrichment was still able to enhance the algae bloom in Laizhou and Bohai Bays. On the other hand, the response of algae growth to phosphorus enrichment was not thus obvious, which suggests that the local phytoplankton dynamics was characterized by the nitrogen limitation. Simulations also show that the nitrogen enrichment is generally accompanied by the phosphorus consumption, so a shift from nitrogen limitation to phosphorus limitation may occur if such a trend continues. (C) 2010 Elsevier Ltd. All rights reserved.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

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

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.