WorldWideScience

Sample records for 2-dimensional numerical modeling

  1. Development of a numerical 2-dimensional beach evolution model

    Baykal, Cüneyt

    2014-01-01

    to compute the nearshore depth-averaged wave-induced current velocities and mean water level changes, a sediment transport model to compute the local total sediment transport rates occurring under the action of wind waves, and a bottom evolution model to compute the bed level changes in time based on......This paper presents the description of a 2-dimensional numerical model constructed for the simulation of beach evolution under the action of wind waves only over the arbitrary land and sea topographies around existing coastal structures and formations. The developed beach evolution numerical model...... is composed of 4 submodels: a nearshore spectral wave transformation model based on an energy balance equation including random wave breaking and diffraction terms to compute the nearshore wave characteristics, a nearshore wave-induced circulation model based on the nonlinear shallow water equations...

  2. 2-dimensional numerical modeling of active magnetic regeneration

    Nielsen, Kaspar Kirstein; Pryds, Nini; Smith, Anders;

    2009-01-01

    Various aspects of numerical modeling of Active Magnetic Regeneration (AMR) are presented. Using a 2-dimensional numerical model for solving the unsteady heat transfer equations for the AMR system, a range of physical effects on both idealized and non-idealized AMR are investigated. The modeled...... system represents a linear, parallel-plate based AMR. The idealized version of the model is able to predict the theoretical performance of AMR in terms of cooling power and temperature span. This is useful to a certain extent, but a model reproducing experiments to a higher degree is desirable. Therefore...... physical effects such as thermal parasitic losses have been included. Furthermore, experimentally found magnetocaloric properties are used when available, since the commonly used mean field model can be too idealized and is not always able to determine the magnetocaloric effect accurately. In the present...

  3. Development of orthogonal 2-dimensional numerical code TFC2D for fluid flow with various turbulence models and numerical schemes

    Park, Ju Yeop; In, Wang Kee; Chun, Tae Hyun; Oh, Dong Seok [Korea Atomic Energy Research Institute, Taejeon (Korea)

    2000-02-01

    The development of orthogonal 2-dimensional numerical code is made. The present code contains 9 kinds of turbulence models that are widely used. They include a standard k-{epsilon} model and 8 kinds of low Reynolds number ones. They also include 6 kinds of numerical schemes including 5 kinds of low order schemes and 1 kind of high order scheme such as QUICK. To verify the present numerical code, pipe flow, channel flow and expansion pipe flow are solved by this code with various options of turbulence models and numerical schemes and the calculated outputs are compared to experimental data. Furthermore, the discretization error that originates from the use of standard k-{epsilon} turbulence model with wall function is much more diminished by introducing a new grid system than a conventional one in the present code. 23 refs., 58 figs., 6 tabs. (Author)

  4. New hidden symmetries in 2-dimensional models

    In an attempt to derive the hidden symmetries for some integrable 2-dimensional models by considering the invariances of the corresponding linearization systems and the Riemann-Hilbert transformations, we arrive at a new ''sub''-algebra of the ordinary Kac-Moody algebra which represents the hidden symmetry for for example the sine-Gordon theory. A similar ''sub''-algebra is found for the Liouville model. These new algebras differ from the ordinary ones in having a different structure according to whether the grading is even or odd. We describe a new systematic way of finding such hidden symmetries from general linearization systems. (orig.)

  5. Antarctic ice sheet GLIMMER model test and its simplified model on 2-dimensional ice flow

    Xueyuan Tang; Zhanhai Zhang; Bo Sun; Yuansheng Li; Na Li; Bangbing Wang; Xiangpei Zhang

    2008-01-01

    The 3-dimensional finite difference thermodynamic coupled model on Antarctic ice sheet, GLIMMER model, is described. An ide-alized ice sheet numerical test was conducted under the EISMINT-I benchmark, and the characteristic curves of ice sheets under steady state were obtained. Based on this, this model was simplified from a 3-dimensional one to 2-dimensional one. Improvement of the dif-ference method and coordinate system was proposed. Evolution of the 2-dimensional ice flow was simulated under coupled temperature field conditions. The results showed that the characteristic curves deriving from the conservation of the mass, momentum and energy agree with the results of ice sheet profile simulated with GLIMMER model and with the theoretical results. The application prospect of the simplified 2-dimensional ice flow model to simulate the relation of age-depth-accumulation in Dome A region was discussed.

  6. Constructive Renormalization of 2-dimensional Grosse-Wulkenhaar Model

    Wang, Zhituo

    2012-01-01

    In this talk we briefly report the recent work on the construction of the 2-dimensional Grosse-Wulkenhaar model with the method of loop vertex expansion. We treat renormalization with this new tool, adapt Nelson's argument and prove Borel summability of the perturbation series. This is the first non-commutative quantum field theory model to be built in a non-perturbative sense.

  7. Damage spreading in 2-dimensional isotropic and anisotropic Bak-Sneppen models

    Bakar, Burhan; Tirnakli, Ugur

    2007-01-01

    We implement the damage spreading technique on 2-dimensional isotropic and anisotropic Bak-Sneppen models. Our extensive numerical simulations show that there exists a power-law sensitivity to the initial conditions at the statistically stationary state (self-organized critical state). Corresponding growth exponent $\\alpha$ for the Hamming distance and the dynamical exponent $z$ are calculated. These values allow us to observe a clear data collapse of the finite size scaling for both versions...

  8. Damage spreading in 2-dimensional isotropic and anisotropic Bak-Sneppen models

    Bakar, B.; Tirnakli, U.

    2008-03-01

    We implement the damage spreading technique on 2-dimensional isotropic and anisotropic Bak-Sneppen models. Our extensive numerical simulations show that there exists a power-law sensitivity to the initial conditions at the statistically stationary state (self-organized critical state). Corresponding growth exponent α for the Hamming distance and the dynamical exponent z are calculated. These values allow us to observe a clear data collapse of the finite size scaling for both versions of the Bak-Sneppen model. Moreover, it is shown that the growth exponent of the distance in the isotropic and anisotropic Bak-Sneppen models is strongly affected by the choice of the transient time.

  9. Numerical models

    Unnikrishnan, A.S.; Manoj, N.T.

    the wetted perimeter and A the area of cross section (excluding mud flats); C = (1.49/n)R1/6, where n is the Manning coefficient. The numerical scheme used by Harleman and Lee (1969) was used to solve the above equations. In this scheme, the continuity... equation is solved at odd grid points to compute eta at the next time step and the momentum equation is solved at even grid points to compute U . The original scheme of Harleman & Lee (1969) was developed for a single channel. For developing a model...

  10. Towards automatic calibration of 2-dimensional flood propagation models

    P. Fabio

    2009-11-01

    Full Text Available Hydraulic models for flood propagation description are an essential tool in many fields, e.g. civil engineering, flood hazard and risk assessments, evaluation of flood control measures, etc. Nowadays there are many models of different complexity regarding the mathematical foundation and spatial dimensions available, and most of them are comparatively easy to operate due to sophisticated tools for model setup and control. However, the calibration of these models is still underdeveloped in contrast to other models like e.g. hydrological models or models used in ecosystem analysis. This has basically two reasons: first, the lack of relevant data against the models can be calibrated, because flood events are very rarely monitored due to the disturbances inflicted by them and the lack of appropriate measuring equipment in place. Secondly, especially the two-dimensional models are computationally very demanding and therefore the use of available sophisticated automatic calibration procedures is restricted in many cases. This study takes a well documented flood event in August 2002 at the Mulde River in Germany as an example and investigates the most appropriate calibration strategy for a full 2-D hyperbolic finite element model. The model independent optimiser PEST, that gives the possibility of automatic calibrations, is used. The application of the parallel version of the optimiser to the model and calibration data showed that a it is possible to use automatic calibration in combination of 2-D hydraulic model, and b equifinality of model parameterisation can also be caused by a too large number of degrees of freedom in the calibration data in contrast to a too simple model setup. In order to improve model calibration and reduce equifinality a method was developed to identify calibration data with likely errors that obstruct model calibration.

  11. Localized Solutions in a 2 Dimensional Landau-Lifshitz Model

    Piette, B.; Zakrzewski, W. J.

    1996-01-01

    We demonstrate the existence of stable time dependent solutions of the Landau-Lifshitz model with a constant external magnetic field. We find such solutions in all topological sectors, including N=0. We discuss some of their properties.

  12. Effective Hydraulic Conductivity Scaling in a 2-Dimensional Geometrical Multifractal Model for Aquifer Heterogeneity

    Gentry, R. W.; Perfect, E.; Sukop, M. C.

    2005-12-01

    Recent analyses of field data suggest that the spatial variation of hydraulic conductivity, K, within an aquifer may be multifractal. We investigated the implications of this finding for the scaling of effective hydraulic conductivity, , by performing numerical simulations of flow in 2-dimensional geometrical multifractal K fields. A theoretical framework for generating such fields is presented based on the parameters of the truncated binomial distribution, TBD. This leads to an approximate analytical expression showing that increases with increasing length scale as a power law, whose exponent, α, is determined by the TBD parameters. Five geometrical multifractal K fields were generated with different minimum length scales. Each domain was discretized using a block center grid consisting of 59,049 uniformly-spaced nodes. A unit cube aquifer was used for the numerical simulations. The boundary conditions were implemented with constant head (unit gradient) parallel planes, and corresponding zero flux planes on the normal axes. A finite difference simulation model based on MODFLOW 2000 was used, and "zone budget" was employed to calculate the flow balance. The discharge into and out of the unit cube was then used to calculate based on Darcy's law. The numerical simulations produced similar increases in with increasing length scale to those predicted by the analytical model. Nonlinear regression analyses yielded estimates of α from the numerical simulations that were within 10% of the analytical value for these fields. These simulations provide a theoretical explanation for effective hydraulic conductivity scaling in terms of multifractals. The advantage of such an approach is that the α-parameter, which controls the degree of scaling, is physically-based and can potentially be estimated from independent measurements.

  13. A 2-dimensional fully analytical model for design of high voltage junction barrier Schottky (JBS) diodes

    Radhakrishnan, Rahul; Zhao, Jian H.

    2011-09-01

    A physics-based closed form analytical model for the reverse leakage current of a high voltage junction barrier Schottky (JBS) diode is developed and shown to agree with experimental results. Maximum electric field "seen" by the Schottky contact is calculated from first principles by a 2-dimensional method as a function of JBS diode design parameters and confirmed by numerical simulations. Considering thermionic emission under image force barrier lowering and quantum mechanical tunneling, electric field at the Schottky contact is then related to reverse current. In combination with previously reported forward current and resistance models, this gives a complete I- V relationship for the JBS diode. A layout of interdigitated stripes of P-N and Schottky contacts at the anode is compared theoretically with a honeycomb layout and the 2-D model is extended to the 3-D honeycomb structure. Although simulation and experimental results from 4H-Silicon Carbide (SiC) diodes are used to validate it, the model itself is applicable to all JBS diodes.

  14. Dynamical Modelling and Controllability Analysis of an Underactuated 2-Dimensional TORA System on a Slope

    Bingtuan Gao

    2014-04-01

    Full Text Available The 2-Dimensional Translational Oscillators with Rotating Actuator (2DTORA is a novel underactuated system which has one actuated rotor and two unactuated translational carts. This paper focuses on dynamical modelling and simulation analysis of the underactuated 2DTORA on a slope. Based on Lagrange equations, the dynamics of the 2DTORA is achieved by selecting a transverse position of a cart, a travelling position of a cart, and the rotor angle as the general coordinates and torque acting on the rotor as the general force. When the slope angle is set to zero, the dynamics of 2DTORA on a slope is reduced to that of 2DTORA on the horizontal plane. Moreover, by eliminating one degree of translational cart motion, the dynamics of 2DTORA is reduced to that of TORA which is a benchmark of underactuated systems. In addition, the equilibrium and controllability of the 2DTORA system on a slop are discussed. Finally, numerical simulations are performed to verify the feasibility of the developed dynamic models.

  15. Investigation of two different anoxia models by 2-dimensional gel electrophoresis

    Wulff, Tune; Jessen, Flemming; Hoffmann, Else Kay

    anoxia obtained by NaN3 is a widely used model for simulating anoxia (Ossum et al., 2004). The effects of anoxia were studied by protein expression analysis using 2-dimensional gel electrophoresis followed by MS/MS. In this way we were able to separate more than 1500 protein spots with an apparent range...

  16. Dynamical analysis and simulation of a 2-dimensional disease model with convex incidence

    Yu, Pei; Zhang, Wenjing; Wahl, Lindi M.

    2016-08-01

    In this paper, a previously developed 2-dimensional disease model is studied, which can be used for both epidemiologic modeling and in-host disease modeling. The main attention of this paper is focused on various dynamical behaviors of the system, including Hopf and generalized Hopf bifurcations which yield bistability and tristability, Bogdanov-Takens bifurcation, and homoclinic bifurcation. It is shown that the Bogdanov-Takens bifurcation and homoclinic bifurcation provide a new mechanism for generating disease recurrence, that is, cycles of remission and relapse such as the viral blips observed in HIV infection.

  17. Numerical Modelling of Streams

    Vestergaard, Kristian

    In recent years there has been a sharp increase in the use of numerical water quality models. Numeric water quality modeling can be divided into three steps: Hydrodynamic modeling for the determination of stream flow and water levels. Modelling of transport and dispersion of a conservative...

  18. 2-dimensional modelling of the steam refomer of the ADAM/EVA II facility

    For the steam reformer of the ADAM/EVA II facility a 2-dimensional computer code is made, which includes local radiation quantities and by-pass circuits. Not only integral sizes (reaction, fission and temperature) and curve states as well as individual tube states are computable by this programme. It is shown that the experimental results are to understand very good with this model. For a steam reformer with NTIW-geometry in the 20 MW range the model was changed. Finally an anchor cooling is proposed to release the thermal-mechanical heavily loaded supporting structure. It is possible to describe quantitatively the cooling effect with the plotted computer code. It shows that the anchor wall temperature falls round about 50 K in a 5% cooling circuit. This could prevent the extension of anchors, as well as prolong the lifetime of the bundle and make possible the competence of license. (orig.)

  19. Development of non-orthogonal and 2-dimensional numerical code TFC2D-BFC for fluid flow

    Park, Ju Yeop; In, Wang Kee; Chun, Tae Hyun; Oh, Dong Seok

    2000-09-01

    The development of algorithm for three dimensional non-orthogonal coordinate system has been made. The algorithm adopts a non-staggered grid system, Cartesian velocity components for independent variables of momentum equations and a SIMPLER algorithm for a pressure correction equation. Except the pressure correction method, the selected grid system and the selected independent variables for momentum equations have been widely used in a commercial code. It is well known that the SIMPLER is superior to the SIMPLE algorithm in the view of convergence rate. Using this algorithm, a two dimensional non-orthogonal numerical code has been completed. The code adopts a structured single square block in a computational domain with a uniform mesh interval. Consequently, any solid body existing in a flow field can be implemented in the numerical code through a blocked-off method which was devised by Patankar.

  20. Performance Of Bathymetric Lidar On Flow Properties Predicted With A 2-Dimensional Hydraulic Model

    Tonina, D.; McKean, J. A.; Wright, C. W.

    2014-12-01

    Increased computer processing speeds and new computational fluid dynamics codes have significantly improved numerical modeling of flow and sediment transport over large domains of streams, up to several kilometers in length. Recent developments in remote sensing technologies have also greatly improved our ability to map the morphology of streams over similar spatial extents. However, limited information is available on whether the remote sensing methods can map channel topography with sufficient accuracy to define the flow boundary necessary for a fluid dynamics model. We assessed the ability of a second generation airborne bathymetric sensor, the Experimental Advanced Airborne Research Lidar (EAARL-B), to support a two dimensional fluid dynamics model of a small morphologically-complex mountain stream. We compared flow model predictions using the lidar bathymetry with those made using a total station field survey of the channel. In this riverscape, results suggest EAARL bathymetric lidar can map channel topography with sufficient accuracy to support a two dimensional computational flow model.

  1. Numerical Transducer Modeling

    Henriquez, Vicente Cutanda

    This thesis describes the development of a numerical model of the propagation of sound waves in fluids with viscous and thermal losses, with application to the simulation of acoustic transducers, in particular condenser microphones for measurement. The theoretical basis is presented, numerical...... tools and implementation techniques are described and performance tests are carried out. The equations that govern the motion of fluids with losses and the corresponding boundary conditions are reduced to a form that is tractable for the Boundary Element Method (BEM) by adopting some hypotheses...... that are allowable in this case: linear variations, absence of flow, harmonic time variation, thermodynamical equilibrium and physical dimensions much larger than the molecular mean free path. A formulation of the BEM is also developed with an improvement designed to cope with the numerical difficulty associated...

  2. Numerical simulations of lower hybrid current drive by the (1+2)-dimensional Fokker-Planck code combined with the ray-tracing code

    The (1+2) dimensional (1-D in real space and 2-D in momentum space) relativistic Fokker-Planck code combined with the ray-tracing code has been newly developed for analyzing the lower hybrid current drive (LHCD) on tokamak plasmas. This numerical code calculates the 2-D MHD equilibrium, ray-tracing, and Fokker-Planck analyses self-consistently. The simulations have been carried out by using the code, and the results were compared with experiments on the lower hybrid current drive in JT-60. As a result, the obtained simulation results agreed with the experimental results on the current drive efficiency and the current density profile (or the plasma internal inductance). (author)

  3. Numerical transducer modelling

    Cutanda, Vicente

    . However, there are several difficulties to be addressed that are derived from the size, internal structure and precision requirements that are characteristic of these devices. One of them, the presence of very close surfaces (e.g. the microphone diaphragm and back-electrode), leads to machine precision......Numerical modelling is of importance for the design, improvement and study of acoustic transducers such as microphones and accelerometers. Techniques like the boundary element method and the finite element method are the most common supplement to the traditional empirical and analytical approaches...

  4. Numerical modeling of creep in metals: numerical results

    A 2-dimensional initial strain direct boundary element method was proposed to numerically model the creep behaviour of metals. The boundary of the body was discretized into quadratic elements and the domain into quadratic quadrilaterals. The variables were also assumed to have a quadratic variation over the elements. Due to the time-dependent nature of creep, the solution was derived over increments of time. Backward Euler method and automatic time incrementation technique for updating the variables were implemented to assure stability and accuracy of results. An algorithm had been developed to implement this method. The results for a square plate under biaxial load were compared to analytical solutions and showed to be in good agreement with errors of 4.17% and 1.91% for effective creep strain and creep strain in the xx-direction, respectively. Parametric study of the effects of varying the initial time step, tolerance range and convergence criteria on the numerical results were also carried out. (Author)

  5. An effective depression filling algorithm for DEM-based 2-dimensional surface flow modelling

    Zhu, D.; Ren, Q.; Xuan, Y.; Y. Chen; Cluckie, I.

    2012-01-01

    The surface runoff process in fluvial/pluvial flood modelling is often simulated employing a two-dimensional (2-D) diffusive wave approximation to described by grid based digital elevation models (DEMs). However, a serious problem of this approach may arise when using a 2-D surface flow model which exchanges flows through adjacent cells, or conventional rink removal algorithms which also allow flow to be exchanged along diagonal directions, due to the existence of artificial depression in ...

  6. An effective depression filling algorithm for DEM-based 2-dimensional surface flow modelling

    D. Zhu

    2012-09-01

    Full Text Available The surface runoff process in fluvial/pluvial flood modelling is often simulated employing a two-dimensional (2-D diffusive wave approximation to described by grid based digital elevation models (DEMs. However, a serious problem of this approach may arise when using a 2-D surface flow model which exchanges flows through adjacent cells, or conventional rink removal algorithms which also allow flow to be exchanged along diagonal directions, due to the existence of artificial depression in DEMs. This study firstly analyses the two types of depressions in DEMs and reviews the current depression filling algorithms with a medium sized basin in South-East England, the Upper Medway Catchment (220 km2 used to demonstrate the depression issue in 2-D surface runoff simulation by MIKE SHE with different DEM resolutions (50 m, 100 m and 200 m. An alternative depression-filling algorithm for 2-D overland flow modelling is developed and evaluated by comparing the simulated flows at the outlet of the catchment. This result suggests that the depression estimates at different grid resolution of DEM highly influences overland flow estimation and the new depression filling algorithm is shown to be effective in tackling this issue when comparing simulations in sink-dominated and sink-free digital elevation models, especially for depressions in relatively flat areas on digital land surface models.

  7. Numerical modeling of economic uncertainty

    Schjær-Jacobsen, Hans

    2007-01-01

    Representation and modeling of economic uncertainty is addressed by different modeling methods, namely stochastic variables and probabilities, interval analysis, and fuzzy numbers, in particular triple estimates. Focusing on discounted cash flow analysis numerical results are presented, comparisons...

  8. Numerical modeling of gas hydrate emplacements in oceanic sediments

    Schnurle, Philippe; Liu, Char-Shine

    2011-01-01

    We have implemented a 2-dimensional numerical model for simulating gas hydrate and free gas accumulation in marine sediments. The starting equations are those of the conservation of the transport of momentum, energy, and mass, as well as those of the thermodynamics of methane hydrate stability and methane solubility in the pore-fluid. These constitutive equations are then integrated into a finite element in space, finite-difference in time scheme. We are then able to examine the formation and...

  9. Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models

    Bakar, Burhan; Tirnakli, Ugur

    2008-01-01

    The conventional Hamming distance measurement captures only the short-time dynamics of the displacement between the uncorrelated random configurations. The minimum difference technique introduced by Tirnakli and Lyra [Int. J. Mod. Phys. C 14, 805 (2003)] is used to study the short-time and long-time dynamics of the two distinct random configurations of the isotropic and anisotropic Bak-Sneppen models on a square lattice. Similar to 1-dimensional case, the time evolution of the displacement is...

  10. Convergence dynamics of 2-dimensional isotropic and anisotropic Bak Sneppen models

    Bakar, Burhan; Tirnakli, Ugur

    2008-09-01

    The conventional Hamming distance measurement captures only short-time dynamics of the displacement between uncorrelated random configurations. The minimum difference technique introduced by Tirnakli and Lyra [U. Tirnakli, M.L. Lyra. Int. J. Mod. Phys. C 14 (2003) 805] is used to study short-time and long-time dynamics of the two distinct random configurations of isotropic and anisotropic Bak-Sneppen models on a square lattice. Similar to a 1-dimensional case, the time evolution of the displacement is intermittent. The scaling behavior of the jump activity rate and waiting time distribution reveal the absence of typical spatial-temporal scales in the mechanism of displacement jumps used to quantify convergence dynamics.

  11. Numerical modeling and analysis of the active magnetic regenerator

    Nielsen, Kaspar Kirstein

    In this thesis the active magnetic regenerator (AMR) is analyzed using various numerical tools and experimental devices. A 2-dimensional transient numerical model of the AMR is developed and implemented and it is used to investigate the in uence of a range of parameters on the performance of the...... investigated using the numerical AMR model. The results show indeed that the performance may be enhanced signicantly and it may thus be concluded that the performance of the AMR is dependent on a vast number of parameters (material composition, magnetic eld source, regenerator geometry, regenerator eciency...... AMR. The model simulates a regenerator made of parallel plates. The operating parameters, such as uid ow rates, thermal utilization, magnetocaloric properties etc. are varied as are geometric properties such as plate and channel thickness, regenerator length and porosity. In this way the performance...

  12. Numerical modeling of parallel-plate based AMR

    In this work we present an improved 2-dimensional numerical model of a parallel-plate based AMR. The model includes heat transfer in fluid and magnetocaloric domains respectively. The domains are coupled via inner thermal boundaries. The MCE is modeled either as an instantaneous change between high...... the direction not resolved through a realistic description of the thermal resistance between localized points in the bed and the ambient. The results show that the additions to the model place numerical modeling of AMR very close to the corresponding experimental results. Thus, the model is verified by...... direct comparison with experiment. This is used as a firm basis for predicting and optimizing performance of a large variety of regenerator configurations in order to study and learn the trends, tendencies and even absolute values of temperature span and cooling powers for the optimal (and buildable...

  13. Numerical models for differential problems

    Quarteroni, Alfio

    2014-01-01

    In this text, we introduce the basic concepts for the numerical modelling of partial differential equations. We consider the classical elliptic, parabolic and hyperbolic linear equations, but also the diffusion, transport, and Navier-Stokes equations, as well as equations representing conservation laws, saddle-point problems and optimal control problems. Furthermore, we provide numerous physical examples which underline such equations. We then analyze numerical solution methods based on finite elements, finite differences, finite volumes, spectral methods and domain decomposition methods, and reduced basis methods. In particular, we discuss the algorithmic and computer implementation aspects and provide a number of easy-to-use programs. The text does not require any previous advanced mathematical knowledge of partial differential equations: the absolutely essential concepts are reported in a preliminary chapter. It is therefore suitable for students of bachelor and master courses in scientific disciplines, an...

  14. Numerical Modelling of Scramjet Combustor

    M. Deepu

    2007-07-01

    Full Text Available Numerical modelling of turbulent-reacting flow field of supersonic combustion ramjet(scramjet combustors are presented. The developed numerical procedure is based on the implicittreatment of chemical source terms by preconditioning and solved along with unstedy turbulentNavier-Stokes equations explicitly. Reaction is modelled using an eight-step hydrogen-airchemistry. Code is validated against a standard wall jet experimental data and is successfullyused to model the turbulent-reacting flow field resulting due to the combustion of hydrogeninjected from diamond-shaped strut and also in the wake region of wedge-shaped strut placedin the heated supersonic airstream. The analysis could demonstrate the effect of interaction ofoblique shock wave with a supersonic stream of hydrogen  in its (fuel-air mixing and reactionfor strut-based scramjet combustors.

  15. Numerical modelling of rapid solidification

    Pryds, Nini; Hattel, Jesper Henri

    1997-01-01

    A mathematical model of the melt spinning process has been developed based on the control-volume finite-difference method. The model avoids some of the limitations of the previous models, for example including the effect of the wheel in the heat how calculations and the temperature dependence of ...... decrease with increasing thermal conductivity of the wheel. The observed increase in the wheel surface temperature suggests the importance of including the wheel in the numerical calculations, especially for a wheel made of a low-conductive material...

  16. Plasma modelling and numerical simulation

    Plasma modelling is an exciting subject in which virtually all physical disciplines are represented. Plasma models combine the electromagnetic, statistical and fluid dynamical theories that have their roots in the 19th century with the modern insights concerning the structure of matter that were developed throughout the 20th century. The present cluster issue consists of 20 invited contributions, which are representative of the state of the art in plasma modelling and numerical simulation. These contributions provide an in-depth discussion of the major theories and modelling and simulation strategies, and their applications to contemporary plasma-based technologies. In this editorial review, we introduce and complement those papers by providing a bird's eye perspective on plasma modelling and discussing the historical context in which it has surfaced. (editorial review)

  17. Numerical Modelling of Overburden Deformations

    J. Barták

    2002-01-01

    Full Text Available This paper focuses on the application and verification of mathematical models of the effect of supporting measures on the reduction of overburden deformations. The study of the behaviour of the models is divided into three parts: reduction of the tunnelling effects on the Minorit monastery by means of a jet-grouting curtain; the behaviour of the Hvížďalka backfilled tunnel and a numerical analysis of the supporting measures affecting the tunnel deformations of the Mrázovka tunnel in Prague.

  18. Comprehensive numerical modelling of tokamaks

    We outline a plan for the development of a comprehensive numerical model of tokamaks. The model would consist of a suite of independent, communicating packages describing the various aspects of tokamak performance (core and edge transport coefficients and profiles, heating, fueling, magnetic configuration, etc.) as well as extensive diagnostics. These codes, which may run on different computers, would be flexibly linked by a user-friendly shell which would allow run-time specification of packages and generation of pre- and post-processing functions, including workstation-based visualization of output. One package in particular, the calculation of core transport coefficients via gyrokinetic particle simulation, will become practical on the scale required for comprehensive modelling only with the advent of teraFLOP computers. Incremental effort at LLNL would be focused on gyrokinetic simulation and development of the shell

  19. Numerical modeling of foam flows

    Liquid foam flows are involved in numerous applications, e.g. food and cosmetics industries, oil extraction, nuclear decontamination. Moreover, their study leads to fundamental knowledge: as it is easier to manipulate and analyse, foam is used as a model material to understand the flow of emulsions, polymers, pastes, or cell aggregates, all of which display both solid and liquid behaviour. Systematic experiments performed by Francois Graner et al. provide precise data that emphasize the non Newtonian properties of the foam. Meanwhile, Pierre Saramito proposed a visco-elasto-plastic continuous tensorial model, akin to predict the behaviour of the foam. The goal of this thesis is to understand this complex behaviour, using these two elements. We have built and validated a resolution algorithm based on a bidimensional finite elements methods. The numerical solutions are in excellent agreement with the spatial distribution of all measured quantities, and confirm the predictive capabilities of the model. The dominant parameters have been identified and we evidenced the fact that the viscous, elastic, and plastic contributions to the flow have to be treated simultaneously in a tensorial formalism. We provide a substantial contribution to the understanding of foams and open the path to realistic simulations of complex VEP flows for industrial applications. (author)

  20. Numerical modeling of leachate production

    A key issue for sustainable development of Pakistan is the management and disposal of increased quantities of waste. The majority of the municipal solid waste is presently disposed off by tipping or dumping on the land surface- a process known as land filling. waste materials disposed off as land filling are progressively decomposed by biochemical and physical reactions , producing a number of By-products. Those of most concerned for environmental protection point of view are landfill gas and leachate. Landfill gas is a mixture of methane, carbon dioxide, oxygen and various trace gases, many of which can be detrimental to human health. Leachate is the liquid pollutant generated during waste decomposition, which has the potential to contaminate surface and ground water resources. The prediction of water resources is a fundamental consideration in managing landfill operations. Landfill sites should be designed and operated so as to control leachate production and hence minimize the risk of surface and groundwater pollution. Whether leachate is to be collect and treated or allowed to discharge to the surrounding soil and water , it may be essential to have estimates of leachate flow and strength. This paper considers the development of an applied numerical model named NUMMOL (Numerical Modeling of Leachate based on the waste balance approach, which simulates the leachate production, movement and distribution within landfill sites. NUMMOL incorporates the most appropriate mathematical models representing the various landfill hydrological processes. The moisture flow through waste layers is modeled using the models derived through experimental investigation. The primary purpose of the model is to estimate the leachate quantity and to assist in the comparison of landfill design alternatives as judged by their water balances. (authors)

  1. Numerical experiments modelling turbulent flows

    Trefilík, Jiří; Kozel, Karel; Příhoda, Jaromír

    Kutná Hora : E D P Sciences, 2014 - (Vít, T.; Dančová, P.; Novotný, P.) ISBN 978-80-260-5375-0. ISSN 2101-6275. - (EPJ Web of Conferences. 67). [Experimental Fluid Mechanics 2013. Kutná Hora (CZ), 19.11.2013-22.11.2013] R&D Projects: GA ČR(CZ) GAP101/10/1329; GA ČR GAP101/12/1271 Institutional support: RVO:61388998 Keywords : numerical simulation * transonic flow * two-equation model Subject RIV: BK - Fluid Dynamics

  2. Numerical Modeling of Shoreline Undulations

    Kærgaard, Kasper Hauberg

    length of the shoreline undulations is determined in the linear regime using a shoreline stability analysis based on the numerical model. The analysis shows that the length of the undulations in the linear regime depends on the incoming wave conditions and on the coastal profile. For larger waves and...... relationship between the shoreline undulations and longshore bars and the relationship between the morphology and the hydrodynamics. In one of the data sets the shoreline undulations are well correlated with undulations on the depth contours between -5 m and +2 m relative to mean sea level. An analysis of the...

  3. NUMERICAL MODELING OF THE RIVER POLLUTION

    M. M. Biliaiev

    2009-04-01

    Full Text Available The 2D numerical model to simulate the pollutant dispersion in rivers is offered. The model is based on the equation of potential flow and the transport model. The results of numerical experiment are presented.

  4. Numerical modeling of water waves

    Lin, Pengzhi

    2008-01-01

    Modelling large-scale wave fields and their interaction with coastal and offshore structures has become much more feasible over the last two decades with increases in computer speeds. Wave modelling can be viewed as an extension of wave theory, a mature and widely published field, applied to practical engineering through the use of computer tools. Information about the various wave models which have been developed is often widely scattered in the literature, and consequently this is one of the first books devoted to wave models and their applications. At the core of the book is an introduction to various types of wave models. For each model, the theoretical assumptions, the application range, and the advantages and limitations are elaborated. The combined use of different wave models from large-scale to local-scale is highlighted with a detailed discussion of the application and matching of boundary conditions. At the same time the book provides a grounding in hydrodynamics, wave theory, and numerical methods...

  5. Numerical methods used in fusion science numerical modeling

    Yagi, M.

    2015-04-01

    The dynamics of burning plasma is very complicated physics, which is dominated by multi-scale and multi-physics phenomena. To understand such phenomena, numerical simulations are indispensable. Fundamentals of numerical methods used in fusion science numerical modeling are briefly discussed in this paper. In addition, the parallelization technique such as open multi processing (OpenMP) and message passing interface (MPI) parallel programing are introduced and the loop-level parallelization is shown as an example.

  6. Numerical Modeling of Ocean Circulation

    Miller, Robert N.

    2007-01-01

    The modelling of ocean circulation is important not only for its own sake, but also in terms of the prediction of weather patterns and the effects of climate change. This book introduces the basic computational techniques necessary for all models of the ocean and atmosphere, and the conditions they must satisfy. It describes the workings of ocean models, the problems that must be solved in their construction, and how to evaluate computational results. Major emphasis is placed on examining ocean models critically, and determining what they do well and what they do poorly. Numerical analysis is introduced as needed, and exercises are included to illustrate major points. Developed from notes for a course taught in physical oceanography at the College of Oceanic and Atmospheric Sciences at Oregon State University, this book is ideal for graduate students of oceanography, geophysics, climatology and atmospheric science, and researchers in oceanography and atmospheric science. Features examples and critical examination of ocean modelling and results Demonstrates the strengths and weaknesses of different approaches Includes exercises to illustrate major points and supplement mathematical and physical details

  7. Numerical modelling of fuel sprays

    Bergstroem, C.

    1999-06-01

    The way the fuel is introduced into the combustion chamber is one of the most important parameters for the power output and the generation of emissions in the combustion of liquid fuels. The interaction between the turbulent gas flow field and the liquid fuel droplets, the vaporisation of them and the mixing of the gaseous fuel with the ambient air that are vital parameters in the combustion process. The use of numerical calculations is an important tool to better understand these complex interacting phenomena. This thesis reports on the numerical modelling of fuel sprays in non-reacting cases using an own developed spray module. The spray module uses the stochastic parcel method to represent the spray. The module was made in such manner that it could by coupled with different gas flow solver. Results obtained from four different gas flow solvers are presented in the thesis, including the use of two different kinds of turbulence models. In the first part the spray module is coupled with a k-{eta} based 2-D cylindrical gas flow solver. A thorough sensitivity analysis was performed on the spray and gas flow solver parameters, such as grid size dependence and sensitivity to initial values of k-{eta}. The results of the spray module were also compared to results from other spray codes, e.g. the well known KIVA code. In the second part of this thesis the spray was injected into a turbulent and fully developed crossflow studied. The spray module was attached to a LES (Large Eddy Simulation) based flow solvers enabling the study of the complex structures and time dependent phenomena involved in spray in crossflows. It was found that the spray performs an oscillatory motion and that the Strouhal number in the wake was about 0.1. Different spray breakup models were evaluated by comparing with experimental results 66 refs, 56 figs

  8. Numerical Modeling of Nanoelectronic Devices

    Klimeck, Gerhard; Oyafuso, Fabiano; Bowen, R. Chris; Boykin, Timothy

    2003-01-01

    Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.

  9. Experimental and numerical analysis of a knee endoprosthesis numerical model

    L. Zach

    2016-07-01

    Full Text Available The aim of this study is to create and verify a numerical model for a Medin Modular orthopedic knee-joint implant by investigating contact pressure, its distribution and contact surfaces. An experiment using Fuji Prescale pressure sensitive films and a finite element analysis (FEA using Abaqus software were carried out. The experimental data were evaluated using a special designed program and were compared with the results of the analysis. The designed evaluation program had been constructed on the basis of results obtained from a supplementary calibration experiment. The applicability of the numerical model for the real endoprosthesis behavior prediction was proven on the basis of their good correlation.

  10. Application of numerical models and codes

    Vyzikas, Thomas

    2014-01-01

    This report indicates the importance of numerical modelling in the modelling process, gradually builds the essential background theory in the fields of fluid mechanics, wave mechanics and numerical modelling, discusses a list of commonly used software and finally recommends which models are more suitable for different engineering applications in a marine renewable energy project.

  11. Numerical modeling of slow shocks

    This paper reviews previous attempt and the present status of efforts to understand the structure of slow shocks by means of time dependent numerical calculations. Studies carried out using MHD or hybrid-kinetic codes have demonstrated qualitative agreement with theory. A number of unresolved issues related to hybrid simulations of the internal shock structure are discussed in some detail. 43 refs., 8 figs

  12. Explicit BCJ numerators of nonlinear sigma model

    Du, Yi-Jian

    2016-01-01

    In this paper, we investigate the color-kinematics duality in nonlinear sigma model (NLSM). We present explicit polynomial expressions for the kinematic numerators (BCJ numerators). The calculation is done separately in two parametrization schemes of the theory using Kawai-Lewellen-Tye relation inspired technique, both lead to polynomial numerators. We summarize the calculation in each case into a set of rules that generates BCJ numerators for all multilplicities. In Cayley parametrization we find the numerator is described by a particularly simple formula solely in terms of momentum kernel.

  13. Numerical Modeling of Microelectrochemical Systems

    Adesokan, Bolaji James

    three independent papers and manuscripts. As a preliminary to the study, we describe a general model for electrochemical systems and study their underlying mechanisms through electroanalytical techniques. We then extend the model to a more realistic model for microelectrochemical systems which......The PhD dissertation is concerned with mathematical modeling and simulation of electrochemical systems. The first three chapters of the thesis consist of the introductory part, the model development chapter and the chapter on the summary of the main results. The remaining three chapters report...... the reactants in the bulk electrolyte that are traveling waves. The first paper presents the mathematical model which describes an electrochemical system and simulates an electroanalytical technique called cyclic voltammetry. The model is governed by a system of advection–diffusion equations with a...

  14. A combined 2-dimensional fuzzy regression model to study effect of climate change on the electrical peak load

    Hamed Shakouri G.

    2010-06-01

    Full Text Available This paper studies the impact of climate change on the electricity consumption by means of a fuzzy regression approach. The climate factors which have been considered in this paper are humidity and temperature, whereas the simultaneous effect of these two climate factors is considered. The impacts of other climate variables, like the wind, with a minor effect on energy consumption are ignored. The innovation which applies in this paper is the division of the year into two parts by using the temperature-day graph in the year. To index the humidity, data of the minimum humidity per day are used. For temperature, the maximum temperature of the first part of the year (warm days and the minimum of the second part (cold days are used. The indicator for the consumption is the daily peak load. The model results show high sensitivity to the temperature but low sensitivity to the humidity. Moreover, it is concluded that the model structure cannot be the same and for the cold par additional variables such as gas consumption should be considered.

  15. Survey of numerical electrostimulation models

    Reilly, J. Patrick

    2016-06-01

    This paper evaluates results of a survey of electrostimulation models of myelinated nerve. Participants were asked to determine thresholds of excitation for 18 cases involving different characteristics of the neuron, the stimulation waveform, and the electrode arrangement. Responses were received from 7 investigators using 10 models. Excitation thresholds differed significantly among these models. For example, with a 2 ms monophasic stimulus pulse and an electrode/fiber distance of 1 cm, thresholds from the least to greatest value differed by a factor of 8.3; with a 5 μs pulse, thresholds differed by the factor 3.8. Significant differences in reported simulations point to the need for experimental validation. Additional efforts are needed to develop computational models for unmyelinated C-fibers, A-delta fibers, CNS neurons, and CNS Synapses.

  16. Numerical Modelling of Jets and Plumes

    Larsen, Torben

    1993-01-01

    An overview on numerical models for prediction of the flow and mixing processes in turbulent jets and plumes is given. The overview is structured to follow an increasing complexity in the physical and numerical principles. The various types of models are briefly mentioned, from the one-dimensiona......An overview on numerical models for prediction of the flow and mixing processes in turbulent jets and plumes is given. The overview is structured to follow an increasing complexity in the physical and numerical principles. The various types of models are briefly mentioned, from the one......-dimensional integral method to the general 3-dimensional solution of the Navier-Stokes equations. Also the predictive capabilities of the models are discussed. The presentation takes the perspective of civil engineering and covers issues like sewage outfalls and cooling water discharges to the sea....

  17. Numerical models of supershell dynamics

    Superbubbles play an important role in determining the state of the ISM in both spiral and irregular galaxies. The authors are modeling supershell dynamics in both homogeneous and stratified atmospheres using ZEUS, a 2-D hydrocode. They find that when a superbubble blows out of a Gaussian atmosphere, the cold, dense shell is not greatly accelerated.In addition, the authors believe that they observe the Vishniac overstability in radiative, decelerating shells

  18. Numerical models of groundwater flow and transport

    This chapter reviews the state-of-the-art in deterministic modeling of groundwater flow and transport processes, which can be used for interpretation of isotope data through groundwater flow analyses. Numerical models which are available for this purpose are described and their applications to complex field problems are discussed. The theoretical bases of deterministic modeling are summarized, and advantages and limitations of numerical models are described. The selection of models for specific applications and their calibration procedures are described, and results of a few illustrative case study type applications are provided. (author). 145 refs, 17 figs, 2 tabs

  19. Numerical 3-D Modelling of Overflows

    Larsen, Torben; Nielsen, L.; Jensen, B.; Christensen, E. D.

    2008-01-01

    -dimensional so-called Volume of Fluid Models (VOF-models) based on the full Navier-Stokes equations (named NS3 and developed by DHI Water & Environment) As a general conclusion, the two numerical models show excellent results when compared with measurements. However, considerable errors occur when......The present study uses laboratory experiments to evaluate the reliability of two types of numerical models of sewers systems: - 1-dimensional model based on the extended Saint-Venant equation including the term for curvature of the water surface (the so-called Boussinesq approximation) - 2- and 3...

  20. Numerical modelling of elastic space tethers

    Kristiansen, Kristian Uldall; Palmer, P. L.; Roberts, R. M.

    2012-01-01

    , the numerical experiments of an orbiting tether system show that bending may introduce significant forces in some regions of phase space. Finally, numerical evidence for the existence of an almost invariant slow manifold of the singularly perturbed, regularised, non-dissipative massive tether model is provided....... It is also shown that on the slow manifold the dynamics of the satellites are well-approximated by the finite dimensional slack-spring model....

  1. Wave Numerical Model for Shallow Water

    徐福敏; 严以新; 张长宽; 宋志尧; 茅丽华

    2000-01-01

    The history of forecasting wind waves by wave energy conservation equation is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave energy conservation models for the simulation of shallow water waves are introduced,with emphasis placed on the SWAN model, which takes use of the most advanced wave research achievements and has been applied to several theoretical and field conditions. The characteristics and applicability of the model, the finite difference numerical scheme of the action balance equation and its source terms computing methods are described in detail. The model has been verified with the propagation refraction numerical experiments for waves propagating in following and opposing currents; finally, the model is applied to the Haian Gulf area to simulate the wave height and wave period field there, and the results are compared with observed data.

  2. Numerical modeling of fluidized-bed combustor

    Sha, W T; Soo, S L

    1977-11-01

    Optimum design of fluidized-bed combustor requires high carbon burn-up, good sulfur retention, minimized sorbent (Ca) utilization, efficient feed distribution and mechanical layout. These parameters are strongly affected by the dynamics of the fluidized bed. The dynamic behavior of fluidized combustor is formulated in terms of multidomain - multiphase mechanics. Fluidization, bubble mechanics, coal combustion, sorbent sulfation, oxidation, solids movement and elutriation, and heat transfer are explicitly taken into account in the proposed numerical model. The model solves conservation equations of mass, momentum and energy coupled with chemical reactions as boundary value problem in space and initial value problem in time. Multi-fluid model and modified implicit multi-field numerical scheme are employed. The objective of this numerical model is for use in engineering design and scaling. Progress to date shows that all necessary relations can be incorporated within the framework of an overall multidomain - multiphase model for deterministic computation. Provisions are made for subsequent refinements of submodels of individual mechanism and improvements of the existing numerical model. These refinements and improvements can be achieved as better understanding of physical phenomena and more experimental data become available. The numerical model outlined in this report is specifically designed for the fluidized-bed combustor; however, it can readily be extended to various coal gasification systems.

  3. Ferrofluids: Modeling, numerical analysis, and scientific computation

    Tomas, Ignacio

    This dissertation presents some developments in the Numerical Analysis of Partial Differential Equations (PDEs) describing the behavior of ferrofluids. The most widely accepted PDE model for ferrofluids is the Micropolar model proposed by R.E. Rosensweig. The Micropolar Navier-Stokes Equations (MNSE) is a subsystem of PDEs within the Rosensweig model. Being a simplified version of the much bigger system of PDEs proposed by Rosensweig, the MNSE are a natural starting point of this thesis. The MNSE couple linear velocity u, angular velocity w, and pressure p. We propose and analyze a first-order semi-implicit fully-discrete scheme for the MNSE, which decouples the computation of the linear and angular velocities, is unconditionally stable and delivers optimal convergence rates under assumptions analogous to those used for the Navier-Stokes equations. Moving onto the much more complex Rosensweig's model, we provide a definition (approximation) for the effective magnetizing field h, and explain the assumptions behind this definition. Unlike previous definitions available in the literature, this new definition is able to accommodate the effect of external magnetic fields. Using this definition we setup the system of PDEs coupling linear velocity u, pressure p, angular velocity w, magnetization m, and magnetic potential ϕ We show that this system is energy-stable and devise a numerical scheme that mimics the same stability property. We prove that solutions of the numerical scheme always exist and, under certain simplifying assumptions, that the discrete solutions converge. A notable outcome of the analysis of the numerical scheme for the Rosensweig's model is the choice of finite element spaces that allow the construction of an energy-stable scheme. Finally, with the lessons learned from Rosensweig's model, we develop a diffuse-interface model describing the behavior of two-phase ferrofluid flows and present an energy-stable numerical scheme for this model. For a

  4. Numerical methods and modelling for engineering

    Khoury, Richard

    2016-01-01

    This textbook provides a step-by-step approach to numerical methods in engineering modelling. The authors provide a consistent treatment of the topic, from the ground up, to reinforce for students that numerical methods are a set of mathematical modelling tools which allow engineers to represent real-world systems and compute features of these systems with a predictable error rate. Each method presented addresses a specific type of problem, namely root-finding, optimization, integral, derivative, initial value problem, or boundary value problem, and each one encompasses a set of algorithms to solve the problem given some information and to a known error bound. The authors demonstrate that after developing a proper model and understanding of the engineering situation they are working on, engineers can break down a model into a set of specific mathematical problems, and then implement the appropriate numerical methods to solve these problems. Uses a “building-block” approach, starting with simpler mathemati...

  5. Numerical modelling of Bose-Einstein correlations

    Utyuzh, O. V.; Wilk, G.; Wlodarczyk, Z.

    2001-01-01

    We propose extension of the algorithm for numerical modelling of Bose-Einstein correlations (BEC), which was presented some time ago in the literature. It is formulated on quantum statistical level for a single event and uses the fact that identical particles subjected to Bose statistics do bunch themselves, in a maximal possible way, in the same cells in phase-space. The bunching effect is in our case obtained in novel way allowing for broad applications and fast numerical calculations. Firs...

  6. Simple numerical techniques for mesoscale polymer models

    E.A. Koopman

    2014-01-01

    There exist various numerical techniques for modeling polymer behavior on small time and length scales. This work introduces some new techniques in this field, and shows novel new combinations of existing techniques. Among the new techniques are multiple new thermostats, a way of modeling polymers i

  7. Amorphous track models: A numerical comparison study

    Greilich, Steffen; Grzanka, L.; Bassler, N.; Andersen, Claus Erik; Jäkel, O.

    2010-01-01

    We present an open-source code library for amorphous track modelling which is suppose to faciliate the application and numerical comparability as well as serve as a frame-work for the implementation of new models. We show an example of using the library indicating the choice of submodels has a...

  8. Numerical FEM modeling in dental implantology

    Roateşi, Iulia; Roateşi, Simona

    2016-06-01

    This paper is devoted to a numerical approach of the stress and displacement calculation of a system made up of dental implant, ceramic crown and surrounding bone. This is the simulation of a clinical situation involving both biological - the bone tissue, and non-biological - the implant and the crown, materials. On the other hand this problem deals with quite fine technical structure details - the threads, tapers, etc with a great impact in masticatory force transmission. Modeling the contact between the implant and the bone tissue is important to a proper bone-implant interface model and implant design. The authors proposed a three-dimensional numerical model to assess the biomechanical behaviour of this complex structure in order to evaluate its stability by determining the risk zones. A comparison between this numerical analysis and clinical cases is performed and a good agreement is obtained.

  9. Fundamentals of Numerical Modelling of Casting Processes

    Pryds, Nini; Thorborg, Jesper; Lipinski, Marek;

    Fundamentals of Numerical Modelling of Casting Processes comprises a thorough presentation of the basic phenomena that need to be addressed in numerical simulation of casting processes. The main philosophy of the book is to present the topics in view of their physical meaning, whenever possible......, rather than relying strictly on mathematical formalism. The book, aimed both at the researcher and the practicing engineer, as well as the student, is naturally divided into four parts. Part I (Chapters 1-3) introduces the fundamentals of modelling in a 1-dimensional framework. Part II (Chapter 4......) presents the most important aspects of solidification theory related to modelling. Part III (Chapter 5) describes the fluid flow phenomena and in part IV (Chapter 6) the stress-strain analysis is addressed. For all parts, both numerical formulations as well as some important analytical solutions...

  10. Numerical modelling of swirling diffusive flames

    Parra-Santos Teresa; Perez Ruben; Szasz Robert Z.; Gutkowski Artur N.; Castro Francisco

    2016-01-01

    Computational Fluid Dynamics has been used to study the mixing and combustion of two confined jets whose setup and operating conditions are those of the benchmark of Roback and Johnson. Numerical model solves 3D transient Navier Stokes for turbulent and reactive flows. Averaged velocity profiles using RNG swirl dominated k-epsilon model have been validated with experimental measurements from other sources for the non reactive case. The combustion model is Probability Density Function. Bearing...

  11. Numerical modelling of Bose-Einstein correlations

    Utyuzh, O V; Wlodarczyk, Z

    2001-01-01

    We propose extension of the algorithm for numerical modelling of Bose-Einstein correlations (BEC), which was presented some time ago in the literature. It is formulated on quantum statistical level for a single event and uses the fact that identical particles subjected to Bose statistics do bunch themselves, in a maximal possible way, in the same cells in phase-space. The bunching effect is in our case obtained in novel way allowing for broad applications and fast numerical calculations. First comparison with $e^+e^-$ annihilations data performed by using simple cascade hadronization model is very encouraging.

  12. High performance computing and numerical modelling

    ,

    2014-01-01

    Numerical methods play an ever more important role in astrophysics. This is especially true in theoretical works, but of course, even in purely observational projects, data analysis without massive use of computational methods has become unthinkable. The key utility of computer simulations comes from their ability to solve complex systems of equations that are either intractable with analytic techniques or only amenable to highly approximative treatments. Simulations are best viewed as a powerful complement to analytic reasoning, and as the method of choice to model systems that feature enormous physical complexity such as star formation in evolving galaxies, the topic of this 43rd Saas Fee Advanced Course. The organizers asked me to lecture about high performance computing and numerical modelling in this winter school, and to specifically cover the basics of numerically treating gravity and hydrodynamics in the context of galaxy evolution. This is still a vast field, and I necessarily had to select a subset ...

  13. Basset force in numerical models of saltation

    Lukerchenko, Nikolay; Dolanský, Jindřich; Vlasák, Pavel

    2012-01-01

    Roč. 60, č. 4 (2012), s. 277-287. ISSN 0042-790X R&D Projects: GA ČR GA103/09/1718 Institutional research plan: CEZ:AV0Z20600510 Keywords : basset force * bed load transport * numerical model * particle-bed collision Subject RIV: BK - Fluid Dynamics Impact factor: 0.653, year: 2012

  14. Numerical modelling of steel arc welding

    Welding is a highly used assembly technique. Welding simulation software would give access to residual stresses and information about the weld's microstructure, in order to evaluate the mechanical resistance of a weld. It would also permit to evaluate the process feasibility when complex geometrical components are to be made, and to optimize the welding sequences in order to minimize defects. This work deals with the numerical modelling of arc welding process of steels. After describing the industrial context and the state of art, the models implemented in TransWeld (software developed at CEMEF) are presented. The set of macroscopic equations is followed by a discussion on their numerical implementation. Then, the theory of re-meshing and our adaptive anisotropic re-meshing strategy are explained. Two welding metal addition techniques are investigated and are compared in terms of the joint size and transient temperature and stresses. The accuracy of the finite element model is evaluated based on experimental results and the results of the analytical solution. Comparative analysis between experimental and numerical results allows the assessment of the ability of the numerical code to predict the thermomechanical and metallurgical response of the welded structure. The models limitations and the phenomena identified during this study are finally discussed and permit to define interesting orientations for future developments. (author)

  15. Numerical modeling of magma-repository interactions

    Bokhove, O.

    2001-01-01

    This report explains the numerical programs behind a comprehensive modeling effort of magma-repository interactions. Magma-repository interactions occur when a magma dike with high-volatile content magma ascends through surrounding rock and encounters a tunnel or drift filled with either a magmatic gas or air at atmospheric pressure. The simplified mathematical model to describe these flow phenomena consists of compressible flow equations of one- or two-dimensional flow in a flow tube with a ...

  16. Numerical modelling in wave energy conversion systems

    El Marjani, A. [Labo. de Turbomachines, Ecole Mohammadia d' Ingenieurs (EMI), Universite Mohammed V Agdal, Av Ibn Sina, B.P. 765 Agdal, Rabat (Morocco); Castro Ruiz, F.; Rodriguez, M.A.; Parra Santos, M.T. [Depto. de Ingenieria Energetica y Fluidomecanica, Escuela Tecnica Superior de Ingenieros Industriales, Universidad de Valladolid, Paseo del Cauce s/n, E-47011 Valladolid (Spain)

    2008-08-15

    This paper deals with a numerical modelling devoted to predict the flow characteristics in the components of an oscillating water column (OWC) system used for the wave energy capture. In the present paper, the flow behaviour is modelled by using the FLUENT code. Two numerical flow models have been elaborated and tested independently in the geometries of an air chamber and a turbine, which is chosen of a radial impulse type. The flow is assumed to be three-dimensional (3D), viscous, turbulent and unsteady. The FLUENT code is used with a solver of the coupled conservation equations of mass, momentum and energy, with an implicit time scheme and with the adoption of the dynamic mesh and the sliding mesh techniques in areas of moving surfaces. Turbulence is modelled with the k-{epsilon} model. The obtained results indicate that the developed models are well suitable to analyse the air flows both in the air chamber and in the turbine. The performances associated with the energy transfer processes have been well predicted. For the turbine, the numerical results of pressure and torque were compared to the experimental ones. Good agreements between these results have been observed. (author)

  17. Lattice Boltzmann Model for Numerical Relativity

    Ilseven, E

    2015-01-01

    In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.

  18. Numerical modeling of fires on gas pipelines

    When natural gas is released through a hole on a high-pressure pipeline, it disperses in the atmosphere as a jet. A jet fire will occur when the leaked gas meets an ignition source. To estimate the dangerous area, the shape and size of the fire must be known. The evolution of the jet fire in air is predicted by using a finite-volume procedure to solve the flow equations. The model is three-dimensional, elliptic and calculated by using a compressibility corrected version of the k - ξ turbulence model, and also includes a probability density function/laminar flamelet model of turbulent non-premixed combustion process. Radiation heat transfer is described using an adaptive version of the discrete transfer method. The model is compared with the experiments about a horizontal jet fire in a wind tunnel in the literature with success. The influence of wind and jet velocity on the fire shape has been investigated. And a correlation based on numerical results for predicting the stoichiometric flame length is proposed. - Research highlights: → We developed a model to predict the evolution of turbulent jet diffusion flames. → Measurements of temperature distributions match well with the numerical predictions. → A correlation has been proposed to predict the stoichiometric flame length. → Buoyancy effects are higher in the numerical results. → The radiative heat loss is bigger in the experimental results.

  19. Numerical analysis of the rebellious voter model

    Swart, Jan M.; Vrbenský, Karel

    2010-01-01

    Roč. 140, č. 5 (2010), s. 873-899. ISSN 0022-4715 R&D Projects: GA ČR GA201/09/1931; GA MŠk 1M0572 Institutional research plan: CEZ:AV0Z10750506 Keywords : rebellious voter model * parity conservation * exactly solvable model * coexistence * interface tightness * cancellative systems * Markov chain Monte Carlo Subject RIV: BA - General Mathematics Impact factor: 1.447, year: 2010 http://library.utia.cas.cz/separaty/2010/SI/swart-numerical analysis of the rebellious voter model.pdf

  20. Economic growth models: symbolic and numerical computations

    P. B. Vasconcelos

    2013-11-01

    Full Text Available Growth economic models play a crucial role in understanding countries development, inter-country macroeconomic relationship and, ultimately, to anticipate the effects on endogenous variables due to political shocks on model parameters or exogenous variables. Yet, build a mathematical model can be a difficult and time consuming task. Symbolic computations can be of great help in the development process. Then, the ability to simulate, under initial assumptions, is, assuredly, a priceless tool for policy makers to take decisions and to adapt them along the time. Robust and efficient solvers are required to allow for reliable and fast answers. This paper deals with an integrated computational approach to economic growth models, based on the exceptional ability of MATLAB's numerical computing and exploring its symbolic computing capabilities. Illustration is made with the Ramsey-Cass-Koopmans model, one of the macroeconomic workhorse models.

  1. A Numerical Model for Atomtronic Circuit Analysis

    Chow, Weng W. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Straatsma, Cameron J. E. [Univ. of Colorado, Boulder, CO (United States); Anderson, Dana Z. [Univ. of Colorado and National Inst. of Standards and Technology, Boulder, CO (United States)

    2015-07-16

    A model for studying atomtronic devices and circuits based on finite-temperature Bose-condensed gases is presented. The approach involves numerically solving equations of motion for atomic populations and coherences, derived using the Bose-Hubbard Hamiltonian and the Heisenberg picture. The resulting cluster expansion is truncated at a level giving balance between physics rigor and numerical demand mitigation. This approach allows parametric studies involving time scales that cover both the rapid population dynamics relevant to nonequilibrium state evolution, as well as the much longer time durations typical for reaching steady-state device operation. This model is demonstrated by studying the evolution of a Bose-condensed gas in the presence of atom injection and extraction in a double-well potential. In this configuration phase locking between condensates in each well of the potential is readily observed, and its influence on the evolution of the system is studied.

  2. Some Experiences with Numerical Modelling of Overflows

    Larsen, Torben; Nielsen, L.; Jensen, B.;

    2007-01-01

    Overflows are commonly applied in storm sewer systems to control flow and water surface level. Therefore overflows play a central role in the control of discharges of pollutants from sewer systems to the environment. The basic hydrodynamic principle of an overflow is the so-called critical flow...... across the edge of the overflow. To ensure critical flow across the edge, the upstream flow must be subcritical whereas the downstream flow is either supercritical or a free jet. Experimentally overflows are well studied. Based on laboratory experiments and Froude number scaling, numerous accurate and...... including the term for curvature of the water surface (the so-called Boussinesq approximation) 2. 2- and 3-dimensional so-called Volume of Fluid Models (VOF-models) based on the full Navier-Stokes equations (named NS3 and developed by DHI Water & Environment) As a general conclusion, the two numerical...

  3. Lattice Boltzmann model for numerical relativity

    Ilseven, E.; Mendoza, M.

    2016-02-01

    In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.

  4. Lattice Boltzmann model for numerical relativity.

    Ilseven, E; Mendoza, M

    2016-02-01

    In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems. PMID:26986435

  5. Numerical modeling in materials science and engineering

    Rappaz, Michel; Deville, Michel

    2003-01-01

    This book introduces the concepts and methodologies related to the modelling of the complex phenomena occurring in materials processing. After a short reminder of conservation laws and constitutive relationships, the authors introduce the main numerical methods: finite differences, finite volumes and finite elements. These techniques are developed in three main chapters of the book that tackle more specific problems: phase transformation, solid mechanics and fluid flow. The two last chapters treat inverse methods to obtain the boundary conditions or the material properties and stochastic methods for microstructural simulation. This book is intended for undergraduate and graduate students in materials science and engineering, mechanical engineering and physics and for engineering professionals or researchers who want to get acquainted with numerical simulation to model and compute materials processing.

  6. Numerical model for atomtronic circuit analysis

    Chow, Weng W; Anderson, Dana Z

    2015-01-01

    A model for studying atomtronic devices and circuits based on finite temperature Bose-condensed gases is presented. The approach involves numerically solving equations of motion for atomic populations and coherences, derived using the Bose-Hubbard Hamiltonian and the Heisenberg picture. The resulting cluster expansion is truncated at a level giving balance between physics rigor and numerical demand mitigation. This approach allows parametric studies involving time scales that cover both the rapid population dynamics relevant to non-equilibrium state evolution, as well as the much longer time durations typical for reaching steady-state device operation. The model is demonstrated by studying the evolution of a Bose-condensed gas in the presence of atom injection and extraction in a double-well potential. In this configuration phase-locking between condensates in each well of the potential is readily observed, and its influence on the evolution of the system is studied.

  7. Numerical Modeling and Investigation of Boiling Phenomena

    Kunkelmann, Christian

    2011-01-01

    The subject of the present thesis is the numerical modeling and investigation of boiling phenomena. The heat transfer during boiling is highly efficient and therefore used for many applications in power generation, process engineering and cooling of high performance electronics. The precise knowledge of particular boiling processes, their relevant parameters and limitations is of utmost importance for an optimized application. Therefore, the fundamentals of boiling heat transfer have been...

  8. Hydrodynamic numerical modelling of Maputo Bay

    Canhanga, S.J.V.

    2003-01-01

    The aim of this work was to understand by using a numerical model, the main characteristics of the Maputo Bay hydrodynamics; and furthermore to assess the hydrodynamics implications on the different processes that can occur in the Bay. The study was initiated with a description of the hydrodynamical characteristics of Maputo Bay, through the analysis of vertical profiles of salinity and temperature, time series of the components of tide currents and water elevation. The wind intensity and ...

  9. Numerical Modeling of Laterally Loaded Piles

    S. T. Kok; B. B.k. Huat

    2008-01-01

    Design of laterally loaded piles due to soil movement relies on a number of theoretical and numerical approaches. However, the magnitude of soil movement is difficult to estimate with reasonable confidence and accuracy. Finite Element Analysis (FEA) offers an excellent alternative to study pile-soil interaction and piles response under lateral loading due to soil movement. This research presents published analytical results and case history modeled in a 2D finite element environment in the ca...

  10. A numerical model of coastal overwash

    Donnelly, Chantal; Hanson, Hans; Larson, Magnus

    2009-01-01

    Overwash, the flow of water and sediment over the crest of a beach, contributes to flooding and the deposition of sand landward of the beach crest. Washover, the sand deposited by overwash, contributes to the sediment budget and migration of barrier islands. The ability to predict the occurrence, location, and thickness of overwash deposits is important for coastal residents, coastal town planners, environmental planners, and engineers alike. In this study, a numerical model that simulates th...

  11. Avoiding numerical pitfalls in social force models.

    Köster, Gerta; Treml, Franz; Gödel, Marion

    2013-06-01

    The social force model of Helbing and Molnár is one of the best known approaches to simulate pedestrian motion, a collective phenomenon with nonlinear dynamics. It is based on the idea that the Newtonian laws of motion mostly carry over to pedestrian motion so that human trajectories can be computed by solving a set of ordinary differential equations for velocity and acceleration. The beauty and simplicity of this ansatz are strong reasons for its wide spread. However, the numerical implementation is not without pitfalls. Oscillations, collisions, and instabilities occur even for very small step sizes. Classic solution ideas from molecular dynamics do not apply to the problem because the system is not Hamiltonian despite its source of inspiration. Looking at the model through the eyes of a mathematician, however, we realize that the right hand side of the differential equation is nondifferentiable and even discontinuous at critical locations. This produces undesirable behavior in the exact solution and, at best, severe loss of accuracy in efficient numerical schemes even in short range simulations. We suggest a very simple mollified version of the social force model that conserves the desired dynamic properties of the original many-body system but elegantly and cost efficiently resolves several of the issues concerning stability and numerical resolution. PMID:23848804

  12. Numerical modeling of explosions for nuclear monitoring

    Stevens, J. L.

    2014-12-01

    Monitoring the Earth for underground nuclear explosions requires a detailed understanding of the explosion source. In this context, "source" refers to the source of seismic waves, and it is generated by the complex nonlinear near-source motion that accompanies the nuclear explosion. In particular, nuclear monitoring requires understanding the transition from the hydrodynamic to elastic regimes, and propagation of waveforms from the source to stations at distances of hundreds to thousands of kilometers. In the transition region, shear strength is critically important, as are changes in shear strength as the shock wave propagates. Numerical modeling using 1D spherically symmetric, 2D axisymmetric and full 3D calculations provides important insights into the seismic source and the waveforms it generates. Important considerations for numerical modeling include emplacement conditions (tamped or in a cavity), source type (chemical or nuclear), material models for strength and strength reduction, and geologic conditions including topography and tectonic stresses in the source region. In addition to calculating the near source ground motion, we propagate the near source solution to regional and teleseismic distances where the observations of seismic signals from nuclear explosions are made. The objectives of nuclear monitoring are detection of seismic events (earthquakes, quarry blasts and other sources in addition to nuclear explosions), accurate location of these events, discrimination of nuclear explosions from other types of sources, and estimation of nuclear explosion yield. Numerical modeling is particularly important for discrimination and yield estimation. Numerical modeling is used to understand unexpected anomalies that occur, such as the large surface waves generated by the three North Korean nuclear tests, which may have been caused by a difference in tectonic stress state between North Korea and other test sites. Another important issue that can be addressed

  13. Advanced Numerical Model for Irradiated Concrete

    Giorla, Alain B. [Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)

    2015-03-01

    In this report, we establish a numerical model for concrete exposed to irradiation to address these three critical points. The model accounts for creep in the cement paste and its coupling with damage, temperature and relative humidity. The shift in failure mode with the loading rate is also properly represented. The numerical model for creep has been validated and calibrated against different experiments in the literature [Wittmann, 1970, Le Roy, 1995]. Results from a simplified model are shown to showcase the ability of numerical homogenization to simulate irradiation effects in concrete. In future works, the complete model will be applied to the analysis of the irradiation experiments of Elleuch et al. [1972] and Kelly et al. [1969]. This requires a careful examination of the experimental environmental conditions as in both cases certain critical information are missing, including the relative humidity history. A sensitivity analysis will be conducted to provide lower and upper bounds of the concrete expansion under irradiation, and check if the scatter in the simulated results matches the one found in experiments. The numerical and experimental results will be compared in terms of expansion and loss of mechanical stiffness and strength. Both effects should be captured accordingly by the model to validate it. Once the model has been validated on these two experiments, it can be applied to simulate concrete from nuclear power plants. To do so, the materials used in these concrete must be as well characterized as possible. The main parameters required are the mechanical properties of each constituent in the concrete (aggregates, cement paste), namely the elastic modulus, the creep properties, the tensile and compressive strength, the thermal expansion coefficient, and the drying shrinkage. These can be either measured experimentally, estimated from the initial composition in the case of cement paste, or back-calculated from mechanical tests on concrete. If some

  14. Efficient numerical integrators for stochastic models

    De Fabritiis, G; Español, P; Coveney, P V

    2006-01-01

    The efficient simulation of models defined in terms of stochastic differential equations (SDEs) depends critically on an efficient integration scheme. In this article, we investigate under which conditions the integration schemes for general SDEs can be derived using the Trotter expansion. It follows that, in the stochastic case, some care is required in splitting the stochastic generator. We test the Trotter integrators on an energy-conserving Brownian model and derive a new numerical scheme for dissipative particle dynamics. We find that the stochastic Trotter scheme provides a mathematically correct and easy-to-use method which should find wide applicability.

  15. Numerical modelling of swirling diffusive flames

    Parra-Santos Teresa

    2016-01-01

    Full Text Available Computational Fluid Dynamics has been used to study the mixing and combustion of two confined jets whose setup and operating conditions are those of the benchmark of Roback and Johnson. Numerical model solves 3D transient Navier Stokes for turbulent and reactive flows. Averaged velocity profiles using RNG swirl dominated k-epsilon model have been validated with experimental measurements from other sources for the non reactive case. The combustion model is Probability Density Function. Bearing in mind the annular jet has swirl number over 0.5, a vortex breakdown appears in the axis of the burner. Besides, the sudden expansion with a ratio of 2 in diameter between nozzle exits and the test chamber produces the boundary layer separation with the corresponding torus shape recirculation. Contrasting the mixing and combustion models, the last one produces the reduction of the vortex breakdown.

  16. Numerical modelling of swirling diffusive flames

    Parra-Santos, Teresa; Perez, Ruben; Szasz, Robert Z.; Gutkowski, Artur N.; Castro, Francisco

    2016-03-01

    Computational Fluid Dynamics has been used to study the mixing and combustion of two confined jets whose setup and operating conditions are those of the benchmark of Roback and Johnson. Numerical model solves 3D transient Navier Stokes for turbulent and reactive flows. Averaged velocity profiles using RNG swirl dominated k-epsilon model have been validated with experimental measurements from other sources for the non reactive case. The combustion model is Probability Density Function. Bearing in mind the annular jet has swirl number over 0.5, a vortex breakdown appears in the axis of the burner. Besides, the sudden expansion with a ratio of 2 in diameter between nozzle exits and the test chamber produces the boundary layer separation with the corresponding torus shape recirculation. Contrasting the mixing and combustion models, the last one produces the reduction of the vortex breakdown.

  17. Numerical Investigation of Lattice Weinberg - Salam Model

    Zubkov, M. A.

    2013-11-01

    Lattice Weinberg - Salam model without fermions for the value of the Weinberg angle θW ~ 30°, and bare fine structure constant around α ˜ (1)/({150)} is investigated numerically. We consider the value of the scalar self coupling corresponding to bare Higgs mass around 150 GeV. We investigate phenomena existing in the vicinity of the phase transition between the physical Higgs phase and the unphysical symmetric phase of the lattice model. This is the region of the phase diagram, where the continuum physics is to be approached. We find the indications that at the energies above 1 TeV nonperturbative phenomena become important in the Weinberg - Salam model.

  18. Numerical modeling of the laser heated solenoid

    A numerical model of the interaction of laser radiation with the magnetically confined plasma in an infinitive solenoid was given. An approximate solution which includes the balance of total pressure instead of momentum equation was also developed. Thus, the time step in computing is not bounded by the characteristics given by the Alfven speed. This approximation approach makes the efficient computing of this problem possible. The results of the approximate solution agree very well with those of the exact solution. They have the same final steady state solution

  19. Numerical time integration for air pollution models

    Verwer, Jan; Hundsdorfer, Willem; Blom, Joke

    1998-01-01

    Due to the large number of chemical species and the three space dimensions, off-the-shelf stiff ODE integrators are not feasible for the numerical time integration of stiff systems of advection-diffusion-reaction equations [ fracpar{c{t + nabla cdot left( vu{u c right) = nabla cdot left( K ,nabla, c right) + R left( c right), c=c(vu{x,t), c in IR^m, vu{x in Omega subset IR^3 ] from the field of air pollution modelling. This has led to the use of special time integration techniques. This paper...

  20. Numerical modeling of atmospheric washout processes

    For the washout of particles from the atmosphere by clouds and rain one has to distinguish between processes which work in the first phase of cloud development, when condensation nuclei build up in saturated air (Nucleation Aerosol Scavenging, NAS) and those processes which work at the following cloud development. In the second case particles are taken off by cloud droplets or by falling rain drops via collision (Collision Aerosol Scavenging, CAS). The physics of both processes is described. For the CAS process a numerical model is presented. The report contains a documentation of the mathematical equations and the computer programs (FORTRAN). (KW)

  1. NUMERICAL MODEL APPLICATION IN ROWING SIMULATOR DESIGN

    Petr Chmátal

    2016-04-01

    Full Text Available The aim of the research was to carry out a hydraulic design of rowing/sculling and paddling simulator. Nowadays there are two main approaches in the simulator design. The first one includes a static water with no artificial movement and counts on specially cut oars to provide the same resistance in the water. The second approach, on the other hand uses pumps or similar devices to force the water to circulate but both of the designs share many problems. Such problems are affecting already built facilities and can be summarized as unrealistic feeling, unwanted turbulent flow and bad velocity profile. Therefore, the goal was to design a new rowing simulator that would provide nature-like conditions for the racers and provide an unmatched experience. In order to accomplish this challenge, it was decided to use in-depth numerical modeling to solve the hydraulic problems. The general measures for the design were taken in accordance with space availability of the simulator ́s housing. The entire research was coordinated with other stages of the construction using BIM. The detailed geometry was designed using a numerical model in Ansys Fluent and parametric auto-optimization tools which led to minimum negative hydraulic phenomena and decreased investment and operational costs due to the decreased hydraulic losses in the system.

  2. Numerical modeling of materials under extreme conditions

    Brown, Eric

    2014-01-01

    The book presents twelve state of the art contributions in the field of numerical modeling of materials subjected to large strain, high strain rates, large pressure and high stress triaxialities, organized into two sections. The first part is focused on high strain rate-high pressures such as those occurring in impact dynamics and shock compression related phenomena, dealing with material response identification, advanced modeling incorporating microstructure and damage, stress waves propagation in solids and structures response under impact. The latter part is focused on large strain-low strain rates applications such as those occurring in technological material processing, dealing with microstructure and texture evolution, material response at elevated temperatures, structural behavior under large strain and multi axial state of stress.

  3. Numerical Model of the DARHT Accelerating Cell

    Hughes, Thomas P; Genoni, Thomas C; Kang, Mike; Prichard, Benjamin A

    2005-01-01

    The DARHT-2 facility at Los Alamos National Laboratory accelerates a 2 microsecond electron beam using a series of inductive accelerating cells. The cell inductance is provided by large Metglas cores, which are driven by a pulse-forming network. The original cell design was susceptible to electrical breakdown near the outer radius of the cores. We developed a numerical model for the magnetic properties of Metglas over the range of dB/dt (magnetization rate) relevant to DARHT. The model was implemented in a radially-resolved circuit code, and in the LSP* electromagnetic code. LSP simulations showed that the field stress distribution across the outer radius of the cores was highly nonuniform. This was subsequently confirmed in experiments at LBNL. The calculated temporal evolution of the electric field stress inside the cores approximately matches experimental measurements. The cells have been redesigned to greatly reduce the field stresses along the outer radius.

  4. Partial Differential Equations Modeling and Numerical Simulation

    Glowinski, Roland

    2008-01-01

    This book is dedicated to Olivier Pironneau. For more than 250 years partial differential equations have been clearly the most important tool available to mankind in order to understand a large variety of phenomena, natural at first and then those originating from human activity and technological development. Mechanics, physics and their engineering applications were the first to benefit from the impact of partial differential equations on modeling and design, but a little less than a century ago the Schrödinger equation was the key opening the door to the application of partial differential equations to quantum chemistry, for small atomic and molecular systems at first, but then for systems of fast growing complexity. Mathematical modeling methods based on partial differential equations form an important part of contemporary science and are widely used in engineering and scientific applications. In this book several experts in this field present their latest results and discuss trends in the numerical analy...

  5. Adaptive numerical algorithms in space weather modeling

    Tóth, Gábor; van der Holst, Bart; Sokolov, Igor V.; De Zeeuw, Darren L.; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Najib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav

    2012-02-01

    Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit

  6. Adaptive numerical algorithms in space weather modeling

    Space weather describes the various processes in the Sun–Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit

  7. DANA: distributed numerical and adaptive modelling framework.

    Rougier, Nicolas P; Fix, Jérémy

    2012-01-01

    DANA is a python framework ( http://dana.loria.fr ) whose computational paradigm is grounded on the notion of a unit that is essentially a set of time dependent values varying under the influence of other units via adaptive weighted connections. The evolution of a unit's value are defined by a set of differential equations expressed in standard mathematical notation which greatly ease their definition. The units are organized into groups that form a model. Each unit can be connected to any other unit (including itself) using a weighted connection. The DANA framework offers a set of core objects needed to design and run such models. The modeler only has to define the equations of a unit as well as the equations governing the training of the connections. The simulation is completely transparent to the modeler and is handled by DANA. This allows DANA to be used for a wide range of numerical and distributed models as long as they fit the proposed framework (e.g. cellular automata, reaction-diffusion system, decentralized neural networks, recurrent neural networks, kernel-based image processing, etc.). PMID:22994650

  8. Submarine sand volcanos: experiments and numerical modelling

    Philippe, P.; Ngoma, J.; Delenne, J.

    2012-12-01

    Fluid overpressure at the bottom of a soil layer may generate fracturation in preferential paths for a cohesive material. But the case of sandy soils is rather different: a significant internal flow is allowed within the material and can potentially induce hydro-mechanical instabilities whose most common example is fluidization. Many works have been devoted to fluidization but very few have the issue of initiation and development of a fluidized zone inside a granular bed, prior entire fluidization of the medium. In this contribution, we report experimental results and numerical simulations on a model system of immersed sand volcanos generated by a localized upward spring of liquid, injected at constant flow-rate at the bottom of a granular layer. Such a localized state of fluidization is relevant for some industrial processes (spouted bed, maintenance of navigable waterways,…) and for several geological issues (kimberlite volcano conduits, fluid venting, oil recovery in sandy soil, More precisely, what is presented here is a comparison between experiments, carried out by direct visualization throughout the medium, and numerical simulations, based on DEM modelling of the grains coupled to resolution of NS equations in the liquid phase (LBM). There is a very good agreement between the experimental phenomenology and the simulation results. When the flow-rate is increased, three regimes are successively observed: static bed, fluidized cavity that does not extend to the top of the layer, and finally fluidization over the entire height of layer that creates a fluidized chimney. A very strong hysteretic effect is present here with an extended range of stability for fluidized cavities when flow-rate is decreased back. This can be interpreted in terms force chains and arches. The influences of grain diameter, layer height and injection width are studied and interpreted using a model previously developed by Zoueshtiagh [1]. Finally, growing rate of the fluidized zone and

  9. Numerical modelling of ion transport in flames

    Han, Jie

    2015-10-20

    This paper presents a modelling framework to compute the diffusivity and mobility of ions in flames. The (n, 6, 4) interaction potential is adopted to model collisions between neutral and charged species. All required parameters in the potential are related to the polarizability of the species pair via semi-empirical formulas, which are derived using the most recently published data or best estimates. The resulting framework permits computation of the transport coefficients of any ion found in a hydrocarbon flame. The accuracy of the proposed method is evaluated by comparing its predictions with experimental data on the mobility of selected ions in single-component neutral gases. Based on this analysis, the value of a model constant available in the literature is modified in order to improve the model\\'s predictions. The newly determined ion transport coefficients are used as part of a previously developed numerical approach to compute the distribution of charged species in a freely propagating premixed lean CH4/O2 flame. Since a significant scatter of polarizability data exists in the literature, the effects of changes in polarizability on ion transport properties and the spatial distribution of ions in flames are explored. Our analysis shows that changes in polarizability propagate with decreasing effect from binary transport coefficients to species number densities. We conclude that the chosen polarizability value has a limited effect on the ion distribution in freely propagating flames. We expect that the modelling framework proposed here will benefit future efforts in modelling the effect of external voltages on flames. Supplemental data for this article can be accessed at http://dx.doi.org/10.1080/13647830.2015.1090018. © 2015 Taylor & Francis.

  10. Numerical Modeling of Ocular Dysfunction in Space

    Nelson, Emily S.; Mulugeta, Lealem; Vera, J.; Myers, J. G.; Raykin, J.; Feola, A. J.; Gleason, R.; Samuels, B.; Ethier, C. R.

    2014-01-01

    Upon introduction to microgravity, the near-loss of hydrostatic pressure causes a marked cephalic (headward) shift of fluid in an astronaut's body. The fluid shift, along with other factors of spaceflight, induces a cascade of interdependent physiological responses which occur at varying time scales. Long-duration missions carry an increased risk for the development of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, a spectrum of ophthalmic changes including posterior globe flattening, choroidal folds, distension of the optic nerve sheath, kinking of the optic nerve and potentially permanent degradation of visual function. In the cases of VIIP found to date, the initial onset of symptoms occurred after several weeks to several months of spaceflight, by which time the gross bodily fluid distribution is well established. We are developing a suite of numerical models to simulate the effects of fluid shift on the cardiovascular, central nervous and ocular systems. These models calculate the modified mean volumes, flow rates and pressures that are characteristic of the altered quasi-homeostatic state in microgravity, including intracranial and intraocular pressures. The results of the lumped models provide initial and boundary data to a 3D finite element biomechanics simulation of the globe, optic nerve head and retrobulbar subarachnoid space. The integrated set of models will be used to investigate the evolution of the biomechanical stress state in the ocular tissues due to long-term exposure to microgravity.

  11. Numerical model of post-DNB film boiling heat transfer

    It is proposed in this paper a physical model for the film boiling heat transfer. The corresponding mathematical descriptions are given in details and the heat transfer characteristic of post-DNB film boiling is analyzed. The numerical model of post-DNB film boiling heat transfer is obtained as the empirical value of the coefficient is determined by the experimental data. The numerical model is compared with the experimental data of different parameters and other numerical models, and the statistical deviations are calculated. The calculating results of the numerical model in this paper show good agreement with the experimental data, and the numerical model in this paper has comprehensive applicability compared with other numerical models. The effects of thermal-hydraulic parameters on the post-DNB film boiling heat transfer have been numerically researched using the numerical model in this paper. The calculating results are as same as the experimental results. (authors)

  12. Modeling Biodegradation and Reactive Transport: Analytical and Numerical Models

    Sun, Y; Glascoe, L

    2005-06-09

    The computational modeling of the biodegradation of contaminated groundwater systems accounting for biochemical reactions coupled to contaminant transport is a valuable tool for both the field engineer/planner with limited computational resources and the expert computational researcher less constrained by time and computer power. There exists several analytical and numerical computer models that have been and are being developed to cover the practical needs put forth by users to fulfill this spectrum of computational demands. Generally, analytical models provide rapid and convenient screening tools running on very limited computational power, while numerical models can provide more detailed information with consequent requirements of greater computational time and effort. While these analytical and numerical computer models can provide accurate and adequate information to produce defensible remediation strategies, decisions based on inadequate modeling output or on over-analysis can have costly and risky consequences. In this chapter we consider both analytical and numerical modeling approaches to biodegradation and reactive transport. Both approaches are discussed and analyzed in terms of achieving bioremediation goals, recognizing that there is always a tradeoff between computational cost and the resolution of simulated systems.

  13. Numerical Non-Equilibrium and Smoothing of Solutions in The Difference Method for Plane 2-Dimensional Adhesive Joints / Nierównowaga Numeryczna i Wygładzanie Rozwiazań w Metodzie Różnicowej Dla Dwuwymiarowych Połączeń Klejowych

    Rapp Piotr

    2016-03-01

    Full Text Available The subject of the paper is related to problems with numerical errors in the finite difference method used to solve equations of the theory of elasticity describing 2- dimensional adhesive joints in the plane stress state. Adhesive joints are described in terms of displacements by four elliptic partial differential equations of the second order with static and kinematic boundary conditions. If adhesive joint is constrained as a statically determinate body and is loaded by a self-equilibrated loading, the finite difference solution is sensitive to kinematic boundary conditions. Displacements computed at the constraints are not exactly zero. Thus, the solution features a numerical error as if the adhesive joint was not in equilibrium. Herein this phenomenon is called numerical non-equilibrium. The disturbances in displacements and stress distributions can be decreased or eliminated by a correction of loading acting on the adhesive joint or by smoothing of solutions based on Dirichlet boundary value problem.

  14. Mathematical models and numerical simulation in electromagnetism

    Bermúdez, Alfredo; Salgado, Pilar

    2014-01-01

    The book represents a basic support for a master course in electromagnetism oriented to numerical simulation. The main goal of the book is that the reader knows the boundary-value problems of partial differential equations that should be solved in order to perform computer simulation of electromagnetic processes. Moreover it includes a part devoted to electric circuit theory  based on ordinary differential equations. The book is mainly oriented to electric engineering applications, going from the general to the specific, namely, from the full Maxwell’s equations to the particular cases of electrostatics, direct current, magnetostatics and eddy currents models. Apart from standard exercises related to analytical calculus, the book includes some others oriented to real-life applications solved with MaxFEM free simulation software.

  15. Numerical linearized MHD model of flapping oscillations

    Korovinskiy, D. B.; Ivanov, I. B.; Semenov, V. S.; Erkaev, N. V.; Kiehas, S. A.

    2016-06-01

    Kink-like magnetotail flapping oscillations in a Harris-like current sheet with earthward growing normal magnetic field component Bz are studied by means of time-dependent 2D linearized MHD numerical simulations. The dispersion relation and two-dimensional eigenfunctions are obtained. The results are compared with analytical estimates of the double-gradient model, which are found to be reliable for configurations with small Bz up to values ˜ 0.05 of the lobe magnetic field. Coupled with previous results, present simulations confirm that the earthward/tailward growth direction of the Bz component acts as a switch between stable/unstable regimes of the flapping mode, while the mode dispersion curve is the same in both cases. It is confirmed that flapping oscillations may be triggered by a simple Gaussian initial perturbation of the Vz velocity.

  16. Numerical modeling capabilities to predict repository performance

    1979-09-01

    This report presents a summary of current numerical modeling capabilities that are applicable to the design and performance evaluation of underground repositories for the storage of nuclear waste. The report includes codes that are available in-house, within Golder Associates and Lawrence Livermore Laboratories; as well as those that are generally available within the industry and universities. The first listing of programs are in-house codes in the subject areas of hydrology, solute transport, thermal and mechanical stress analysis, and structural geology. The second listing of programs are divided by subject into the following categories: site selection, structural geology, mine structural design, mine ventilation, hydrology, and mine design/construction/operation. These programs are not specifically designed for use in the design and evaluation of an underground repository for nuclear waste; but several or most of them may be so used.

  17. Numerical modeling capabilities to predict repository performance

    This report presents a summary of current numerical modeling capabilities that are applicable to the design and performance evaluation of underground repositories for the storage of nuclear waste. The report includes codes that are available in-house, within Golder Associates and Lawrence Livermore Laboratories; as well as those that are generally available within the industry and universities. The first listing of programs are in-house codes in the subject areas of hydrology, solute transport, thermal and mechanical stress analysis, and structural geology. The second listing of programs are divided by subject into the following categories: site selection, structural geology, mine structural design, mine ventilation, hydrology, and mine design/construction/operation. These programs are not specifically designed for use in the design and evaluation of an underground repository for nuclear waste; but several or most of them may be so used

  18. Numerical Modeling of Suspension HVOF Spray

    Jadidi, M.; Moghtadernejad, S.; Dolatabadi, A.

    2016-02-01

    A three-dimensional two-way coupled Eulerian-Lagrangian scheme is used to simulate suspension high-velocity oxy-fuel spraying process. The mass, momentum, energy, and species equations are solved together with the realizable k-ɛ turbulence model to simulate the gas phase. Suspension is assumed to be a mixture of solid particles [mullite powder (3Al2O3·2SiO2)], ethanol, and ethylene glycol. The process involves premixed combustion of oxygen-propylene, and non-premixed combustion of oxygen-ethanol and oxygen-ethylene glycol. One-step global reaction is used for each mentioned reaction together with eddy dissipation model to compute the reaction rate. To simulate the droplet breakup, Taylor Analogy Breakup model is applied. After the completion of droplet breakup, and solvent evaporation/combustion, the solid suspended particles are tracked through the domain to determine the characteristics of the coating particles. Numerical simulations are validated against the experimental results in the literature for the same operating conditions. Seven or possibly eight shock diamonds are captured outside the nozzle. In addition, a good agreement between the predicted particle temperature, velocity, and diameter, and the experiment is obtained. It is shown that as the standoff distance increases, the particle temperature and velocity reduce. Furthermore, a correlation is proposed to determine the spray cross-sectional diameter and estimate the particle trajectories as a function of standoff distance.

  19. Numerical Modelling of Flow and Settling in Secondary Settling Tanks

    Dahl, Claus Poulsen

    This thesis discusses the development of a numerical model for the simulation of secondary settling tanks. In the first part, the status on the development of numerical models for settling tanks and a discussion of the current design practice are presented. A study of the existing numerical models...

  20. Numerical modeling of subaqueous sand dune morphodynamics

    Doré, Arnaud; Bonneton, Philippe; Marieu, Vincent; Garlan, Thierry

    2016-03-01

    The morphodynamic evolution of subaqueous sand dunes is investigated, using a 2-D Reynolds-averaged Navier-Stokes numerical model. A laboratory experiment where dunes are generated under stationary unidirectional flow conditions is used as a reference case. The model reproduces the evolution of the erodible bed until a state of equilibrium is reached. In particular, the simulation exhibits the different stages of the bed evolution, e.g., the incipient ripple generation, the nonlinear bed form growing phase, and the dune field equilibrium phase. The results show good agreement in terms of dune geometrical dimensions and time to equilibrium. After the emergence of the first ripple field, the bed growth is driven by cascading merging sequences between bed forms of different heights. A sequence extracted from the simulation shows how the downstream bed form is first eroded before merging with the upstream bed form. Superimposed bed forms emerge on the dune stoss sides during the simulation. An analysis of the results shows that they emerge downstream of a slight deflection on the dune profile. The deflection arises due to a modification of the sediment flux gradient consecutive to a reduction in the turbulence relaxation length while the upstream bed form height decreases. As they migrate, superimposed bed forms grow on the dune stoss side and eventually provoke the degeneration of the dune crest. Cascading merging sequences and superimposed bed forms dynamics both influence the dune field evolution and size and therefore play a fundamental role in the dune field self-organization process.

  1. A numerical model of aerosol scavenging

    Using a three-dimensional numerical cloud/smoke-plume model, we have simulated the burning of a large, mid-latitude city following a nuclear exchange. The model includes 18 dynamic and microphysical equations that predict the fire-driven airflow, cloud processes, and smoke-cloud interactions. In the simulation, the intense heating from the burning city produces a firestorm with updraft velocities exceeding 60 m/s. Within 15 minutes of ignition, the smoke plume penetrates the tropopause. The updraft triggers a cumulonimbus cloud that produces significant quantities of ice, snow, and hail. These solid hydrometeors, as well as cloud droplets and rain, interact with the smoke particles from the fire. At the end of the one-hour simulation, over 20% of the smoke is in slowly falling snowflakes. If the snow reaches the ground before the flakes completely sublimate (or melt and then evaporate), then only approximately 50% of the smoke will survive the scavenging processes and remain in the atmosphere to affect the global climate

  2. Understanding Etna flank instability through numerical models

    Apuani, Tiziana; Corazzato, Claudia; Merri, Andrea; Tibaldi, Alessandro

    2013-02-01

    As many active volcanoes, Mount Etna shows clear evidence of flank instability, and different mechanisms were suggested to explain this flank dynamics, based on the recorded deformation pattern and character. Shallow and deep deformations, mainly associated with both eruptive and seismic events, are concentrated along recognised fracture and fault systems, mobilising the eastern and south-eastern flank of the volcano. Several interacting causes were postulated to control the phenomenon, including gravity force, magma ascent along the feeding system, and a very complex local and/or regional tectonic activity. Nevertheless, the complexity of such dynamics is still an open subject of research and being the volcano flanks heavily urbanised, the comprehension of the gravitative dynamics is a major issue for public safety and civil protection. The present research explores the effects of the main geological features (in particular the role of the subetnean clays, interposed between the Apennine-Maghrebian flysch and the volcanic products) and the role of weakness zones, identified by fracture and fault systems, on the slope instability process. The effects of magma intrusions are also investigated. The problem is addressed by integrating field data, laboratory tests and numerical modelling. A bi- and tri-dimensional stress-strain analysis was performed by a finite difference numerical code (FLAC and FLAC3D), mainly aimed at evaluating the relationship among geological features, volcano-tectonic structures and magmatic activity in controlling the deformation processes. The analyses are well supported by dedicated structural-mechanical field surveys, which allowed to estimate the rock mass strength and deformability parameters. To take into account the uncertainties which inevitably occur in a so complicated model, many efforts were done in performing a sensitivity analysis along a WNW-ESE section crossing the volcano summit and the Valle del Bove depression. This was

  3. Numerical Modelling of Sediment Transport in Combined Sewer Systems

    Schlütter, Flemming

    A conceptual sediment transport model has been developed. Through a case study a comparison with other numerical models is performed.......A conceptual sediment transport model has been developed. Through a case study a comparison with other numerical models is performed....

  4. Precise numerical modeling of next generation multimode fiber based links

    Maksymiuk, L.; Stepniak, G.

    2015-12-01

    In order to numerically model modern multimode fiber based links we are required to take into account modal and chromatic dispersion, profile dispersion and spectral dependent coupling. In this paper we propose a complete numerical model which not only is precise but also versatile. Additionally to the detailed mathematical description of the model we provide also a bunch of numerical calculations performed with the use of the model.

  5. Numerical modelling of nearshore wave transformation

    Chandramohan, P.; Nayak, B.U.; SanilKumar, V.

    A software has been developed for numerical refraction study based on finite amplitude wave theories. Wave attenuation due to shoaling, bottom friction, bottom percolation and viscous dissipation has also been incorporated. The software...

  6. Numerical models for high beta magnetohydrodynamic flow

    The fundamentals of numerical magnetohydrodynamics for highly conducting, high-beta plasmas are outlined. The discussions emphasize the physical properties of the flow, and how elementary concepts in numerical analysis can be applied to the construction of finite difference approximations that capture these features. The linear and nonlinear stability of explicit and implicit differencing in time is examined, the origin and effect of numerical diffusion in the calculation of convective transport is described, and a technique for maintaining solenoidality in the magnetic field is developed. Many of the points are illustrated by numerical examples. The techniques described are applicable to the time-dependent, high-beta flows normally encountered in magnetically confined plasmas, plasma switches, and space and astrophysical plasmas. 40 refs

  7. Numerical model of post-DNB transition boiling heat transfer

    In this paper a physical model for the transition boiling heat transfer is proposed. The corresponding mathematical descriptions are given in detail and the heat transfer characteristics of post-DNB transition boiling is analyzed. The numerical model of post-DNB transition boiling heat transfer is obtained as the empirical value of the coefficient is determined by the experimental data. The numerical model is compared with the experimental data of different parameters and other numerical models, and the statistical deviations are calculated. The calculating results of the numerical model in this paper show good agreement with the experimental data and the numerical model in this paper is with good applicability compared with other numerical models. (authors)

  8. Numerical Modeling of Earthquake Dynamic Rupture : Requirements for Realistic Modeling

    Fukuyama, Eiichi

    2003-01-01

    I propose a strategy to make a numerical computation applicable to the realistic modeling of an earthquake dynamic rupture process.To do this, it is important to introduce any observables into the simulation as initial and boundary conditions.As an initial condition, distribution of total stress before the dynamic rupture, and as boundary conditions, fault constitutive relation and geometry of the fault are necessary.The initial stress distribution would be obtained by both in-situ measuremen...

  9. Speckle Patterns and 2-Dimensional Brownian Motion

    We present the results of a Monte Carlo simulation of Brownian Motion on a 2-dimensional lattice with nearest-neighbor interactions described by a linear model. These nearest-neighbor interactions lead to a spatial variance structure on the lattice. The resulting Brownian pattern fluctuates in value from point to point in a manner characteristic of a stationary stochastic process. The value at a lattice point is interpreted as an intensity level. The difference in values in neighboring cells produces a fluctuating intensity pattern on the lattice. Changing the size of the mesh changes the relative size of the speckles. Increasing the mesh size tends to average out the intensity in the direction of the mean of the stationary process. (Author)

  10. Masonry constructions mechanical models and numerical applications

    Lucchesi, Massimiliano; Padovani, Cristina

    2008-01-01

    Numerical methods for the structural analysis of masonry constructions can be of great value in assessing the safety of artistically important masonry buildings and optimizing potential operations of maintenance and strengthening in terms of their cost-effectiveness, architectural impact and static effectiveness. This monograph firstly provides a detailed description of the constitutive equation of masonry-like materials, clearly setting out its most important features. It then goes on to provide a numerical procedure to solve the equilibrium problem of masonry solids. A large portion of the w

  11. Large scale experiments as a tool for numerical model development

    Kirkegaard, Jens; Hansen, Erik Asp; Fuchs, Jesper;

    2003-01-01

    improvement of the reliability of physical model results. This paper demonstrates by examples that numerical modelling benefits in various ways from experimental studies (in large and small laboratory facilities). The examples range from very general hydrodynamic descriptions of wave phenomena to specific......Experimental modelling is an important tool for study of hydrodynamic phenomena. The applicability of experiments can be expanded by the use of numerical models and experiments are important for documentation of the validity of numerical tools. In other cases numerical tools can be applied for...

  12. Conceptual and Numerical Models for UZ Flow and Transport

    The purpose of this Analysis/Model Report (AMR) is to document the conceptual and numerical models used for modeling of unsaturated zone (UZ) fluid (water and air) flow and solute transport processes. This is in accordance with ''AMR Development Plan for U0030 Conceptual and Numerical Models for Unsaturated Zone (UZ) Flow and Transport Processes, Rev 00''. The conceptual and numerical modeling approaches described in this AMR are used for models of UZ flow and transport in fractured, unsaturated rock under ambient and thermal conditions, which are documented in separate AMRs. This AMR supports the UZ Flow and Transport Process Model Report (PMR), the Near Field Environment PMR, and the following models: Calibrated Properties Model; UZ Flow Models and Submodels; Mountain-Scale Coupled Processes Model; Thermal-Hydrologic-Chemical (THC) Seepage Model; Drift Scale Test (DST) THC Model; Seepage Model for Performance Assessment (PA); and UZ Radionuclide Transport Models

  13. Numerical model of compressible gas flow in soil pollution control

    2002-01-01

    Based on the theory of fluid dynamics in porous media, a numerical model of gas flow in unsaturated zone is developed with the consideration of gas density change due to variation of air pressure. This model is characterized of its wider range of availability. The accuracy of this numerical model is analyzed through comparison with modeling results by previous model with presumption of little pressure variation and the validity of this numerical model is shown. Thus it provides basis for the designing and management of landfill gas control system or soil vapor ex.action system in soil pollution control.

  14. NUMERICAL MODELING OF FINE SEDIMENT PHYSICAL PROCESSES.

    Schoellhamer, David H.

    1985-01-01

    Fine sediment in channels, rivers, estuaries, and coastal waters undergo several physical processes including flocculation, floc disruption, deposition, bed consolidation, and resuspension. This paper presents a conceptual model and reviews mathematical models of these physical processes. Several general fine sediment models that simulate some of these processes are reviewed. These general models do not directly simulate flocculation and floc disruption, but the conceptual model and existing functions are shown to adequately model these two processes for one set of laboratory data.

  15. Numerical Modeling in Geodynamics: Success, Failure and Perspective

    Ismail-Zadeh, A.

    2005-12-01

    A real success in numerical modeling of dynamics of the Earth can be achieved only by multidisciplinary research teams of experts in geodynamics, applied and pure mathematics, and computer science. The success in numerical modeling is based on the following basic, but simple, rules. (i) People need simplicity most, but they understand intricacies best (B. Pasternak, writer). Start from a simple numerical model, which describes basic physical laws by a set of mathematical equations, and move then to a complex model. Never start from a complex model, because you cannot understand the contribution of each term of the equations to the modeled geophysical phenomenon. (ii) Study the numerical methods behind your computer code. Otherwise it becomes difficult to distinguish true and erroneous solutions to the geodynamic problem, especially when your problem is complex enough. (iii) Test your model versus analytical and asymptotic solutions, simple 2D and 3D model examples. Develop benchmark analysis of different numerical codes and compare numerical results with laboratory experiments. Remember that the numerical tool you employ is not perfect, and there are small bugs in every computer code. Therefore the testing is the most important part of your numerical modeling. (iv) Prove (if possible) or learn relevant statements concerning the existence, uniqueness and stability of the solution to the mathematical and discrete problems. Otherwise you can solve an improperly-posed problem, and the results of the modeling will be far from the true solution of your model problem. (v) Try to analyze numerical models of a geological phenomenon using as less as possible tuning model variables. Already two tuning variables give enough possibilities to constrain your model well enough with respect to observations. The data fitting sometimes is quite attractive and can take you far from a principal aim of your numerical modeling: to understand geophysical phenomena. (vi) If the number of

  16. Exercises in 80223 Numerical Modelling of Thermal Processing of Materials

    Frandsen, Jens Ole

    , guidelines are given on how to write the report which has to be handed in at the end of the course. The exercise book is a updated version of the exercise book from 1999. The exercise book is used in the course 42224 'Numerical Process Modelling' which earlier was called 80223 'Numerical Modelling of Thermal...

  17. Numerical Modelling of Electromagnetic Field in a Tornado

    Pavel Fiala

    2008-01-01

    Full Text Available This study deals with the numerical model of both the physical and the chemical processes in the tornado. Within the paper, a basic theoretical model and a numerical solution are presented. We prepared numerical models based on the combined finite element method (FEM and the finite volume method (FVM. The model joins the magnetic, electric and current fields, the flow field and a chemical nonlinear ion model. The results were obtained by means of the FEM/FVM as a main application in ANSYS software.

  18. Mathematical modelling and numerical simulation of casting processes

    Hattel, Jesper Henri

    1998-01-01

    The control volume method applied to numerical modelling of castning. Analytical solutions based on the error function.Riemann-temperature. Modelling of release of latent heat with the enthalpy method....

  19. Regression Models for 2-Dimensional Cartesian Coordinates Prediction: A Case study at University of Mines and Technology (UMaT, Tarkwa-Ghana.

    Yao Yevenyo Ziggah

    2013-01-01

    Full Text Available The aim of this research is to study and analyze statistical models applicable in bringing out arelationship between global coordinates and cartesian planimetric coordinates of some known controlstations in the University of Mines and Technology (UMaT campus. To achieve the aims of this research,the Global Position System (GPS latitudes and longitudes of selected control stations with knowncartesian planimetric coordinates were determined using the Handheld GPS receiver at different epoch(morning and evening. Linear Regression analysis was then conducted to establish the correlationbetween global and cartesian planimetric coordinates of the selected control stations and regressionmodels generated to show the results. The correlation coefficient r, a t-test for non -zero slope, t-test oncorrelation coefficient, graphical residual analysis, test of normality, comparing model predictions toobserved data, were used to evaluate and check the adequacy of the models. The obtained resultsindicated that the proposed linear regression models are suitable for predictions at 95% confidenceinterval and do not violate any of the statistical assumptions of a linear model. However, the proposedregression models for the evening observation gave better prediction accuracy than the morning. Acomputer programming algorithm and a designed interface was created for the proposed regressionmodels established using Microsoft C++ standard edition 6.0, thus making it easier in applying themodels in making cartesian planimetric coordinates prediction at different epoch at UMaT.

  20. Software Simplifies the Sharing of Numerical Models

    2014-01-01

    To ease the sharing of climate models with university students, Goddard Space Flight Center awarded SBIR funding to Reston, Virginia-based Parabon Computation Inc., a company that specializes in cloud computing. The firm developed a software program capable of running climate models over the Internet, and also created an online environment for people to collaborate on developing such models.

  1. Economic growth models: symbolic and numerical computations

    Vasconcelos, P. B.

    2013-01-01

    Growth economic models play a crucial role in understanding countries development, inter-country macroeconomic relationship and, ultimately, to anticipate the effects on endogenous variables due to political shocks on model parameters or exogenous variables. Yet, build a mathematical model can be a difficult and time consuming task. Symbolic computations can be of great help in the development process. Then, the ability to simulate, under initial assumptions, is, assuredly, a priceless tool f...

  2. A numerical model of the VKS experiment

    Gissinger, Christophe

    2009-01-01

    We present numerical simulations of the magnetic field generated by the flow of liquid sodium driven by two counter-rotating impellers (VKS experiment). Using a dynamo kinematic code in cylindrical geometry, it is shown that different magnetic modes can be generated depending on the flow configuration. While the time averaged axisymmetric mean flow generates an equatorial dipole, our simulations show that an axial field of either dipolar or quadrupolar symmetry can be generated by taking into account non-axisymmetric components of the flow. Moreover, we show that by breaking a symmetry of the flow, the magnetic field becomes oscillatory. This leads to reversals of the axial dipole polarity, involving a competition with the quadrupolar component.

  3. Numerical MHD Codes for Modeling Astrophysical Flows

    Koldoba, A V; Lii, P S; Comins, M L; Dyda, S; Romanova, M M; Lovelace, R V E

    2015-01-01

    We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.

  4. A numerical model for a room fire

    A simple room fire experiment has been used as a development case for numerical fire simulation based on computational fluid dynamics. Experimental values (Steckler et al. 1982) are used to verify simulation. A limited method to use a flee flow boundary condition has been described (flow is assumed to be oblique to the boundary). Description of buoyancy terms, atmospheric pressure and additional terms of turbulence due to buoyancy are also given. Fire source is described as a predefined heat source. Comparison to the experiments shows that flow velocity at the room door, which is the only opening to the room, can be predicted with a good and temperature with a moderate accuracy. (au) (3 refs., 13 figs., 2 tabs.)

  5. Numerical MHD codes for modeling astrophysical flows

    Koldoba, A. V.; Ustyugova, G. V.; Lii, P. S.; Comins, M. L.; Dyda, S.; Romanova, M. M.; Lovelace, R. V. E.

    2016-05-01

    We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.

  6. Topological 2-Dimensional Quantum Mechanics

    Dasnières de Veigy, A; Veigy, Alain Dasnieres de; Ouvry, Stephane

    1993-01-01

    We define a Chern- Simons Lagrangian for a system of planar particles topologically interacting at a distance. The anyon model appears as a particular case where all the particles are identical. We propose exact N-body eigenstates, set up a perturbative algorithm, discuss the case where some particles are fixed on a lattice, and also consider curved manifolds. PACS numbers: 05.30.-d, 11.10.-z

  7. Numerical model of Ca(OH)

    Koster, T.; Peelen, W.; Larbi, J.; Rooij, M. de; Polder, R.

    2010-01-01

    A mathematical model is being developed to describe a repair method in concrete, called cathodic protection (CP). The model is in principle also useful to describe electrodeposition in concrete, e.g. the process of re-precipitation of Ca(OH)2 invoked by an electrical current. In CP, the c

  8. A numerical model for multigroup radiation hydrodynamics

    We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M1) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well-established kinetic model through Marshak wave tests with frequency-dependent opacities. Very good agreement between the multigroup M1 and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.

  9. Traffic Flow Models and Their Numerical Solutions

    Jin, Wenlong

    2004-01-01

    In this thesis, Riemann problems and Godunov methods are developed for higher order traffic flow models. A rigorous analysis of the first order traffic flow model of inhomogeneous road is presented. A two-level simulation framework of network vehicular traffic is proposed as a Godunov-type finite difference system with the supply-demand method as an alternative of Riemann solver.

  10. Numerical Modelling of Wave Run-Up

    Ramirez, Jorge Robert Rodriguez; Frigaard, Peter; Andersen, Thomas Lykke;

    2011-01-01

    Wave loads are important in problems related to offshore structure, such as wave run-up, slamming. The computation of such wave problems are carried out by CFD models. This paper presents one model, NS3, which solve 3D Navier-Stokes equations and use Volume of Fluid (VOF) method to treat the free...

  11. Characterisation of the Bahía Blanca estuary by data analysis and numerical modelling

    Campuzano, Francisco Javier; Pierini, Jorge O.; Leitão, Paulo C.; Gómez, Eduardo A.; Neves, Ramiro J.

    2014-01-01

    The Bahía Blanca estuary is a complex system of channels and tidal flats where the most important deep water harbour system of Argentina is located. The main goal of the present work was to obtain a hydrodynamic conceptual model for the Bahía Blanca coastal area. For this reason, a combined analysis of observed data and numerical modelling has been performed for the whole area. The gained knowledge on the system hydrodynamics could aid in the decision support for navigation security, waste water discharges management, sediment dredging and rejection operations among other applications. Due to the Bahía Blanca coastal vast area, hydrodynamic observations are scarce and located near the populated areas. In order to describe the hydrodynamics of such a complex and large system, the analysed tidal and current data from different periods have been completed through numerical modelling. Data analysis served to determine the main processes governing the Bahía Blanca hydrodynamics, to characterise the area using general descriptors, to provide inputs for the numerical model and to aid in evaluating its performance. In addition, a 2-dimensional application was set up using the MOHID water modelling system for the Bahía Blanca estuary. This application aimed to gain a better understanding of the system dynamics, to explain and test the consistency of the observed data and to reproduce the processes taking place. Model results were in good agreement with the analysed data and served to confirm an inconsistency found on the sea level observations. The combination of both methodologies served to further describe the hydrodynamic processes governing this coastal area and also to obtain a conceptual model for the water and property circulation in the Bahía Blanca estuary.

  12. Deterministic combination of numerical and physical coastal wave models

    Zhang, H.W.; Schäffer, Hemming Andreas; Jakobsen, K.P.

    2007-01-01

    the interface between the numerical and physical models. The link between numerical and physical models is given by an ad hoc unified wave generation theory which is devised in the study. This wave generation theory accounts for linear dispersion and shallow water non-linearity. Local wave phenomena......A deterministic combination of numerical and physical models for coastal waves is developed. In the combined model, a Boussinesq model MIKE 21 BW is applied for the numerical wave computations. A piston-type 2D or 3D wavemaker and the associated control system with active wave absorption provides...... (evanescent modes) near the wavemaker are taken into account. With this approach, the data transfer between the two models is thus on a deterministic level with detailed wave information transmitted along the wavemaker....

  13. Numerical Modeling and Experimental Testing of a Wave Energy Converter

    Zurkinden, Andrew Stephen; Kramer, Morten; Ferri, Francesco;

    numerical values for comparison with the experimental test results which were carried out in the same time. It is for this reason why Chapter 4 does consist exclusively of numerical values. Experimental values and measured time series of wave elevations have been used throughout the report in order to a......) validate the numerical model and b) preform stochastic analysis. The latter technique is introduced in order to optimize the control parameters of the power take off system....

  14. The influence of numerical models on determining the drag coefficient

    Dobeš Josef

    2014-03-01

    Full Text Available The paper deals with numerical modelling of body aerodynamic drag coefficient in the transition from laminar to turbulent flow regimes, where the selection of a suitable numerical model is problematic. On the basic problem of flow around a simple body – sphere selected computational models are tested. The values obtained by numerical simulations of drag coefficients of each model are compared with the graph of dependency of the drag coefficient vs. Reynolds number for a sphere. Next the dependency of Strouhal number vs. Reynolds number is evaluated, where the vortex shedding frequency values for given speed are obtained numerically and experimentally and then the values are compared for each numerical model and experiment. The aim is to specify trends for the selection of appropriate numerical model for flow around bodies problem in which the precise description of the flow field around the obstacle is used to define the acoustic noise source. Numerical modelling is performed by finite volume method using CFD code.

  15. Dependent Risk Modelling and Ruin Probability: Numerical Computation and Applications

    Zhao, Shouqi

    2014-01-01

    In this thesis, we are concerned with the finite-time ruin probabilities in two alternative dependent risk models, the insurance risk model and the dual risk model, including the numerical evaluation of the explicit expressions for these quantities and the application of the probabilistic results obtained. We first investigate the numerical properties of the formulas for the finite-time ruin probability derived by Ignatov and Kaishev (2000, 2004) and Ignatov et al. (2001) for a generalized in...

  16. Numerical Calculation of Model Rocket Trajectories.

    Keeports, David

    1990-01-01

    Discussed is the use of model rocketry to teach the principles of Newtonian mechanics. Included are forces involved; calculations for vertical launches; two-dimensional trajectories; and variations in mass, drag, and launch angle. (CW)

  17. Computational numerical modelling of plasma focus

    Several models for calculation of the dynamics of Plasma Focus have been developed. All of them begin from the same physic principle: the current sheet run down the anode length, ionizing and collecting the gas that finds in its way.This is known as snow-plow model.Concerning pinch's compression, a MHD model is proposed.The plasma is treated as a fluid , particularly as a high ionized gas.However, there are not many models that, taking into account thermal equilibrium inside the plasma, make approximated calculations of the maximum temperatures reached in the pinch.Besides, there are no models which use those temperatures to estimate the termofusion neutron yield for the Deuterium or Deuterium-Tritium gas filled cases.In the PLADEMA network (Dense Magnetized Plasmas) a code was developed with the objective of describe the plasma focus dynamics, in a conceptual engineering stage.The codes calculates the principal variables (currents, time to focus, etc) and estimates the neutron yield in Deuterium-filled plasma focus devices.It can be affirmed that the code's experimental validation, in its axial and radial stages, was very successfully. However, it was accepted that the compression stage should be formulated again, to find a solution for a large variation of a parameter related with velocity profiles for the particles trapped inside the pinch.The objectives of this work can be stated in the next way : - Check the compression's model hypothesis. Develop a new model .- Implement the new model in the code. Compare results against experimental data of Plasma Focus devices from all around the world

  18. Modeling, Analysis, and Numerics in Electrohydrodynamics

    Schmuck, Markus

    2008-01-01

    The main subject of this thesis is to analyze the incompressible Navier-Stokes-Nernst-Planck-Poisson system for bounded domains. Such a system is used as a model in electrohydrodynamics or physicochemical models. First, we verify existence of weak and strong solutions. Moreover, we are able to characterize the weak solutions by an energy and an entropy law. The concentrations in the Nernst-Planck equations additionally are non-negative and bounded. These results motivate to construct conv...

  19. Evaluation of wave runup predictions from numerical and parametric models

    Stockdon, Hilary F.; Thompson, David M.; Plant, Nathaniel G.; Long, Joseph W.

    2014-01-01

    Wave runup during storms is a primary driver of coastal evolution, including shoreline and dune erosion and barrier island overwash. Runup and its components, setup and swash, can be predicted from a parameterized model that was developed by comparing runup observations to offshore wave height, wave period, and local beach slope. Because observations during extreme storms are often unavailable, a numerical model is used to simulate the storm-driven runup to compare to the parameterized model and then develop an approach to improve the accuracy of the parameterization. Numerically simulated and parameterized runup were compared to observations to evaluate model accuracies. The analysis demonstrated that setup was accurately predicted by both the parameterized model and numerical simulations. Infragravity swash heights were most accurately predicted by the parameterized model. The numerical model suffered from bias and gain errors that depended on whether a one-dimensional or two-dimensional spatial domain was used. Nonetheless, all of the predictions were significantly correlated to the observations, implying that the systematic errors can be corrected. The numerical simulations did not resolve the incident-band swash motions, as expected, and the parameterized model performed best at predicting incident-band swash heights. An assimilated prediction using a weighted average of the parameterized model and the numerical simulations resulted in a reduction in prediction error variance. Finally, the numerical simulations were extended to include storm conditions that have not been previously observed. These results indicated that the parameterized predictions of setup may need modification for extreme conditions; numerical simulations can be used to extend the validity of the parameterized predictions of infragravity swash; and numerical simulations systematically underpredict incident swash, which is relatively unimportant under extreme conditions.

  20. Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling

    Gunzburger, Max [Florida State Univ., Tallahassee, FL (United States)

    2015-02-17

    We have treated the modeling, analysis, numerical analysis, and algorithmic development for nonlocal models of diffusion and mechanics. Variational formulations were developed and finite element methods were developed based on those formulations for both steady state and time dependent problems. Obstacle problems and optimization problems for the nonlocal models were also treated and connections made with fractional derivative models.

  1. Elements of Constitutive Modelling and Numerical Analysis of Frictional Soils

    Jakobsen, Kim Parsberg

    This thesis deals with elements of elasto-plastic constitutive modelling and numerical analysis of frictional soils. The thesis is based on a number of scientific papers and reports in which central characteristics of soil behaviour and applied numerical techniques are considered. The development...

  2. Numerical Considerations for Lagrangian Stochastic Dispersion Models: Eliminating Rogue Trajectories, and the Importance of Numerical Accuracy

    Bailey, Brian N.

    2016-07-01

    When Lagrangian stochastic models for turbulent dispersion are applied to complex atmospheric flows, some type of ad hoc intervention is almost always necessary to eliminate unphysical behaviour in the numerical solution. Here we discuss numerical strategies for solving the non-linear Langevin-based particle velocity evolution equation that eliminate such unphysical behaviour in both Reynolds-averaged and large-eddy simulation applications. Extremely large or `rogue' particle velocities are caused when the numerical integration scheme becomes unstable. Such instabilities can be eliminated by using a sufficiently small integration timestep, or in cases where the required timestep is unrealistically small, an unconditionally stable implicit integration scheme can be used. When the generalized anisotropic turbulence model is used, it is critical that the input velocity covariance tensor be realizable, otherwise unphysical behaviour can become problematic regardless of the integration scheme or size of the timestep. A method is presented to ensure realizability, and thus eliminate such behaviour. It was also found that the numerical accuracy of the integration scheme determined the degree to which the second law of thermodynamics or `well-mixed condition' was satisfied. Perhaps more importantly, it also determined the degree to which modelled Eulerian particle velocity statistics matched the specified Eulerian distributions (which is the ultimate goal of the numerical solution). It is recommended that future models be verified by not only checking the well-mixed condition, but perhaps more importantly by checking that computed Eulerian statistics match the Eulerian statistics specified as inputs.

  3. Numerical modeling of a cutting torch

    A two-dimensional turbulent model of a low-current intensity (30 A) cutting plasma torch was developed using the local thermodynamic equilibrium approximation. A good agreement was found between measured and modelled results of plasma temperature and velocity, the latter of which has not been previously reported in the literature for high energy density torches. The cutting performance was also studied in terms of the heat flux to the work-piece and the value of the force exerted by the impinging plasma jet.

  4. Amorphous track models: a numerical comparison study

    Greilich, Steffen; Grzanka, Leszek; Hahn, Ute; Kiderlen, Markus; Bassler, Niels; Andersen, Claus E.; Jäkel, Oliver

    carbon ion treatment at the particle facility HIT in Heidelberg. Apparent differences between the LEM and the Katz model are the way how interactions of individual particle tracks and how extended targets are handled. Complex scenarios, however, can mask the actual effect of these differences. Here, we...

  5. A numerical reference model for themomechanical subduction

    Quinquis, Matthieu; Chemia, Zurab; Tosi, Nicola; Buiter, Susanne; Dolejš, David

    2010-01-01

    devolatilization reactions. Our reference model represents ocean-ocean convergence and describes initial geometries and lithological stratification for a three-layered subducting slab and overriding plate along with their respective flow laws and chemical composition. It also includes kinematic and thermal...

  6. Numerical modelling of instantaneous plate tectonics

    Minster, J. B.; Haines, E.; Jordan, T. H.; Molnar, P.

    1974-01-01

    Assuming lithospheric plates to be rigid, 68 spreading rates, 62 fracture zones trends, and 106 earthquake slip vectors are systematically inverted to obtain a self-consistent model of instantaneous relative motions for eleven major plates. The inverse problem is linearized and solved iteratively by a maximum-likelihood procedure. Because the uncertainties in the data are small, Gaussian statistics are shown to be adequate. The use of a linear theory permits (1) the calculation of the uncertainties in the various angular velocity vectors caused by uncertainties in the data, and (2) quantitative examination of the distribution of information within the data set. The existence of a self-consistent model satisfying all the data is strong justification of the rigid plate assumption. Slow movement between North and South America is shown to be resolvable.

  7. Analytical and Numerical Modeling for Flexible Pipes

    WANG Wei; CHEN Geng

    2011-01-01

    The unbonded flexible pipe of eight layers,in which all the layers except the carcass layer are assumed to have isotropic properties,has been analyzed.Specifically,the carcass layer shows the orthotropic characteristics.The effective elastic moduli of the carcass layer have been developed in terms of the influence of deformation to stiffness.With consideration of the effective elastic moduli,the structure can be properly analyzed.Also the relative movements of tendons and relative displacements of wires in helical armour layer have been investigated.A three-dimensional nonlinear finite element model has been presented to predict the response of flexible pipes under axial force and torque.Further,the friction and contact of interlayer have been considered.Comparison between the finite element model and experimental results obtained in literature has been given and discussed,which might provide practical and technical support for the application of unbonded flexible pipes.

  8. Multidimensional numerical modeling of heat exchanges

    A comprehensive, multidimensional, thermal-hydraulic model is developed for the analysis of shell-and-tube heat exchanges for liquid metal services. For the shellside fluid, the conservation equations of mass, momentum and energy for continuum fluids are modified using the concept of porosity, surface premeability and distributed resistance to account for the blockage effects due to the presence of heat transfer tubes, flow baffles/shrouds, the support plates, etc. On the tubeside, heat transfer tubes are connected in parallel between the inlet and outlet plenums, and tubeside flow distribution is calculated based on the plenum-to-plenum pressure difference being equal for all tubes. It is assumed that the fluid remains single-phased on the shell side and may undergo phase-change in the tube side, thereby simulating the conditions of Liquid Metal Fast Breeder Reaction (LMFBR) intermediate heat exchanges (IHX) and steam generators (SG). The analytical model predictions are compared with three sets of test data (one for IHX and two for SG) and favorable results are obtained, thus providing a limited validation of the model

  9. Numerical Modelling of the Mining Induced Horizontal Displacement

    Tajduś, Krzysztof

    2015-12-01

    The paper presents results of numerical calculations and modeling of mining-induced surface deformation based on Finite Element Method (FEM). Applying the numerical method discussed to calculations allows us to assume a larger number of factors, such as rock mass structure, fracture network, rock properties, etc., which essentially affect the results obtained. On the basis of an elastic transversely isotropic model, an analysis of horizontal displacement distribution and surface subsidence was carried out for two sample regions of mines. The results of numerical calculations were later compared with the measured values. Such an analysis proved that the applied numerical model properly described distribution and values of subsidence and slope of subsidence trough, though there were serious differences in the values of calculated horizontal displacement, especially in areas of far influence range. In order to improve the matching, the influence of boundary conditions of the model on the value of calculated horizontal displacement was analyzed. The results are presented in graphs.

  10. Summary of Numerical Modeling for Underground Nuclear Test Monitoring Symposium

    This document contains the Proceedings of the Numerical Modeling for Underground Nuclear Test Monitoring Symposium held in Durango, Colorado on March 23-25, 1993. The symposium was sponsored by the Office of Arms Control and Nonproliferation of the United States Department of Energy and hosted by the Source Region Program of Los Alamos National Laboratory. The purpose of the meeting was to discuss state-of-the-art advances in numerical simulations of nuclear explosion phenomenology for the purpose of test ban monitoring. Another goal of the symposium was to promote discussion between seismologists and explosion source-code calculators. Presentation topics include the following: numerical model fits to data, measurement and characterization of material response models, applications of modeling to monitoring problems, explosion source phenomenology, numerical simulations and seismic sources

  11. Mathematical and numerical foundations of turbulence models and applications

    Chacón Rebollo, Tomás

    2014-01-01

    With applications to climate, technology, and industry, the modeling and numerical simulation of turbulent flows are rich with history and modern relevance. The complexity of the problems that arise in the study of turbulence requires tools from various scientific disciplines, including mathematics, physics, engineering, and computer science. Authored by two experts in the area with a long history of collaboration, this monograph provides a current, detailed look at several turbulence models from both the theoretical and numerical perspectives. The k-epsilon, large-eddy simulation, and other models are rigorously derived and their performance is analyzed using benchmark simulations for real-world turbulent flows. Mathematical and Numerical Foundations of Turbulence Models and Applications is an ideal reference for students in applied mathematics and engineering, as well as researchers in mathematical and numerical fluid dynamics. It is also a valuable resource for advanced graduate students in fluid dynamics,...

  12. Terrane accretion: Insights from numerical modelling

    Vogt, Katharina; Gerya, Taras

    2016-04-01

    The oceanic crust is not homogenous, but contains significantly thicker crust than norm, i.e. extinct arcs, spreading ridges, detached continental fragments, volcanic piles or oceanic swells. These (crustal) fragments may collide with continental crust and form accretionary complexes, contributing to its growth. We analyse this process using a thermo-mechanical computer model (i2vis) of an ocean-continent subduction zone. In this model the oceanic plate can bend spontaneously under the control of visco-plastic rheologies. It moreover incorporates effects such as mineralogical phase changes, fluid release and consumption, partial melting and melt extraction. Based on our 2-D experiments we suggest that the lithospheric buoyancy of the downgoing slab and the rheological strength of crustal material may result in a variety of accretionary processes. In addition to terrane subduction, we are able to identify three distinct modes of terrane accretion: frontal accretion, basal accretion and underplating plateaus. We show that crustal fragments may dock onto continental crust and cease subduction, be scrapped off the downgoing plate, or subduct to greater depth prior to slab break off and subsequent exhumation. Direct consequences of these processes include slab break off, subduction zone transference, structural reworking, formation of high-pressure terranes, partial melting and crustal growth.

  13. Comparisons of numerical modelling of the Selective Laser Melting

    Van Belle, Laurent; Vansteenkiste, Guillaume; Boyer, Jean-Claude

    2012-01-01

    International audience Selective laser melting (SLM) first developed for rapid prototyping (RP) is now used for rapid manufacturing of parts with inner complex shapes that cannot be made by more conventional routes. For example, production of injection moulds with cooling channels is of special interest. In this paper, a numerical model of SLM process was investigated to simulate the genesis of residual stresses. The proposed numerical modelling is based upon a double meshing with a multi-...

  14. Complexities in coastal sediment transport studies by numerical modeling

    Dandayudapani, I.; M. Murali

    2013-01-01

    Marine environmental studies related to erosion, accretion, pollution transport, dredge disposal, location of seawater intake, effluent disposal, etc., involve sediment transport studies. Numerical models use set of well linked mathematical equations arrived based on scientific principles as all natural phenomena are governed by certain rules which can be explained by scientific principles. Efficiency of numerical modeling greatly depends on quality of input parameters. When input parameters ...

  15. Numerical modeling in photonic crystals integrated technology: the COPERNICUS Project

    Malaguti, Stefania; Armaroli, Andrea; Bellanca, Gaetano; Trillo, Stefano; Kaunga-Nyirenda, Simeon; Lim, Jun; Larkins, Eric; Kristensen, Philip Trøst; Yvind, Kresten; Mørk, Jesper; Dumeige, Yannick; Gay, Mathilde; Colman, Pierre; Combrie, Sylvain; De Rossi, Alfredo

    Photonic crystals will play a fundamental role in the future of optical communications. The relevance of the numerical modeling for the success of this technology is assessed by using some examples concerning the experience of the COPERNICUS Project.......Photonic crystals will play a fundamental role in the future of optical communications. The relevance of the numerical modeling for the success of this technology is assessed by using some examples concerning the experience of the COPERNICUS Project....

  16. Numerical modelling of concentrated leak erosion during Hole Erosion Tests

    Mercier, F.; S. Bonelli; Golay, F.; Anselmet, F; Philippe, P.; Borghi, R.

    2015-01-01

    This study focuses on the numerical modelling of concentrated leak erosion of a cohesive soil by a turbulent flow in axisymmetrical geometry, with application to the Hole Erosion Test (HET). The numerical model is based on adaptive remeshing of the water/soil interface to ensure accurate description of the mechanical phenomena occurring near the soil/water interface. The erosion law governing the interface motion is based on two erosion parameters: the critical shear stress and the erosion co...

  17. Numerical simulation of a contractivity based multiscale cancer invasion model

    Kolbe, Niklas; Lukacova-Medvidova, Maria; Sfakianakis, Nikolaos; Wiebe, Bettina

    2016-01-01

    We present a problem-suited numerical method for a particularly challenging cancer invasion model. This model is a multiscale haptotaxis advection-reaction-diffusion system that describes the macroscopic dynamics of two types of cancer cells coupled with microscopic dynamics of the cells adhesion on the extracellular matrix. The difficulties to overcome arises from the non-constant advection and diffusion coefficients, a time delay term, as well as stiff reaction terms. Our numerical method i...

  18. Numerical modelling of composite laminates with through-thickness-reinforcements

    Grassi, Marcello

    2004-01-01

    The main objective of the present research study was to develop numerical models to investigate the mechanical properties and effectiveness of z-fibre reinforced laminates. A survey of relevant literature on through-thickness reinforcements (TTR) was undertaken and z-fibre pinning was chosen as the main topic of study. The development of numerical tools was mainly based on the finite element (FE) method and was carried out at different model scale levels. At a micro-mechanical level of ana...

  19. Numerical solution of dynamic equilibrium models under Poisson uncertainty

    Posch, Olaf; Trimborn, Timo

    2013-01-01

    We propose a simple and powerful numerical algorithm to compute the transition process in continuous-time dynamic equilibrium models with rare events. In this paper we transform the dynamic system of stochastic differential equations into a system of functional differential equations of the retar...... solution to Lucas' endogenous growth model under Poisson uncertainty are used to compute the exact numerical error. We show how (potential) catastrophic events such as rare natural disasters substantially affect the economic decisions of households....

  20. The Numerical Modeling of Transient Regimes of Diesel Generator Sets

    Cristian Roman

    2010-07-01

    Full Text Available This paper deals with the numerical modeling of a diesel generator set used as amain energy source in isolated areas and as a back-up energy source in the case ofrenewable energy systems. The numerical models are developed using a Matlab/Simulinksoftware package and they prove to be a powerful tool for the computer aided design ofcomplex hybrid power systems. Several operation regimes of the equipment are studied.The numerical study is completed with experimental measurements on a Kipor type dieselelectricgenerator set.

  1. Multidimensional numerical modeling of heat exchangers

    Sha, W. T.; Yang, C. I.; Kao, T. T.; Cho, S. M.

    A comprehensive, multidimensional, thermal-hydraulic model is developed for the analysis of shell-and-tube heat exchangers for liquid-metal services. For the shellside fluid, the conservation equations of mass, momentum, and energy for continuum fluids are modified using the concept of porosity, surface permeability and distributed resistance to account for the blockage effects due to the presence of heat-transfer tubes, flow baffles/shrouds, the support plates, etc. On the tubeside, the heat-transfer tubes are connected in parallel between the inlet and outlet plenums, and tubeside flow distribution is calculated based on the plenum-to-plenum pressure difference being equal for all tubes. It is assumed that the fluid remains single-phase on the shell side and may undergo phase-change on the tube side, thereby simulating the conditions of Liquid Metal Fast Breeder Reactor (LMFBR) intermediate heat exchangers (IHX) and steam generators (SG).

  2. Multidimensional numerical modeling of heat exchangers

    A comprehensive, multidimensional, thermal-hydraulic model is developed for the analysis of shell-and-tube heat exchangers for liquid-metal services. For the shellside fluid, the conservation equations of mass, momentum, and energy for continuum fluids are modified using the concept of porosity, surface permeability and distributed resistance to account for the blockage effects due to the presence of heat-transfer tubes, flow baffles/shrouds, the support plates, etc. On the tubeside, the heat-transfer tubes are connected in parallel between the inlet and outlet plenums, and tubeside flow distribution is calculated based on the plenum-to-plenum pressure difference being equal for all tubes. It is assumed that the fluid remains single-phase on the shell side and may undergo phase-change on the tube side, thereby simulating the conditions of Liquid Metal Fast Breeder Reactor (LMFBR) intermediate heat exchangers (IHX) and steam generators (SG)

  3. A modular approach to numerical human body modeling

    Forbes, P.A.; Griotto, G.; Rooij, L. van

    2007-01-01

    The choice of a human body model for a simulated automotive impact scenario must take into account both accurate model response and computational efficiency as key factors. This study presents a "modular numerical human body modeling" approach which allows the creation of a customized human body mod

  4. Numerical Modelling of Induction Heating Process for Testing Plant

    HRENIUC Ruslan Ovidiu; TIURBE Cristian

    2012-01-01

    This paper presents the numerical modelling of electromagnetic and thermal fields in the induction heating of a 2 kW testing plant, made by the author. The purpose of this modelling is to analyze thedevelopment of heating process, to notice any changes required to increase its efficiency. Modelling is performed by means of FLUX software package.

  5. Numerical model for learning concepts of streamflow simulation

    DeLong, L.L.

    1993-01-01

    Numerical models are useful for demonstrating principles of open-channel flow. Such models can allow experimentation with cause-and-effect relations, testing concepts of physics and numerical techniques. Four PT is a numerical model written primarily as a teaching supplement for a course in one-dimensional stream-flow modeling. Four PT options particularly useful in training include selection of governing equations, boundary-value perturbation, and user-programmable constraint equations. The model can simulate non-trivial concepts such as flow in complex interconnected channel networks, meandering channels with variable effective flow lengths, hydraulic structures defined by unique three-parameter relations, and density-driven flow.The model is coded in FORTRAN 77, and data encapsulation is used extensively to simplify maintenance and modification and to enhance the use of Four PT modules by other programs and programmers.

  6. Stepped spillway optimization through numerical and physical modeling

    Hamed Sarkardeh, Morteza Marosi, Raza Roshan

    2015-01-01

    Full Text Available The spillway is among the most important structures of a dam. It is importance for the spillway to be designed properly and passes flood flow safely with more energy dissipation. The zone which ogee spillway crest and stepped chute profile are joined with each other is important in design view. In the present study, a physical model as well as a numerical model was employed on a case study of stepped spillway to modify the transitional zone and improve flow pattern over the spillway. Many alternatives were examined and optimized. Finally, the performance of the selected alternative was checked for different flow conditions, air entrainment and energy dissipation. To simulate the turbulence phenomenon, RNG model and for free surface VOF model was selected in the numerical model. Results of the numerical and physical models were compared and good agreement concluded in flow conditions and energy dissipation.

  7. Numerical modeling of shoreline undulations part 1: Constant wave climate

    Kærgaard, Kasper Hauberg; Fredsøe, Jørgen

    2013-01-01

    This paper presents a numerical study of the non-linear development of alongshore undulations up to fully developed quasi-steady equilibrium. A numerical model which describes the longshore sediment transport along arbitrarily shaped shorelines is applied, based on a spectral wave model, a depth...... integrated flow model, a wave-phase resolving sediment transport description and a one-line shoreline model.First the length of the shoreline undulations is determined in the linear regime using a stability analysis. Next the further evolution from the linear to the fully non-linear regime is described. In...

  8. Numerical simulations of a reduced model for blood coagulation

    Pavlova, Jevgenija; Fasano, Antonio; Sequeira, Adélia

    2016-04-01

    In this work, the three-dimensional numerical resolution of a complex mathematical model for the blood coagulation process is presented. The model was illustrated in Fasano et al. (Clin Hemorheol Microcirc 51:1-14, 2012), Pavlova et al. (Theor Biol 380:367-379, 2015). It incorporates the action of the biochemical and cellular components of blood as well as the effects of the flow. The model is characterized by a reduction in the biochemical network and considers the impact of the blood slip at the vessel wall. Numerical results showing the capacity of the model to predict different perturbations in the hemostatic system are discussed.

  9. Quantitative comparisons of numerical models of brittle wedge dynamics

    Buiter, Susanne

    2010-05-01

    Numerical and laboratory models are often used to investigate the evolution of deformation processes at various scales in crust and lithosphere. In both approaches, the freedom in choice of simulation method, materials and their properties, and deformation laws could affect model outcomes. To assess the role of modelling method and to quantify the variability among models, we have performed a comparison of laboratory and numerical experiments. Here, we present results of 11 numerical codes, which use finite element, finite difference and distinct element techniques. We present three experiments that describe shortening of a sand-like, brittle wedge. The material properties of the numerical ‘sand', the model set-up and the boundary conditions are strictly prescribed and follow the analogue setup as closely as possible. Our first experiment translates a non-accreting wedge with a stable surface slope of 20 degrees. In agreement with critical wedge theory, all models maintain the same surface slope and do not deform. This experiment serves as a reference that allows for testing against analytical solutions for taper angle, root-mean-square velocity and gravitational rate of work. The next two experiments investigate an unstable wedge in a sandbox-like setup, which deforms by inward translation of a mobile wall. The models accommodate shortening by formation of forward and backward shear zones. We compare surface slope, rate of dissipation of energy, root-mean-square velocity, and the location, dip angle and spacing of shear zones. We show that we successfully simulate sandbox-style brittle behaviour using different numerical modelling techniques and that we obtain the same styles of deformation behaviour in numerical and laboratory experiments at similar levels of variability. The GeoMod2008 Numerical Team: Markus Albertz, Michelle Cooke, Tony Crook, David Egholm, Susan Ellis, Taras Gerya, Luke Hodkinson, Boris Kaus, Walter Landry, Bertrand Maillot, Yury Mishin

  10. Numerical Modeling and Mechanical Analysis of Flexible Risers

    J. Y. Li

    2015-01-01

    Full Text Available ABAQUS is used to create a detailed finite element model for a 10-layer unbonded flexible riser to simulate the riser’s mechanical behavior under three load conditions: tension force and internal and external pressure. It presents a technique to create detailed finite element model and to analyze flexible risers. In FEM model, all layers are modeled separately with contact interfaces; interaction between steel trips in certain layers has been considered as well. FEM model considering contact interaction, geometric nonlinearity, and friction has been employed to accurately simulate the structural behavior of riser. The model includes the main features of the riser geometry with very little simplifying assumptions. The model was solved using a fully explicit time-integration scheme implemented in a parallel environment on an eight-processor cluster and 24 G memory computer. There is a very good agreement obtained from numerical and analytical comparisons, which validates the use of numerical model here. The results from the numerical simulation show that the numerical model takes into account various details of the riser. It has been shown that the detailed finite element model can be used to predict riser’s mechanics behavior under various load cases and bound conditions.

  11. A Semi-implicit Numerical Scheme for a Two-dimensional, Three-field Thermo-Hydraulic Modeling

    Hwang, Moonkyu; Jeong, Jaejoon

    2007-07-15

    The behavior of two-phase flow is modeled, depending on the purpose, by either homogeneous model, drift flux model, or separated flow model, Among these model, in the separated flow model, the behavior of each flow phase is modeled by its own governing equation, together with the interphase models which describe the thermal and mechanical interactions between the phases involved. In this study, a semi-implicit numerical scheme for two-dimensional, transient, two-fluid, three-field is derived. The work is an extension to the previous study for the staggered, semi-implicit numerical scheme in one-dimensional geometry (KAERI/TR-3239/2006). The two-dimensional extension is performed by specifying a relevant governing equation set and applying the related finite differencing method. The procedure for employing the semi-implicit scheme is also described in detail. Verifications are performed for a 2-dimensional vertical plate for a single-phase and two-phase flows. The calculations verify the mass and energy conservations. The symmetric flow behavior, for the verification problem, also confirms the momentum conservation of the numerical scheme.

  12. 2D-Modelling of pellet injection in the poloidal plane: results of numerical tests

    A time-dependent two-dimensional resistive MHD code is being developed for computing the expansion of pellet-produced clouds in the poloidal plane. The various components of the code complex are being tested by means of simplified model calculations. In the present paper, results pertaining to the expansion and drift of a high density plasmoid in a magnetically confined homogeneous background plasma are reported. The confining (toroidal) field may be uniform or has a prescribed gradient. In this test phase, the 2-dimensional code is ran in a one-dimensional mode: it is assumed that all changes are restricted to the 'x' direction, which represents in our case the radial direction in the poloidal plane. The full set of time-dependent resistive MHD equations consisting of the conservation equations for mass, momentum, and energy, and supplemented by Maxwell's equations, a number of rate equations (ionization rate, etc.), and equations describing diffusive transport processes (internal energy, magnetic field) is solved by applying a second order Godunov numerical scheme and the GMRES method. The numerical scheme is based on a Riemann solver with Roe's approximation. In this analysis, three cases are considered in detail: a) the magnetic field strength is zero; b) an initially homogeneous magnetic field of B = 2 Tesla is applied over the whole domain; c) a spatially varying magnetic field with a gradient of the order of 1 tesla/m is applied over the width of the plasmoid

  13. Numerical modelling of river morphodynamics: Latest developments and remaining challenges

    Siviglia, Annunziato; Crosato, Alessandra

    2016-07-01

    Numerical morphodynamic models provide scientific frameworks for advancing our understanding of river systems. The research on involved topics is an important and socially relevant undertaking regarding our environment. Nowadays numerical models are used for different purposes, from answering questions about basic morphodynamic research to managing complex river engineering problems. Due to increasing computer power and the development of advanced numerical techniques, morphodynamic models are now more and more used to predict the bed patterns evolution to a broad spectrum of spatial and temporal scales. The development and the success of application of such models are based upon a wide range of disciplines from applied mathematics for the numerical solution of the equations to geomorphology for the physical interpretation of the results. In this light we organized this special issue (SI) soliciting multidisciplinary contributions which encompass any aspect needed for the development and applications of such models. Most of the papers in the SI stem from contributions to session HS9.5/GM7.11 on numerical modelling and experiments in river morphodynamics at the European Geosciences Union (EGU) General Assembly held in Vienna, April 27th to May 2nd 2014.

  14. Numerical simulations of bubbly flows using an averaged equations' model

    The paper discusses a numerical method for solving a two phase flow model based on the interpenetrating continua hypothesis. The model incorporates terms to account for the effects of virtual mass force, different pressures for the two phases and the viscous dissipation. Our numerical scheme extends the incremental projection scheme for the incompressible Navier-Stokes equation toward the multiphase flows. An optimal stability is obtained by slightly modifying the Galerkin formulation. The stabilized Galerkin technique we used is based on a two-level hierarchical decomposition of the approximation space. Numerical simulations of the three-dimensional bubbly flows in a periodic domain are presented. These simulations are compared with experiments. The stability of this flow with respect to 3D perturbations is studied numerically and a discussion of the results is presented. (author)

  15. Numerical modelling of low temperature radio-frequency hydrogen plasmas

    Zorat, Roberto

    2003-01-01

    In this thesis low temperature hydrogen radio frequency (rf) plasmas discharges are modelled numerically by using both a global model and a Particle- In-Cell (PIC) simulation. Such plasmas are of interest because of their industrial applications and for the development of negative ion sources for fusion plasmas. The global model technique was adapted and then implemented to model rf inductively coupled hydrogen plasma discharges created in the DENISE experiment, with particular attention to t...

  16. Numerical modelling of tunnel construction in anisotropic foliated soft rock

    Markovič, Jernej

    2009-01-01

    The present work focuses on the influence on tunnelling in the anisotropic foliated soft rock. The excavation initiates stress redistribution around an opening and thus causes the deformation to occur. The numerical problem of the tunnel excavation was modelled in the Plaxis 2D code using different soil constitutive models for modelling the rock mass behaviour. A parametric study was performed to obtain the model response to alteration of the rock mass parameters. The analysis was divided int...

  17. Numerical Model of Radical Photopolymerization Based on Interdiffusion

    Shuhei Yoshida; Yosuke Takahata; Shuma Horiuchi; Hiroyuki Kurata; Manabu Yamamoto

    2014-01-01

    An accurate reaction model is required to analyze the characteristics of photopolymers. For this purpose, we propose a numerical model for radical photopolymerization. In the proposed model, elementary reactions such as initiation, propagation, and termination are considered, and we assume interdiffusion for each component in the material. We analyzed the diffraction characteristics of a radical photopolymer based on the proposed interdiffusion model with the beam propagation method. Moreover...

  18. Numerical Modeling of Electromagnetic Field Effects on the Human Body

    Zuzana Psenakova

    2006-01-01

    Full Text Available Interactions of electromagnetic field (EMF with environment and with tissue of human beings are still under discussion and many research teams are investigating it. The human simulation models are used for biomedical research in a lot of areas, where it is advantage to replace real human body (tissue by the numerical model. Biological effects of EMF are one of the areas, where numerical models are used with many advantages. On the other side, this research is very specific and it is always quite hard to simulate realistic human tissue. This paper deals with different possibilities of numerical modelling of electromagnetic field effects on the human body (especially calculation of the specific absorption rate (SAR distribution in human body and thermal effect.

  19. Numerical modeling of surf beat generated by moving breakpoint

    2009-01-01

    As an important hydrodynamic phenomenon in the nearshore zone, the cross-shore surf beat is numerically studied in this paper with a fully nonlinear Boussinesq-type model, which resolves the primary wave motion as well as the long waves. Compared with the classical Boussinesq equations, the equations adopted here allow for improved linear dispersion characteristics. Wave breaking and run-up in the swash zone are included in the numerical model. Mutual interactions between short waves and long waves are inherent in the model. The numerical study of long waves is based on bichromatic wave groups with a wide range of mean frequencies, group frequencies and modulation rates. The cross-shore variation in the amplitudes of short waves and long waves is investigated. The model results are compared with laboratory experiments from the literature and good agreement is found.

  20. Numerical modeling to investigate slopes and mass flow phenomena

    Heinz Konietzky; Lei NIE; Youhong SUN

    2006-01-01

    An overview is given about up-to-date techniques for slope stability and deformation analysis as well as mass flow phenomena simulation. The paper concentrates on a few aspects in respect to the use of numerical modeling techniques, especially in relation to the shear strength reduction techniques, discontinuum modeling, probabilistic concepts, the combination of GIS and numerical modeling as well as sophisticated hydro-mechanical coupling with time-dependent material behavior. At present these topics are preferred topics of scientific and technical research.

  1. Dimensionless numerical model for simulation of active magnetic regenerator refrigerator

    Sarlah, A.; Poredos, A. [Faculty of Mechanical Engineering, University of Ljubljana, Askerceva 6, SI-1000 Ljubljana (Slovenia)

    2010-09-15

    In order to obtain a better reliability, consistency and accuracy of results obtained with a numerical simulation of an AMRR (active magnetic regenerator refrigerator), a dimensionless numerical model was developed, which can equally be used for determination of regenerator's heat transfer coefficient and simulation of passive heat regenerators or AMRR operation. Regenerator's heat transfer coefficient {alpha}{sub f}, is a crucial input parameter in the simulation of AMRR operation and has a primal effect on the outcome of a solution. This paper deals with a derived dimensionless model and discusses errors involved when using different models for heat transfer coefficient and AMRR operation simulation. (author)

  2. Numerical modelling of the jet nozzle enrichment process

    A numerical model was developed for the simulation of the isotopic enrichment produced by the jet nozzle process. The flow was considered stationary and under ideal gas conditions. The model calculates, for any position of the skimmer piece: (a) values of radial mass concentration profiles for each isotopic species and (b) values of elementary separation effect (Σ sub(A)) and uranium cut (theta). The comparison of the numerical results obtained with the experimental values given in the literature proves the validity of the present work as an initial step in the modelling of the process. (Author)

  3. Hysteresis model of magnetostrictive actuators and its numerical realization

    TANG Zhi-feng; LV Fu-zai; XIANG Zhan-qin

    2007-01-01

    This paper presents two numerical realization of Preisach model by Density Function Method (DFM) and F Function Method (FFM) for a giant magnetostrictive actuator (GMA). Experiment and simulation showed that FFM is better than DFM for predicting precision of hysteresis loops. Lagrange bilinear interpolation algorithm is used in Preisach numerical realization to enhance prediction performance. A set of hysteresis loops and higher order reversal curves are predicted and experimentally verified. The good agreement between the measured and predicted curves shows that the classical Preisach model is effective for modelling the quasi-static hysteresis of the GMA.

  4. On Numerical Modeling of Corporate Strategic Debt Service

    Li, Quan

    2009-01-01

    The point(s) which divide the state variables space into different phases and give rise to different debt service flows based on game theory can be determined either analytically for perpetual debt, or numerically for debt of finite maturity. The fact that debt service flow can vary through its life span and that such a change can be modeled is of importance to credit risk management. Numerical implementation of the strategic debt service theory poses many challenging aspects. The aim of this...

  5. Numerical modeling in electroporation-based biomedical applications

    Pavšelj, Nataša; Miklavčič, Damijan

    2015-01-01

    Background. Numerous experiments have to be performed before a biomedical application is put to practical use in clinical environment. As a complementary work to in vitro, in vivo and medical experiments, we can use analytical and numerical models to represent, as realistically as possible, real biological phenomena of, in our case, electroporation. In this way we canevaluate different electrical parameters in advance, such as pulse amplitude, duration, number of pulses, or different electrod...

  6. Numerical modeling in electroporation-based biomedical applications:

    Miklavčič, Damijan; Pavšelj, Nataša

    2008-01-01

    Background. Numerous experiments have to be performed before a biomedical application is put to practical use in clinical environment. As a complementary work to in vitro, in vivo and medical experiments, we can use analytical and numerical models to represent, as realistically as possible, real biological phenomena of, in our case, electroporation. In this way we canevaluate different electrical parameters in advance, such as pulse amplitude, duration, number of pulses, or different electrod...

  7. Local Duality for 2-Dimensional Local Ring

    Belgacem Draouil

    2008-11-01

    We prove a local duality for some schemes associated to a 2-dimensional complete local ring whose residue field is an -dimensional local field in the sense of Kato–Parshin. Our results generalize the Saito works in the case =0 and are applied to study the Bloch–Ogus complex for such rings in various cases.

  8. Numerical human model for impact and seating comfort

    Hoof, J.F.A.M. van; Lange, R. de; Verver, M.M.

    2003-01-01

    This paper presents a detailed numerical model of the human body that can be used to evaluate both safety and comfort aspects of vehicle interiors. The model is based on a combination of rigid body and finite element techniques to provide an optimal combination of computational efficiency and accura

  9. A Climate System Model, Numerical Simulation and Climate Predictability

    ZENG Qingcun; WANG Huijun; LIN Zhaohui; ZHOU Guangqing; YU Yongqiang

    2007-01-01

    @@ The implementation of the project has lasted for more than 20 years. As a result, the following key innovative achievements have been obtained, ranging from the basic theory of climate dynamics, numerical model development and its related computational theory to the dynamical climate prediction using the climate system models:

  10. PROBABILITY MODELS FOR OBTAINING NON-NUMERICAL DATA

    Orlov A. I.

    2015-01-01

    Full Text Available The statistics of objects of non-numerical nature (statistics of non-numerical objects, non-numerical data statistics, non-numeric statistics is the area of mathematical statistics, devoted to the analysis methods of non-numeric data. Basis of applying the results of mathematical statistics are probabilistic-statistical models of real phenomena and processes, the most important (and often only which are models for obtaining data. The simplest example of a model for obtaining data is the model of the sample as a set of independent identically distributed random variables. In this article we have considered the basic probabilistic models for obtaining non-numeric data. Namely, the models of dichotomous data, results of paired comparisons, binary relations, ranks, the objects of general nature. We have discussed the various options of probabilistic models and their practical use. For example, the basic probabilistic model of dichotomous data - Bernoulli vector (Lucian i.e. final sequence of independent Bernoulli trials, for which the probabilities of success may be different. The mathematical tools of solutions of various statistical problems associated with the Bernoulli vectors are useful for the analysis of random tolerances; random sets with independent elements; in processing the results of independent pairwise comparisons; statistical methods for analyzing the accuracy and stability of technological processes; in the analysis and synthesis of statistical quality control plans (for dichotomous characteristics; the processing of marketing and sociological questionnaires (with closed questions like "yes" - "no"; the processing of socio-psychological and medical data, in particular, the responses to psychological tests such as MMPI (used in particular in the problems of human resource management, and analysis of topographic maps (used for the analysis and prediction of the affected areas for technological disasters, distributing corrosion

  11. Entropic analysis of evolving 2-dimensional scalar fields: Field-dependent MFM datasets and non-local random-switch hysteresis

    Different means of applying numerical techniques are proposed, typically used for pattern recognition purposes and for artificial texture classification and segmentation, to study the evolution of magnetization at the nano- and micro-scale in ferromagnetic samples. This framework is applied both to the Magnetic Force Microscope (MFM) dataset and to synthetic binary patterns of 2-dimensional spins, based on a nucleation and growth model for the hysteresis. In general 2-dimensional gray-scale frames carry information characterized by a spatial correlation between each pixel. Multiple datasets as the one analyzed, either 256-bit-valued magnetic field dependent MFM matrices or 2-bit-valued synthetic field dependent scatter matrices, are arranged as multidimensional arrays and their 2-dimensional entropy is computed.

  12. Ensemble-type numerical uncertainty information from single model integrations

    Rauser, Florian, E-mail: florian.rauser@mpimet.mpg.de; Marotzke, Jochem; Korn, Peter

    2015-07-01

    We suggest an algorithm that quantifies the discretization error of time-dependent physical quantities of interest (goals) for numerical models of geophysical fluid dynamics. The goal discretization error is estimated using a sum of weighted local discretization errors. The key feature of our algorithm is that these local discretization errors are interpreted as realizations of a random process. The random process is determined by the model and the flow state. From a class of local error random processes we select a suitable specific random process by integrating the model over a short time interval at different resolutions. The weights of the influences of the local discretization errors on the goal are modeled as goal sensitivities, which are calculated via automatic differentiation. The integration of the weighted realizations of local error random processes yields a posterior ensemble of goal approximations from a single run of the numerical model. From the posterior ensemble we derive the uncertainty information of the goal discretization error. This algorithm bypasses the requirement of detailed knowledge about the models discretization to generate numerical error estimates. The algorithm is evaluated for the spherical shallow-water equations. For two standard test cases we successfully estimate the error of regional potential energy, track its evolution, and compare it to standard ensemble techniques. The posterior ensemble shares linear-error-growth properties with ensembles of multiple model integrations when comparably perturbed. The posterior ensemble numerical error estimates are of comparable size as those of a stochastic physics ensemble.

  13. On the numerical modeling of shear banding around an excavation

    A benchmark has been proposed by the research group of MoMas concerning the numerical prediction of the extend and evolution of the Excavation Damaged Zones (EDZ) around the nuclear waste deep geological disposals. The final objective is the evaluation of the relevance of the numerical codes used by the participants to model the different tests of the benchmark with the same constitutive equation. The different scientific teams participating are listed in a companion paper by Chavant C and Fernandes R. In the proposed benchmark exercises, thermal aspects and creep effects are ignored. The modelling focuses on the coupling between mechanical behaviour of the host formation and the hydraulic diffusion (in saturated conditions) during the excavation process. The proposed constitutive laws should model two in-situ observations: the degradation of the rock mechanical properties and the induced increase of medium permeability. The mechanical constitutive model is an associated Drucker Prager elastoplastic model with a softening due to the degradation of the cohesion related with the deviatoric plastic strain. Softening models which have been extensively used in the past are now known to generate numerical problems especially when localized ruptures modes (like shear bands) appear in a computation. Numerical modelling of such kind of localized shear zones suffers by a pathological sensitivity to the finite element discretization. Around a tunnel, this latter mesh sensitivity leads to undefined position and number of shear bands. In a coupled hydro-mechanical simulation, this renders more difficult the modelling of the increase of the permeability in the damage zone (EDZ) and more specifically in the localized bands. We present first the numerical computations obtained with the proposed model and the proposed boundary Conditions. (authors)

  14. Numerical Modeling of Electrodynamic Aggregation of Magnetized Nanodust

    Neverov, V. S.; Kukushkin, A. B.

    2011-01-01

    The recent results of applying the parallel numerical code SELFAS-3 to modelling of electrodynamic aggregation of magnetized nanodust are presented. The modelling describes evolution of a many-body system of basic blocks which are taken as strongly magnetized thin rods (i.e., one-dimensional static magnetic dipoles), with electric conductivity and static electric charge, screened with its own plasma sheath. The code provides continuous modelling of the following stages of evolution: (i) align...

  15. Numerical modeling of complex heat transfer phenomena in cooling applications

    Hou, Xiaofei

    2015-01-01

    Multiphase and multicomponent flows are frequently encountered in the cooling applications due to combined heat transfer and phase change phenomena. Two-fluid and homogeneous mixture models are chosen to numerically study these flows in the cooling phenomena. Therefore this work is divided in two main parts. In the first part, a two-fluid model algorithm for free surface flows is presented. The two fluid model is usually used as a tool to simulate dispersed flow. With its extension, it may al...

  16. Numerical modeling for heatsink emissions in power electronics

    J. Kulanayagam; Hagmann, J. H.; Schenke, S.; K. F. Hoffmann; S. Dickmann

    2012-01-01

    The parasitic coupling between power semiconductors and the heat sink is responsible for noise current in Switching Mode Power Supply (SMPS) systems. In this paper, the variations in the radiation characteristics of heatsinks are investigated with respect to their geometries by use of numerical models. Analyses are facilitated by using a mopole antenna as an EMI receiver and by using simplified heatsink models as EMI transmitters to model the heatsink radiated emissions. In ...

  17. Ocean wave prediction using numerical and neural network models

    Mandal, S.; Prabaharan, N.

    formula- tions [1]. Subsequently, many first and second generation numerical models were developed during last few decades and presently, a few including third generation operational wave models exist for wave forecasting. With the rapid de- velopments... not be fully analyzed till late 1970’s, was the non-linear wave- wave interaction [2]. The third generation operational wave models use the components of the source function without any prior restrictions on the spectral shape [3]. Even though the third...

  18. Numerical Models of Sewage Dispersion and Statistica Bathing Water Standards

    Petersen, Ole; Larsen, Torben

    1991-01-01

    As bathing water standards usually are founded in statistical methods, the numerical models used in outfall design should reflect this. A statistical approach, where stochastic variations in source strength and bacterial disappearance is incorporated into a numerical dilution model is presented. It...... is demonstrated for a specific outfall how the method can be used to estimate the bathing water quality. The ambition with the paper has been to demonstrate how stochastic variations in a simple manner can be included in the analysis of water quality....

  19. Feedbacks Between Numerical and Analytical Models in Hydrogeology

    Zlotnik, V. A.; Cardenas, M. B.; Toundykov, D.; Cohn, S.

    2012-12-01

    Hydrogeology is a relatively young discipline which combines elements of Earth science and engineering. Mature fundamental disciplines (e.g., physics, chemistry, fluid mechanics) have centuries-long history of mathematical modeling even prior to discovery of Darcy's law. Thus, in hydrogeology, relatively few classic analytical models (such those by Theis, Polubarinova-Kochina, Philip, Toth, Henry, Dagan, Neuman) were developed by the early 1970's. The advent of computers and practical demands refocused mathematical models towards numerical techniques. With more diverse but less mathematically-oriented training, most hydrogeologists shifted from analytical methods to use of standardized computational software. Spatial variability in internal properties and external boundary conditions and geometry, and the added complexity of chemical and biological processes will remain major challenges for analytical modeling. Possibly, analytical techniques will play a subordinate role to numerical approaches in many applications. On the other hand, the rise of analytical element modeling of groundwater flow is a strong alternative to numerical models when data demand and computational efficiency is considered. The hallmark of analytical models - transparency and accuracy - will remain indispensable for scientific exploration of complex phenomena and for benchmarking numerical models. Therefore, there will always be feedbacks and complementarities between numerical and analytical techniques, as well as a certain ideological schism among various views to modeling. We illustrate the idea of feedbacks by reviewing evolution of Joszef Toth's analytical model of gravity driven flow systems. Toth's (1963) approach was to reduce the flow domain to a rectangle which allowed for closed-form solution of the governing equations. Succeeding numerical finite-element models by Freeze and Witherspoon (1966-1968) explored the effects of geometry and heterogeneity on regional groundwater flow

  20. Non-stationary iterative methods for solving macroeconomic numeric models

    Bogdan OANCEA

    2006-01-01

    Full Text Available Macroeconometric modeling was influenced by the development of new and efficient computational techniques. Rational Expectations models, a particular class of macroeconometric models, give raise to very large systems of equations, the solution of which requires heavy computations. Therefore, such models are an interesting testing ground for the numerical methods addressed in this research. The most difficult problem is to obtain the solution of the linear system that arises during the Newton step. As an alternative to the direct methods, we propose non-stationary iterative methods, also called Krylov methods, to solve these models. Numerical experiments conducted by authors confirm the interesting features of these methods: low computational complexity and storage requirements.

  1. Numerical Simulation and Cold Modeling experiments on Centrifugal Casting

    Keerthiprasad, Kestur Sadashivaiah; Murali, Mysore Seetharam; Mukunda, Pudukottah Gopaliengar; Majumdar, Sekhar

    2011-02-01

    In a centrifugal casting process, the fluid flow eventually determines the quality and characteristics of the final product. It is difficult to study the fluid behavior here because of the opaque nature of melt and mold. In the current investigation, numerical simulations of the flow field and visualization experiments on cold models have been carried out for a centrifugal casting system using horizontal molds and fluids of different viscosities to study the effect of different process variables on the flow pattern. The effects of the thickness of the cylindrical fluid annulus formed inside the mold and the effects of fluid viscosity, diameter, and rotational speed of the mold on the hollow fluid cylinder formation process have been investigated. The numerical simulation results are compared with corresponding data obtained from the cold modeling experiments. The influence of rotational speed in a real-life centrifugal casting system has also been studied using an aluminum-silicon alloy. Cylinders of different thicknesses are cast at different rotational speeds, and the flow patterns observed visually in the actual castings are found to be similar to those recorded in the corresponding cold modeling experiments. Reasonable agreement is observed between the results of numerical simulation and the results of cold modeling experiments with different fluids. The visualization study on the hollow cylinders produced in an actual centrifugal casting process also confirm the conclusions arrived at from the cold modeling experiments and numerical simulation in a qualitative sense.

  2. Development, validation and application of numerical space environment models

    Honkonen, Ilja

    2013-10-01

    Currently the majority of space-based assets are located inside the Earth's magnetosphere where they must endure the effects of the near-Earth space environment, i.e. space weather, which is driven by the supersonic flow of plasma from the Sun. Space weather refers to the day-to-day changes in the temperature, magnetic field and other parameters of the near-Earth space, similarly to ordinary weather which refers to changes in the atmosphere above ground level. Space weather can also cause adverse effects on the ground, for example, by inducing large direct currents in power transmission systems. The performance of computers has been growing exponentially for many decades and as a result the importance of numerical modeling in science has also increased rapidly. Numerical modeling is especially important in space plasma physics because there are no in-situ observations of space plasmas outside of the heliosphere and it is not feasible to study all aspects of space plasmas in a terrestrial laboratory. With the increasing number of computational cores in supercomputers, the parallel performance of numerical models on distributed memory hardware is also becoming crucial. This thesis consists of an introduction, four peer reviewed articles and describes the process of developing numerical space environment/weather models and the use of such models to study the near-Earth space. A complete model development chain is presented starting from initial planning and design to distributed memory parallelization and optimization, and finally testing, verification and validation of numerical models. A grid library that provides good parallel scalability on distributed memory hardware and several novel features, the distributed cartesian cell-refinable grid (DCCRG), is designed and developed. DCCRG is presently used in two numerical space weather models being developed at the Finnish Meteorological Institute. The first global magnetospheric test particle simulation based on the

  3. Frame Design and Reality of Numerical Model for Sculptured Part Machining

    2000-01-01

    The importance of the numerical model for sculptured part machining based on virtual environment is introduced. Meanwhile, the general frame of the numerical model is proposed, and the techniques of developing the numerical model are discussed in detail.

  4. Physical and numerical modeling of Joule-heated melters

    The Joule-heated ceramic-lined melter is an integral part of the high level waste immobilization process under development by the US Department of Energy. Scaleup and design of this waste glass melting furnace requires an understanding of the relationships between melting cavity design parameters and the furnace performance characteristics such as mixing, heat transfer, and electrical requirements. Developing empirical models of these relationships through actual melter testing with numerous designs would be a very costly and time consuming task. Additionally, the Pacific Northwest Laboratory (PNL) has been developing numerical models that simulate a Joule-heated melter for analyzing melter performance. This report documents the method used and results of this modeling effort. Numerical modeling results are compared with the more conventional, physical modeling results to validate the approach. Also included are the results of numerically simulating an operating research melter at PNL. Physical Joule-heated melters modeling results used for qualiying the simulation capabilities of the melter code included: (1) a melter with a single pair of electrodes and (2) a melter with a dual pair (two pairs) of electrodes. The physical model of the melter having two electrode pairs utilized a configuration with primary and secondary electrodes. The principal melter parameters (the ratio of power applied to each electrode pair, modeling fluid depth, electrode spacing) were varied in nine tests of the physical model during FY85. Code predictions were made for five of these tests. Voltage drops, temperature field data, and electric field data varied in their agreement with the physical modeling results, but in general were judged acceptable. 14 refs., 79 figs., 17 tabs

  5. Mathematical modeling and numerical simulation of Czochralski Crystal Growth

    Jaervinen, J.; Nieminen, R. [Center for Scientific Computing, Espoo (Finland)

    1996-12-31

    A detailed mathematical model and numerical simulation tools based on the SUPG Finite Element Method for the Czochralski crystal growth has been developed. In this presentation the mathematical modeling and numerical simulation of the melt flow and the temperature distribution in a rotationally symmetric crystal growth environment is investigated. The temperature distribution and the position of the free boundary between the solid and liquid phases are solved by using the Enthalpy method. Heat inside of the Czochralski furnace is transferred by radiation, conduction and convection. The melt flow is governed by the incompressible Navier-Stokes equations coupled with the enthalpy equation. The melt flow is numerically demonstrated and the temperature distribution in the whole Czochralski furnace. (author)

  6. Numerical strategy for model correction using physical constraints

    He, Yanyan; Xiu, Dongbin

    2016-05-01

    In this paper we present a strategy for correcting model deficiency using observational data. We first present the model correction in a general form, involving both external correction and internal correction. The model correction problem is then parameterized and casted into an optimization problem, from which the parameters are determined. More importantly, we discuss the incorporation of physical constraints from the underlying physical problem. Several representative examples are presented, where the physical constraints take very different forms. Numerical tests demonstrate that the physics constrained model correction is an effective way to address model-form uncertainty.

  7. Integrating Numerical Groundwater Modeling Results With Geographic Information Systems

    Witkowski, M. S.; Robinson, B. A.; Linger, S. P.

    2001-12-01

    Many different types of data are used to create numerical models of flow and transport of groundwater in the vadose zone. Results from water balance studies, infiltration models, hydrologic properties, and digital elevation models (DEMs) are examples of such data. Because input data comes in a variety of formats, for consistency the data need to be assembled in a coherent fashion on a single platform. Through the use of a geographic information system (GIS), all data sources can effectively be integrated on one platform to store, retrieve, query, and display data. In our vadoze zone modeling studies in support of Los Alamos National Laboratory's Environmental Restoration Project, we employ a GIS comprised of a Raid storage device, an Oracle database, ESRI's spatial database engine (SDE), ArcView GIS, and custom GIS tools for three-dimensional (3D) analysis. We store traditional GIS data, such as, contours, historical building footprints, and study area locations, as points, lines, and polygons with attributes. Numerical flow and transport model results from the Finite Element Heat and Mass Transfer Code (FEHM) are stored as points with attributes, such as fluid saturation, or pressure, or contaminant concentration at a given location. We overlay traditional types of GIS data with numerical model results, thereby allowing us to better build conceptual models and perform spatial analyses. We have also developed specialized analysis tools to assist in the data and model analysis process. This approach provides an integrated framework for performing tasks such as comparing the model to data and understanding the relationship of model predictions to existing contaminant source locations and water supply wells. Our process of integrating GIS and numerical modeling results allows us to answer a wide variety of questions about our conceptual model design: - Which set of locations should be identified as contaminant sources based on known historical building operations

  8. Numerical investigation of quantumfield model of strong-coupling binucleon

    The quantumfield binucleon model for the case of the nucleon spot interaction with the scalar and pseudoscalar mesons fields is considered. For the strong coupling limit the nonlinear equations describing two nucleons in the meson field are developed. Methods of numerical investigation and obtained results are discussed. 10 refs., 3 figs., 3 tabs

  9. Numerical modelling and experimental assessment of concrete spalling in fire

    Shamalta, M.; Breunese, A.; Peelen, W.; Fellinger, J.

    2005-01-01

    In this paper, the phenomenon of spalling of concrete in fire has been studied using a numerical model. Spalling is the violent or non-violent breaking off of layers or pieces of concrete when it is exposed to high temperatures as experienced in fires. The types and mechanisms of spalling have been

  10. Numerical modeling, calibration, and validation of an ultrasonic separator

    Cappon, H.J.; Keesman, K.J.

    2013-01-01

    Our overall goal is to apply acoustic separation technology for the recovery of valuable particulate matter from wastewater in industry. Such large-scale separator systems require detailed design and evaluation to optimize the system performance at the earliest stage possible. Numerical models can f

  11. Design of advanced industrial furnaces using numerical modeling method

    Dong, Wei

    2000-01-01

    This doctoral thesis describes the fundamentals ofmathematical modeling for the industrial furnaces and boilersand presents the results from the numerical simulations of sometypical applications in advanced industrial furnaces andboilers. The main objective of this thesis work is to employcomputational fluid dynamics (CFD) technology as an effectivecomputer simulation tool to study and develop the newcombustion concepts, phenomena and processes in advancedindustrial furnaces and boilers. The ...

  12. Development of numerical Grids for UZ Flow and Transport Modeling

    P. Dobson

    2004-08-31

    This report describes the methods used to develop numerical grids of the unsaturated hydrogeologic system beneath Yucca Mountain, Nevada. Numerical grid generation is an integral part of the development of the unsaturated zone (UZ) flow and transport model, a complex, three-dimensional (3-D) model of Yucca Mountain. This revision contains changes made to improve the clarity of the description of grid generation. The numerical grids, developed using current geologic, hydrogeologic, and mineralogic data, provide the necessary framework to: (1) develop calibrated hydrogeologic property sets and flow fields, (2) test conceptual hypotheses of flow and transport, and (3) predict flow and transport behavior under a variety of climatic and thermal-loading conditions. The technical scope, content, and management for the current revision of this report are described in the planning document ''Technical Work Plan for: Unsaturated Zone Flow Analysis and Model Report Integration'' (BSC 2004 [DIRS 169654], Section 2). Grids generated and documented in this report supersede those documented in Revision 00 of this report, ''Development of Numerical Grids for UZ Flow and Transport Modeling'' (BSC 2001 [DIRS 159356]). The grids presented in this report are the same as those developed in Revision 01 (BSC 2003 [DIRS 160109]); however, the documentation of the development of the grids in Revision 02 has been updated to address technical inconsistencies and achieve greater transparency, readability, and traceability. The constraints, assumptions, and limitations associated with this report are discussed in the appropriate sections that follow.

  13. Numerical model of human cardiovascular system-Korotkoff's sound simulation

    Maršík, František; Převorovská, Světlana; Štembera, Vítězslav

    Graz, Rakousko: Departments of Mathematics and Physiology University of Graz, 2003 - (Kappel, F.; Batzel, J.; Fink, M.; Schneditz, D.). s. 24 [Workshop on Cardiovascular, Respiratory and Metabolic Control Modeling. 11.06.2003-14.06.2003, Graz] Institutional research plan: CEZ:AV0Z2076919 Keywords : cardiovascular system * Korotkoff's sound * numerical simulation Subject RIV: ED - Physiology

  14. Numerical Modeling of Electromagnetic Field in the Biological Cell

    Vlachová-Hutová, E.; Kříž, P.; Gescheidtová, E.; Bartušek, Karel

    Cambridge: The Electromagnetics Academy, 2014, s. 1890-1894. ISBN 978-1-934142-28-8. [PIERS 2014. Progress In Electromagnetics Research Symposium /35./. Guangzhou (CN), 25.08.2014-28.08.2014] R&D Projects: GA ČR GAP102/12/1104 Institutional support: RVO:68081731 Keywords : electromagnetic field * iological ell * numerical modeling Subject RIV: JA - Electronics ; Optoelectronics, Electrical Engineering

  15. Numerical Modeling of Accuracy of Air Ion Field Measurement

    Bartušek, Karel; Fiala, P.; Bachorec, T.; Kadlecová, E.

    Cambridge : The Electromagnetic Academy, 2007, s. 578-581. ISBN 978-1-934142-00-4. [Progress in Electromagnetics Research Symposium - PIERS 2007. Beijing (CN), 26.03.2007-20.03.2007] Institutional research plan: CEZ:AV0Z20650511 Keywords : air ion * numerical modeling Subject RIV: JA - Electronics ; Optoelectronics, Electrical Engineering

  16. Numerical modeling of secondary side thermohydraulics of horizontal steam generator

    Melikhov, V.I.; Melikhov, O.I.; Nigmatulin, B.I. [Research and Engineering Centre of LWR Nuclear Plants Safety, Moscow (Russian Federation)

    1995-12-31

    A mathematical model for the transient three-dimensional secondary side thermal hydraulics of the horizontal steam generator has been developed. The calculations of the steam generator PGV-1000 and PGV-4 nominal regimes and comparison of numerical and experimental results have been carried out. 7 refs.

  17. Numeral eddy current sensor modelling based on genetic neural network

    Yu A-Long

    2008-01-01

    This paper presents a method used to the numeral eddy current sensor modelling based on the genetic neural network to settle its nonlinear problem. The principle and algorithms of genetic neural network are introduced. In this method, the nonlinear model parameters of the numeral eddy current sensor are optimized by genetic neural network (GNN) according to measurement data. So the method remains both the global searching ability of genetic algorithm and the good local searching ability of neural network. The nonlinear model has the advantages of strong robustness,on-line modelling and high precision.The maximum nonlinearity error can be reduced to 0.037% by using GNN.However, the maximum nonlinearity error is 0.075% using the least square method.

  18. Physicochemical and numerical modeling of electrokinetics in inhomogenous matrices

    Paz-Garcia, Juan Manuel

    into porous solid matrices of different kinds. These techniques are typically denoted as electrokinetic treatments. In these kind of electrochemically-induced transport processes, the driving force is related the concentration gradients and the unbalanced in ionic charge produced by the electrochemical...... the corroboration of the assumptions and, eventually, the development of prediction tools for optimizing the different treatments. In accordance to this, the primary objectives proposed for the present project “Physicochemical and numerical modeling of electrokinetics in inhomogeneous matrices” are: 1.to formulate...... processes. Furthermore, the physicochemical model is described in a generalized manner, so it can be used to simulate a wide range of electrokinetic treatments for different inhomogeneous matrices. Anumericalmodelhasbeenbuiltforthemathematicalsolutionoftheproposed physicochemical model. The numerical model...

  19. 2D numerical modelling of meandering channel formation

    XIAO, Y.; ZHOU, G.; YANG, F. S.

    2016-03-01

    A 2D depth-averaged model for hydrodynamic sediment transport and river morphological adjustment was established. The sediment transport submodel takes into account the influence of non-uniform sediment with bed surface armoring and considers the impact of secondary flow in the direction of bed-load transport and transverse slope of the river bed. The bank erosion submodel incorporates a simple simulation method for updating bank geometry during either degradational or aggradational bed evolution. Comparison of the results obtained by the extended model with experimental and field data, and numerical predictions validate that the proposed model can simulate grain sorting in river bends and duplicate the characteristics of meandering river and its development. The results illustrate that by using its control factors, the improved numerical model can be applied to simulate channel evolution under different scenarios and improve understanding of patterning processes.

  20. Accounting for Errors in Model Analysis Theory: A Numerical Approach

    Sommer, Steven R.; Lindell, Rebecca S.

    2004-09-01

    By studying the patterns of a group of individuals' responses to a series of multiple-choice questions, researchers can utilize Model Analysis Theory to create a probability distribution of mental models for a student population. The eigenanalysis of this distribution yields information about what mental models the students possess, as well as how consistently they utilize said mental models. Although the theory considers the probabilistic distribution to be fundamental, there exists opportunities for random errors to occur. In this paper we will discuss a numerical approach for mathematically accounting for these random errors. As an example of this methodology, analysis of data obtained from the Lunar Phases Concept Inventory will be presented. Limitations and applicability of this numerical approach will be discussed.

  1. A reference dataset for verifying numerical electrophysiological heart models

    Paetsch Ingo

    2011-01-01

    Full Text Available Abstract Background The evaluation, verification and comparison of different numerical heart models are difficult without a commonly available database that could be utilized as a reference. Our aim was to compile an exemplary dataset. Methods The following methods were employed: Magnetic Resonance Imaging (MRI of heart and torso, Body Surface Potential Maps (BSPM and MagnetoCardioGraphy (MCG maps. The latter were recorded simultaneously from the same individuals a few hours after the MRI sessions. Results A training dataset is made publicly available; datasets for blind testing will remain undisclosed. Conclusions While the MRI data may provide a common input that can be applied to different numerical heart models, the verification and comparison of different models can be performed by comparing the measured biosignals with forward calculated signals from the models.

  2. A numerical framework for modelling floating wind turbines

    Vire, Axelle; Xiang, Jiansheng; Piggott, Matthew; Latham, John-Paul; Pain, Christopher

    2012-11-01

    This work couples a fluid/ocean- and a solid- dynamics model in order to numerically study fluid-structure interactions. The fully non-linear Navier-Stokes and solid-dynamics equations are solved on two distinct finite-element and unstructured grids. The interplay between fluid and solid is represented through a penalty force in the momentum balances of each material. The present algorithm is novel in that it spatially conserves the discrete penalty force, when exchanging it between both models, independently of the mesh resolution and of the shape-function orders in each model. This numerical framework targets the modelling of offshore floating wind turbines. Results will be shown for the flow past a moving pile and an actuator-disk representation of a turbine. This research is supported by the European Union Seventh Framework Programme (grant agreement PIEF-GA-2010-272437).

  3. Numeral eddy current sensor modelling based on genetic neural network

    This paper presents a method used to the numeral eddy current sensor modelling based on the genetic neural network to settle its nonlinear problem. The principle and algorithms of genetic neural network are introduced. In this method, the nonlinear model parameters of the numeral eddy current sensor are optimized by genetic neural network (GNN) according to measurement data. So the method remains both the global searching ability of genetic algorithm and the good local searching ability of neural network. The nonlinear model has the advantages of strong robustness, on-line modelling and high precision. The maximum nonlinearity error can be reduced to 0.037% by using GNN. However, the maximum nonlinearity error is 0.075% using the least square method

  4. A comprehensive numerical model of wind-blown sand

    Kok, Jasper F

    2009-01-01

    Wind-blown sand, or "saltation", ejects dust aerosols into the atmosphere, creates sand dunes, and erodes geological features. We present a comprehensive numerical model of steady-state saltation that, in contrast to most previous studies, can simulate saltation over mixed soils. Our model simulates the motion of saltating particles due to gravity, fluid drag, particle spin, fluid shear, and turbulence. Moreover, the model explicitly accounts for the retardation of the wind due to drag from saltating particles. We also developed a physically-based parameterization of the ejection of surface particles by impacting saltating particles which matches experimental results. Our numerical model is the first to reproduce measurements of the wind shear velocity at the impact threshold (i.e., the lowest shear velocity for which saltation is possible) and of the aerodynamic roughness length in saltation. It also correctly predicts a wide range of other saltation processes, including profiles of the wind speed and partic...

  5. Climate system model, numerical simulation and climate predictability

    2010-01-01

    @@ Thanks to its work of past more than 20 years,a research team led by Prof.ZENG Qingcun and Prof.WANG Huijun from the CAS Institute of Atmospheric Physics (IAP) has scored innovative achievements in their studies of basic theory of climate dynamics,numerical model development,its related computational theory,and the dynamical climate prediction using the climate system models.Their work received a second prize of the National Award for Natural Sciences in 2005.

  6. Numerical modeling of tool-rock interaction in laminated formations

    Ouyang, Hsien-Ju

    1990-01-01

    Based on a review of the literature, a numerical model of the chip formation process in laminated rock formations has been developed. The model uses a finite element approach to simulate the anisotropic behavior of laminated rock formations. . The finite element program has been developed with the assumption of two-dimensional plane strain. Anisotropic elements and dynamic loading are used to represent the actual penetration process of a bit tooth. An iteration method, using a...

  7. Interaction of tide and salinity barrier: Limitation of numerical model

    Suphat Vongvisessomjai; Phairot Chatanantavet2; Patchanok Srivihok

    2008-01-01

    Nowadays, the study of interaction of the tide and the salinity barrier in an estuarine area is usually accomplished vianumerical modeling, due to the speed and convenience of modern computers. However, numerical models provide littleinsight with respect to the fundamental physical mechanisms involved. In this study, it is found that all existing numericalmodels work satisfactorily when the barrier is located at some distance far from upstream and downstream boundary conditions.Results are co...

  8. Numerical Detection of Ergodicity Breaking in a Glass Model

    Sasaki, Munetaka; Hukushima, Koji

    2016-07-01

    We present a numerical method of directly detecting ergodicity breaking in glassy systems. To examine the validity of the proposed method, we applied it to the Biroli-Mézard glass model on a regular random graph. The obtained results clearly indicate that the model exhibits a dynamical transition with ergodicity breaking at an occupation density, that is consistent with the prediction obtained by the cavity method. The present method is applicable to glassy systems in finite dimensions.

  9. Numerical modelling of compressible viscous flow in turbine cascades

    Louda, P.; Kozel, K.; Příhoda, Jaromír

    Vol. 490. Bristol : IOP Publishing, 2014. ISSN 1742-6588. [International Conference on Mathematical Modeling in Physical Sciences 2013 /2./. Praha (CZ), 01.09.2013-05.09.2013] R&D Projects: GA ČR(CZ) GAP101/10/1329; GA ČR GAP101/12/1271 Institutional support: RVO:61388998 Keywords : numerical simulation * turbulence modelling * turbomachinery flow Subject RIV: BK - Fluid Dynamics

  10. Numerical optimisation for model evaluation in combustion kinetics

    Fischer, Marc; Jiang, Xi

    2015-01-01

    Numerical optimisation related to the estimation of kinetic parameters and model evaluation is playing an increasing role in combustion as well as in other areas of applied energy research. The present work aims at presenting the current probability-based approaches along applications to real problems of combustion chemical kinetics. The main methods related to model and parameter evaluation have been explicated. An in-house program for the systematic adjustment of kinetic parameters to exper...