Two-Stream Approximation to the Radiative Transfer Equation:A New Improvement and Comparative Accuracy with Existing Methods

Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

Quadrature-based moment methods for the population balance equation: An algorithm review

The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source terms. The equation is therefore quite complex and difficult to solve analytically or numerically. The quadrature-based moment methods(QBMMs) are a class of methods that solve the PBE by converting the transport equation of the number density function(NDF) into moment transport equations. The unknown source terms are closed by numerical quadrature. Over the years, many QBMMs have been developed for different problems, such as the quadrature method of moments(QMOM), direct quadrature method of moments(DQMOM),extended quadrature method of moments(EQMOM), conditional quadrature method of moments(CQMOM),extended conditional quadrature method of moments(ECQMOM) and hyperbolic quadrature method of moments(Hy QMOM). In this paper, we present a comprehensive algorithm review of these QBMMs. The mathematical equations for spatially homogeneous systems with first-order point processes and second-order point processes are derived in detail. The algorithms are further extended to the inhomogeneous system for multiphase flows, in which the computational fluid dynamics(CFD) can be coupled with the PBE. The physical limitations and the challenging numerical problems of these QBMMs are discussed. Possible solutions are also summarized.

An analysis of the radar backscatter from oil-covered sea surfaces using moment method and Monte-Carlo simulation: preliminary results

An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face is numerically generated with an ocean waveheight spectrum for a given wind velocity. A two-layered medium is then generated by adding a thin oil layer on the simulated rough sea surface. The electric fields backscattered from the sea surface with two-layered medium are computed with the method of moments （MoM）, and the backscattering coefficients are statistically obtained with N independent samples for each oil-spilled surface using the Monte-Carlo technique for various conditions of surface roughness, oil-layer thickness, frequency, polarization and incidence angle. The numerical simulation results are compared with theoretical models for clean sea surfaces and SAR images of an off-spilled sea surface caused by the Hebei （Hebei province, China） Spirit oil tanker in 2007. Further, conditions for better oil spill extraction are sought by the numerical simulation on the effects of wind speed and oil-layer thickness at different inci- dence angles on the backscattering coefficients.

Quadrature Method of Moments for Nanoparticle Coagulation and Diffusion in the Planar Impinging Jet Flow

A computational model combining large .eddy simulation with quadrature moment method was em-ployed to study nanoparticle evolution in a confined impinging jet. The investigated particle size is limited in the transient regime, and the particle collision kernel was obtained by using the theory of flux matching. The simulation was validated by comparing it with the experimental results. The numerical results show coherent structure acts to dominate particle number intensity, size and polydispersity distributions, and it also induce particle-laden iet to be diluted by .the ambient.The evolution of particle dynarnics in.the impinging jet flow are strongly related to the Rey-nolds number and nozzle-to-plate distance, and their relationships were analyzed.

RANDOM SINGULAR INTEGRAL OF RANDOM PROCESS WITH SECOND ORDER MOMENT

This paper discussses the random singular integral of random process with second order moment, establishes the concepts of the random singular integral and proves that it's a linear bounded operator of space Ha(L)(m, s). Then Plemelj formula and some other properties for random singular integral are proved.

Effects of slip, viscous dissipation and Joule heating on the MHD flow and heat transfer of a second grade fluid past a radially stretching sheet

The flow and heat transfer of an electrically conducting non-Newtonian second grade fluid due to a radially stretching surface with partial slip is considered. The partial slip is controlled by a dimensionless slip factor, which varies between zero （total adhesion） and infinity （full slip）. Suitable similarity transformations are used to reduce the resulting highly nonlinear partial differential equations into ordinary differential equations. The issue of paucity of boundary conditions is addressed and an effective numerical scheme is adopted to solve the obtained differential equations even without augmenting any extra boundary conditions. The important findings in this communication are the combined effects of the partial slip, magnetic interaction parameter and the second grade fluid parameter on the velocity and temperature fields. It is interesting to find that the slip increases the momentum and thermal boundary layer thickness. As the slip increases in magnitude, permitting more fluid to slip past the sheet, the skin friction coefficient decreases in magnitude and approaches zero for higher values of the slip parameter, i.e., the fluid behaves as though it were inviscid. The presence of a magnetic field has also substantial effects on velocity and temperature fields.

A NEW MOMENT METHOD FOR THE FAST AND ACCURATE ANALYSIS OF NORMAL MODE HELICAL ANTENNAS

In this letter, a new moment method using helical segments is presented to model Normal Mode Helical Antenna (NMHA). Using this method, the NMHA can be modeled by a few segments. The current distributions and radiation patterns of some NMHAs are calculated.A comparison is made between results obtained using this helical segment algorithm and a linear segment algorithm, and the results of the two algorithms agree fairly well. When calculating the impedance matrix [Z], all the elements of the matrix can be obtained by only calculating a few elements with the application of the symmetric and periodic characteristics of the NMHA.Therefore, the CPU time and the memory storage are significantly reduced, with the accuracy and speed enhanced.

Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

APPLICATION OF THE MOMENT METHODS TO ANALYSIS OF X-RAY DIFFRACTION LINE PROFILE FOR PA1010-BMI SYSTEM

Pure X-ray diffraction profiles have been analysed for polyamide 1010 and PA1O1O-BMI system by means of multipeak fitting resolution of X-ray diffraction. The methods of variance and fourth moment have been applied to determine the particle size and strain values for the paracrystalline materials. The results indicated that both variance and fourth moment of X-ray diffraction line profile yielded approximately the same values of the particle size and the strain. The particle sizes of (100) reflection have been found to decrease with increasing BMI content, whereas the strain values increased.

Structure and Performance of Fe-N_x-C Catalyst for Oxygen Reduction Reaction Prepared by Vacuum Casting Method and the Second Pyrolysis

Affordable non-precious metal（NPM） catalysts played a vital role in the wide application of polymer electrolyte membrane fuel cells（PEMFC）. In current work, a facile vacuum casting reacting method based on vacuum casting was introduced to prepare Fe-N_x-C oxygen reduction reaction（ORR） catalysts with high efficient in acid medium. The catalysts were prepared with ammonium ferrous sulfate hexahydrate（AFS） and 1,10-phenanthroline monohydrate utilizing homemade mesoporous silica template. The heat treatment and its influence on structure and performance were systematically evaluated to achieve superior ORR performance and some clues were found. And 850 ℃ was found to be the best temperature for the first and second pyrolysis. The linear sweep voltammetry（LSV） results showed that there were only 18 mV slightly negative shifts of half-wave potential（E_（1/2）） of the optimal catalyst（749 mV） compared with the commercial Pt/C（20 μg·Pt·cm^-2）. Besides, I850 R also showed better electrochemical stability and methanol-tolerance than that of Pt/C. All evidences proved that our vacuum casting reacting strategy and heat treatment process were prospective for the future R＆D of high performance Fe-N_x-C ORR catalysts.

A novel method for integrating chromatographic fingerprint analytical units of Chinese materia medica:the matching frequency statistical moment method

Objective To facilitate the quality evaluation suitable for the unique characteristics of Chinese materia medica(CMM)by developing and implementing a novel approach known as the matching frequency statistical moment(MFSM)method.Methods This study established the MFSM method.To demonstrate its effectiveness,we applied this novel approach to analyze Danxi Granules(丹膝颗粒,DXG)and its constituent herbal materials.To begin with,the ultra-performance liquid chromatography(UPLC)was applied to obtain the chromatographic fingerprints of DXG and its constituent herbal materi-als.Next,the MFSM was leveraged to compress and integrate them into a new fingerprint with fewer analytical units.Then,we characterized the properties and variability of both the original and integrated fingerprints by calculating total quantum statistical moment(TQSM)parameters,information entropy and information amount,along with their relative standard deviation(RSD).Finally,we compared the TQSM parameters,information entropy and infor-mation amount,and their RSD between the traditional and novel fingerprints to validate the new analytical method.Results The chromatographic peaks of DXG and its 12 raw herbal materials were divided and integrated into peak families by the MFSM method.Before integration,the ranges of the peak number,three TQSM parameters,information entropy and information amount for each peak or peak family of UPLC fingerprints of DXG and its 12 raw herbal materials were 95.07−209.73,9390−183064μv·s,5.928−21.33 min,22.62−106.69 min^(2),4.230−6.539,and 50530−974186μv·s,respectively.After integration,the ranges of these parameters were 10.00−88.00,9390−183064μv·s,5.951−22.02 min,22.27−104.73 min^(2),2.223−5.277,and 38159−807200μv·s,respectively.Correspondingly,the RSD of all the aforementioned pa-rameters before integration were 2.12%−9.15%,6.04%−49.78%,1.15%−23.10%,3.97%−25.79%,1.49%−19.86%,and 6.64%−51.20%,respectively.However,after integration,they changed to 0.00%,6.04%−49.87%,1.73%−23.02%,3.84%−26.85%,1.17%−16.54%,and 6.40%−48.59%,respectively.The results demonstrated that in the newly integrated fingerprint,the analytical units of constituent herbal materials,information entropy and information amount were significantly reduced(P<0.05),while the TQSM parameters remained unchanged(P>0.05).Additionally,the RSD of the TQSM parameters,information entropy,and information amount didn’t show significant difference before and after integration(P>0.05),but the RSD of the number and area of the integrated analytical units significantly decreased(P<0.05).Conclusion The MFSM method could reduce the analytical units of constituent herbal mate-rials while maintain the properties and variability from their original fingerprint.Thus,it could serve as a feasible and reliable tool to reduce difficulties in analyzing multi-compo-nents within CMMs and facilitating the evaluation of their quality.

The b value spectrum and η value by the moment method

The earthquake size distribution is generally considered to obey the Gutenberg Richter (GR) law. We have introduced the concept of the b value spectrum based on the moment method to investigate the deviation of the actual magnitude distribution of earthquakes from this law. This enables us to describe characteristic features of the magnitude frequency distribution of earthquakes. We found also a simple relation between the η value and the b value spectrum. Analysis using this scheme showed that the actual size distributions of earthquakes have large variations from case to case and sometimes deviate considerably from the widely assumed the GR formula.

QM/MM Study of the Second Harmonic Generation and Two-photon Absorption Properties of a Fluorescent Proteins-Dreiklang

A new reversibly switchable fluorescent protein（RSFP）, namely Dreiklang, exhibits prominent feature that the wavelengths for switching and fluorescence are decoupled due to its great different structures between bright and dark states. This feature might also induce some nonlinear optic（NLO） properties changing as switching between two states, which might promote new method of biological science. We employ the QM/MM method to simulate the structures of different states, and study their second harmonic generation（SHG） and two-photon absorption（TPA） properties. And we found different states of Dreiklang have different SHG and TPA responses. The SHG and TPA properties of Dreiklang are correlated to particularly geometrical structures of different states, especially the centrosymmetric or nocentrosymmetric π-stacking structures which are formed by chromophore and beside residue Tyr203, so the SHG and TPA responses could be changed as the light induces switching among different states of Dreiklang. This work would prospectively guide the application of Dreiklang on the NLO technology, and help the development of new RSFP with special NLO function.

Microscopic Mechanism of the ω Variation in Moments of Inertia for the Yrast Superdeformed Bands 194Tl(1a, 1b)

The variation in moments of inertia (J(1) and J(2)) with rotational frequency for the superdeformed bandsin odd-odd nuclei, 194Tl(la,lb), is investigated by using the particle-number conserving method for treating the pairinginteraction (monopole and quadrupole). The observed variations of J(1) and J(2) with ω are reproduced quite well inthe calculation and the contributions from each major shell are clearly displayed.

Explicit Iterative Methods of Second Order and Approximate Inverse Preconditioners for Solving Complex Computational Problems

Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is discussed. The second order iterative methods behave quite similar to first order methods and the development of efficient preconditioners for solving the original linear system is a decisive factor for making the second order iterative methods superior to the first order iterative methods. Adaptive preconditioned Conjugate Gradient methods using explicit approximate preconditioners for solving efficiently large sparse systems of algebraic equations are also presented. The generalized Approximate Inverse Matrix techniques can be efficiently used in conjunction with explicit iterative schemes leading to effective composite semi-direct solution methods for solving large linear systems of algebraic equations.

Theoretical study of the hyperfine interaction constants, Landé g-factors, and electric quadrupole moments for the low-lying states of the 61Niq+(q=11,12,14,and 15) ions

Highly charged nickel ions have been suggested as candidates for the ultra-precise optical clock, meanwhile the relevant experiment has been carried out. In the framework of the multiconfiguration Dirac–Hartree–Fock(MCDHF)method, we calculated the hyperfine interaction constants, the Landég-factors, and the electric quadrupole moments for the low-lying states in the 61Ni11+,61Ni12+,61Ni14+, and61Ni15+ ions. These states are clock states of the selected clock transitions in highly charged nickel ions(see Fig. 1). Based on discussing the effects of the electron correlations, the Breit interaction, and quantum electrodynamics(QED) effect on these physical quantities, reasonable uncertainties were obtained for our calculated results. In addition, the electric quadrupole frequency shifts and the Zeeman frequency shifts of the clock transitions concerned were analyzed.

Research on the Related Factors of the Second Molar Dislocation and Orthodontic Erect Method

The second molar dislocation is more common clinically.To investigate the related factors of the second permanent molar dislocation,and provide reference for the clinical diagnosis and treatment of orthodontics.From the current clinical research,the clinical methods of orthodontic erect secondary molars are also diverse and clinical.The narrower first molar alveolar arch width,smaller ANB angle,and crowded maxillary posterior segment arch are the factors that cause the maxillary second permanent molar dislocation.The narrow alveolar arch width,the smaller SNB angle,the larger ANB angle,and the crowded lower mandibular arch are the factors leading to the dislocation of the mandibular second permanent molar.In addition,for the second mandibular molar malposition,it is particularly important to select the corrective treatment plan.It is especially important to improve the treatment.

SECOND-MOMENT CLOSURE FOR MODELLING THE NEAR-WALL TURBULENCE IN 3D MEAN FLOWS

A second-moment closure for the near-wall turbulence is proposed. The limiting behaviour of this closure near a wall is consistent with that of the exact Reynolds-stress transport equations, and it converts asymptotically into a high- Reynolds-number closure remote from the wall. The closure is applied to a pressure- driven 3D transient channel flow. The predicted results are in fair agreement with the DNS data.

TDDFT-SOS Studies on the Polarizabilities and Second-order Hyperpolarizabilities of Zn_3O_3 II-VI Semiconductor Clusters

The geometrical structure of semiconductor clusters including Zn3O3 was optimized by the DFT B3LYP method. With the same basis sets, dipole moments, polarizabilities and secondorder hyperpolarizabilities have been calculated and compared with the results obtained by TDDFT B3LYP method combined with sum-over-state （SOS） formula. The calculation results indicate that the dipole moments of the ground state depend on the atom radius and electronegative differences between elements and are their balance point as well. The polarizabilities of the clusters accord with the rule of the corresponding energy transformation from ground to excited state. The results predict an increase of second-order hyperpolarizabilities with increasing the distances between atoms in the clusters as well as a decrease of the polarizabilities and second-order hyperpolarizabilities in the same serial of semiconductor clusters with increasing the dipole moments of the ground states. The changes of dipole moments in ground states are inconsistent with transition moments. Spatial structure, charge transfer and other factors play an important role in composing the transition moments.