Rheology and Microfluidic flows of a cationic surfactant and hydrotropic salt mixture

Thesis (Ph.D.)--University of Washington, 2016-08

Upplevelser av "Employer Branding" : En kvalitativ fallstudie av anställda hos BA

Syftet med denna studie var att undersöka hur medarbetarna på BA (fiktivt namn) upplever organisationens interna Employer Branding-arbete samt att belysa de aspekter som eventuellt skiljer strategi mot upplevelse. Det gjordes ett målinriktat urval av organisation medan valet av deltagare var ett slumpmässigt urval vilket resulterade i åtta respondenter (n= 8). Studien var en kvalitativ fallstudie med ett psykologiskt angreppssätt och innehöll en ostrukturerad intervju som tillsammans med dokumentation utgjorde intervjuguiden för den primära datainsamlingen genom åtta semistrukturerade intervjuer med medarbetarna. Fördelningen av urvalet var fyra kvinnor och fyra män med varierade yrkesroller inom en bemanningsdivision. Både den ostrukturerade intervjun och de semistrukturerade intervjuerna analyserades med hjälp av analysmetoden tematisk analys. Studiens resultat visade att BA:s interna Employer Branding och medarbetarnas upplevelser i stora drag överensstämde. Värderingar och interna karriärmöjligheter var centrala aspekter som både organisationen och respondenterna belyste. 

Physical and chemical kinetics of structural build-up of cement suspensions

Abstract : The structural build-up of fresh cement-based materials has a great impact on their structural performance after casting. Accordingly, the mixture design should be tailored to adapt the kinetics of build-up given the application on hand. The rate of structural build-up of cement-based suspensions at rest is a complex phenomenon affected by both physical and chemical structuration processes. The structuration kinetics are strongly dependent on the mixture’s composition, testing parameters, as well as the shear history. Accurate measurements of build-up rely on the efficiency of the applied pre-shear regime to achieve an initial well-dispersed state as well as the applied stress during the liquid-solid transition. Studying the physical and chemical mechanisms of build-up of cement suspensions at rest can enhance the fundamental understanding of this phenomenon. This can, therefore, allow a better control of the rheological and time-dependent properties of cement-based materials. The research focused on the use of dynamic rheology in investigating the kinetics of structural build-up of fresh cement pastes. The research program was conducted in three different phases. The first phase was devoted to evaluating the dispersing efficiency of various disruptive shear techniques. The investigated shearing profiles included rotational, oscillatory, and combination of both. The initial and final states of suspension’s structure, before and after disruption, were determined by applying a small-amplitude oscillatory shear (SAOS). The difference between the viscoelastic values before and after disruption was used to express the degree of dispersion. An efficient technique to disperse concentrated cement suspensions was developed. The second phase aimed to establish a rheometric approach to dissociate and monitor the individual physical and chemical mechanisms of build-up of cement paste. In this regard, the non-destructive dynamic rheometry was used to investigate the evolutions of both storage modulus and phase angle of inert calcium carbonate and cement suspensions. Two independent build-up indices were proposed. The structural build-up of various cement suspensions made with different cement contents, silica fume replacement percentages, and high-range water reducer dosages was evaluated using the proposed indices. These indices were then compared to the well-known thixotropic index (Athix.). Furthermore, the proposed indices were correlated to the decay in lateral pressure determined for various cement pastes cast in a pressure column. The proposed pre-shearing protocol and build-up indices (phases 1 and 2) were then used to investigate the effect of mixture’s parameters on the kinetics of structural build-up in phase 3. The investigated mixture’s parameters included cement content and fineness, alkali sulfate content, and temperature of cement suspension. Zeta potential, calorimetric, spectrometric measurements were performed to explore the corresponding microstructural changes in cement suspensions, such as inter-particle cohesion, rate of Brownian flocculation, and nucleation rate. A model linking the build-up indices and the microstructural characteristics was developed to predict the build-up behaviour of cement-based suspensions The obtained results showed that oscillatory shear may have a greater effect on dispersing concentrated cement suspension than the rotational shear. Furthermore, the increase in induced shear strain was found to enhance the breakdown of suspension’s structure until a critical point, after which thickening effects dominate. An effective dispersing method is then proposed. This consists of applying a rotational shear around the transitional value between the linear and non-linear variations of the apparent viscosity with shear rate, followed by an oscillatory shear at the crossover shear strain and high angular frequency of 100 rad/s. Investigating the evolutions of viscoelastic properties of inert calcite-based and cement suspensions and allowed establishing two independent build-up indices. The first one (the percolation time) can represent the rest time needed to form the elastic network. On the other hand, the second one (rigidification rate) can describe the increase in stress-bearing capacity of formed network due to cement hydration. In addition, results showed that combining the percolation time and the rigidification rate can provide deeper insight into the structuration process of cement suspensions. Furthermore, these indices were found to be well-correlated to the decay in the lateral pressure of cement suspensions. The variations of proposed build-up indices with mixture’s parameters showed that the percolation time is most likely controlled by the frequency of Brownian collisions, distance between dispersed particles, and intensity of cohesion between cement particles. On the other hand, a higher rigidification rate can be secured by increasing the number of contact points per unit volume of paste, nucleation rate of cement hydrates, and intensity of inter-particle cohesion.

An experimental investigation by towing tank on VIV of a long flexible cylinder for deepwater riser application

In this study, the dynamic response of a vertical flexible cylinder vibrating at low mode numbers with combined x-y motion was investigated in a towing tank. The uniform flow was simulated by towing the flexible cylinder along the tank in still water; therefore, the turbulence intensity of the free flow was negligible in obtaining more reliable results. A lower branch of dominant frequencies with micro vibration amplitude was found in both cross-flow and in-line directions. This justifiable discrepancy was likely caused by an initial lock-in. The maximum attainable amplitude, modal analysis and x-y trajectory in cross-flow and in-line directions are reported here and compared with previous literature, along with some good agreements and different observations that were obtained from the study. Drag and lift coefficients are also evaluated by making use of a generalized integral transform technique approach, yielding an alternative method to study fluid force acting upon a flexible cylinder.

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

Three dimensional modeling of contaminants transportation (oil spill) in Caspian Sea

Nowadays, Caspian Sea is in focus of more attentions than past because of its individualistic as the biggest lake in the world and the existing of very large oil and gas resources within it. Very large scale of oil pollution caused by development of oil exploration and excavation activities not only make problem for coastal facilities but also make severe damage on environment. In the first stage of this research, the location and quality of oil resources in offshore and onshore have been determined and then affected depletion factors on oil spill such as evaporation, emulsification, dissolution, sedimentation and so on have been studied. In second stage, sea hydrodynamics model is offered and tested by determination of governing hydrodynamic equations on sea currents and on pollution transportation in sea surface and by finding out main parameters in these equations such as Coriolis, bottom friction, wind and etc. this model has been calculated by using cell vertex finite volume method in an unstructured mesh domain. According to checked model; sea currents of Caspian Sea in different seasons of the year have been determined and in final stage different scenarios of oil spill movement in Caspian sea on various conditions have been investigated by modeling of three dimensional oil spill movement on surface (affected by sea currents) and on depth (affected by buoyancy, drag and gravity forces) by applying main above mentioned depletion factors.

Large eddy simulation of shock boundary layer interaction control using micro-vortex generators

The performance of supersonic engine inlets and external aerodynamic surfaces can be critically affected by shock wave / boundary layer interactions (SBLIs), whose severe adverse pressure gradients can cause boundary layer separation. Currently such problems are avoided primarily through the use of boundary layer bleed/suction which can be a source of significant performance degradation. This study investigates a novel type of flow control device called micro-vortex generators (µVGs) which may offer similar control benefits without the bleed penalties. µVGs have the ability to alter the near-wall structure of compressible turbulent boundary layers to provide increased mixing of high speed fluid which improves the boundary layer health when subjected to flow disturbance. Due to their small size,µVGs are embedded in the boundary layer which provide reduced drag compared to the traditional vortex generators while they are cost-effective, physically robust and do not require a power source. To examine the potential of µVGs, a detailed experimental and computational study of micro-ramps in a supersonic boundary layer at Mach 3 subjected to an oblique shock was undertaken. The experiments employed a flat plate boundary layer with an impinging oblique shock with downstream total pressure measurements. The moderate Reynolds number of 3,800 based on displacement thickness allowed the computations to use Large Eddy Simulations without the subgrid stress model (LES-nSGS). The LES predictions indicated that the shock changes the structure of the turbulent eddies and the primary vortices generated from the micro-ramp. Furthermore, they generally reproduced the experimentally obtained mean velocity profiles, unlike similarly-resolved RANS computations. The experiments and the LES results indicate that the micro-ramps, whose height is h≈0.5δ, can significantly reduce boundary layer thickness and improve downstream boundary layer health as measured by the incompressible shape factor, H. Regions directly behind the ramp centerline tended to have increased boundary layer thickness indicating the significant three-dimensionality of the flow field. Compared to baseline sizes, smaller micro-ramps yielded improved total pressure recovery. Moving the smaller ramps closer to the shock interaction also reduced the displacement thickness and the separated area. This effect is attributed to decreased wave drag and the closer proximity of the vortex pairs to the wall. In the second part of the study, various types of µVGs are investigated including micro-ramps and micro-vanes. The results showed that vortices generated from µVGs can partially eliminate shock induced flow separation and can continue to entrain high momentum flux for boundary layer recovery downstream. The micro-ramps resulted in thinner downstream displacement thickness in comparison to the micro-vanes. However, the strength of the streamwise vorticity for the micro-ramps decayed faster due to dissipation especially after the shock interaction. In addition, the close spanwise distance between each vortex for the ramp geometry causes the vortex cores to move upwards from the wall due to induced upwash effects. Micro-vanes, on the other hand, yielded an increased spanwise spacing of the streamwise vortices at the point of formation. This resulted in streamwise vortices staying closer to the wall with less circulation decay, and the reduction in overall flow separation is attributed to these effects. Two hybrid concepts, named “thick-vane” and “split-ramp”, were also studied where the former is a vane with side supports and the latter has a uniform spacing along the centerline of the baseline ramp. These geometries behaved similar to the micro-vanes in terms of the streamwise vorticity and the ability to reduce flow separation, but are more physically robust than the thin vanes. Next, Mach number effect on flow past the micro-ramps (h~0.5δ) are examined in a supersonic boundary layer at M=1.4, 2.2 and 3.0, but with no shock waves present. The LES results indicate that micro-ramps have a greater impact at lower Mach number near the device but its influence decays faster than that for the higher Mach number cases. This may be due to the additional dissipation caused by the primary vortices with smaller effective diameter at the lower Mach number such that their coherency is easily lost causing the streamwise vorticity and the turbulent kinetic energy to decay quickly. The normal distance between the vortex core and the wall had similar growth indicating weak correlation with the Mach number; however, the spanwise distance between the two counter-rotating cores further increases with lower Mach number. Finally, various µVGs which include micro-ramp, split-ramp and a new hybrid concept “ramped-vane” are investigated under normal shock conditions at Mach number of 1.3. In particular, the ramped-vane was studied extensively by varying its size, interior spacing of the device and streamwise position respect to the shock. The ramped-vane provided increased vorticity compared to the micro-ramp and the split-ramp. This significantly reduced the separation length downstream of the device centerline where a larger ramped-vane with increased trailing edge gap yielded a fully attached flow at the centerline of separation region. The results from coarse-resolution LES studies show that the larger ramped-vane provided the most reductions in the turbulent kinetic energy and pressure fluctuation compared to other devices downstream of the shock. Additional benefits include negligible drag while the reductions in displacement thickness and shape factor were seen compared to other devices. Increased wall shear stress and pressure recovery were found with the larger ramped-vane in the baseline resolution LES studies which also gave decreased amplitudes of the pressure fluctuations downstream of the shock.

Investigating wetting characteristics on microstructured surfaces for superhydrophobicity and metal microcasting

The engineering of liquid behavior on surfaces is important for infrastructure, transportation, manufacturing, and sensing. Surfaces can be rendered superhydrophobic by microstructuring, and superhydrophobic devices could lead to practical corrosion inhibition, self-cleaning, fluid flow control, and surface drag reduction. To more fully understand how liquid interacts with microstructured surfaces, this dissertation introduces a direct method for determining droplet solid-liquid-vapor interfacial geometry on microstructured surfaces. The technique performs metrology on molten metal droplets deposited onto microstructured surfaces and then frozen. Unlike other techniques, this visualization technique can be used on large areas of curved and opaque microstructured surfaces to determine contact line. This dissertation also presents measurements and models for how curvature and flexing of microstructured polymers affects hydrophobicity. Increasing curvature of microstructured surfaces leads to decreased slide angle for liquid droplets suspended on the surface asperities. For a surface with regularly spaced asperities, as curvature becomes more positive, droplets suspended on the tops of asperities are suspended on fewer asperities. Curvature affects superhydrophobicity because microscopic curvature changes solid-liquid interaction, pitch is altered, and curvature changes the shape of the three phase contact line. This dissertation presents a model of droplet interactions with curved microstructured surfaces that can be used to design microstructure geometries that maintain the suspension of a droplet when curved surfaces are covered with microstructured polymers. Controlling droplet dynamics could improve microfluidic devices and the shedding of liquids from expensive equipment, preventing corrosion and detrimental performance. This dissertation demonstrates redirection of dynamic droplet spray with anisotropic microstructures. Superhydrophobic microstructured surfaces can be economically fabricated using metal embossing masters, so this dissertation describes casting-based microfabrication of metal microstructures and nanostructures. Low melting temperature metal was cast into flexible silicone molds which were themselves cast from microfabricated silicon templates. The flexibility of the silicone mold permits casting of curved surfaces, which this dissertation demonstrates by fabricating a cylindrical metal roller with microstructures. The metal microstructures can be in turn used as a reusable molding tool. This dissertation also describes an industrial investment casting process to produce aluminum molds having integrated microstructures. Unlike conventional micromolding tools, the aluminum mold was large and had complex curved surfaces. The aluminum was cast into curved microstructured ceramic molds which were themselves cast from curved microstructured rubber. Many structures were successfully cast into the aluminum with excellent replication fidelity, including circular, square, and triangular holes. This dissertation demonstrates molding of large, curved surfaces having surface microstructures using the aluminum mold. This work contributes a more full understanding of the phenomenon of superhydrophobicity and techniques for the economic fabrication of superhydrophobic microstructures.

Tidal Energy ”Round Robin” Tests Comparisons between towing tank and circulating tank results

One key step of the industrial development of a tidal energy device is the testing of scale prototype devices within a controlled laboratory environment. At present, there is no available experimental protocol which addresses in a quantitative manner the differences which can be expected between results obtained from the different types of facilities currently employed for this type of testing. As a consequence, where differences between results are found it has been difficult to confirm the extent to which these differences relate to the device performance or to the test facility type. In the present study, a comparative ”Round Robin” testing programme has been conducted as part of the EC FP VII MaRINET program in order to evaluate the impact of different experimental facilities on the test results. The aim of the trials was to test the same model tidal turbine in four different test facilities to explore the sensitivity of the results to the choice of facility. The facilities comprised two towing tanks, of very different size, and two circulating water channels. Performance assessments in terms of torque, drag and inflow speed showed very similar results in all facilities. However, expected differences between the different tank types (circulating and towing) were observed in the fluctuations of torque and drag measurements. The main facility parameters which can influence the behaviour of the turbine were identified; in particular the effect of blockage was shown to be significant in cases yielding for high thrust coefficients, even at relatively small blockage ratios.

Free-flight experiments in LISA Pathfinder

The LISA Path finder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Path finder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Path finder system that are not present in the full LISA con figuration. While LISA Path finder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this 'suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

Effect of Geometric Uncertainties on the Aerodynamic Characteristic of Offshore Wind Turbine Blades

Offshore wind turbines operate in a complex unsteady flow environment which causes unsteady aerodynamic loads. The unsteady flow environment is characterized by a high degree of uncertainty. In addition, geometry variations and material imperfections also cause uncertainties in the design process. Probabilistic design methods consider these uncertainties in order to reach acceptable reliability and safety levels for offshore wind turbines. Variations of the rotor blade geometry influence the aerodynamic loads which also affect the reliability of other wind turbine components. Therefore, the present paper is dealing with geometric uncertainties of the rotor blades. These can arise from manufacturing tolerances and operational wear of the blades. First, the effect of geometry variations of wind turbine airfoils on the lift and drag coefficients are investigated using a Latin hypercube sampling. Then, the resulting effects on the performance and the blade loads of an offshore wind turbine are analyzed. The variations of the airfoil geometry lead to a significant scatter of the lift and drag coefficients which also affects the damage-equivalent flapwise bending moments. In contrast to that, the effects on the power and the annual energy production are almost negligible with regard to the assumptions made.

Modelling the relationship between biomechanics and performance of young sprinting swimmers

The aim of this study was to compute a swimming performance confirmatory model based on biomechanical parameters. The sample included 100 young swimmers (overall: 12.3 ± 0.74 years; 49 boys: 12.5 ± 0.76 years; 51 girls: 12.2 ± 0.71 years; both genders in Tanner stages 1-2 by self-report) participating on a regular basis in regional and national-level events. The 100 m freestyle event was chosen as the performance indicator. Anthropometric (arm span), strength (throwing velocity), power output (power to overcome drag), kinematic (swimming velocity) and efficiency (propelling efficiency) parameters were measured and included in the model. The path-flow analysis procedure was used to design and compute the model. The anthropometric parameter (arm span) was excluded in the final model, increasing its goodness-of-fit. The final model included the throw velocity, power output, swimming velocity and propelling efficiency. All links were significant between the parameters included, but the throw velocity-power output. The final model was explained by 69% presenting a reasonable adjustment (model's goodness-of-fit; x(2)/df = 3.89). This model shows that strength and power output parameters do play a mediator and meaningful role in the young swimmers' performance.

The influence of nano-scale second-phase particles on deformation of fine grained calcite mylonites

Grey and white carbonate mylonites were collected along thrust planes of the Helvetic Alps. They are characterised by very small grain sizes and non-random grain shape (SPO) and crystallographic preferred orientation (CPO). Presumably they deformed in the field of grain size sensitive flow by recrystallisation accommodated intracrystalline deformation in combination with granular flow. Both mylonites show a similar mean grain size, but in the grey mylonites the grain size range is larger, the grain shapes are more elongate and the dynamically recrystallised calcite grains are more often twinned. Grey mylonites have an oblique CPO, while the CPO in white mylonites is symmetric with respect to the shear plane. Combustion analysis and TEM investigations revealed that grey mylonites contain a higher amount of highly structured kerogens with particle sizes of a few tens of nanometers, which are finely dispersed at the grain boundaries. During deformation of the rock, nano-scale particles reduced the migration velocity of grain boundaries by Zener drag resulting in slower recrystallisation rates of the calcite aggregate. In the grey mylonites, more strain increments were accommodated by individual grains before they became refreshed by dynamic recrystallisation than in white mylonites, where grain boundary migration was less hindered and recrystallisation cycles were faster. Consequently, grey mylonites represent ‘deformation’ microfabrics while white mylonites are characterised by ‘recrystallisation’ microfabrics. Field geologists must utilise this different deformation behavior when applying the obliquity in CPO and SPO of the respective mylonites as reliable shear sense indicators.

U-Pb geochronology and Nd isotope contributions to the interpretation of a peculiar ring massif: the Santa Eulália plutonic complex (SW Iberia, Portugal)

The Santa Eulalia plutonic complex (SEPC) is a late-Variscan granitic body placed in the Ossa-Morena Zone. The host rocks of the complex belong to metamorphic formations from Proterozoic to Lower Paleozoic. The SEPC is a ring massif (ca. 400 km2 area) composed by two main granitic facies with different colours and textures. From the rim to the core, there is (i) a peripheral pink medium- to coarse-grained granite (G0 group) involving large elongated masses of mafic and intermediate rocks, from gabbros to granodiorites (M group), and (ii) a central gray medium-grained granite (G1 group). The mafic to intermediate rocks (M group) are metaluminous and show wide compositions: 3.34–13.51 wt% MgO; 0.70–7.20 ppm Th; 0.84–1.06 (Eu/Eu*)N (Eu* calculated between Sm and Tb); 0.23–0.97 (Nb/Nb*)N (Nb* calculated between Th and La). Although involving the M-type bodies and forming the outer ring, the G0 granites are the most differentiated magmatic rocks of the SEPC, with a transitional character between metaluminous and peraluminous: 0.00–0.62 wt% MgO; 15.00–56.00 ppm Th; and 0.19–0.42 (Eu/Eu*)N ; 0.08–0.19 (Nb/Nb*)N [1][2]. The G1 group is composed by monzonitic granites with a dominant peraluminous character and represents the most homogeneous compositional group of the SEPC: 0.65–1.02 wt% MgO; 13.00–16.95 ppm Th; 0.57–0.70 (Eu/Eu*)N ; 0.14–0.16 (Nb/Nb*)N . According to the SiO2 vs. (Na2O+K2O–CaO) relationships, the M and G1 groups predominantly fall in the calc-alkaline field, while the G0 group is essencially alkali-calcic; on the basis of the SiO2 vs. FeOt/(FeOt+MgO) correlation, SEPC should be considered as a magnesian plutonic association [3]. New geochronological data (U-Pb on zircons) slightly correct the age of the SEPC, previously obtained by other methods (290 Ma, [4]). They provide ages of 306  2 Ma for the M group, 305  6 Ma for the G1 group, and 301  4 Ma for the G0 group, which confirm the late-Variscan character of the SEPC, indicating however a faintly older emplacement, during the Upper Carboniferous. Recent whole-rock isotopic data show that the Rb-Sr system suffered significant post-magmatic disturbance, but reveal a consistent set of Sm-Nd results valuable in the approach to the magmatic sources of this massif: M group (2.9 < Ndi < +1.8); G1 group (5.8 < Ndi < 4.6); G0 group (2.2 < Ndi < 0.8). These geochemical data suggest a petrogenetic model for the SEPC explained by a magmatic event developed in two stages. Initially, magmas derived from long-term depleted mantle sources (Ndi < +1.8 in M group) were extracted to the crust promoting its partial melting and extensive mixing and/or AFC magmatic evolution, thereby generating the G1 granites (Ndi < 4.6). Subsequently, a later extraction of similar primary magmas in the same place or nearby, could have caused partial melting of some intermediate facies (e.g. diorites) of the M group, followed by magmatic differentiation processes, mainly fractional crystallization, able to produce residual liquids compositionally close to the G0 granites (Ndi < 0.8). The kinetic energy associated with the structurally controlled (cauldron subsidence type?) motion of the G0 liquids to the periphery, would have been strong enough to drag up M group blocks as those occurring inside the G0 granitic ring.

An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken.