Aerostructural Optimization with Tubular Spar

With aerodynamic- and structural-only analyses done, we now examine an aerostructural design problem. The construction of the problem follows the same logic as outlined in Aerodynamic Optimization, though with some added details. For example, we use an AerostructPoint group instead of an AeroGroup because it contains the additional components needed for aerostructural optimization. Additionally, we have more variable connections due to the more complex problem formulation.

For more details about mesh_dict and surface in the following script, see Mesh and Surface Dictionaries.

import numpy as np

from openaerostruct.geometry.utils import generate_mesh
from openaerostruct.integration.aerostruct_groups import AerostructGeometry, AerostructPoint
from openaerostruct.utils.constants import grav_constant

import openmdao.api as om

# Create a dictionary to store options about the surface
mesh_dict = {"num_y": 5, "num_x": 2, "wing_type": "CRM", "symmetry": True, "num_twist_cp": 5}

mesh, twist_cp = generate_mesh(mesh_dict)

surface = {
    # Wing definition
    "name": "wing",  # name of the surface
    "symmetry": True,  # if true, model one half of wing
    # reflected across the plane y = 0
    "S_ref_type": "wetted",  # how we compute the wing area,
    # can be 'wetted' or 'projected'
    "fem_model_type": "tube",
    "thickness_cp": np.array([0.1, 0.2, 0.3]),
    "twist_cp": twist_cp,
    "mesh": mesh,
    # Aerodynamic performance of the lifting surface at
    # an angle of attack of 0 (alpha=0).
    # These CL0 and CD0 values are added to the CL and CD
    # obtained from aerodynamic analysis of the surface to get
    # the total CL and CD.
    # These CL0 and CD0 values do not vary wrt alpha.
    "CL0": 0.0,  # CL of the surface at alpha=0
    "CD0": 0.015,  # CD of the surface at alpha=0
    # Airfoil properties for viscous drag calculation
    "k_lam": 0.05,  # percentage of chord with laminar
    # flow, used for viscous drag
    "t_over_c_cp": np.array([0.15]),  # thickness over chord ratio (NACA0015)
    "c_max_t": 0.303,  # chordwise location of maximum (NACA0015)
    # thickness
    "with_viscous": True,
    "with_wave": False,  # if true, compute wave drag
    # Structural values are based on aluminum 7075
    "E": 70.0e9,  # [Pa] Young's modulus of the spar
    "G": 30.0e9,  # [Pa] shear modulus of the spar
    "yield": 500.0e6 / 2.5,  # [Pa] yield stress divided by 2.5 for limiting case
    "mrho": 3.0e3,  # [kg/m^3] material density
    "fem_origin": 0.35,  # normalized chordwise location of the spar
    "wing_weight_ratio": 2.0,
    "struct_weight_relief": False,  # True to add the weight of the structure to the loads on the structure
    "distributed_fuel_weight": False,
    # Constraints
    "exact_failure_constraint": False,  # if false, use KS function
}

# Create the problem and assign the model group
prob = om.Problem()

# Add problem information as an independent variables component
indep_var_comp = om.IndepVarComp()
indep_var_comp.add_output("v", val=248.136, units="m/s")
indep_var_comp.add_output("alpha", val=5.0, units="deg")
indep_var_comp.add_output("Mach_number", val=0.84)
indep_var_comp.add_output("re", val=1.0e6, units="1/m")
indep_var_comp.add_output("rho", val=0.38, units="kg/m**3")
indep_var_comp.add_output("CT", val=grav_constant * 17.0e-6, units="1/s")
indep_var_comp.add_output("R", val=11.165e6, units="m")
indep_var_comp.add_output("W0", val=0.4 * 3e5, units="kg")
indep_var_comp.add_output("speed_of_sound", val=295.4, units="m/s")
indep_var_comp.add_output("load_factor", val=1.0)
indep_var_comp.add_output("empty_cg", val=np.zeros((3)), units="m")

prob.model.add_subsystem("prob_vars", indep_var_comp, promotes=["*"])

aerostruct_group = AerostructGeometry(surface=surface)

name = "wing"

# Add tmp_group to the problem with the name of the surface.
prob.model.add_subsystem(name, aerostruct_group)

point_name = "AS_point_0"

# Create the aero point group and add it to the model
AS_point = AerostructPoint(surfaces=[surface])

prob.model.add_subsystem(
    point_name,
    AS_point,
    promotes_inputs=[
        "v",
        "alpha",
        "Mach_number",
        "re",
        "rho",
        "CT",
        "R",
        "W0",
        "speed_of_sound",
        "empty_cg",
        "load_factor",
    ],
)

com_name = point_name + "." + name + "_perf"
prob.model.connect(
    name + ".local_stiff_transformed", point_name + ".coupled." + name + ".local_stiff_transformed"
)
prob.model.connect(name + ".nodes", point_name + ".coupled." + name + ".nodes")

# Connect aerodyamic mesh to coupled group mesh
prob.model.connect(name + ".mesh", point_name + ".coupled." + name + ".mesh")

# Connect performance calculation variables
prob.model.connect(name + ".radius", com_name + ".radius")
prob.model.connect(name + ".thickness", com_name + ".thickness")
prob.model.connect(name + ".nodes", com_name + ".nodes")
prob.model.connect(name + ".cg_location", point_name + "." + "total_perf." + name + "_cg_location")
prob.model.connect(name + ".structural_mass", point_name + "." + "total_perf." + name + "_structural_mass")
prob.model.connect(name + ".t_over_c", com_name + ".t_over_c")

prob.driver = om.ScipyOptimizeDriver()
prob.driver.options["tol"] = 1e-9

# Setup problem and add design variables, constraint, and objective
prob.model.add_design_var("wing.twist_cp", lower=-10.0, upper=15.0)
prob.model.add_design_var("wing.thickness_cp", lower=0.01, upper=0.5, scaler=1e2)
prob.model.add_constraint("AS_point_0.wing_perf.failure", upper=0.0)
prob.model.add_constraint("AS_point_0.wing_perf.thickness_intersects", upper=0.0)

# Add design variables, constraisnt, and objective on the problem
prob.model.add_design_var("alpha", lower=-10.0, upper=10.0)
prob.model.add_constraint("AS_point_0.L_equals_W", equals=0.0)
prob.model.add_objective("AS_point_0.fuelburn", scaler=1e-5)

# Set up the problem
prob.setup(check=True)

# Inserting a small unit test here. Verify that beta is correctly promoted in an Aerostruct
# group.
print(prob["AS_point_0.beta"])
[0.]
prob.run_driver()
INFO: checking out_of_order
INFO: checking system
INFO: checking solvers
INFO: checking dup_inputs
INFO: checking missing_recorders
WARNING: The Problem has no recorder of any kind attached
INFO: checking unserializable_options
INFO: checking comp_has_no_outputs
INFO: checking auto_ivc_warnings

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 6213058.65 1
NL: NLBGS 2 ; 8416718.24 1.35468192
NL: NLBGS 3 ; 56428.312 0.00908221138
NL: NLBGS 4 ; 721.304619 0.000116094932
NL: NLBGS 5 ; 40.9168761 6.58562528e-06
NL: NLBGS 6 ; 1.70725765 2.74785375e-07
NL: NLBGS 7 ; 0.0154420698 2.48542154e-09
NL: NLBGS 8 ; 0.000698482795 1.12421729e-10
NL: NLBGS 9 ; 1.95871709e-05 3.15258104e-12
NL: NLBGS 10 ; 3.34545036e-07 5.38454657e-14
NL: NLBGS 11 ; 1.19197691e-08 1.91850259e-15
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 1.47723283e-09 1
NL: NLBGS Converged
/home/docs/checkouts/readthedocs.org/user_builds/mdolab-openaerostruct/envs/latest/lib/python3.7/site-packages/openmdao/core/total_jac.py:1564: DerivativesWarning:Design variables [('wing.thickness_cp', inds=[1])] have no impact on the constraints or objective.

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 5253321.36 1
NL: NLBGS 2 ; 7051639.76 1.34232027
NL: NLBGS 3 ; 98662.5092 0.0187809773
NL: NLBGS 4 ; 4531.17799 0.000862535847
NL: NLBGS 5 ; 416.907275 7.93607027e-05
NL: NLBGS 6 ; 24.6364461 4.68968952e-06
NL: NLBGS 7 ; 0.772258013 1.47003764e-07
NL: NLBGS 8 ; 0.0311794168 5.9351817e-09
NL: NLBGS 9 ; 0.00274431129 5.22395472e-10
NL: NLBGS 10 ; 0.000137658354 2.6204061e-11
NL: NLBGS 11 ; 4.11694255e-06 7.83683743e-13
NL: NLBGS 12 ; 2.61999232e-07 4.98730639e-14
NL: NLBGS 13 ; 1.93778988e-08 3.68869472e-15
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 554400.711 1
NL: NLBGS 2 ; 747332.392 1.34800042
NL: NLBGS 3 ; 9097.73465 0.0164100343
NL: NLBGS 4 ; 374.587226 0.000675661518
NL: NLBGS 5 ; 40.6952397 7.34040179e-05
NL: NLBGS 6 ; 2.42650774 4.3768121e-06
NL: NLBGS 7 ; 0.0840023587 1.51519212e-07
NL: NLBGS 8 ; 0.003495648 6.30527331e-09
NL: NLBGS 9 ; 0.000366601733 6.61257689e-10
NL: NLBGS 10 ; 9.53570814e-06 1.72000287e-11
NL: NLBGS 11 ; 4.00465125e-07 7.22338765e-13
NL: NLBGS 12 ; 4.99686502e-08 9.01309274e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 23267.4518 1
NL: NLBGS 2 ; 31803.7322 1.36687646
NL: NLBGS 3 ; 209.647674 0.00901034096
NL: NLBGS 4 ; 21.8347374 0.000938424095
NL: NLBGS 5 ; 1.03444069 4.44587013e-05
NL: NLBGS 6 ; 0.0433189361 1.86178256e-06
NL: NLBGS 7 ; 0.00205742708 8.84251141e-08
NL: NLBGS 8 ; 0.000155316022 6.67524845e-09
NL: NLBGS 9 ; 1.12761168e-05 4.84630498e-10
NL: NLBGS 10 ; 3.33782993e-07 1.43454898e-11
NL: NLBGS 11 ; 4.16059913e-08 1.78816278e-12
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 112527.344 1
NL: NLBGS 2 ; 153899.605 1.36766406
NL: NLBGS 3 ; 1031.82158 0.00916951859
NL: NLBGS 4 ; 106.40532 0.000945595232
NL: NLBGS 5 ; 5.36538666 4.76807365e-05
NL: NLBGS 6 ; 0.212013197 1.8841038e-06
NL: NLBGS 7 ; 0.00970379611 8.62350051e-08
NL: NLBGS 8 ; 0.000793447976 7.05115704e-09
NL: NLBGS 9 ; 5.40483884e-05 4.80313374e-10
NL: NLBGS 10 ; 1.6395825e-06 1.45705251e-11
NL: NLBGS 11 ; 7.51214764e-08 6.67584194e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 555127.211 1
NL: NLBGS 2 ; 759242.514 1.36769104
NL: NLBGS 3 ; 5190.27945 0.00934971182
NL: NLBGS 4 ; 552.34134 0.000994981563
NL: NLBGS 5 ; 30.5536703 5.50390427e-05
NL: NLBGS 6 ; 1.17922689 2.12424624e-06
NL: NLBGS 7 ; 0.0543909741 9.79792974e-08
NL: NLBGS 8 ; 0.00462198994 8.3260014e-09
NL: NLBGS 9 ; 0.00032994448 5.94358326e-10
NL: NLBGS 10 ; 1.01113849e-05 1.82145366e-11
NL: NLBGS 11 ; 4.68180948e-07 8.43375966e-13
NL: NLBGS 12 ; 9.19200824e-08 1.65583817e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2603416.5 1
NL: NLBGS 2 ; 3561784.87 1.3681195
NL: NLBGS 3 ; 30061.8416 0.0115470735
NL: NLBGS 4 ; 3782.95725 0.00145307416
NL: NLBGS 5 ; 322.651361 0.000123933824
NL: NLBGS 6 ; 12.5767994 4.83088259e-06
NL: NLBGS 7 ; 0.588512419 2.26053887e-07
NL: NLBGS 8 ; 0.0709349311 2.72468624e-08
NL: NLBGS 9 ; 0.00732973015 2.81542741e-09
NL: NLBGS 10 ; 0.000240812838 9.24987755e-11
NL: NLBGS 11 ; 1.21578784e-05 4.66997057e-12
NL: NLBGS 12 ; 1.5230382e-06 5.85015192e-13
NL: NLBGS 13 ; 1.32494346e-07 5.08924893e-14
NL: NLBGS 14 ; 3.57281452e-08 1.3723561e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 3800983.23 1
NL: NLBGS 2 ; 5223509.11 1.37425208
NL: NLBGS 3 ; 268006.91 0.0705098901
NL: NLBGS 4 ; 77774.5661 0.0204616967
NL: NLBGS 5 ; 29849.9993 0.00785323101
NL: NLBGS 6 ; 2503.6912 0.000658695671
NL: NLBGS 7 ; 393.167565 0.000103438385
NL: NLBGS 8 ; 203.198551 5.34594705e-05
NL: NLBGS 9 ; 68.1665394 1.79339227e-05
NL: NLBGS 10 ; 4.15488617 1.09310826e-06
NL: NLBGS 11 ; 1.72638619 4.54194636e-07
NL: NLBGS 12 ; 0.593375661 1.56111097e-07
NL: NLBGS 13 ; 0.0506676832 1.33301517e-08
NL: NLBGS 14 ; 0.0111271631 2.92744336e-09
NL: NLBGS 15 ; 0.0055479688 1.45961412e-09
NL: NLBGS 16 ; 0.000499827847 1.31499619e-10
NL: NLBGS 17 ; 8.8176119e-05 2.31982394e-11
NL: NLBGS 18 ; 3.45938217e-05 9.10128239e-12
NL: NLBGS 19 ; 1.2257601e-05 3.22485007e-12
NL: NLBGS 20 ; 9.24876353e-07 2.43325555e-13
NL: NLBGS 21 ; 2.8981733e-07 7.62479897e-14
NL: NLBGS 22 ; 1.44794907e-07 3.80940662e-14
NL: NLBGS 23 ; 3.89579258e-08 1.02494338e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 85763.5713 1
NL: NLBGS 2 ; 115322.151 1.34465192
NL: NLBGS 3 ; 8834.73708 0.1030127
NL: NLBGS 4 ; 1581.59145 0.0184412965
NL: NLBGS 5 ; 245.537846 0.002862962
NL: NLBGS 6 ; 146.688504 0.00171038242
NL: NLBGS 7 ; 14.8036867 0.000172610427
NL: NLBGS 8 ; 3.24443038 3.78299356e-05
NL: NLBGS 9 ; 0.954055 1.11242452e-05
NL: NLBGS 10 ; 0.330240913 3.85059657e-06
NL: NLBGS 11 ; 0.02694311 3.14155644e-07
NL: NLBGS 12 ; 0.00420865333 4.90727388e-08
NL: NLBGS 13 ; 0.00242689306 2.82974813e-08
NL: NLBGS 14 ; 0.000801460542 9.34499963e-09
NL: NLBGS 15 ; 4.12706405e-05 4.81214108e-10
NL: NLBGS 16 ; 1.78939762e-05 2.08643086e-10
NL: NLBGS 17 ; 7.53116097e-06 8.78130523e-11
NL: NLBGS 18 ; 6.34216338e-07 7.39493853e-12
NL: NLBGS 19 ; 9.61868704e-08 1.12153527e-12
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 47022.0122 1
NL: NLBGS 2 ; 64395.9755 1.36948575
NL: NLBGS 3 ; 2818.94618 0.0599495013
NL: NLBGS 4 ; 730.501481 0.0155353088
NL: NLBGS 5 ; 209.206581 0.00444912012
NL: NLBGS 6 ; 18.553908 0.000394579201
NL: NLBGS 7 ; 2.38654413 5.07537644e-05
NL: NLBGS 8 ; 1.23296155 2.6220944e-05
NL: NLBGS 9 ; 0.685129037 1.45703896e-05
NL: NLBGS 10 ; 0.0432214071 9.19173916e-07
NL: NLBGS 11 ; 0.0117293793 2.49444435e-07
NL: NLBGS 12 ; 0.00443859272 9.43939339e-08
NL: NLBGS 13 ; 0.000393827218 8.37537994e-09
NL: NLBGS 14 ; 5.3737418e-05 1.14281409e-09
NL: NLBGS 15 ; 2.57043615e-05 5.4664529e-10
NL: NLBGS 16 ; 1.33524835e-05 2.83962401e-10
NL: NLBGS 17 ; 9.16258031e-07 1.94857257e-11
NL: NLBGS 18 ; 2.55386389e-07 5.43120928e-12
NL: NLBGS 19 ; 1.18728281e-07 2.5249511e-12
NL: NLBGS 20 ; 2.58683749e-08 5.50133303e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 231370.302 1
NL: NLBGS 2 ; 316821.103 1.36932484
NL: NLBGS 3 ; 14153.856 0.0611740393
NL: NLBGS 4 ; 3675.67184 0.0158865325
NL: NLBGS 5 ; 1043.26241 0.00450905931
NL: NLBGS 6 ; 93.5424958 0.000404297764
NL: NLBGS 7 ; 11.8401746 5.11741331e-05
NL: NLBGS 8 ; 6.04614456 2.61318955e-05
NL: NLBGS 9 ; 3.52623097 1.52406378e-05
NL: NLBGS 10 ; 0.229766667 9.93068968e-07
NL: NLBGS 11 ; 0.0595427883 2.57348448e-07
NL: NLBGS 12 ; 0.0226253437 9.77884518e-08
NL: NLBGS 13 ; 0.00190636206 8.23944144e-09
NL: NLBGS 14 ; 0.000258031525 1.11523183e-09
NL: NLBGS 15 ; 0.000123858276 5.35324866e-10
NL: NLBGS 16 ; 7.74749711e-05 3.34852703e-10
NL: NLBGS 17 ; 5.47437943e-06 2.36606832e-11
NL: NLBGS 18 ; 1.34510138e-06 5.81363023e-12
NL: NLBGS 19 ; 5.3048015e-07 2.29277545e-12
NL: NLBGS 20 ; 5.52633918e-08 2.38852572e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 1106121.73 1
NL: NLBGS 2 ; 1513563.38 1.36835155
NL: NLBGS 3 ; 74301.5364 0.0671730194
NL: NLBGS 4 ; 19720.9318 0.0178288982
NL: NLBGS 5 ; 5294.38173 0.00478643678
NL: NLBGS 6 ; 513.160087 0.000463927317
NL: NLBGS 7 ; 60.6480332 5.48294384e-05
NL: NLBGS 8 ; 28.4821347 2.57495481e-05
NL: NLBGS 9 ; 20.888531 1.88844776e-05
NL: NLBGS 10 ; 1.63338853 1.47668063e-06
NL: NLBGS 11 ; 0.359665283 3.25158863e-07
NL: NLBGS 12 ; 0.132307108 1.19613515e-07
NL: NLBGS 13 ; 0.0206409622 1.86606606e-08
NL: NLBGS 14 ; 0.00205446269 1.85735677e-09
NL: NLBGS 15 ; 0.000696830521 6.29976342e-10
NL: NLBGS 16 ; 0.00052394578 4.73678227e-10
NL: NLBGS 17 ; 5.39323413e-05 4.87580524e-11
NL: NLBGS 18 ; 1.40994131e-05 1.27467102e-11
NL: NLBGS 19 ; 3.97091452e-06 3.5899435e-12
NL: NLBGS 20 ; 4.8771558e-07 4.40923966e-13
NL: NLBGS 21 ; 1.04857736e-07 9.47976458e-14
NL: NLBGS 22 ; 5.33599631e-08 4.82405884e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 3068391.77 1
NL: NLBGS 2 ; 4174991.59 1.36064489
NL: NLBGS 3 ; 919552.839 0.299685603
NL: NLBGS 4 ; 91157.6081 0.0297085949
NL: NLBGS 5 ; 17979.2009 0.00585948674
NL: NLBGS 6 ; 3062.2674 0.000998004045
NL: NLBGS 7 ; 770.560272 0.000251128386
NL: NLBGS 8 ; 201.18644 6.556739e-05
NL: NLBGS 9 ; 45.2357673 1.47425005e-05
NL: NLBGS 10 ; 11.9128743 3.88244891e-06
NL: NLBGS 11 ; 0.38204926 1.24511239e-07
NL: NLBGS 12 ; 0.156114229 5.08781931e-08
NL: NLBGS 13 ; 0.0565311508 1.84237069e-08
NL: NLBGS 14 ; 0.0146690284 4.78068951e-09
NL: NLBGS 15 ; 0.00118201881 3.85224215e-10
NL: NLBGS 16 ; 0.000230273785 7.50470612e-11
NL: NLBGS 17 ; 4.92244644e-05 1.60424314e-11
NL: NLBGS 18 ; 1.18672933e-05 3.86759391e-12
NL: NLBGS 19 ; 2.19522373e-06 7.1543137e-13
NL: NLBGS 20 ; 1.04929553e-06 3.41969218e-13
NL: NLBGS 21 ; 8.06583109e-08 2.62868359e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2305451.69 1
NL: NLBGS 2 ; 3082351.38 1.33698372
NL: NLBGS 3 ; 247095.803 0.107178911
NL: NLBGS 4 ; 55764.319 0.0241880232
NL: NLBGS 5 ; 5465.22203 0.00237056454
NL: NLBGS 6 ; 2405.56307 0.00104342376
NL: NLBGS 7 ; 395.985943 0.000171760677
NL: NLBGS 8 ; 39.8323922 1.72774786e-05
NL: NLBGS 9 ; 11.2179315 4.86582803e-06
NL: NLBGS 10 ; 5.06076293 2.19512859e-06
NL: NLBGS 11 ; 0.440589326 1.91107594e-07
NL: NLBGS 12 ; 0.0619956605 2.68908955e-08
NL: NLBGS 13 ; 0.0239070247 1.0369779e-08
NL: NLBGS 14 ; 0.0057981577 2.51497687e-09
NL: NLBGS 15 ; 0.000481658958 2.08921731e-10
NL: NLBGS 16 ; 0.000123106489 5.3397991e-11
NL: NLBGS 17 ; 5.42358638e-05 2.35250489e-11
NL: NLBGS 18 ; 7.68720137e-06 3.33435803e-12
NL: NLBGS 19 ; 1.10288768e-06 4.7838247e-13
NL: NLBGS 20 ; 2.12015275e-07 9.19625753e-14
NL: NLBGS 21 ; 3.80273801e-08 1.64945465e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2859960.77 1
NL: NLBGS 2 ; 3882062.97 1.3573833
NL: NLBGS 3 ; 812630.627 0.284140481
NL: NLBGS 4 ; 80718.1361 0.0282235117
NL: NLBGS 5 ; 23210.6752 0.00811573203
NL: NLBGS 6 ; 3265.70269 0.00114186975
NL: NLBGS 7 ; 815.152446 0.000285022248
NL: NLBGS 8 ; 262.595522 9.18178756e-05
NL: NLBGS 9 ; 49.8559573 1.7432392e-05
NL: NLBGS 10 ; 13.1924189 4.61279717e-06
NL: NLBGS 11 ; 2.5202635 8.81223103e-07
NL: NLBGS 12 ; 0.858626606 3.00223212e-07
NL: NLBGS 13 ; 0.186476732 6.52025491e-08
NL: NLBGS 14 ; 0.0705080694 2.46535093e-08
NL: NLBGS 15 ; 0.0114059125 3.988136e-09
NL: NLBGS 16 ; 0.00324398568 1.13427629e-09
NL: NLBGS 17 ; 0.000304655891 1.065245e-10
NL: NLBGS 18 ; 0.000105351063 3.68365412e-11
NL: NLBGS 19 ; 2.84565515e-05 9.94997969e-12
NL: NLBGS 20 ; 1.0288675e-05 3.59748815e-12
NL: NLBGS 21 ; 1.74030653e-06 6.08507135e-13
NL: NLBGS 22 ; 5.69003678e-07 1.98955064e-13
NL: NLBGS 23 ; 1.94250154e-07 6.79205662e-14
NL: NLBGS 24 ; 4.50923418e-08 1.57667694e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2276873.28 1
NL: NLBGS 2 ; 3044342.38 1.3370715
NL: NLBGS 3 ; 241293.463 0.105975798
NL: NLBGS 4 ; 55402.0349 0.0243325069
NL: NLBGS 5 ; 5282.43844 0.00232004059
NL: NLBGS 6 ; 2157.14604 0.000947415939
NL: NLBGS 7 ; 338.916912 0.000148851899
NL: NLBGS 8 ; 34.6354621 1.52118532e-05
NL: NLBGS 9 ; 9.0352003 3.96824908e-06
NL: NLBGS 10 ; 4.17389659 1.83317035e-06
NL: NLBGS 11 ; 0.372718289 1.63697423e-07
NL: NLBGS 12 ; 0.0553799689 2.43228156e-08
NL: NLBGS 13 ; 0.0188031565 8.25832366e-09
NL: NLBGS 14 ; 0.00391935329 1.72137524e-09
NL: NLBGS 15 ; 0.000357433243 1.56984249e-10
NL: NLBGS 16 ; 6.70702704e-05 2.94571819e-11
NL: NLBGS 17 ; 3.90295595e-05 1.71417355e-11
NL: NLBGS 18 ; 5.67636133e-06 2.49305105e-12
NL: NLBGS 19 ; 8.35047206e-07 3.66751726e-13
NL: NLBGS 20 ; 2.04797757e-07 8.9946928e-14
NL: NLBGS 21 ; 3.17179419e-08 1.39304819e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2626505.07 1
NL: NLBGS 2 ; 3572167.52 1.36004592
NL: NLBGS 3 ; 802162.556 0.305410625
NL: NLBGS 4 ; 68960.0813 0.0262554533
NL: NLBGS 5 ; 20444.846 0.0077840497
NL: NLBGS 6 ; 3008.37342 0.00114539029
NL: NLBGS 7 ; 783.229778 0.000298202271
NL: NLBGS 8 ; 240.931193 9.17307166e-05
NL: NLBGS 9 ; 51.5504072 1.9626997e-05
NL: NLBGS 10 ; 11.0241999 4.19728862e-06
NL: NLBGS 11 ; 1.96543193 7.48306921e-07
NL: NLBGS 12 ; 0.590938836 2.24990556e-07
NL: NLBGS 13 ; 0.158835422 6.04740588e-08
NL: NLBGS 14 ; 0.0351033545 1.33650435e-08
NL: NLBGS 15 ; 0.00918848358 3.49836887e-09
NL: NLBGS 16 ; 0.00290379317 1.10557303e-09
NL: NLBGS 17 ; 0.000550426602 2.09566167e-10
NL: NLBGS 18 ; 6.72340719e-05 2.55983027e-11
NL: NLBGS 19 ; 9.66668126e-06 3.68043502e-12
NL: NLBGS 20 ; 2.02226244e-06 7.69944234e-13
NL: NLBGS 21 ; 3.00814653e-07 1.14530391e-13
NL: NLBGS 22 ; 7.74291212e-08 2.94799054e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2085385.17 1
NL: NLBGS 2 ; 2783773.3 1.33489647
NL: NLBGS 3 ; 230629.67 0.11059332
NL: NLBGS 4 ; 52110.8099 0.024988578
NL: NLBGS 5 ; 5203.58137 0.00249526152
NL: NLBGS 6 ; 2058.43478 0.000987076542
NL: NLBGS 7 ; 247.685072 0.000118771858
NL: NLBGS 8 ; 29.7956705 1.428785e-05
NL: NLBGS 9 ; 7.99616442 3.83438251e-06
NL: NLBGS 10 ; 3.71236407 1.78018149e-06
NL: NLBGS 11 ; 0.334953456 1.60619468e-07
NL: NLBGS 12 ; 0.0445727123 2.13738512e-08
NL: NLBGS 13 ; 0.0160632693 7.70278296e-09
NL: NLBGS 14 ; 0.00374070853 1.79377344e-09
NL: NLBGS 15 ; 0.000339411166 1.62757063e-10
NL: NLBGS 16 ; 6.6763467e-05 3.20149333e-11
NL: NLBGS 17 ; 3.30488696e-05 1.58478491e-11
NL: NLBGS 18 ; 5.32801611e-06 2.55493143e-12
NL: NLBGS 19 ; 7.05090674e-07 3.38110525e-13
NL: NLBGS 20 ; 1.10811099e-07 5.31369939e-14
NL: NLBGS 21 ; 2.2228853e-08 1.06593512e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 1415054.57 1
NL: NLBGS 2 ; 1898318.13 1.34151585
NL: NLBGS 3 ; 234771.779 0.165910053
NL: NLBGS 4 ; 27859.4813 0.0196879201
NL: NLBGS 5 ; 6265.65293 0.00442785251
NL: NLBGS 6 ; 714.956177 0.000505249898
NL: NLBGS 7 ; 116.873225 8.25927332e-05
NL: NLBGS 8 ; 21.7636815 1.53801005e-05
NL: NLBGS 9 ; 3.20071956 2.26190539e-06
NL: NLBGS 10 ; 0.642310106 4.53911898e-07
NL: NLBGS 11 ; 0.126459866 8.93674837e-08
NL: NLBGS 12 ; 0.0324580023 2.2937633e-08
NL: NLBGS 13 ; 0.00394279429 2.7863196e-09
NL: NLBGS 14 ; 0.0011216791 7.92675508e-10
NL: NLBGS 15 ; 0.000157340558 1.11190452e-10
NL: NLBGS 16 ; 3.09106563e-05 2.18441443e-11
NL: NLBGS 17 ; 4.36675577e-06 3.08592747e-12
NL: NLBGS 18 ; 1.36998158e-06 9.68147527e-13
NL: NLBGS 19 ; 2.16775413e-07 1.53192264e-13
NL: NLBGS 20 ; 5.39051072e-08 3.8094013e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 1165600.31 1
NL: NLBGS 2 ; 1558200.88 1.33682264
NL: NLBGS 3 ; 116396.531 0.0998597291
NL: NLBGS 4 ; 26266.1245 0.0225344179
NL: NLBGS 5 ; 2470.19159 0.0021192441
NL: NLBGS 6 ; 822.679351 0.000705798843
NL: NLBGS 7 ; 93.9625319 8.06129948e-05
NL: NLBGS 8 ; 12.2024005 1.04687691e-05
NL: NLBGS 9 ; 2.12587581 1.82384629e-06
NL: NLBGS 10 ; 1.04521284 8.96716337e-07
NL: NLBGS 11 ; 0.104942063 9.0032631e-08
NL: NLBGS 12 ; 0.017771513 1.52466611e-08
NL: NLBGS 13 ; 0.00455797591 3.91041068e-09
NL: NLBGS 14 ; 0.000667069278 5.72296757e-10
NL: NLBGS 15 ; 8.33127573e-05 7.14762654e-11
NL: NLBGS 16 ; 9.35484409e-06 8.0257735e-12
NL: NLBGS 17 ; 2.3659869e-06 2.02984409e-12
NL: NLBGS 18 ; 1.19303797e-06 1.02353951e-12
NL: NLBGS 19 ; 1.13005107e-07 9.69501347e-14
NL: NLBGS 20 ; 2.19255228e-08 1.88104984e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 748023.952 1
NL: NLBGS 2 ; 1002770.96 1.34055996
NL: NLBGS 3 ; 83232.0379 0.111269215
NL: NLBGS 4 ; 14241.6623 0.0190390458
NL: NLBGS 5 ; 2492.33963 0.00333189816
NL: NLBGS 6 ; 386.605401 0.000516835591
NL: NLBGS 7 ; 47.8104682 6.39156916e-05
NL: NLBGS 8 ; 8.29756009 1.10926396e-05
NL: NLBGS 9 ; 0.97099336 1.2980779e-06
NL: NLBGS 10 ; 0.191491733 2.55996794e-07
NL: NLBGS 11 ; 0.0222199454 2.97048582e-08
NL: NLBGS 12 ; 0.00339850791 4.5433143e-09
NL: NLBGS 13 ; 0.000455699589 6.09204542e-10
NL: NLBGS 14 ; 0.000100309366 1.34099137e-10
NL: NLBGS 15 ; 1.22514472e-05 1.63784157e-11
NL: NLBGS 16 ; 1.27900441e-06 1.7098442e-12
NL: NLBGS 17 ; 2.7643403e-07 3.69552378e-13
NL: NLBGS 18 ; 1.04512043e-07 1.39717509e-13
NL: NLBGS 19 ; 2.97014554e-08 3.97065567e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 502281.116 1
NL: NLBGS 2 ; 671686.401 1.33727186
NL: NLBGS 3 ; 48848.6336 0.0972535739
NL: NLBGS 4 ; 10644.5873 0.0211924896
NL: NLBGS 5 ; 1111.67374 0.00221325012
NL: NLBGS 6 ; 295.585999 0.000588487183
NL: NLBGS 7 ; 27.9909159 5.57275897e-05
NL: NLBGS 8 ; 4.48555313 8.93036387e-06
NL: NLBGS 9 ; 0.659979327 1.31396405e-06
NL: NLBGS 10 ; 0.219026175 4.3606293e-07
NL: NLBGS 11 ; 0.0229194024 4.56306273e-08
NL: NLBGS 12 ; 0.00225306379 4.48566294e-09
NL: NLBGS 13 ; 0.000743041291 1.47933352e-09
NL: NLBGS 14 ; 0.000228732535 4.55387487e-10
NL: NLBGS 15 ; 2.29514359e-05 4.56944033e-11
NL: NLBGS 16 ; 2.44931861e-06 4.87638999e-12
NL: NLBGS 17 ; 1.18397052e-06 2.357187e-12
NL: NLBGS 18 ; 1.30871505e-07 2.60554301e-13
NL: NLBGS 19 ; 1.59690533e-08 3.17930593e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 190569.711 1
NL: NLBGS 2 ; 254981.881 1.33799794
NL: NLBGS 3 ; 29410.1424 0.154327476
NL: NLBGS 4 ; 2737.29815 0.014363763
NL: NLBGS 5 ; 570.148064 0.00299180842
NL: NLBGS 6 ; 127.881523 0.000671048526
NL: NLBGS 7 ; 14.0462654 7.37067047e-05
NL: NLBGS 8 ; 1.15632639 6.06773439e-06
NL: NLBGS 9 ; 0.365274236 1.91674865e-06
NL: NLBGS 10 ; 0.0829301083 4.35169409e-07
NL: NLBGS 11 ; 0.0081141095 4.25781698e-08
NL: NLBGS 12 ; 0.00128519638 6.74396982e-09
NL: NLBGS 13 ; 0.000223871767 1.17474999e-09
NL: NLBGS 14 ; 2.90163966e-05 1.52261325e-10
NL: NLBGS 15 ; 3.76322809e-06 1.9747252e-11
NL: NLBGS 16 ; 9.92640591e-07 5.20880568e-12
NL: NLBGS 17 ; 1.16504625e-07 6.11349123e-13
NL: NLBGS 18 ; 1.81297207e-08 9.51343244e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 205063.806 1
NL: NLBGS 2 ; 282899.296 1.37956718
NL: NLBGS 3 ; 15550.9158 0.0758345223
NL: NLBGS 4 ; 2881.26278 0.0140505671
NL: NLBGS 5 ; 379.003053 0.00184822012
NL: NLBGS 6 ; 40.1512905 0.000195799011
NL: NLBGS 7 ; 3.43435421 1.67477346e-05
NL: NLBGS 8 ; 0.684766356 3.33928433e-06
NL: NLBGS 9 ; 0.168578374 8.22077659e-07
NL: NLBGS 10 ; 0.0158528357 7.73068443e-08
NL: NLBGS 11 ; 0.00161305867 7.86613055e-09
NL: NLBGS 12 ; 0.000315245504 1.53730446e-09
NL: NLBGS 13 ; 4.11380337e-05 2.00610895e-10
NL: NLBGS 14 ; 2.42305225e-06 1.18160893e-11
NL: NLBGS 15 ; 7.80616566e-07 3.80670085e-12
NL: NLBGS 16 ; 8.78510191e-08 4.28408214e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 167833.485 1
NL: NLBGS 2 ; 233270.014 1.38988959
NL: NLBGS 3 ; 17903.8465 0.106676248
NL: NLBGS 4 ; 2136.05555 0.012727231
NL: NLBGS 5 ; 243.875058 0.00145307748
NL: NLBGS 6 ; 26.4634176 0.000157676625
NL: NLBGS 7 ; 1.90213317 1.13334545e-05
NL: NLBGS 8 ; 0.299359241 1.78366814e-06
NL: NLBGS 9 ; 0.069726559 4.15450819e-07
NL: NLBGS 10 ; 0.00714786694 4.25890395e-08
NL: NLBGS 11 ; 0.000431959866 2.57374066e-09
NL: NLBGS 12 ; 0.000118148218 7.03960941e-10
NL: NLBGS 13 ; 1.46540352e-05 8.73129406e-11
NL: NLBGS 14 ; 6.75309439e-07 4.02368715e-12
NL: NLBGS 15 ; 2.14111392e-07 1.27573703e-12
NL: NLBGS 16 ; 3.03270787e-08 1.80697426e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 236679.875 1
NL: NLBGS 2 ; 329164.226 1.39075714
NL: NLBGS 3 ; 25938.0054 0.10959109
NL: NLBGS 4 ; 2768.50203 0.0116972431
NL: NLBGS 5 ; 257.570067 0.00108826349
NL: NLBGS 6 ; 27.5474691 0.000116391261
NL: NLBGS 7 ; 2.01997173 8.53461552e-06
NL: NLBGS 8 ; 0.22989372 9.71327706e-07
NL: NLBGS 9 ; 0.0438032484 1.85073819e-07
NL: NLBGS 10 ; 0.00446660443 1.88719232e-08
NL: NLBGS 11 ; 0.000231935746 9.79955504e-10
NL: NLBGS 12 ; 6.03531361e-05 2.5499902e-10
NL: NLBGS 13 ; 8.10404421e-06 3.42405294e-11
NL: NLBGS 14 ; 3.78151009e-07 1.597732e-12
NL: NLBGS 15 ; 7.48045364e-08 3.16057867e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 226068.888 1
NL: NLBGS 2 ; 314521.931 1.39126589
NL: NLBGS 3 ; 25964.2458 0.114851035
NL: NLBGS 4 ; 2554.64066 0.0113002753
NL: NLBGS 5 ; 205.453523 0.00090880937
NL: NLBGS 6 ; 21.6755974 9.58804976e-05
NL: NLBGS 7 ; 1.59578059 7.05882446e-06
NL: NLBGS 8 ; 0.145715356 6.44561741e-07
NL: NLBGS 9 ; 0.0238555351 1.055233e-07
NL: NLBGS 10 ; 0.00237694342 1.05142439e-08
NL: NLBGS 11 ; 0.000127120068 5.62306777e-10
NL: NLBGS 12 ; 2.83774211e-05 1.25525549e-10
NL: NLBGS 13 ; 3.86822315e-06 1.71108161e-11
NL: NLBGS 14 ; 1.93505211e-07 8.55956843e-13
NL: NLBGS 15 ; 1.88309802e-08 8.32975309e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 267216.9 1
NL: NLBGS 2 ; 371885.917 1.39170059
NL: NLBGS 3 ; 33744.8917 0.126282775
NL: NLBGS 4 ; 2994.02542 0.0112044763
NL: NLBGS 5 ; 215.887071 0.000807909496
NL: NLBGS 6 ; 22.2152977 8.31358261e-05
NL: NLBGS 7 ; 1.61611875 6.04796611e-06
NL: NLBGS 8 ; 0.123431948 4.61916698e-07
NL: NLBGS 9 ; 0.0174454892 6.52858753e-08
NL: NLBGS 10 ; 0.00166440816 6.22867853e-09
NL: NLBGS 11 ; 0.00010051137 3.76141516e-10
NL: NLBGS 12 ; 1.77969848e-05 6.66012698e-11
NL: NLBGS 13 ; 2.22661963e-06 8.3326303e-12
NL: NLBGS 14 ; 1.1085212e-07 4.14839481e-13
NL: NLBGS 15 ; 2.79256982e-08 1.04505734e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 376124.656 1
NL: NLBGS 2 ; 523210.394 1.39105582
NL: NLBGS 3 ; 55489.4207 0.147529336
NL: NLBGS 4 ; 4078.77844 0.0108442198
NL: NLBGS 5 ; 280.397766 0.000745491585
NL: NLBGS 6 ; 27.3810413 7.27977834e-05
NL: NLBGS 7 ; 1.92131539 5.10818782e-06
NL: NLBGS 8 ; 0.123898093 3.29406996e-07
NL: NLBGS 9 ; 0.0131114646 3.48593594e-08
NL: NLBGS 10 ; 0.00101941297 2.71030615e-09
NL: NLBGS 11 ; 7.36209105e-05 1.95735401e-10
NL: NLBGS 12 ; 6.53615259e-06 1.73776233e-11
NL: NLBGS 13 ; 4.58935179e-07 1.22016776e-12
NL: NLBGS 14 ; 5.94721986e-08 1.58118319e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 608192.872 1
NL: NLBGS 2 ; 843297.883 1.38656325
NL: NLBGS 3 ; 112600.801 0.185139954
NL: NLBGS 4 ; 5437.66782 0.00894069638
NL: NLBGS 5 ; 471.460619 0.000775182743
NL: NLBGS 6 ; 44.1050623 7.25182164e-05
NL: NLBGS 7 ; 3.05806984 5.02812509e-06
NL: NLBGS 8 ; 0.317653096 5.22290067e-07
NL: NLBGS 9 ; 0.0462342074 7.60189891e-08
NL: NLBGS 10 ; 0.00482365055 7.93111983e-09
NL: NLBGS 11 ; 0.000527996116 8.6813927e-10
NL: NLBGS 12 ; 7.52390554e-05 1.23709203e-10
NL: NLBGS 13 ; 2.33265731e-06 3.83539073e-12
NL: NLBGS 14 ; 2.80092364e-07 4.60532138e-13
NL: NLBGS 15 ; 3.52883314e-08 5.80216129e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 785411.511 1
NL: NLBGS 2 ; 1074251.18 1.36775583
NL: NLBGS 3 ; 164977.678 0.210052535
NL: NLBGS 4 ; 35876.6448 0.0456787865
NL: NLBGS 5 ; 7075.35162 0.00900846438
NL: NLBGS 6 ; 636.178748 0.000809994174
NL: NLBGS 7 ; 164.163424 0.000209015811
NL: NLBGS 8 ; 99.9361301 0.00012724047
NL: NLBGS 9 ; 16.2225666 2.06548623e-05
NL: NLBGS 10 ; 1.58637744 2.01980416e-06
NL: NLBGS 11 ; 0.644996156 8.21220656e-07
NL: NLBGS 12 ; 0.108043967 1.37563513e-07
NL: NLBGS 13 ; 0.0264509448 3.36778166e-08
NL: NLBGS 14 ; 0.00355528068 4.52664702e-09
NL: NLBGS 15 ; 0.0016942092 2.15709749e-09
NL: NLBGS 16 ; 0.000393552113 5.01077597e-10
NL: NLBGS 17 ; 7.16933475e-05 9.12812538e-11
NL: NLBGS 18 ; 7.53039227e-06 9.58783029e-12
NL: NLBGS 19 ; 1.54965296e-06 1.9730459e-12
NL: NLBGS 20 ; 4.6397939e-07 5.90746868e-13
NL: NLBGS 21 ; 6.75145998e-08 8.59607975e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 161710.975 1
NL: NLBGS 2 ; 225205.432 1.39264161
NL: NLBGS 3 ; 25446.863 0.157360149
NL: NLBGS 4 ; 3338.23543 0.0206432212
NL: NLBGS 5 ; 914.400249 0.00565453427
NL: NLBGS 6 ; 103.06815 0.000637360266
NL: NLBGS 7 ; 23.8026686 0.00014719266
NL: NLBGS 8 ; 2.96193892 1.83162517e-05
NL: NLBGS 9 ; 0.771479692 4.77073182e-06
NL: NLBGS 10 ; 0.112592933 6.96260307e-07
NL: NLBGS 11 ; 0.00966251544 5.97517605e-08
NL: NLBGS 12 ; 0.00361039029 2.23261921e-08
NL: NLBGS 13 ; 0.00167634843 1.03663244e-08
NL: NLBGS 14 ; 0.000263006451 1.62639828e-09
NL: NLBGS 15 ; 3.31929214e-05 2.05260784e-10
NL: NLBGS 16 ; 1.17461246e-05 7.2636533e-11
NL: NLBGS 17 ; 2.90480749e-06 1.79629583e-11
NL: NLBGS 18 ; 3.52882568e-07 2.1821807e-12
NL: NLBGS 19 ; 1.04172315e-07 6.44188284e-13
NL: NLBGS 20 ; 1.1892248e-08 7.35401418e-14
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 198172.536 1
NL: NLBGS 2 ; 276080.745 1.39313322
NL: NLBGS 3 ; 36153.6273 0.182435104
NL: NLBGS 4 ; 4316.1075 0.0217795441
NL: NLBGS 5 ; 1311.81257 0.00661954778
NL: NLBGS 6 ; 156.528689 0.000789860653
NL: NLBGS 7 ; 33.6692723 0.000169898781
NL: NLBGS 8 ; 3.56216624 1.79750752e-05
NL: NLBGS 9 ; 1.09381798 5.51952353e-06
NL: NLBGS 10 ; 0.556221534 2.80675387e-06
NL: NLBGS 11 ; 0.0735723943 3.7125424e-07
NL: NLBGS 12 ; 0.0105649349 5.33118016e-08
NL: NLBGS 13 ; 0.00329497112 1.66267798e-08
NL: NLBGS 14 ; 0.000587907942 2.96664691e-09
NL: NLBGS 15 ; 3.56364724e-05 1.79825485e-10
NL: NLBGS 16 ; 1.65051416e-05 8.32867253e-11
NL: NLBGS 17 ; 1.0999262e-05 5.55034625e-11
NL: NLBGS 18 ; 1.54328807e-06 7.78759812e-12
NL: NLBGS 19 ; 1.80288185e-07 9.09753631e-13
NL: NLBGS 20 ; 7.64522306e-08 3.85786204e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 12128.0352 1
NL: NLBGS 2 ; 16876.8002 1.39155271
NL: NLBGS 3 ; 2353.80129 0.194079359
NL: NLBGS 4 ; 269.539313 0.0222244831
NL: NLBGS 5 ; 82.5324065 0.00680509295
NL: NLBGS 6 ; 12.4151744 0.00102367566
NL: NLBGS 7 ; 2.24962336 0.000185489515
NL: NLBGS 8 ; 0.332773665 2.74383823e-05
NL: NLBGS 9 ; 0.0790721281 6.51978059e-06
NL: NLBGS 10 ; 0.0339098742 2.79599076e-06
NL: NLBGS 11 ; 0.00439348427 3.62258537e-07
NL: NLBGS 12 ; 0.000707718 5.83538876e-08
NL: NLBGS 13 ; 0.000199226608 1.64269484e-08
NL: NLBGS 14 ; 3.48952392e-05 2.87723764e-09
NL: NLBGS 15 ; 2.24516219e-06 1.85121676e-10
NL: NLBGS 16 ; 1.25402546e-06 1.03398897e-10
NL: NLBGS 17 ; 4.87612543e-07 4.02054032e-11
NL: NLBGS 18 ; 6.33820826e-08 5.22608005e-12
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 402.994881 1
NL: NLBGS 2 ; 560.768764 1.39150344
NL: NLBGS 3 ; 70.1487165 0.174068505
NL: NLBGS 4 ; 9.62057244 0.0238726914
NL: NLBGS 5 ; 2.76141879 0.00685224284
NL: NLBGS 6 ; 0.343123843 0.000851434744
NL: NLBGS 7 ; 0.0683439211 0.000169590048
NL: NLBGS 8 ; 0.0143729312 3.56652948e-05
NL: NLBGS 9 ; 0.00176770152 4.38641185e-06
NL: NLBGS 10 ; 0.000568593559 1.41092005e-06
NL: NLBGS 11 ; 8.25730495e-05 2.04898507e-07
NL: NLBGS 12 ; 6.4866007e-06 1.60959878e-08
NL: NLBGS 13 ; 2.872908e-06 7.12889453e-09
NL: NLBGS 14 ; 1.39906346e-06 3.47166559e-09
NL: NLBGS 15 ; 2.20919129e-07 5.48193387e-10
NL: NLBGS 16 ; 3.03256508e-08 7.52507096e-11
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 358.06785 1
NL: NLBGS 2 ; 494.935081 1.38223826
NL: NLBGS 3 ; 49.2032137 0.137413101
NL: NLBGS 4 ; 15.941309 0.0445203583
NL: NLBGS 5 ; 2.44724567 0.00683458645
NL: NLBGS 6 ; 0.502053219 0.00140211756
NL: NLBGS 7 ; 0.0724275979 0.00020227339
NL: NLBGS 8 ; 0.0272758149 7.61749901e-05
NL: NLBGS 9 ; 0.00407333666 1.13758793e-05
NL: NLBGS 10 ; 0.000366795747 1.02437498e-06
NL: NLBGS 11 ; 0.000183428522 5.12273086e-07
NL: NLBGS 12 ; 3.0982012e-05 8.65255341e-08
NL: NLBGS 13 ; 4.91611129e-06 1.37295523e-08
NL: NLBGS 14 ; 1.70105421e-06 4.75064772e-09
NL: NLBGS 15 ; 2.67851103e-07 7.48045665e-10
NL: NLBGS 16 ; 1.87502719e-08 5.23651367e-11
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 597.801042 1
NL: NLBGS 2 ; 811.101901 1.35680911
NL: NLBGS 3 ; 165.087411 0.276157783
NL: NLBGS 4 ; 56.2269285 0.0940562571
NL: NLBGS 5 ; 17.2274127 0.0288179704
NL: NLBGS 6 ; 2.45585872 0.00410815397
NL: NLBGS 7 ; 0.366749414 0.000613497448
NL: NLBGS 8 ; 0.136923557 0.000229045364
NL: NLBGS 9 ; 0.0263428763 4.40662936e-05
NL: NLBGS 10 ; 0.00241229517 4.03528098e-06
NL: NLBGS 11 ; 0.0010212935 1.70841706e-06
NL: NLBGS 12 ; 0.000202901701 3.39413429e-07
NL: NLBGS 13 ; 3.75068714e-05 6.27413951e-08
NL: NLBGS 14 ; 6.52381568e-06 1.09130216e-08
NL: NLBGS 15 ; 2.72762039e-06 4.56275616e-09
NL: NLBGS 16 ; 6.61197726e-07 1.1060498e-09
NL: NLBGS 17 ; 1.12117081e-07 1.87549157e-10
NL: NLBGS 18 ; 2.30907774e-08 3.86261912e-11
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2627.04709 1
NL: NLBGS 2 ; 3569.87543 1.35889282
NL: NLBGS 3 ; 863.214544 0.328587388
NL: NLBGS 4 ; 223.899681 0.0852286515
NL: NLBGS 5 ; 82.4987147 0.0314035919
NL: NLBGS 6 ; 12.4118817 0.0047246514
NL: NLBGS 7 ; 1.41324082 0.00053795793
NL: NLBGS 8 ; 0.740687133 0.000281946652
NL: NLBGS 9 ; 0.122913993 4.67878908e-05
NL: NLBGS 10 ; 0.0104685521 3.98491223e-06
NL: NLBGS 11 ; 0.00516959002 1.96783302e-06
NL: NLBGS 12 ; 0.000815057002 3.10255954e-07
NL: NLBGS 13 ; 9.23834896e-05 3.51662861e-08
NL: NLBGS 14 ; 4.19672264e-05 1.59750567e-08
NL: NLBGS 15 ; 6.02262031e-06 2.29254372e-09
NL: NLBGS 16 ; 5.93946021e-07 2.26088837e-10
NL: NLBGS 17 ; 2.97630234e-07 1.13294594e-10
NL: NLBGS 18 ; 6.69583023e-08 2.54880479e-11
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 12987.0497 1
NL: NLBGS 2 ; 17667.4909 1.36039296
NL: NLBGS 3 ; 4331.44526 0.333520342
NL: NLBGS 4 ; 1067.11778 0.0821678364
NL: NLBGS 5 ; 409.343143 0.0315193329
NL: NLBGS 6 ; 62.1176191 0.00478304314
NL: NLBGS 7 ; 6.5229629 0.000502266723
NL: NLBGS 8 ; 3.49469227 0.000269090544
NL: NLBGS 9 ; 0.536107907 4.12801921e-05
NL: NLBGS 10 ; 0.0444787626 3.42485503e-06
NL: NLBGS 11 ; 0.0264832142 2.03920172e-06
NL: NLBGS 12 ; 0.00428715206 3.30109775e-07
NL: NLBGS 13 ; 0.000415152896 3.19666826e-08
NL: NLBGS 14 ; 0.000219919733 1.69337716e-08
NL: NLBGS 15 ; 3.73736384e-05 2.87776202e-09
NL: NLBGS 16 ; 2.46440957e-06 1.89759e-10
NL: NLBGS 17 ; 1.3392646e-06 1.03123082e-10
NL: NLBGS 18 ; 4.37205364e-07 3.36647178e-11
NL: NLBGS 19 ; 5.91667029e-08 4.55582324e-12
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 63434.4715 1
NL: NLBGS 2 ; 86296.0887 1.36039738
NL: NLBGS 3 ; 21133.7378 0.33315857
NL: NLBGS 4 ; 5197.75953 0.0819390374
NL: NLBGS 5 ; 1996.03514 0.0314660955
NL: NLBGS 6 ; 302.165235 0.00476342323
NL: NLBGS 7 ; 31.4598159 0.000495941956
NL: NLBGS 8 ; 16.8219477 0.00026518622
NL: NLBGS 9 ; 2.53114735 3.99017646e-05
NL: NLBGS 10 ; 0.216082381 3.40638734e-06
NL: NLBGS 11 ; 0.129564086 2.04248704e-06
NL: NLBGS 12 ; 0.0206507343 3.25544358e-07
NL: NLBGS 13 ; 0.00199843769 3.150397e-08
NL: NLBGS 14 ; 0.00107798406 1.69936635e-08
NL: NLBGS 15 ; 0.000191592499 3.02032151e-09
NL: NLBGS 16 ; 1.31296806e-05 2.06980216e-10
NL: NLBGS 17 ; 6.95993503e-06 1.097185e-10
NL: NLBGS 18 ; 1.73899508e-06 2.74140391e-11
NL: NLBGS 19 ; 2.27163553e-07 3.58107426e-12
NL: NLBGS 20 ; 6.92997804e-08 1.09246249e-12
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 282068.848 1
NL: NLBGS 2 ; 383394.674 1.35922374
NL: NLBGS 3 ; 92568.0451 0.328175358
NL: NLBGS 4 ; 23455.5405 0.083155374
NL: NLBGS 5 ; 8833.7906 0.0313178526
NL: NLBGS 6 ; 1314.23207 0.00465925989
NL: NLBGS 7 ; 139.218686 0.000493562787
NL: NLBGS 8 ; 73.2970899 0.000259855317
NL: NLBGS 9 ; 10.7455073 3.80953351e-05
NL: NLBGS 10 ; 0.954167756 3.38274774e-06
NL: NLBGS 11 ; 0.555507131 1.96940263e-06
NL: NLBGS 12 ; 0.089481212 3.17231813e-07
NL: NLBGS 13 ; 0.00931942322 3.30395338e-08
NL: NLBGS 14 ; 0.00485995442 1.72296744e-08
NL: NLBGS 15 ; 0.000933080995 3.30799024e-09
NL: NLBGS 16 ; 6.67848144e-05 2.36767778e-10
NL: NLBGS 17 ; 3.00887303e-05 1.06671582e-10
NL: NLBGS 18 ; 6.17422243e-06 2.18890618e-11
NL: NLBGS 19 ; 9.44458661e-07 3.34832673e-12
NL: NLBGS 20 ; 2.4823254e-07 8.80042381e-13
NL: NLBGS 21 ; 3.66572244e-08 1.29958429e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 247637.821 1
NL: NLBGS 2 ; 336460.36 1.35867921
NL: NLBGS 3 ; 80715.9055 0.325943368
NL: NLBGS 4 ; 20795.8173 0.0839767413
NL: NLBGS 5 ; 7841.35671 0.0316646168
NL: NLBGS 6 ; 1151.31405 0.00464918504
NL: NLBGS 7 ; 121.531624 0.00049076358
NL: NLBGS 8 ; 63.7306242 0.000257354163
NL: NLBGS 9 ; 9.15653542 3.69755127e-05
NL: NLBGS 10 ; 0.893820346 3.60938544e-06
NL: NLBGS 11 ; 0.494211728 1.99570375e-06
NL: NLBGS 12 ; 0.0756009918 3.05288552e-07
NL: NLBGS 13 ; 0.00780553917 3.151998e-08
NL: NLBGS 14 ; 0.00421065754 1.70032894e-08
NL: NLBGS 15 ; 0.000723533461 2.92174053e-09
NL: NLBGS 16 ; 5.57011197e-05 2.24929776e-10
NL: NLBGS 17 ; 2.80304798e-05 1.13191433e-10
NL: NLBGS 18 ; 5.16998657e-06 2.08772091e-11
NL: NLBGS 19 ; 6.62908309e-07 2.67692675e-12
NL: NLBGS 20 ; 2.20991706e-07 8.92398847e-13
NL: NLBGS 21 ; 3.92646791e-08 1.58556875e-13
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 2741.68026 1
NL: NLBGS 2 ; 3755.06501 1.36962178
NL: NLBGS 3 ; 939.923361 0.34282749
NL: NLBGS 4 ; 155.856395 0.0568470354
NL: NLBGS 5 ; 83.427472 0.0304293222
NL: NLBGS 6 ; 12.6636334 0.00461893154
NL: NLBGS 7 ; 1.36323855 0.0004972274
NL: NLBGS 8 ; 0.535225605 0.000195218098
NL: NLBGS 9 ; 0.13776091 5.02468914e-05
NL: NLBGS 10 ; 0.0240064586 8.75611169e-06
NL: NLBGS 11 ; 0.00341768921 1.24656739e-06
NL: NLBGS 12 ; 0.00153534796 5.60002558e-07
NL: NLBGS 13 ; 0.000219949989 8.02245223e-08
NL: NLBGS 14 ; 2.23166464e-05 8.13976988e-09
NL: NLBGS 15 ; 9.64896112e-06 3.51936047e-09
NL: NLBGS 16 ; 1.62861719e-06 5.94021563e-10
NL: NLBGS 17 ; 3.12998243e-07 1.14162927e-10
NL: NLBGS 18 ; 6.54395567e-08 2.3868413e-11
NL: NLBGS Converged

==================
AS_point_0.coupled
==================
NL: NLBGS 1 ; 267.481945 1
NL: NLBGS 2 ; 373.566684 1.39660523
NL: NLBGS 3 ; 56.6413585 0.211757688
NL: NLBGS 4 ; 6.18836755 0.0231356458
NL: NLBGS 5 ; 2.18620181 0.00817326867
NL: NLBGS 6 ; 0.286694653 0.00107182805
NL: NLBGS 7 ; 0.0340437936 0.000127275108
NL: NLBGS 8 ; 0.0121480556 4.54163574e-05
NL: NLBGS 9 ; 0.00439650227 1.64366319e-05
NL: NLBGS 10 ; 0.000622350496 2.32670095e-06
NL: NLBGS 11 ; 0.000104372258 3.90203002e-07
NL: NLBGS 12 ; 2.8714497e-05 1.0735116e-07
NL: NLBGS 13 ; 5.98056732e-06 2.23587701e-08
NL: NLBGS 14 ; 9.1077044e-07 3.40497912e-09
NL: NLBGS 15 ; 2.82472617e-07 1.05604368e-09
NL: NLBGS 16 ; 4.2818422e-08 1.60079672e-10
NL: NLBGS Converged
Optimization terminated successfully    (Exit mode 0)
            Current function value: [0.92369793]
            Iterations: 39
            Function evaluations: 44
            Gradient evaluations: 39
Optimization Complete
-----------------------------------
print(prob["AS_point_0.fuelburn"][0])
92369.79279575823