Numerical Simulation of Reactive Extrusion Processes for Activated Anionic Polymerization

In order to deal with the complicated relationships among the variables of the reactive extrusion process for activated anionic polymerization, a three-dimensional equivalent model of closely intermeshing co-rotating twin screw extruders was established. Then the numerical computation expressions of the monomer concentration, the monomer conversion, the average molecular weight and the fluid viscosity were deduced, and the numerical simulation of the reactive extrusion process of Styrene was carried out. At last, our simulated results were compared with Michaeli's simulated results and experimental results. (C) 2007 Elsevier B.V. All rights reserved

Numerical simulation of coil-helix transition processes of gelatin

Gelatin is widely used in food, pharmaceutical, and photographic industries due to the coil-helix transition, whereas the structural inhomogeneity considerably affects its essential properties closely connecting with the industrial applications. The spatially structural inhomogeneity of the gelatin caused by the uneven and unstable temperature field is analyzed by the finite element method during the cooling-induced coil-helix transition process. The helix conversion and the crosslinking density as functions of time and spatial grid are calculated by the incremental method. A length distribution density function is introduced to describe the continuous length distributions of two kinds of triple helices.

Numerical simulation of reactive extrusion processes of PA6

To analyze the complicated relationships among the variables during the reactive extrusion process of polyamide 6 (PA6), and then control the chemical reaction and the material structures, the process of continuous polymerization of caprolactam into PA6 in a closely intermeshing co-rotating twin screw extruder was simulated by means of the finite volume method, and the influences of three key processing parameters on the reactive extrusion process were discussed. The simulated results of an example were in good agreement with the experimental results.

Numerical simulation of mesoscopic mechanical behaviors of gradual multi-fiber-reinforced polymer matrix composites

To simulate the deformation and the fracture of gradual multi-fiber-reinforced matrix composites, a numerical simulation method for the mesoscopic mechanical behaviors was developed on the basis of the finite element and the Monte Carlo methods. The results indicate that the strength of a composite increases if the variability of statistical fiber strengths is decreased.

On the relationships between the radial tidal current field and the radial sand ridges in the southern Yellow Sea: a numerical simulation

The ori-in of the radial sand ridges (RSRs) in the southern Yellow Sea has been a controversial problem since they were discovered in the early 1960s. To resolve the problem, two key questions need to be answered: (1) was the radial tidal current field in the RSR area generated by the submarine topography, or (2) did it exist before the RSRs occurred? In this study, the M-2 tide and tidal currents in the RSR area were simulated with a two-dimensional tidal model using a flat bottom and a shelving slope topography, the results being then compared with the field data. It is demonstrated that the radial tidal current field in the southern Yellow Sea is independent of bottom topography, and may thus be the controlling factor generating the RSRs. The radial tidal current field probably existed before the RSRs were formed.

Numerical simulation of tides and tidal currents in Liaodong Bay with POM

A three-dimensional ocean circulation model, called Princeton Ocean Model (POM), is employed to simulate tides and tidal currents in Liaodong Bay. The nested grid technique is adopted to improve the computation precision. Computed harmonic constants of M-1, M-2 tides at five tidal gauge stations and surface elevations at two oil platforms are compared with those observed, and show good agreements with them. Based on the calculated results, the co-amplitude and co-phase tag lines of nil and M-2 tidal constituents, the residual current field of M-2 constituent, tidal form, tidal Current ellipse and the moving style of tidal current are given. It is found that diurnal tidal constituents have no amphidromic point whereas semi-diurnal constituents have one in the region of interest. Meanwhile, some meaningful results are concluded and presented, which are conducive to a thorough knowledge of the characteristics of tides and tidal currents in the Liaodong Bay.

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 Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the source region, the seafloor lifts to a designated height by a generation function. The numerical tests show that file linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations eta(max)(0) are carried out by both the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on eta(max)(0) are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences between the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the eta(max)(0) near-linearly varies with the wave amplitudes of the surface waves, and the eta(max)(0) has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, aid these differences are significantly affected by The wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.

Numerical simulation and validation of ocean waves measured by an Along-Track Interferometric Synthetic Aperture Radar

A new nonlinear integral transform of ocean wave spectra into Along-Track Interferometric Synthetic Aperture Radar (ATI-SAR) image spectra is described. ATI-SAR phase image spectra are calculated for various sea states and radar configurations based on the nonlinear integral transform. The numerical simulations show that the slant range to velocity ratio (R/V), significant wave height to ocean wavelength ratio (H-s/lambda), the baseline (2B) and incident angle (theta) affect ATI-SAR imaging. The ATI-SAR imaging theory is validated by means of Two X-band, HH-polarized ATI-SAR phase images of ocean waves and eight C-band, HH-polarized ATI-SAR phase image spectra of ocean waves. It is shown that ATI-SAR phase image spectra are in agreement with those calculated by forward mapping in situ directional wave spectra collected simultaneously with available ATI-SAR observations. ATI-SAR spectral correlation coefficients between observed and simulated are greater than 0.6 and are not sensitive to the degree of nonlinearity. However, the ATI-SAR phase image spectral turns towards the range direction, even if the real ocean wave direction is 30 degrees. It is also shown that the ATI-SAR imaging mechanism is significantly affected by the degree of velocity bunching nonlinearity, especially for high values of R/V and H-s/lambda.

Numerical simulation on the process of saltwater intrusion and its impact on the suspended sediment concentration in the Changjiang (Yangtze) estuary

To study the relationship between sediment transportation and saltwater intrusion in the Changjiang (Yangtze) estuary, a three-dimensional numerical model for temperature, salinity, velocity field, and suspended sediment concentration was established based on the ECOMSED model. Using this model, sediment transportation in the flood season of 2005 was simulated for the Changjiang estuary. A comparison between simulated results and observation data for the tidal level, flow velocity and direction, salinity and suspended sediment concentration indicated that they were consistent in overall. Based on model verification, the simulation of saltwater intrusion and its effect on sediment in the Changjiang estuary was analyzed in detail. The saltwater intrusion in the estuary including the formation, evolution, and disappearance of saltwater wedge and the induced vertical circulation were reproduced, and the crucial impact of the wedge on cohesive and non-cohesive suspended sediment distribution and transportation were successfully simulated. The result shows that near the salinity front, the simulated concentrations of both cohesive and non-cohesive suspended sediment at the surface layer had a strong relationship with the simulated velocity, especially when considering a 1-hour lag. However, in the bottom layer, there was no obvious correlation between them, because the saltwater wedge and its inducing vertical circulation may have resuspended loose sediment on the bed, thus forming a high-concentration area near the bottom even if the velocity near the bottom was very low during the transition phase from flood to ebb.

Electro-Osmotic Flow and Mixing in Heterogeneous Microchannels

Analytical and numerical studies of secondary electro-osmotic flow EOF and its mixing in microchannels with heterogeneous zeta potentials are carried out in the present work. The secondary EOFs are analyzed by solving the Stokes equation with heterogeneous slip velocity boundary conditions. The analytical results obtained are compared with the direct numerical simulation of the Navier-Stokes equations. The secondary EOFs could transport scalar in larger areas and increase the scalar gradients, which significantly improve the mixing rate of scalars. It is shown that the heterogeneous zeta potentials could generate complex flow patterns and be used to enhance scalar mixing.

定常压力梯度边界层相干结构的直接数值模拟

根据流动稳定性理论，将不稳定波的一个周期作为相干结构的初值，采用直接数值模拟方法对有压力梯度湍流边界层中相干结构的演化进行了研究，得出其各种特性的变化与实验凤测到的结果一致。

Mapping Closure Approximation for Reactive Scalars in Random Flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows.

Coherent Structures in Transition of a Flat-Plate Boundary Layer at Ma=0.7

Direct numerical simulation (DNS) of a spatially evolving flat-plate boundary layer transition process at free stream Mach number 0.7 is performed. Tollmien-Schlichting (T-S) waves are added on the inlet boundary as the disturbances before transition. Typical coherent structures in the transition process are investigated based on the second invariant of velocity gradient tensor. The instantaneous shear stress and the mean velocity profile in the transition region are studied. In our view, the fact that the peak value of shear stress in the stress concentration area increases and exceeds a threshold value during the later stage of the transition process plays an important role in the laminar breakdown process.