Complete solution of the columbus problem for large-disturbance cases

The Columbus problem has been rigorously solved by Lyapunov's direct approach to the continuous system in gencral cases of large disturbance and the theory has proved to be in strict consistency with Kelvin's experiments.

The non-axisymmetrical configuration of a large-scale solar and stellar magnetic field

in the corona, consisting of an eruptive prominence and/or a magnetic flux region (loop or arcade, or blob) in front of the prominence. Ahead of the piston, there is a compressed flow, which produces a shock front. This high-density region corresponds to the bright feature of the transient. Behind the piston, there is a rarefaction region, which corresponds to the dark feature of the transient. Therefore, both the bright and dark features of the transient may be explained at the same time by the dynamical process of the moving piston.

Control of Large-size SiC Growth Process Using Modeling tool

Large size bulk silicon carbide (SiC) crystals are commonly grown by the physical vapor transport (PVT) method. The PVT growth of SiC crystals involves sublimation and condensation, chemical reactions, stoichiometry, mass transport, induced thermal stress, as well as defect and micropipes generation and propagation. The quality and polytype of as-grown SiC crystals are related to the temperature distribution inside the growth chamber during the growth process, it is critical to predict the temperature distribution from the measured temperatures outside the crucible by pyrometers. A radio-frequency induction-heating furnace was used for the growth of large-size SiC crystals by the PVT method in the present study. Modeling and simulation have been used to develop the SiC growth process and to improve the SiC crystal quality. Parameters such as the temperature measured at the top of crucible, temperature measured at the bottom of the crucible, and inert gas pressure are used to control the SiC growth process. By measuring the temperatures at the top and bottom of the crucible, the temperatures inside the crucible were predicted with the help of modeling tool. SiC crystals of 6H polytype were obtained and characterized by the Raman scattering spectroscopy and SEM, and crystals of few millimeter size grown inside the crucible were found without micropipes. Expansion of the crystals were also performed with the help of modeling and simulation.

Research on Bearing Capacity Calculation Method of Large Diameter Cylinder in Soft Clay

A large diameter cylinder inserted in soils is a new type of engineering structures used in offshore and port engineering. The mechanism of its bearing capacity and the analysis of its stability are important to its design and applications. In this paper, the finite element method is used to analyze the reacting forces of the soft soil foundation on the structure under the wave action. A simplified method is proposed, based on the plastic limit method, for the safety and stability analysis. Our analysis shows that the assumptions made in this paper and the mechanism used are reasonable, and the results obtained are appropriate. The calculation method is very efficient and can be used to evaluate main parameters of the structure in its preliminary designs.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Parametric studies of tension leg platform with large amplitude motions

Tension leg platform (TLP) is an important kind of working station for deep water exploration and development in ocean, whose dynamic responses deserve a serious thought. It is shown that for severe sea state, the effects of nonlinearities induced by large displacements of TLP may be noteworthy, and then employment of small displacements model should be restrained. In such situation, large amplitude motion model may be an appropriate alternative. The numerical experiments are performed to study the differences of dynamic responses between the two models. It is shown that for most cases, differences between results of the two models are significant. The variances of the differences vs. the wave period are the most remarkable, and that of the differences vs. wave heading angle are also apparent.

An extension of the two dimensional JKR theory to the case with a large contact width

In the Hertz and JKR theories, parabolic assumptions for the rounded profiles of the sphere or cylinder are adopted under the condition that the contact radius (width) should be very small compared to the radius of the sphere or cylinder. However, a large contact radius (width) is often found in experiments even under a zero external loading. We aim at extending the plane strain JKR theory to the case with a large contact width. The relation between the external loading and the contact width is given. Solutions for the Hertz, JKR and rounded-profile cases are compared and analyzed. It is found that when the ratio of a/R is approximately larger than about 0.4, the parabolic assumptions in the Hertz and JKR theories are no longer valid and the exact rounded profile function should be used.

A kinematic subgrid scale model for large-eddy simulation of turbulence generated sound

In the hybrid approach of large-eddy simulation (LES) and Lighthill’s acoustic analogy for turbulence-generated sound, the turbulence source fields are obtained using an LES and the turbulence-generated sound at far fields is calculated from Lighthill’s acoustic analogy. As only the velocity fields at resolved scales are available from the LES, the Lighthill stress tensor, serving as a source term in Lighthill’s acoustic equation, has to be evaluated from the resolved velocity fields. As a result, the contribution from the unresolved velocity fields is missing in the conventional LES. The sound of missing scales is shown to be important and hence needs to be modeled. The present study proposes a kinematic subgrid-scale (SGS) model which recasts the unresolved velocity fields into Lighthill’s stress tensors. A kinematic simulation is used to construct the unresolved velocity fields with the imposed temporal statistics, which is consistent with the random sweeping hypothesis. The kinematic SGS model is used to calculate sound power spectra from isotropic turbulence and yields an improved result: the missing portion of the sound power spectra is approximately recovered in the LES.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Energetic protons emitted from coulomb explosion dynamics of large-sized hydrogen clusters driven by an ultrashort intense laser pulse

Protons with very high kinetic energy of about 10keV and the saturation effect of proton energy for laser intensity have been observed in the interaction of an ultrashort intense laser pulse with large-sized hydrogen clusters. Including the cluster-size distribution as well as the laser-intensity distribution on the focus spot, the theoretical calculations based on a simplified Coulomb explosion model have been compared with our experimental measurements, which are in good agreement with each other.

Continuous radial superresolution filters with large focal depth

A new set of continuous superresolution filters is proposed which exhibits a radial superresolution performance with an extended depth of focus in an optical system by properly choosing the design parameters. Numerical simulation results of the performance parameters of the superresolution gain, the radial central core size, the Strehl ratio, the side-lobe factor and the depth of focus with different design parameters for the optimized patterns are displayed. We also give a design example for this kind of filter characterized by a birefringent element inserted between two parallel polarizers. This kind of filter would be useful in fields such as optical data storage systems.