Moving boundaries due to distributed sources in a slab

Analytical and numerical solutions have been obtained for some moving boundary problems associated with Joule heating and distributed absorption of oxygen in tissues. Several questions have been examined which are concerned with the solutions of classical formulation of sharp melting front model and the classical enthalpy formulation in which solid, liquid and mushy regions are present. Thermal properties and heat sources in the solid and liquid regions have been taken as unequal. The short-time analytical solutions presented here provide useful information. An effective numerical scheme has been proposed which is accurate and simple.

Dynamics of crystal size distributions with size-dependent rates

The growth and dissolution dynamics of nonequilibrium crystal size distributions (CSDs) can be determined by solving the governing population balance equations (PBEs) representing reversible addition or dissociation. New PBEs are considered that intrinsically incorporate growth dispersion and yield complete CSDs. We present two approaches to solving the PBEs, a moment method and a numerical scheme. The results of the numerical scheme agree with the moment technique, which can be solved exactly when powers on mass-dependent growth and dissolution rate coefficients are either zero or one. The numerical scheme is more general and can be applied when the powers of the rate coefficients are non-integers or greater than unity. The influence of the size dependent rates on the time variation of the CSDs indicates that as equilibrium is approached, the CSDs become narrow when the exponent on the growth rate is less than the exponent on the dissolution rate. If the exponent on the growth rate is greater than the exponent on the dissolution rate, then the polydispersity continues to broaden. The computation method applies for crystals large enough that interfacial stability issues, such as ripening, can be neglected. (C) 2002 Elsevier Science B.V. All rights reserved.

Model exact low-lying states and spin dynamics in ferric wheels: Fe-6 to Fe-12

Using an efficient numerical scheme that exploits spatial symmetries and spin parity, we have obtained the exact low-lying eigenstates of exchange Hamiltonians for ferric wheels up to Fe-12. The largest calculation involves the Fe-12 ring which spans a Hilbert space dimension of about 145x10(6) for the M-S=0 subspace. Our calculated gaps from the singlet ground state to the excited triplet state agree well with the experimentally measured values. Study of the static structure factor shows that the ground state is spontaneously dimerized for ferric wheels. The spin states of ferric wheels can be viewed as quantized states of a rigid rotor with the gap between the ground and first excited states defining the inverse of the moment of inertia. We have studied the quantum dynamics of Fe-10 as a representative of ferric wheels. We use the low-lying states of Fe-10 to solve exactly the time-dependent Schrodinger equation and find the magnetization of the molecule in the presence of an alternating magnetic field at zero temperature. We observe a nontrivial oscillation of the magnetization which is dependent on the amplitude of the ac field. We have also studied the torque response of Fe-12 as a function of a magnetic field, which clearly shows spin-state crossover.

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian-Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier-Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant. (C) 2012 Elsevier Inc. All rights reserved.

Operator-splitting finite element algorithms for computations of high-dimensional parabolic problems

An operator-splitting finite element method for solving high-dimensional parabolic equations is presented. The stability and the error estimates are derived for the proposed numerical scheme. Furthermore, two variants of fully-practical operator-splitting finite element algorithms based on the quadrature points and the nodal points, respectively, are presented. Both the quadrature and the nodal point based operator-splitting algorithms are validated using a three-dimensional (3D) test problem. The numerical results obtained with the full 3D computations and the operator-split 2D + 1D computations are found to be in a good agreement with the analytical solution. Further, the optimal order of convergence is obtained in both variants of the operator-splitting algorithms. (C) 2012 Elsevier Inc. All rights reserved.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

Large eddy simulations for fan-OGV broadband noise prediction

The work in this paper forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. Here we focus on LES of the flow field near a fan blade trailing edge. The first part of the paper aims to evaluate LES suitability for predicting the near-field velocity field for a blunt NACA-0012 airfoil at moderate Reynolds numbers (2× 10 5 and 4× 10 5). Preliminary computations of turbulent mean and root-mean-square velocities, as well as energy spectra at the trailing edge, are compared with those from a recent experiment.1 The second part of the paper describes preliminary progress on an LES calculation of the fan wakes on a fan rig. 2 The CFD code uses a mixed element unstructured mesh with a median dual control volume. A wall-adapting local eddy-viscosity sub-grid scale model is employed. A very small amount of numerical dissipation is added in the numerical scheme to keep the compressible solver stable. Further results for the fan turbulentmean and RMS velocity, and especially the aeroacoustics field will be presented at a later stage. Copyright © 2008 by Qinling LI, Nigel Peake & Mark Savill.

Studies on terrestrial interface processes in arid areas, China

In the present paper, we have elucidated the importance of energy and water cycling in arid areas to investigate global climate and local economics. Then, we were concerned with the physical arguments as how to stratify the soil, and the stability of the numerical scheme in the mathematical model for predicting temperature variation and water motion. Furthermore, we discuss the methods to estimate evaporation in arid areas. Numerical simulation of energy and water cycling at the Acsu Observatory, CAS, Xinjiang province and Shapuotou Observatory, CAS, Ningxia Province are conducted as case studies. The results show that the laws of terrestrial processes are rather typical in these arid areas. Planting drought-endurable trees can alleviate unfavourable conditions to a certain extent. (C) 1997 Academic Press Limited.

Simulation of shock-induced bubble collapse with application to vascular injury in shockwave lithotripsy

Shockwave lithotripsy is a noninvasive medical procedure wherein shockwaves are repeatedly focused at the location of kidney stones in order to pulverize them. Stone comminution is thought to be the product of two mechanisms: the propagation of stress waves within the stone and cavitation erosion. However, the latter mechanism has also been implicated in vascular injury. In the present work, shock-induced bubble collapse is studied in order to understand the role that it might play in inducing vascular injury. A high-order accurate, shock- and interface-capturing numerical scheme is developed to simulate the three-dimensional collapse of the bubble in both the free-field and inside a vessel phantom. The primary contributions of the numerical study are the characterization of the shock-bubble and shock-bubble-vessel interactions across a large parameter space that includes clinical shockwave lithotripsy pressure amplitudes, problem geometry and tissue viscoelasticity, and the subsequent correlation of these interactions to vascular injury. Specifically, measurements of the vessel wall pressures and displacements, as well as the finite strains in the fluid surrounding the bubble, are utilized with available experiments in tissue to evaluate damage potential. Estimates are made of the smallest injurious bubbles in the microvasculature during both the collapse and jetting phases of the bubble's life cycle. The present results suggest that bubbles larger than 1 μm in diameter could rupture blood vessels under clinical SWL conditions.

O pré-aquecimento do universo: da ressonância paramétrica à turbulência

A inflação consegue dar conta de uma série de problemas do Modelo padrão da Cosmologia, preservando ainda o sucesso do modelo do Big Bang. Na sua versão mais simples, a inflação é controlada por um campo escalar, o ínflaton, que faz com que o universo se expanda exponencialmente. Após, o ínflaton decai e ocorre o reaquecimento do universo. Contudo, alguns autores apontam a existência de uma fase intermediária, chamada de pré-aquecimento. O decaimento do ínflaton possui uma rica dinâmica não-linear. No primeiro estágio, a ressonância paramétrica promove o crescimento exponencial de alguns modos do ínflaton. Isto altera a dinâmica do modo homogêneo do ínflaton, promovendo uma reestruturação das cartas de ressonâncias da equação de movimento dos modos perturbativos. Desta forma, ocorre a transferência de energia para estes modos, até que o universo termaliza. Esta transferência de energia é típica de um sistema turbulento. Por se tratar de uma evolução não-linear, torna-se conveniente a implementação computacional de métodos numéricos. Neste contexto, os métodos espectrais têm se mostrado uma excelente ferramenta para estudar este tipo de sistema. Esta dissertação apresenta os resultados do esquema numérico desenvolvido para o modelo com potencial quártico, que será a base para os demais estudos a serem desenvolvidos. Como mostrado, este esquema é extremamente preciso e eficiente.

Formação de nanopadrões em superfícies por sputtering iônico: Estudo numérico da equação anisotrópica amortecida de Kuramoto-Sivashinsky.

Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência.

Generalized phenomenological model for the viscoelasticity of idealized asphalts

A generalized theory for the viscoelastic behavior of idealized bituminous mixtures (asphalts) is presented. The mathematical model incorporates strain rate and temperature dependency as well as nonmonotonic loading and unloading with shape recovery. The stiffening effect of the aggregate is included. The model is of phenomenological nature. It can be calibrated using a relatively limited set of experimental parameters, obtainable by uniaxial tests. It is shown that the mathematical model can be represented as a special nonlinear form of the Burgers model. This facilitates the derivation of numerical algorithms for solving the constitutive equations. A numerical scheme is implemented in a user material subroutine (UMAT) in the finite-element analysis (FEA) code ABAQUS. Simulation results are compared with uniaxial and indentation tests on an idealized asphalt mix. © 2014 American Society of Civil Engineers.

青藏高原壳幔形变与层圈耦合数值模拟

Based on multi-principle (such as structures, tectonics and kinematics) exploratory data and related results of continental dynamics in the Tibetan plateau, the author reconstructed the geological-geophysical model of lithospherical structure and tectonic deformation, and the kinetics boundary conditions for the model. Then, the author used the numerical scheme of Fast Lagrangian Analysis of Continua (FLAC), to stimulate the possible process of the stress field and deformational field in the Tibetan plateau and its adjacent area, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. With the above-mentioned results, the author discussed the relationship between crustal movement in shallow layer and the deformational process in interior layers, and its possible dynamic constraints in deep. At the end of the paper, an integrative model has been put forward to explain the outline images of crust-mantle deformation and coupling in the Tibetan Plateau. (1) The characteristics of crust-mantle structure of the Tibetan plateau have been shown to be very complex, and vertical and horizontal difference is significant. The general characteristics of crust-mantle of the Tibetan plateau may be that it's layering in depth direction, and shows blocking from south to north and belting from east to west, mainly according to the results of about 20 seismic sections, such as wide-angle seismic profiles, CMP, seismic tomography and so on. (2) The crust had shortened about 2200km, while the shortening is different for different block from south to north in the Tibetan plateau. It is about 11.5mm/a in Himalayan block, about 9.0mm/a in Lhas-Gangdese block, about 7.0mm/a in Qiangtang block and Songpan-Ganzi-Kekexili block, about 8.0mm/a in Kunlun-Qaidam, and about ll.Omm/a in Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. Which - in demonstrates the shortening rate decreases from south to north, but this rate increases near the north edge of the Tibetan plateau. The crust thickening rate is about 0.4mm/a in the whole Tibetan plateau; and this rate is about 0.5mm/a in Himalayan block, about 0.4mm/a in Lhas-Gangdese block, about 0.3mm/a in Qiangtang block, about 0.2mm/a in Songpan-Ganzi-Kekexili block and about O.lmm/a in Kunlun-Qaidam-Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. This implies that the thickening rate decreases in the blocks of the Tibetan plateau. From south to north, the displacement of eastern boundary in the Tibetan plateau is about 37mm/a in Himalayan block, about 45mm/a in Lhas-Gangdese block, about 47mm/a in Qiangtang block, about 43mm/a in Songpan-Ganzi-Kekexili block, and about 35mm/a in Kunlun-Qaidam-Qilian block, since the collision-matching between the Indian continent and Eurasia continent had happened about 50Ma ago. This implies that the rate of eastward displacement is biggest in the middle of plateau, and decreases to both sides. The transition of S-N compression stress field in Tibetan Plateau, since about 28Ma+ ago, may be caused by two reasons: On one hand, the movement direction of Eurasia continent changed from northward to southward about 28Ma± ago in the northern plateau. On the other hand, the front belt that is located between India continent's and Eurasia continent's convergence-collision, had moved southward to high Himalayan from Indus-Brahmaputra suture almost at the same time in southern plateau. Affected by the stress field, the earlier tectonics rotated clockwise, NE and NW conjugate strike-slip faults developed, and the SN rift formed. This indicated that the EW movement started. The ratio between upper crust and lower crust of different blocks from south to north in the Tibetan plateau during the process of deformation are as following: about 3.5~5:1 in Himalayan block, about 1~5: 3-4 (which is about 1:3o--4 in south and about 4~5:3 in north) in Lhas-Gangdese block, about 1:3~447mm/a in these blocks: Which is located to the north of Banggong-nujiang suture.

A domain decomposition algorithm for inverse welding problems

A mathematical model and a numerical scheme for the inverse determination of heat sources generated by means of a welding process is presented in this paper. The accuracy of the heat source retrieval is discussed.