Numerical investigation of thermocapillary migration of the drop for large Marangoni numbers

An axisymmetric model is adopted to simulate the problem of unsteady drop thermocapillary motion for large Marangoni numbers. Front tracking methods are used in the investigation. It is found that the non-dimensional drop migration velocity will decrease with increasing Marangoni number. This agrees well with the experimental results obtained from the 4th Shen-Zhou space ship. In the meanwhile, this is also the first time for numerical simulations to verify the experimental phenomenon under large Marangoni numbers.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Experimental Study and Numerical Simulation of Cellular Structures and Mach Reflection of Gaseous Detonation Wave

In this paper the Deflagration to Detonation Transition (DDT) process of gaseous H-2-O-2 mixture and Mach reflection of gaseous detonation wave on a wedge have been conducted experimentally. The cellular pattern of DDT process and Mach reflection were obtained from experiments with wedge angle theta = 10(0) similar to 40(0) and initial pressure of gaseous mixture 16kPa similar to 26.7kPa. The 2-D numerical simulations of DDT process and Mach reflection of detonation wave were performed by using the simplified ZND model and improved space-time conservation element and solution element (CE/SE) method. The numerical cellular structures were compared with the cellular patterns of soot track. Compared results were shown that it is satisfactory. The characteristic comparisons on Mach reflection of air shock wave and detonation wave were carried also out and their differences were given.

A numerical investigation on the flow physics of koi carp (Cyprinus carpio koi) burst-and-coast swimming

Abstract—Burst-and-coast is the most common locomotion type in freely routine swimming of koi carps (Cyprinus carpio koi), which consists of a burst phase and a coast phase in each cycle and mostly leads to a straight-line trajectory. Combining with the tracking experiment, the flow physics of koi carp’s burst-andcoast swimming is investigated using a novel integrated CFD method solving the body-fluid interaction problem. The dynamical equations of a deforming body are formulated. Following that, the loose-coupled equations of the body dynamics and the fluid dynamics are numerically solved with the integrated method. The two burst modes, MT (Multiple Tail-beat) and HT (Half Tail-beat), which have been reported by the experiments, are investigated by numerical simulations in this paper. The body kinematics is predicted and the flow physics is visualized, which are in good agreement with the corresponding experiments. Furthermore, the optimization on the energy cost and several critical control mechanisms in burst-and-coast swimming of koi carps are explored, by varying the parameters in its selfpropelled swimming. In this paper, energetics is measured by the two mechanical quantities, total output power CP and Froude efficiency Fr. Results and discussion show that from the standpoint of mechanical energy, burst-and-coast swimming does not actually save energy comparing with steady swimming at the same average speed, in that frequently changing of speed leads to decrease of efficiency.

LES and hybrid LES/RANS simulations for conjugate heat transfer over a matrix of cubes

Large Eddy Simulation (LES) and a novel k -l based hybrid LES/RANS approach have been applied to simulate a conjugate heat transfer problem involving flow over a matrix of surface mounted cubes. In order to assess the capability and reliability of the newly developed k -l based hybrid LES/RANS, numerical results are compared with new LES and existing RANS results. Comparisons include mean velocity profiles, Reynolds stresses and conjugate heat transfer. As well as for hybrid LES/RANS validation purposes, the LES results are used to gain insights into the complex flow physics and heat transfer mechanisms. Numerical simulations show that the hybrid LES/RANS approach is effective. Mean and instantaneous fluid temperatures adjacent to the cube surface are found to strongly correlate with flow structure. Although the LES captures more mean velocity field complexities, broadly time averaged wake temperature fields are found similar for the LES and hybrid LES/RANS. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Numerical study of the impact of process variations on the motional resistance of weakly coupled MEMS resonators

This paper presents a numerical study of the impact of process-induced variations on the achievable motional resistance Rx of one-dimensional, cyclic and cross-coupled architectures of electrostatically transduced MEMS resonators operating in the 250 kHz range. Monte Carlo numerical simulations which accounted for up to 0.75% variation in critical resonator feature sizes were initiated on 1, 2, 3, 4, 5 and 9 coupled MEMS resonators for three distinct coupling architectures. Improvements of 100X in the spread of Rx and 2.7X in mean achievable Rx are reported for the case of 9 resonators when implemented in the cross-coupled topology, as opposed to the traditional one-dimensional chain. © 2013 IEEE.

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.

The design of curved optical waveguides : analytical and numerical analysis

A model for understanding the formation and propagation of modes in curved optical waveguides is developed. A numerical method for the calculation of curved waveguide mode profiles and propagation constants in two dimensional waveguides is developed, implemented and tested. A numerical method for the analysis of propagation of modes in three dimensional curved optical waveguides is developed, implemented and tested. A technique for the design of curved waveguides with reduced transition loss is presented. A scheme for drawing these new waveguides and ensuring that they have constant width is also provided. Claims about the waveguide design technique are substantiated through numerical simulations.

Cross-correlation frequency resolved optical gating analysis of broadband continuum generation in photonic crystal fiber: Simulations and experiments

Numerical simulations are used to study the temporal and spectral characteristics of broadband supercontinua generated in photonic crystal fiber. In particular, the simulations are used to follow the evolution with propagation distance of the temporal intensity, the spectrum, and the cross-correlation frequency resolved optical gating (XFROG) trace. The simulations allow several important physical processes responsible for supercontinuum generation to be identified and, moreover, illustrate how the XFROG trace provides an intuitive means of interpreting correlated temporal and spectral features of the supercontinuum. Good qualitative agreement with preliminary XFROG measurements is observed. © 2002 Optical Society of America.

Numerical experiments of some Krylov subspace methods for black oil model

This paper discusses preconditioned Krylov subspace methods for solving large scale linear systems that originate from oil reservoir numerical simulations. Two types of preconditioners, one being based on an incomplete LU decomposition and the other being based on iterative algorithms, are used together in a combination strategy in order to achieve an adaptive and efficient preconditioner. Numerical tests show that different Krylov subspace methods combining with appropriate preconditioners are able to achieve optimal performance.

Propagating Wave Phenomena Detected in Observations and Simulations of the Lower Solar Atmosphere

We present high-cadence observations and simulations of the solar photosphere, obtained using the Rapid Oscillations in the Solar Atmosphere imaging system and the MuRAM magnetohydrodynamic (MHD) code, respectively. Each data set demonstrates a wealth of magnetoacoustic oscillatory behavior, visible as periodic intensity fluctuations with periods in the range 110–600 s. Almost no propagating waves with periods less than 140 s and 110 s are detected in the observational and simulated data sets, respectively. High concentrations of power are found in highly magnetized regions, such as magnetic bright points and intergranular lanes. Radiative diagnostics of the photospheric simulations replicate our observational results, confirming that the current breed of MHD simulations are able to accurately represent the lower solar atmosphere. All observed oscillations are generated as a result of naturally occurring magnetoconvective processes, with no specific input driver present. Using contribution functions extracted from our numerical simulations, we estimate minimum G-band and 4170 Å continuum formation heights of 100 km and 25 km, respectively. Detected magnetoacoustic oscillations exhibit a dominant phase delay of −8◦ between the G-band and 4170 Å continuum observations, suggesting the presence of upwardly propagating waves.More than 73% of MBPs (73% from observations and 96% from simulations) display upwardly propagating wave phenomena, suggesting the abundant nature of oscillatory behavior detected higher in the solar atmosphere may be traced back to magnetoconvective processes occurring in the upper layers of the Sun’s convection zone.

Numerical and Experimental Investigation of Aircraft Panel Deformations During Riveting Process

In collaboration with Airbus-UK, the dimensional growth of aircraft panels while being riveted with stiffeners is investigated. Small panels are used in this investigation. The stiffeners have been fastened to the panels with rivets and it has been observed that during this operation the panels expand in the longitudinal and transverse directions. It has been observed that the growth is variable and the challenge is to control the riveting process to minimize this variability. In this investigation, the assembly of the small panels and longitudinal stiffeners has been simulated using static stress and nonlinear explicit finite element models. The models have been validated against a limited set of experimental measurements; it was found that more accurate predictions of the riveting process are achieved using explicit finite element models. Yet, the static stress finite element model is more time efficient, and more practical to simulate hundreds of rivets and the stochastic nature of the process. Furthermore, through a series of numerical simulations and probabilistic analyses, the manufacturing process control parameters that influence panel growth have been identified. Alternative fastening approaches were examined and it was found that dimensional growth can be controlled by changing the design of the dies used for forming the rivets.

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Vlasov-kinetic computer simulations of electrostatic waves in dusty plasmas: an overview of recent results

A series of numerical simulations based on a recurrence-free Vlasov kinetic model using kinetic phase point trajectories are presented. Electron-ion plasmas and three-component (electron-ion-dust) dusty or complex plasmas are considered, via independent simulations. Considering all plasma components modeled through a kinetic approach, the linear and nonlinear behavior of ion-acoustic excitations is investigated. Maxwellian and kappa-type (superthermal) distribution functions are assumed, as initial conditions, in separate simulations for the sake of comparison. The focus is on the parametric dependence of ion-acoustic waves on the electron-to-ion temperature ratio and on the dust concentration. © 2014 EDP Sciences, SIF, Springer-Verlag Berlin Heidelberg.