Exercise sheet 3 pdf

Exercises and solutions in PDF

An HMV mapping study of a pulsating jet system

A new system for the generation of hydrodynamic modulated voltammetry (HMV) is presented. This system consists of an oscillating jet produced through the mechanical vibration of a large diaphragm. The structure of the cell is such that a relatively small vibration is transferred to a large fluid flow at the jet outlet. Positioning of an electrode (Pt, 0.5 mm or 25 mu m diameter) over the exit of this jet enables the detection of the modulated flow of liquid. While this flow creates modest mass transfer rates (time averaged similar to 0.015 cm s(-1)) it can also be used to create a HMV system where a 'lock-in' approach is adopted to investigate the redox chemistry in question. This is demonstrated for the Fe(CN)(6)(3-/4-) redox system. Here 'lock-in' to the modulated hydrodynamic signal is achieved through the deployment of bespoke software. The apparatus and procedure is shown to produce a simple and efficient way to obtain the desired signal. In addition the spatial variation of the HMV signal, phase correction and time averaged current with respect to the jet orifice is presented. (C) 2008 Elsevier B.V. All rights reserved.

Approximate Riemann solutions of the two-dimensional shallow-water equations

A finite-difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow-water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gasdynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. An extension to the two-dimensional equations with source terms, is included. The scheme is applied to a dam-break problem with cylindrical symmetry.

Efficient shock capturing for isentropic flows using arithmetic averaging

A shock capturing scheme is presented for the equations of isentropic flow based on upwind differencing applied to a locally linearized set of Riemann problems. This includes the two-dimensional shallow water equations using the familiar gas dynamics analogy. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency, leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver where the computational expense can be prohibitive. The scheme is applied to a two-dimensional dam-break problem and the approximate solution compares well with those given by other authors.

An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

Results of an aqueous source term model for a radiological risk assessment of the Drigg LLW Site, U.K.

A radionuclide source term model has been developed which simulates the biogeochemical evolution of the Drigg low level waste (LLW) disposal site. The DRINK (DRIgg Near field Kinetic) model provides data regarding radionuclide concentrations in groundwater over a period of 100,000 years, which are used as input to assessment calculations for a groundwater pathway. The DRINK model also provides input to human intrusion and gaseous assessment calculations through simulation of the solid radionuclide inventory. These calculations are being used to support the Drigg post closure safety case. The DRINK model considers the coupled interaction of the effects of fluid flow, microbiology, corrosion, chemical reaction, sorption and radioactive decay. It represents the first direct use of a mechanistic reaction-transport model in risk assessment calculations.

State estimation using model order reduction for unstable systems

The problem of state estimation occurs in many applications of fluid flow. For example, to produce a reliable weather forecast it is essential to find the best possible estimate of the true state of the atmosphere. To find this best estimate a nonlinear least squares problem has to be solved subject to dynamical system constraints. Usually this is solved iteratively by an approximate Gauss–Newton method where the underlying discrete linear system is in general unstable. In this paper we propose a new method for deriving low order approximations to the problem based on a recently developed model reduction method for unstable systems. To illustrate the theoretical results, numerical experiments are performed using a two-dimensional Eady model – a simple model of baroclinic instability, which is the dominant mechanism for the growth of storms at mid-latitudes. It is a suitable test model to show the benefit that may be obtained by using model reduction techniques to approximate unstable systems within the state estimation problem.

Similarity solutions describing the melting of a mushy layer

A model of the melting of a mushy region in the absence of fluid flow is presented. Similarity solutions are obtained which are used to describe melting from a hot plate with and without the generation of a completely molten region. These solutions are extended to describe the melting of a mushy region in contact with a hot liquid. A significant feature of melting mushy regions is that the phase change occurs internally by dissolution. Our solutions for melting of a mushy region are used to investigate this internal phase change and are compared with the classical Neumann solutions for melting of a pure substance.

The non-conservation of potential vorticity by a dynamical core compared with the effects of parametrized physical processes

Numerical models of the atmosphere combine a dynamical core, which approximates solutions to the adiabatic, frictionless governing equations for fluid dynamics, with tendencies arising from the parametrization of other physical processes. Since potential vorticity (PV) is conserved following fluid flow in adiabatic, frictionless circumstances, it is possible to isolate the effects of non-conservative processes by accumulating PV changes in an air-mass relative framework. This “PV tracer technique” is used to accumulate separately the effects on PV of each of the different non-conservative processes represented in a numerical model of the atmosphere. Dynamical cores are not exactly conservative because they introduce, explicitly or implicitly, some level of dissipation and adjustment of prognostic model variables which acts to modify PV. Here, the PV tracers technique is extended to diagnose the cumulative effect of the non-conservation of PV by a dynamical core and its characteristics relative to the PV modification by parametrized physical processes. Quantification using the Met Office Unified Model reveals that the magnitude of the non-conservation of PV by the dynamical core is comparable to those from physical processes. Moreover, the residual of the PV budget, when tracing the effects of the dynamical core and physical processes, is at least an order of magnitude smaller than the PV tracers associated with the most active physical processes. The implication of this work is that the non-conservation of PV by a dynamical core can be assessed in case studies with a full suite of physics parametrizations and directly compared with the PV modification by parametrized physical processes. The nonconservation of PV by the dynamical core is shown to move the position of the extratropical tropopause while the parametrized physical processes have a lesser effect at the tropopause level.

Osteocyte Density in the Peri-Implant Bone of Immediately Loaded and Submerged Dental Implants

Background: The role of osteocytes in bone structure and function remains partially unresolved. Their participation in mechanotransduction, i.e., the conversion of a physical stimulus into a cellular response, has been hypothesized. The present study was an evaluation of the osteocyte density in the peri-implant bone of immediately loaded and submerged dental implants. Methods: Fourteen male patients were included in the study; all of them were partially edentulous and needed a posterior mandibular restoration. Implants were inserted in these areas; half of the sample was loaded immediately (included in a fixed provisional prosthesis on the same day as implant surgery), whereas the other half was left to heal submerged. Fourteen implants (seven immediately loaded and seven unloaded) were retrieved with a trephine after a healing period of 8 weeks. The specimens were treated to obtain thin ground sections, and histomorphometry was used to evaluate the osteocyte index in the peri-implant bone. Results: A higher and statistically significant number of osteocytes was found in the peri-implant bone around immediately loaded implants (P=0.0081). A correlation between the percentage of bone-implant contact and osteocyte density was found for immediately loaded implants (P=0.0480) but not for submerged implants (P=0.2667). Conclusion: The higher number of osteocytes in the peri-implant bone around immediately loaded implants could be related to the functional adaptation required by the loading stimulus, which also explains the hypothesized involvement of the osteocytes in the maintenance of the bone matrix. J Periodontol 2009;80:499-504.

Numerical simulation of the liquid phase in SnO(2) thin film deposition by sol-gel-dip-coating

The fluid flow of the liquid phase in the sol-gel-dip-coating process for SnO(2) thin film deposition is numerically simulated. This calculation yields useful information on the velocity distribution close to the substrate, where the film is deposited. The fluid modeling is done by assuming Newtonian behavior, since the linear relation between shear stress and velocity gradient is observed. Besides, very low viscosities are used. The fluid governing equations are the Navier-Stokes in the two dimensional form, discretized by the finite difference technique. Results of optical transmittance and X-ray diffraction on films obtained from colloidal suspensions with regular viscosity, confirm the substrate base as the thickest part of the film, as inferred from the numerical simulation. In addition, as the viscosity increases, the fluid acquires more uniform velocity distribution close to the substrate, leading to more homogenous and uniform films.

Stability of numerical schemes on staggered grids

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthemore, we consider the cases when both explicit and implicit approximations of the boundary conditions arc employed. Why we choose to do this is clearly motivated and arises front solving fluid flow equations with free surfaces when the Reynolds number can be very small. in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t =n delta t rather than t = (n + 1)delta t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar. thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally, stable. Copyright (C) 2008 John Wiley & Sons, Ltd.

An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.