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.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.

A Generalized Soft-Sphere for Monte Carlo Simulation

A new collision model, called the generalized soft-sphere (GSS) model, is introduced. It has the same total cross section as the generalized hard-sphere model [Phys. Fluids A 5, 738 (1993)], whereas the deflection angle is calculated by the soft-sphere scattering model [Phys. Fluids A 3, 2459 (1991)]. In virtue of a two-term formula given to fit the numerical solutions of the collision integrals for the Lennard-Jones (6-12) potential and for the Stockmayer potential, the parameters involved in the GSS model are determined explicitly that may fully reproduce the transport coefficients under these potentials. Coefficients of viscosity, self-diffusion and diffusion for both polar and nonpolar molecules given by the GSS model and experiment are in excellent agreement over a wide range of temperature from low to high.

A Numerical Study of Indentation Using Indenters of Different Geometry

Finite element simulation of the Berkovich, Vickers, Knoop, and cone indenters was carried out for the indentation of elastic-plastic material. To fix the semiapex angle of the cone, several rules of equivalence were used and examined. Despite the asymmetry and differences in the stress and strain fields, it was established that for the Berkovich and Vickers indenters, the load-displacement relation can closely be simulated by a single cone indenter having a semiapex angle equal to 70.3degrees in accordance with the rule of the volume equivalence. On the other hand, none of the rules is applicable to the Knoop indenter owing to its great asymmetry. The finite element method developed here is also applicable to layered or gradient materials with slight modifications.